Python实现点到直线距离、直线交点与夹角计算
本文介绍使用Python进行几何计算的三种常见场景:点到直线距离、两条直线交点坐标及夹角。通过向量法和一般方程法推导公式,并结合NumPy库实现代码。内容涵盖数学原理、参数定义及完整Python实现示例,适用于需要处理几何数据的开发项目。
前言
在实际开发项目中,经常涉及点到直线距离的计算、两条直线交点坐标的求解以及两条直线夹角的计算等需求。
一、点到直线距离计算
由于在项目中获取点的坐标相对容易,采用向量法进行所有数学计算最为清晰直观。以下将使用向量法推导点到直线的距离公式。

import numpy as np
array_longi = np.array([x2-x1, y2-y1])
array_trans = np.array([x2-line_start_x, y2-line_start_y])
# 用向量计算点到直线的举例
array_temp = (float(array_trans .dot(array_longi)) / array_longi.dot(array_longi))
array_temp = array_longi.dot(array_temp)
distance = np.sqrt((array_trans - array_temp).dot(array_trans - array_temp )) # 距离
二、两条直线交点坐标计算
使用一般方程法进行计算:
直线的一般方程表示为 F(x) = ax + by + c = 0。假设直线上有两个已知点 (x0, y0) 和 (x1, y1),则可以推导出系数 a = y0 – y1,b = x1 – x0,c = x0y1 – x1y0。
因此,可以将两条直线分别表示为以下形式:
- F0(x) = a0*x + b0*y + c0 = 0
- F1(x) = a1*x + b1*y + c1 = 0
两条直线的交点必须同时满足以下方程组:
a0*x + b0*y + c0 = a1*x + b1*y + c1
由此可以推导出交点坐标公式:
- x = (b0*c1 – b1*c0) / D
- y = (a1*c0 – a0*c1) / D
- 其中 D = a0*b1 – a1*b0。当 D 为 0 时,表示两条直线平行,无交点。
上述公式实际上是联立方程组的叉积应用,具体表示为:
- F0(x) = a0*x + b0*y + c0 = 0
- F1(x) = a1*x + b1*y + c1 = 0
i j k
a0 b0 c0
a1 b1 c1
class Point:
x = 0
y = 0
def __init__(self, x=0, y=0):
self.x = x
self.y = y
class Line:
def __init__(self, p1, p2):
self.p1 = p1
self.p2 = p2
def GetLinePara(line):
line.a = line.p1.y - line.p2.y
line.b = line.p2.x - line.p1.x
line.c = line.p1.x * line.p2.y - ine.p2.x * line.p1.y
def GetCrossPoint(l1, l2):
GetLinePara(l1)
GetLinePara(l2)
d = l1.a * l2.c - l2.a * l1.b
p = Point()
p.x = (l1.b * l2.c - l2.b * l1.c) * 1.0 /d
p.y = (l1.c * l2.a - l2.c * l1.a) * 1.0 /d
p1 = Point(1, 1)
p2 = Point(3, 3)
line1 = Line(p1, p2)
p3 = Point(2, 3.1)
p = Point(3.1, 2)
line2 = Line(p3, p4)
Pc = GetCrossPoint(line1, line2)
print(Pc.x, Pc.y)
三、两条直线夹角计算
利用向量法求两条直线的夹角。其基本思路与求点到直线距离类似,主要利用余弦定理进行计算。

import numpy as np
def GetCrossAngle(l1, l2):
arr_0 = np.array([(l1.p2.x - l1.p1.x), (l1.p2.y - l1.p1.y)])
arr_1 = np.array([(l2.p2.x - l2.p1.x), (l2.p2.y - l2.p1.y)])
cos_value = (float(arr_0.dot(arr_1)) / (np.sqrt(arr_0.dot(arr_0)) * np.sqrt(arr_1.dot(arr_1))))
return np.arccos(cos_value) * (180 / np.pi)
angle = GetCrossAngle(line1, line2) # 计算得到的角度
参考来源:https://blog.csdn.net/panfengzjz/article/details/80377501
到此这篇关于Python计算点到直线距离、直线间交点夹角的文章就介绍到这了,更多相关python计算点到直线距离内容请搜索本站以前的文章或继续浏览下面的相关文章希望大家以后多多支持本站!
本站声明
本文内容仅供参考,不代表本站观点。转载或引用请注明出处。如发现内容问题,请通过联系我们页面反馈。













