AlphaShape: incorrect geometry

26 次查看（过去 30 天）
Irene Yang
Irene Yang 2019-3-12
评论： Steven Lord 2019-3-12
Hi,
I have a 3D CAD file which I am importing into Matlab via stlread.
I then use alphaShape using the points generated from stlread of the geometry. However, the alphaShape is not accurate to my geometry and adds a whole bunch of triangles that do not match my geometry.
This is what it is supposed to be:
trisurf.PNG
However, this is what I get instead:
alphaShape.PNG
How can I fix this? Does anyone have an idea on how I could possibly go in and edit the triangulation (is this even feasible?)? Or alternatively, another function in matlab that can allow me the same thing as alphaShape? Should I edit the mesh further? I have tried everything that I can think of!
Thanks.
  4 个评论
显示 2更早的评论
Jan
Jan 2019-3-12
编辑：Jan 2019-3-12
alphaShape is used to define a surface based on a point cloud. As Sean has written already, you do have a surface already, so there is no need for alphaShape.
But your problem seems to be, that you used inapropriate values for 'RegionThreshold' or 'HoleThreshold' or 'Alpha' radius. Which values did you use? Please post your code.
Irene Yang
Irene Yang 2019-3-12
编辑：Irene Yang 2019-3-12
Hi Jan,
My primary goal of using alphaShape is not to obtain a surface but to use the function inshape. However, perhaps I am able to use inshape without creating a surface using alphashape but perhaps using trisurf or trimesh should be sufficient?
My code is very basic:
model=stlread('model.STL');
shp = alphaShape(model.Points(:,1),model.Points(:,2),model.Points(:,3));
plot(shp)
Therefore, I have been using the default values for 'RegionThreshold' (which is 0) or 'HoleThreshold' (which is 0) or 'Alpha' radius (which is the critical radius of 2.5). However, I have tried all alpharadiuses of 0-3 in 0.1 incremements but this still doesnt work as the function alphaShape does not deal well with concave regions and adds in all the additional surfaces.
Thanks

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Sean de Wolski
Sean de Wolski 2019-3-12
How about pointlocation?
  2 个评论
Irene Yang
Irene Yang 2019-3-12
编辑：Irene Yang 2019-3-12
Oh this seems quite promising. I will have a look to see if it can be done in 3D as opposed to the example which is in 2D. Thanks for the suggestion!
Steven Lord
Steven Lord 2019-3-12
See the "Enclosing Tetrahedra" example.
x = gallery('uniformdata',[20 1],0);
y = gallery('uniformdata',[20 1],1);
z = gallery('uniformdata',[20 1],2);
TR = delaunayTriangulation(x,y,z);
P = [0.7 0.6 0.3; 0.5 0.5 0.5];
[ID,B] = pointLocation(TR,P)
I'll plot the delaunayTriangulation using tetramesh and turn off the face colors.
h = tetramesh(TR, 'FaceColor', 'none');
Now I'll add the two points inside.
hold on;
plot3(P(:, 1), P(:, 2), P(:, 3), 'ro', 'MarkerFaceColor', 'r')
Finally I'll color the tetrahedron containing the first point in blue and make it somewhat transparent.
set(h(ID(1)), 'FaceColor', 'b', 'FaceAlpha', 0.25)
If you rotate the axes around you should be able to convince yourself that one of the red dots is inside the blue shape. It looks like it's very close to the surface of the shape but inside.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Bounding Regions 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by