Python计算点到直线距离、直线间交点夹角
作者:''竹先森゜ 发布时间:2022-09-05 10:27:04
标签: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
来源:https://www.cnblogs.com/zhuminghui/p/15719831.html
0
投稿
猜你喜欢
- 我参与了IE7的开发过程,看到了在IE浏览器中形形色色使用MSXML的方法。显然有一些东西困扰着开发者:MSXML“混乱”的版本以及如何创建
- 我们做WEB程序开发过程中,经常会为了表格Table中的cellSpacing,cellpadding等样式的设置发愁,经常没有办法可想了,
- 基于很多原因,有不少朋友的空间不支持FSO,当大量的页面生成静态页面时,占用大量的空间,修改样式时,不方便,所以动态页面还是有动态页面的好处
- 一. XSL入门1.XSL---XML的样式表 HTML网页使用预先确定的标识(tags),这就是说所有的标记都有明确的含义,例如<p
- 在 Python 中,函数可以通过以下语法定义和使用:def function_name(parameter1, parameter2, .
- 用python的matplotlib画图时,往往需要加图例说明。如果不设置任何参数,默认是加到图像的内侧的最佳位置。import matpl
- 本文实例为大家分享了python使用Matplotlib画条形图的具体代码,供大家参考,具体内容如下数据中国的四个直辖市分别为北京市、上海市
- 简单介绍正则表达式并不是Python的一部分。正则表达式是用于处理字符串的强大工具,拥有自己独特的语法以及一个独立的处理引擎,效率上可能不如
- CONVERT将某种数据类型的表达式显式转换为另一种数据类型。由于某些需求经常用到取日期格式的不同.现以下可在SQL Server中将日期格
- 关于《回访确认》这件事情,在我blog和5Gme的回复里普遍出现这么几个疑问:1、下了订单却放弃支付的心理因素是什么?2、回访电话可以促进订
- 有时候要通过asp代码在数据库中创建表和列,下面的就是这些操作的函数。1.检测表是否存在tbName 检测的表的名称dbTp 数据库的类型1
- 本文将想大家简单介绍一下XML HttpRequst对象基础方法,希望通过本文能够使大家对其有一个初步的了解readyState一共有5个可
- 导语说到童年爱玩的电脑游戏,你会想到什么?最国民的莫过于金山打字通,接着是扫雷、红心大战,而红极一时的单机游戏当属《大富翁》。嘻嘻 打字游戏
- 首先让我们首先看一个例子。我们知道在*NIX下面,uname命令是查看系统的发行版。可以写这样一个Fabric脚本:from fabric.
- 1. 创建微信公众号首先,你需要注册一个微信公众号。访问微信公众平台,使用你的微信账号登录,并按照提示创建一个新的公众号。2. 开通微信公众
- 这最近在PJ的function库里看到的这个函数,感觉思路差了点,不过相对比较完美,只是闭合标签时的顺序问题,呵呵 修改一下数组arrTag
- 函数input()的工作原理函数input()让程序暂停运行,等待用户输入一些文本。获取用户输入后,Python将其存储在一个变量中,以方便
- 1 环境操作系统:Windows 10Python版本:3.9.0Google Chrome 87.0.4280.88
- 1、yield,将函数变为 generator (生成器)例如:斐波那契数列def fib(num): a, b, c = 1,
- 关于Ajax在使用中要使浏览器产生前进后退的方法,网上比较多的方法有两种:一是采用hash值的方式,这是我们在地图preview版中使用的方