The problem of finding intersection point of triangle mesh

Question

Vivaldi 2019-11-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/492316-the-problem-of-finding-intersection-point-of-triangle-mesh

评论： Vivaldi 2019-11-25

face and vertex.zip

Now I have drawn a triangle mesh. (The function to generate image I use are as follows. )I want to ask for the intersection point, as shown in the figure(Red line with the mesh). But I don't know how to solve it.

PS:You can define the expression of a line at will,like "x=5".

If you want, I have some expression here.

trimesh(face', vertex(1,:), vertex(2,:), vertex(3,:),'LineWidth',1,'EdgeColor','k');
hold off;
vertex0=vertex';
V=vertex0(:,1:2);
x1=V(:,1);
y1=V(:,2);
z1=repmat(-20,1380,1);
tr = triangulation(face',x1,y1,z1); %x1.y1.z1 is a set of three coordinates of vertices
trimesh(tr);
axis equal;
colorbar('vert');
brighten(-0.1);axis on; % 'face' is the index list of triangle patch vertices;'vertex' is the index list of triangle vertex coordinates.

2 个评论
显示无隐藏无

darova 2019-11-21

What about in-built function such as polyxpoly?

Vivaldi 2019-11-24

Thanks for your answer. But I dont know how this function should I use. Cloud you write a function？Sorry to bother you.

请先登录，再进行评论。

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

Answer 1

darova 2019-11-24

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/492316-the-problem-of-finding-intersection-point-of-triangle-mesh#answer_403040

编辑：darova 2019-11-24

在 MATLAB Online 中打开

To check every edge:

for j = 1:size(face,2)
    i1 = 1;
    for i2 = [2 3 1]
            % checking edges [1 2], [2 3], [3 1]
        xx = [V(face(i2,j),1) V(face(i1,j),1)];
        yy = [V(face(i2,j),2) V(face(i1,j),2)];
            % intersection of edge and line x = [20 -150], y = [-40 -40]
        [xc,yc] = polyxpoly(xx,yy,[20 -150],[-40 -40]);    
        i1 = i1+1;
            % if intersection exists
        if ~isempty(xc)
            plot(xc,yc,'.r')
        end
    end
end