python 直线求交_Python中与numpy的直线相交

我有一个相对简单的问题,我知道答案,但我似乎找不到使用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

任何帮助或指导,在正确的方向将不胜感激。在

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
假设两条直线分别为 L1 和 L2,它们分别由两个点 A1、A2 和 B1、B2 定义。为了求出它们的交点,我们需要先计算出这两条直线的方向向量和法向量。 L1 的方向向量为 (A2 - A1),即从 A1 到 A2 的向量。同样地,L2 的方向向量为 (B2 - B1)。将这两个向量叉乘,得到它们的法向量 N: N = (A2 - A1) × (B2 - B1) 接下来,我们需要计算出这两条直线在平面上的投影。我们可以选择其一条直线作为基准线,将另一条直线投影到这条直线所在的平面上。这个平面的法向量就是 N。假设我们选择 L1 作为基准线,投影点为 P1: P1 = A1 + dot(B1 - A1, N) / dot(A2 - A1, N) * (A2 - A1) 其,dot 表示向量的点积。同样地,我们可以计算出 L2 在 L1 所在平面上的投影点 P2: P2 = B1 + dot(A1 - B1, N) / dot(B2 - B1, N) * (B2 - B1) 最后,我们只需要计算出 P1 和 P2 的点,即为这两条直线交点: P = (P1 + P2) / 2 完整的 Python 代码如下: ```python import numpy as np def line_intersection_3d(A1, A2, B1, B2): # 计算方向向量和法向量 dir1 = A2 - A1 dir2 = B2 - B1 N = np.cross(dir1, dir2) # 计算 L2 在 L1 所在平面上的投影点 P1 P1 = A1 + np.dot(B1 - A1, N) / np.dot(dir1, N) * dir1 # 计算 L1 在 L2 所在平面上的投影点 P2 P2 = B1 + np.dot(A1 - B1, N) / np.dot(dir2, N) * dir2 # 计算交点 P P = (P1 + P2) / 2 return P ``` 其,A1、A2、B1 和 B2 分别是三维空间的四个点,可以用 numpy 的数组表示。函数返回的是交点 P,也是一个 numpy 数组。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值