一种解决方案是使用在
this tutorial for finding the intersection point of two lines in 2-D中导出的等式(更新:这是因特网存档链接,因为该站点不再存在).您可以先创建两个矩阵:一个用于保存线端点的x坐标,另一个用于保持y坐标.
x = [0 0; 6 6]; %# Starting points in first row, ending points in second row
y = [0 6; 6 0];
然后可以将来自上述源的等式编码如下:
dx = diff(x); %# Take the differences down each column
dy = diff(y);
den = dx(1)*dy(2)-dy(1)*dx(2); %# Precompute the denominator
ua = (dx(2)*(y(1)-y(3))-dy(2)*(x(1)-x(3)))/den;
ub = (dx(1)*(y(1)-y(3))-dy(1)*(x(1)-x(3)))/den;
现在,您可以计算两条线的交点:
xi = x(1)+ua*dx(1);
yi = y(1)+ua*dy(1);
对于问题中的示例,上面的代码给出了xi = 3和yi = 3,如预期的那样.如果要检查交点是否位于线的端点之间(即它们是有限线段),则只需检查值ua和ub是否介于0和1之间:
isInSegment = all(([ua ub] >= 0) & ([ua ub] <= 1));
我上面链接的教程还有几点:
>如果分母den为0,那么两条线是平行的.>如果ua和ub方程的分母和分子为0,则两条线重合.