Python計(jì)算點(diǎn)到直線距離、直線間交點(diǎn)夾角
前言
項(xiàng)目中會(huì)有點(diǎn)到直線距離計(jì)算、兩條直線交點(diǎn)坐標(biāo)計(jì)算、兩條直線夾角計(jì)算的需求。
一、點(diǎn)到直線距離計(jì)算
由于項(xiàng)目中得到點(diǎn)的坐標(biāo)最容易,因此采用向量法進(jìn)行所有的數(shù)學(xué)計(jì)算最清晰明了。點(diǎn)到直線距離就采用向量法推導(dǎo)。

import numpy as np array_longi = np.array([x2-x1, y2-y1]) array_trans = np.array([x2-line_start_x, y2-line_start_y]) # 用向量計(jì)算點(diǎn)到直線的舉例 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 )) # 距離
二、兩條直線交點(diǎn)坐標(biāo)計(jì)算
一般方程法:
直線的一般方程為 F(x) = ax + by + c = 0。假設(shè)直線的兩個(gè)點(diǎn)為(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
那么兩條直線的交點(diǎn)應(yīng)該滿足
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時(shí),表示兩直線平行)
二者實(shí)際上就是連立方程組的叉積應(yīng)用
- 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)
三、兩條直線夾角計(jì)算
利用向量法求兩條直線夾角。大致思路與求點(diǎn)到直線距離類似,也是利用余弦定理。

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) # 計(jì)算得到的角度
參考:https://blog.csdn.net/panfengzjz/article/details/80377501
到此這篇關(guān)于Python計(jì)算點(diǎn)到直線距離、直線間交點(diǎn)夾角的文章就介紹到這了,更多相關(guān)python計(jì)算點(diǎn)到直線距離內(nèi)容請(qǐng)搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家!
相關(guān)文章
Python基礎(chǔ)第三方模塊requests openpyxl
這篇文章主要為大家介紹了Python基礎(chǔ)第三方模塊requests openpyxl使用示例詳解,有需要的朋友可以借鑒參考下,希望能夠有所幫助,祝大家多多進(jìn)步,早日升職加薪2022-11-11
Python遍歷指定文件夾下的所有文件名的方法小結(jié)
當(dāng)需要遍歷指定文件夾下的所有文件名時(shí),Python提供了多種方法來(lái)實(shí)現(xiàn)這個(gè)任務(wù),本文將介紹如何使用Python來(lái)完成這一任務(wù),有需要的小伙伴可以參考下2024-01-01
Python 多進(jìn)程和數(shù)據(jù)傳遞的理解
這篇文章主要介紹了Python 多進(jìn)程和數(shù)據(jù)傳遞的理解的相關(guān)資料,希望通過(guò)本文能幫助到大家,需要的朋友可以參考下2017-10-10
利于python腳本編寫可視化nmap和masscan的方法
這篇文章主要介紹了利于python腳本編寫可視化nmap和masscan的方法,本文給大家介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值,需要的朋友可以參考下2020-12-12
python django集成cas驗(yàn)證系統(tǒng)
cas是什么東西就不多說(shuō)了,簡(jiǎn)而言之就是單點(diǎn)登陸系統(tǒng),一處登陸,全網(wǎng)有權(quán)限的系統(tǒng)均可以訪問(wèn)2014-07-07

