Python计算点到直线距离、直线间交点夹角

 更新时间:2021年12月22日 16:03:27   作者:''竹先森゜  
这篇文章主要介绍了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计算点到直线距离内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家!

相关文章

  • 通过5个例子让你学会Pandas中的字符串过滤

    通过5个例子让你学会Pandas中的字符串过滤

    毋庸置疑Pandas是使用最广泛的Python库之一,它提供了许多功能和方法来执行有效的数据处理和数据分析,下面这篇文章主要给大家介绍了关于如何通过5个例子让你学会Pandas中字符串过滤的相关资料,需要的朋友可以参考下
    2022-08-08
  • Python创建一个自定义视频播放器的实现

    Python创建一个自定义视频播放器的实现

    本文主要介绍了Python创建一个自定义视频播放器的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2023-02-02
  • Python函数式编程之返回函数实例详解

    Python函数式编程之返回函数实例详解

    函数式编程的一个特点就是,允许把函数本身作为参数传入另一个函数,还允许返回一个函数,下面这篇文章主要给大家介绍了关于Python函数式编程之返回函数的相关资料,需要的朋友可以参考下
    2022-09-09
  • 利用Python3分析sitemap.xml并抓取导出全站链接详解

    利用Python3分析sitemap.xml并抓取导出全站链接详解

    因为最近更换了网址,所以需要在百度站长平台提交网址,不管是主动推送还是手动提交,前提都是要整理网站的链接,手动添加太麻烦,于是就想写个脚本直接抓取全站链接并导出,本文详细介绍的是实现的方法及过程,需要的朋友们一起来看看吧。
    2017-07-07
  • PyCharm更换pip源、模块安装以及PyCharm依赖包导入导出功能

    PyCharm更换pip源、模块安装以及PyCharm依赖包导入导出功能

    这篇文章主要给大家介绍了关于PyCharm更换pip源、模块安装以及PyCharm依赖包导入导出功能的相关资料,我们在使用pycharm的时候,pycharm中的虚拟环境依赖包需要导出成一个文件,需要的朋友可以参考下
    2023-11-11
  • python 爬取微信文章

    python 爬取微信文章

    本文给大家分享的是使用python通过搜狗入口,爬取微信文章的小程序,非常的简单实用,有需要的小伙伴可以参考下
    2016-01-01
  • Python retrying 重试机制详解

    Python retrying 重试机制详解

    这篇文章主要为大家介绍了Python retrying 重试机制,具有一定的参考价值,感兴趣的小伙伴们可以参考一下,希望能够给你带来帮助
    2021-11-11
  • python3处理word文档实例分析

    python3处理word文档实例分析

    在本篇文章里小编给大家分享了一篇关于python3处理word文档实例分析内容,有兴趣的朋友们可以参考下。
    2020-12-12
  • python3 cvs将数据读取为字典的方法

    python3 cvs将数据读取为字典的方法

    今天小编就为大家分享一篇python3 cvs将数据读取为字典的方法,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2018-12-12
  • Python中关于print和return的区别

    Python中关于print和return的区别

    这篇文章主要介绍了Python中关于print和return的区别,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教
    2024-02-02

最新评论