我有一个相对简单的问题,我知道答案,但我似乎找不到使用Python和Numpy的正确实现。想法是,我有两条线,我需要找到虚拟交点(我使用了https://www.youtube.com/watch?v=kCyoaidiXAU&t=313s中的示例)。在
两条直线的形式都是r=r0+t*V,r0是位置向量(直线经过的点),t是变量,V是方向向量。方向向量V可以简单地通过直线的两个点找到,例如V=A-B
有了这一点,就可以把这条线表述为:
L1=r0(已知点)+t(未知变量)*V(方向向量)
现在,我可以很容易地手动找到t,但是我不知道如何告诉Python。。我试过了numpy.linalg.求解,但这给了我一个矩阵,而我需要一个单一的值。在
例如:# Line A and Line B that intersect somewhere
A = LineString([(4., 0.), (4., -3.)])
B = LineString([(6., 2.), (10., 2.)])
# direction vectors for line A and B
v1 = (A[0].x - A[1].x, A[0].y, A[1].y) # --> [0,3]
v2 = (B[0].x - B[1].x, B[0].y, B[1].y) # --> [-4,0]
L1 = A[1] + x * v1
L2 = B[1] + y * v2
现在,我将通过解L1得到x来手动解决这个问题:
^{pr2}$
但我不知道如何告诉numpy/python如何求解x和y
任何帮助或指导,在正确的方向将不胜感激。在