块状和直线相交

Question

如何使用numpy计算两条线段之间的交点？

在代码中，我有segment1 = ((x1,y1),(x2,y2))和segment2 = ((x1,y1),(x2,y2))。注意segment1不等于segment2。因此，在我的代码中，我还一直在计算斜率和y截距，如果可以避免，那将是很好的选择，但我不知道该怎么做。

我一直在用我在Python中编写的函数使用Cramer规则，但我想找到一种更快的方法。

Answer 1

直接从http://www.cs.mun.ca/~rod/2500/notes/numpy-arrays/numpy-
arrays.html
失窃

#
# line segment intersection using vectors
# see Computer Graphics by F.S. Hill
#
from numpy import *
def perp( a ) :
    b = empty_like(a)
    b[0] = -a[1]
    b[1] = a[0]
    return b

# line segment a given by endpoints a1, a2
# line segment b given by endpoints b1, b2
# return 
def seg_intersect(a1,a2, b1,b2) :
    da = a2-a1
    db = b2-b1
    dp = a1-b1
    dap = perp(da)
    denom = dot( dap, db)
    num = dot( dap, dp )
    return (num / denom.astype(float))*db + b1

p1 = array( [0.0, 0.0] )
p2 = array( [1.0, 0.0] )

p3 = array( [4.0, -5.0] )
p4 = array( [4.0, 2.0] )

print seg_intersect( p1,p2, p3,p4)

p1 = array( [2.0, 2.0] )
p2 = array( [4.0, 3.0] )

p3 = array( [6.0, 0.0] )
p4 = array( [6.0, 3.0] )

print seg_intersect( p1,p2, p3,p4)