isSelfIntersecting function gives apparently incorrect result

2023 6 23
1 个回答
回答已采纳
更新时间：2023 6 28
10 次查看（30 天）

0 个投票

I'm trying to use the isSelfIntersecting function from the Lidar Toolbox to check whether or not a surfaceMesh object is self intersecting. In most cases it works fine, but in certain situations it returns true when the surface is not self intersecting.
The problem seems to arise in situations where all vertices in the mesh are contained, or nearly contained, by one plane (as in the example below). I wondered if it might be an issue with leaving the face normals undefined, but I get the same issue even if I define them when calling surfaceMesh (not included in example for readability).
Is this an issue with the way I'm defining the mesh and using the function? Or is this a known issue with isSelfIntersecting?
Many thanks for any help!
% define the vertices for the mesh
vertices = [0 0.0785926524819716 0.0150000000000000; ...
0.00223394141607314 0.0787352254127876 0.0140425952680485; ...
0.00428764473102951 0.0788687035475108 0.0131624353567936; ...
0 0.0822253462122803 0.0150000000000000; ...
0.00224235572430482 0.0822979726488426 0.0140389888997977; ...
0.00430286261362254 0.0823659662099516 0.0131559129739490; ...
0 0.0860000000000000 0.0150000000000000; ...
0.00225109241982507 0.0860000000000000 0.0140352443565181; ...
0.00431866355685131 0.0860000000000000 0.0131491406913786];
%define the faces for the mesh
faces = [1 2 5;
1 4 5;
2 3 6;
2 5 6;
4 5 8;
4 7 8;
5 6 9;
5 8 9];
%create the surface mesh object
mesh = surfaceMesh(vertices,faces);
%plot the mesh and verify visually that it does not self intersect
%(it should be an inclined plane sloping downward along the x-axis)
patch('vertices',vertices,'faces',faces,'EdgeColor',[1 .5 0],'FaceColor',[.5 1 .5]);
view(3)
%call isSelfintersecting, returns true despite mesh appearing to be
%non-self intersecting
isSelfIntersecting(mesh)

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

 采纳的回答

colordepth
colordepth 2023-6-27

0 个投票

Hi Ed,
No, you have defined the mesh correctly and are using isSelfIntersection just fine. This is likely due to floating point errors. On inputting your vertices data, it appears slightly differently in MATLAB workspace variables due to the precision limitations of double.
If you look closely, you can see some differences from your original input. For example, in the 2nd column 1st row, there's a diff of 4e-16. Unfortunately, this seems to be a limitation with floating point computation. Hope this helps!

3 个评论
显示 1更早的评论 隐藏 1更早的评论

Ed Wheatcroft
Ed Wheatcroft 2023-6-27
Hi Deep,
Thanks for your answer, and for finding this discrepancy.
I tried reducing the precision of vertices using:
decPlaces = 10
vertices = round(vertices,decPlaces);
That fixed the issue so long as I kept decPlaces at or below 10, so this would seem to confirm your theory.
I'm a bit confused as to how this came about in the first place - the mesh in the example is a very small piece of a much larger mesh which exhibited this same issue, which was in turn generated by another MATLAB script. I pasted this piece directly out of the workspace to give a minumum working example, so I'm wondering where the extra decimal place of precison appeared from in vertices(1,2) for example.
colordepth
colordepth 2023-6-27
I am not aware why that would happen. This might require more information about your system setup and these scripts. Perhaps the output was copied from a computer with a different hardware platform? You may create a new question with more information about these scripts so it can be identified why this is happening.
Ed Wheatcroft
Ed Wheatcroft 2023-6-28
编辑：Ed Wheatcroft 2023-6-28
Ok, no worries :) I actually found another configuration which gives an apparently erroneous result while I was playing around with this. This time, the mesh is clearly self intersecting, but isSelfIntersecting returns false.
In this new example the entries in vertices should be exactly representable using double (or indeed single) precision, so I wonder if this is even a separate issue.
Either way, it seems like isSelfIntersecting isn't totally reliable.
vertices = [0 0 0;
1 0 0;
0 1 0;
0 0 1;
0.5 0.5 1;
0.5 0.5 -1];
faces = [1 2 3;
1 2 4;
1 4 5;
1 5 6];
mesh = surfaceMesh(vertices,faces);
patch('vertices',vertices,'faces',faces,'EdgeColor',[1 .5 0],'FaceColor',[.5 1 .5]);
view(3)
isSelfIntersecting(mesh)

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Labeling, Segmentation, and Detection 的更多信息

产品

版本

R2023a

标签

Community Treasure Hunt

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

Start Hunting!

Translated by