pointLocation return only one ID if a point lies in more than one triangle

Question

Matteo Nicoletta 2022-3-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1668304-pointlocation-return-only-one-id-if-a-point-lies-in-more-than-one-triangle

编辑： VINAYAK LUHA 2023-9-27

Hi everyone,

this is my code in which I've two triangles and one point and I want to check in which triangle my point lies (I've to use pointLocation for more complex case like a mesh of STL file).

clear all
points=[0 0;1 0; 0 1; 1 1];
tri=[1 2 3; 1 2 4];
TR=triangulation(tri,points);
triplot(TR);
check=[0.6 0.2];
hold on
plot(check(1,1),check(1,2),'o');
ID = pointLocation(TR,check);

I'd like to be returned both the ID of the triangles in which the point lies. How can I do?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

VINAYAK LUHA 2023-9-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1668304-pointlocation-return-only-one-id-if-a-point-lies-in-more-than-one-triangle#answer_1319832

编辑：VINAYAK LUHA 2023-9-27

在 MATLAB Online 中打开

Hi Matteo,

I understand that you want to find the IDs of all triangles which enclose a point.The “pointLocation” function returns the triangle ID of only the first enclosing triangle as defined in the “Triangulation connectivity”.

To get triangleIDs of all the enclosing triangles, you can do a linear search over all the triangles and register the ones which enclose the given point.

The code snippet to achieve this can be seen below:

points = [0 0; 1 0; 0 1; 1 1];
tri = [1 2 3; 1 2 4];
check = [0.6 0.2];
plot(check(1,1), check(1,2), 'o');
hold on
enclosingTriangleIDs = [];
for i = 1:size(tri, 1)
    TR = triangulation(tri(i,:), points);
    ID = pointLocation(TR, check);
    if ~isnan(ID) 
        enclosingTriangleIDs = [enclosingTriangleIDs, i];
    end
    triplot(TR);
end
disp("IDs of triangles enclosing the point:");
disp(enclosingTriangleIDs);
hold off