python 两直线距离 垂线中点法

最新推荐文章于 2024-07-04 16:04:39 发布
最新推荐文章于 2024-07-04 16:04:39 发布
阅读量1.9k 收藏 9
点赞数 1
文章标签： python vim sublime text pycharm
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/tough__man/article/details/119460074
版权

需求：求两条直线的距离

方法：通过两条直线上的两个点分别向另一条直线做垂线，两个垂点的中心点连接作为距离。

如下图：黑线表示待求的两条直线间的距离，红线表示右侧直线向左侧直线的垂线，蓝线表示左侧直线向右侧直线的垂线，紫线表示两条直线上垂点中点的连线（代表距离）。

实现代码：

import math
import cv2
def intersectionLines(a1, a2, b1, b2):
    return (
          ((a1[0] * a2[1]  -  a1[1] * a2[0]) * (b1[0] - b2[0]) -
           (b1[0] * b2[1]  -  b1[1] * b2[0]) * (a1[0] - a2[0])) /
          ((a1[0] - a2[0]) * (b1[1] - b2[1]) -
           (a1[1] - a2[1]) * (b1[0] - b2[0])),
          ((a1[0] * a2[1]  -  a1[1] * a2[0]) * (b1[1] - b2[1]) -
           (b1[0] * b2[1]  -  b1[1] * b2[0]) * (a1[1] - a2[1])) /
          ((a1[0] - a2[0]) * (b1[1] - b2[1]) -
           (a1[1] - a2[1]) * (b1[0] - b2[0]))
          )
def p_to_l(point,l_p1,l_p2):
    if l_p1[0]-l_p2[0]==0:
        crop=(l_p1[0],point[1])
        # print("crop",crop)
    elif l_p1[1]-l_p2[1]==0:
        crop=(point[0],l_p1[1])
    else:
        k_l=(l_p1[1]-l_p2[1])/(l_p1[0]-l_p2[0])
        # print(k_l)
        k=-1/k_l
        # print(k)
        b=point[1]-k*point[0]
        # print(b)
        x=l_p1[0]
        # print(x)
        y=k*x+b
        point2=(x,y)
        # print(point)
        # print(point2)
        # print(l_p1,l_p2)
        crop=intersectionLines(point, point2, l_p1, l_p2)
    return crop
def p_to_distance(p1,p2):
    return math.sqrt((p1[0]-p2[0])**2+(p1[1]-p2[1])**2)


def to_overlap_pts(a1, a2, b1,b2):
    if a1[1]>a2[1]:
        l1_hight=a1
        l1_low=a2
    else:
        l1_hight=a2
        l1_low=a1
    if b1[1]>b2[1]:
        l2_hight=b1
        l2_low=b2
    else:
        l2_hight=b2
        l2_low=b1
    

    l1_p1, l1_p2, l2_p1, l2_p2 = (0., 0.), (0., 0.), (0., 0.), (0., 0.)
    point1_state=0
    point2_state=0
    #直线没有重合部分l1
    l1_l=p_to_l(l1_low,l2_low,l2_hight)

    l1_h=p_to_l(l1_hight,l2_low,l2_hight)
    # print("l1",l1_h)
    # print("l2_low",l2_low)
    # print("l2_hight",l2_hight)
    list=[l2_low[0],l2_hight[0]]
    l=min(list)
    r=max(list)
    if l<=l1_h[0]<=r:
        point1_state=1
    if l<=l1_l[0]<=r:
        point2_state=1
    print("state",point1_state,point2_state)
    if point1_state==0 and point2_state==0:
        print("框选错误，请重新框选")
    else:
        l1_b=0
        l1_t=0
        l2_b=0
        l2_t=0
        #左边直线
        if l2_low[0]<l2_hight[0]:
            l=l2_low[0]
            r=l2_hight[0]
        else:
            r = l2_low[0]
            l = l2_hight[0]

        cro_p1 = p_to_l(l1_low, l2_low, l2_hight)
        if l<=cro_p1[0]<=r:
            p1=l1_low
            p2=cro_p1
            cro_p1 = (cro_p1[0], cro_p1[1])
            l1_b=1
        else:
            pass
        cro_p2 = p_to_l(l1_hight, l2_low, l2_hight)
        if l<=cro_p2[0]<=r:

            p1=l1_hight
            p2=cro_p2
            cro_p2 = (cro_p2[0], cro_p2[1])
            l1_t=1
        else:
            pass

        #右边直线
        if l1_low[0]<l1_hight[0]:
            l=l1_low[0]
            r=l1_hight[0]
        else:
            r= l1_low[0]
            l = l1_hight[0]
        cro_p3 = p_to_l(l2_low, l1_low, l1_hight)
        if l<=cro_p3[0]<=r:
            p1=l2_low
            p2=cro_p3
            cro_p3 = (cro_p3[0], cro_p3[1])
            l2_b=1
        else:
            pass
        cro_p4 = p_to_l(l2_hight, l1_low, l1_hight)
        if l<=cro_p4[0]<=r:
            p1=l2_hight
            p2=cro_p4
            cro_p4 = (cro_p4[0], cro_p4[1])
            l2_t=1
        else:
            pass
        #右直线投影到左直线
        if l1_b+l1_t==0 and l2_b+l2_t==2:
            l1_p1=cro_p3
            l1_p2=cro_p4
            l2_p1=l2_low
            l2_p2=l2_hight
       #左直线投影到右直线
        elif l1_b+l1_t==2 and l2_b+l2_t==0:
            l1_p1=l1_low
            l1_p2=l1_hight
            l2_p1=cro_p1
            l2_p2=cro_p2
        #右直线投影到作直线，左直线一个点投影到右直线
        elif l1_b + l1_t == 1 and l2_b + l2_t == 2:
            #左直线上点
            if l1_t == 1:
                l2_p1 = l2_low
                l2_p2 = cro_p2
            #左直线下
            else:
                l2_p1 = cro_p1
                l2_p2 = l2_hight
            l1_p1 = cro_p3
            l1_p2 = cro_p4
        #左直线投影到右直线，右直线一个点投影到左直线
        elif l1_b + l1_t == 2 and l2_b + l2_t == 1:
            #右直线上点
            if l1_t == 1:
                l1_p1 = l1_low
                l1_p2 = cro_p4
            #右直线下点
            else:
                l1_p1 = cro_p3
                l1_p2 = l1_hight
            l2_p1 = cro_p1
            l2_p2 = cro_p2
        elif l1_b == 1 and l2_t == 1:
            l1_p1 = l1_low
            l1_p2 = cro_p4
            l2_p1 = cro_p1
            l2_p2 = l2_hight
        elif l1_t == 1 and l2_b == 1:
            l1_p1 = cro_p3
            l1_p2 = l1_hight
            l2_p1 = l2_low
            l2_p2 = cro_p2

    return l1_p1, l1_p2, l2_p1, l2_p2

def coincidence_part(a1, a2, b1,b2):
    # img=cv2.imread("1030_1_undistorted.png")

    if a1[1]>a2[1]:
        l1_hight=a1
        l1_low=a2
    else:
        l1_hight=a2
        l1_low=a1
    if b1[1]>b2[1]:
        l2_hight=b1
        l2_low=b2
    else:
        l2_hight=b2
        l2_low=b1
    #直线没有重合部分l1

    if l1_hight[1]<=l2_low[1] or l1_low[1]>=l2_hight[1]:
        print("框选错误，请重新框选")
    else:
        l1_b=0
        l1_t=0
        l2_b=0
        l2_t=0
        #左边直线
        if l2_low[0]<l2_hight[0]:
            l=l2_low[0]
            r=l2_hight[0]
        else:
            r = l2_low[0]
            l = l2_hight[0]

        cro_p1 = p_to_l(l1_low, l2_low, l2_hight)
        if l<=cro_p1[0]<=r:
            p1=l1_low
            p2=cro_p1
            cro_p1 = (cro_p1[0], cro_p1[1])
            l1_b=1
        else:
            pass
        cro_p2 = p_to_l(l1_hight, l2_low, l2_hight)
        if l<=cro_p2[0]<=r:

            p1=l1_hight
            p2=cro_p2
            cro_p2 = (cro_p2[0], cro_p2[1])
            l1_t=1
        else:
            pass

        #右边直线
        if l1_low[0]<l1_hight[0]:
            l=l1_low[0]
            r=l1_hight[0]
        else:
            r= l1_low[0]
            l = l1_hight[0]
        cro_p3 = p_to_l(l2_low, l1_low, l1_hight)
        if l<=cro_p3[0]<=r:
            p1=l2_low
            p2=cro_p3
            cro_p3 = (cro_p3[0], cro_p3[1])
            l2_b=1
        else:
            pass
        cro_p4 = p_to_l(l2_hight, l1_low, l1_hight)
        if l<=cro_p4[0]<=r:
            p1=l2_hight
            p2=cro_p4
            cro_p4 = (cro_p4[0], cro_p4[1])
            l2_t=1
        else:
            pass
        #右直线投影到左直线
        if l1_b+l1_t==0 and l2_b+l2_t==2:
            l1_p1=cro_p3
            l1_p2=cro_p4
            l2_p1=l2_low
            l2_p2=l2_hight
       #左直线投影到右直线
        elif l1_b+l1_t==2 and l2_b+l2_t==0:
            l1_p1=l1_low
            l1_p2=l1_hight
            l2_p1=cro_p1
            l2_p2=cro_p2
        #右直线投影到作直线，左直线一个点投影到右直线
        elif l1_b + l1_t == 1 and l2_b + l2_t == 2:
            #左直线上点
            if l1_t == 1:
                l2_p1 = l2_low
                l2_p2 = cro_p2
            #左直线下
            else:
                l2_p1 = cro_p1
                l2_p2 = l2_hight
            l1_p1 = cro_p3
            l1_p2 = cro_p4
        #左直线投影到右直线，右直线一个点投影到左直线
        elif l1_b + l1_t == 2 and l2_b + l2_t == 1:
            #右直线上点
            if l1_t == 1:
                l1_p1 = l1_low
                l1_p2 = cro_p4
            #右直线下点
            else:
                l1_p1 = cro_p3
                l1_p2 = l1_hight
            l2_p1 = cro_p1
            l2_p2 = cro_p2
        elif l1_b == 1 and l2_t == 1:
            l1_p1 = l1_low
            l1_p2 = cro_p4
            l2_p1 = cro_p1
            l2_p2 = l2_hight
        elif l1_t == 1 and l2_b == 1:
            l1_p1 = cro_p3
            l1_p2 = l1_hight
            l2_p1 = l2_low
            l2_p2 = cro_p2


        # print(l1_b, l1_t, l2_b, l2_t)
        # print(l1_p1,l1_p2,l2_p1,l2_p2)
        l1_c = ((l1_p1[0] + l1_p2[0])/ 2, (l1_p1[1] + l1_p2[1]) / 2)
        l1_to_l2=p_to_l(l1_c,l2_p1,l2_p2)
        l1_to_l2=(l1_to_l2[0],l1_to_l2[1])

        l2_c = ((l2_p1[0] + l2_p2[0]) / 2, (l2_p1[1] + l2_p2[1]) / 2)
        # print(l2_c,l1_p1,l1_p2)

        l2_to_l1=p_to_l(l2_c,l1_p1,l1_p2)
        l2_to_l1 = (int(l2_to_l1[0]), int(l2_to_l1[1]))

        l1=p_to_distance(l1_c,l1_to_l2)
        l2=p_to_distance(l2_c,l2_to_l1)
        # print(l1,l2)
        distance=(l1+l2)/2

        # print(distance)
        return distance
# coincidence_part((200,100),(200,600),(800,100),(800,621))
# coincidence_part((100,100),(600,100),(200,400),(500,400))
tough__man
关注
  • 1
    点赞
  • 踩
  • 9
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
已知俩点求俩点之间的直线，俩点间的中垂线，俩条直线的交点
dz4543的博客
06-20 7727
参考博客 https://www.cnblogs.com/DHUtoBUAA/p/8057056.html 用的时候才发现这些基础知识还真忘了不少。 俩点之间的直线 直线的一般表达式是：Ax+By+C=0。如果一直俩点（x1,y2）、(x2,y2).那么对于过着俩点的一般式是这样： A=y2-y1 B=x1-x2 C=X2×Y1-X1×Y2 所以用Python代码来表示是...
利用python计算两个空间直线的公垂线距离并计算和两条直线之间的交点
hunter206的博客
03-31 3264
利用python计算两个空间直线的公垂线距离并计算和两条直线之间的交点
Python求两点之间的直线距离（两种方法）
热门推荐
Ethan的博客
08-13 8万+
方法一： #导入math包 import math #定义点的函数 class Point: def __init__(self,x=0,y=0): self.x=x self.y=y def getx(self): return self.x def gety(self): return self.y ...
利用 python 计算点到直线距离、直线间交点夹角
panfengzjz的博客
05-19 1万+
近期在用 OpenCV 时需要涉及到一些点与直线关系。这都是很早以前的中学知识了，无奈只好重新捡起来。本次项目主要需要用到以下三种计算1. 点到直线距离2. 两条直线交点3. 两条直线夹角1. 点到直线距离由于项目中得到点的坐标最容易，因此采用向量法进行所有的数学计算最清晰明了。点到直线距离就采用向量法推导。 推导过程： 代码实现：array_longi = np.array([x2-x1, y...
python通过两点经纬度计算两点间直线距离
09-29
python通过两点经纬度计算两点间直线距离
Python求两点之间的直线距离(2种实现方法)
09-19
今天小编就为大家分享一篇Python求两点之间的直线距离(2种实现方法)，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
Python求平面内点到直线距离的实现
09-18
首先，我们需要理解点到直线距离的基本数学原理。点到直线的距离可以使用点到直线的投影来计算。假设直线的方程是 Ax + By + C = 0，其中A、B是直线的斜率参数，C是截距。给定点P(x, y)，点到直线的距离d可以通过...
中点法画直线的程序
02-06
中点法是一种在计算机图形学中用于绘制直线的算法，特别是在二维平面上。这个程序的扩展版能够处理各种斜率的直线，使得它更加通用。本文将深入探讨中点法的原理、实现细节以及如何扩展以支持不同斜率的直线。 ### ...
最小最大距离法（python实现）
04-30
最小最大距离法（python实现），用于实现样本聚类，含数据示例
垂距限值法_抽稀_垂距限值法_python抽稀_
09-30
垂距限值抽稀算法python实现。空间数据抽稀梳理。
python绘制直线的垂线_p countour中的Python垂直线
weixin_39797780的博客
12-08 1217
在二维绘图中，我通常使用l, = pylab.plot([10, 10], [-1000,1000], color="g", lw=0.5)在x位置10在垂直轴上绘制范围从-1000到1000的垂直线。在我想在等高线图上也这样做如何在等高线图中绘制垂直线？在我特别使用下面的代码# Blue-White-Red colorbar for plots with negative and positiv...
python实现两线段最短距离并实现线段分簇
DETACH-MENT的博客
10-17 2080
思路：线段取n个点，计算两线段线段不同点之间的距离。 import matplotlib.pylab as plt import numpy as np from numpy import array, int32 from typing import List from itertools import combinations def show_lines_data(lines: List, label: str): plt.xlim(200, 650) plt.ylim(160,
cvtcolor python opencv_二值分析 | OpenCV + skimage如何提取中心线
weixin_39572764的博客
11-24 538
点击上方蓝字关注我们微信公众号：OpenCV学堂关注获取更多计算机视觉与深度学习知识问题前几天有个人问了我一个问题，问题是这样的，他有如下的一张二值图像：怎么得到白色Blob中心线，他希望的效果如下：显然OpenCV中常见的轮廓分析无法获得上面的中心红色线段，本质上这个问题是如何提取二值对象的骨架，提取骨架的方法在OpenCV的扩展模块中，另外skimage包也支持图像的骨架提取。这里就...
python：实现两条线之间的最短距离(附完整源码)
希望我的博客，能帮上你解决学习中工作中所遇到的问题
03-17 710
python：实现两条线之间的最短距离(附完整源码)
python公式代码记录（四）-求点到直线的垂直距离
LuckyMon的博客
10-29 1297
已知一点p（a，b），已知p1（a1，b1），p2（a2，b2）在同一直线，求p到p1，p2的垂直距离。 area = 1/2 * h * c3 h=2*area/c3 代码如下 // An highlighted block # 计算点到直线的距离 def point_to_string(a,b,a1,b1,a2,b2): c1=base_distance(a,b,a1,b1) c2=base_distance(a,b,a2,b2) c3=base_distance(a1,b1,a
python-opencv直线检测+直线拟合+计算直线间距离
weixin_41921564的博客
06-13 1万+
python-opencv直线检测+直线拟合+计算偏移量 计算边界直线与x坐标轴的距离=首先canny边缘检测——>直线检测——>直线拟合——>计算与x坐标轴的距离。 直接上代码： import cv2 #opencv读取的格式是BGR import numpy as np vc = cv2.VideoCapture('D:/Desktop/duizhong/fangzhen.mp4') if vc.isOpened(): #判断能否读取 open, frame =
利用python计算各类距离
光尘的博客
05-27 5151
Jaccard相似度 1、当两个集合元素个数相同，则直接调包 from numpy import * import scipy.spatial.distance as dist # 导入scipy距离公式 matV = mat([[1,1,0,1,0,1,0,0,1],[0,1,1,0,0,0,1,1,1]]) print ("dist.jaccard:", dist.pdist(matV,'jaccard')) 2、当集合元素个数不同 def correlation(set_a,set_b
Python-类03（求两点之间直线距离）
Python学习笔记
04-01 2万+
#导入math包 import math #定义点的函数 class Point: def __init__(self,x=0,y=0): self.x=x self.y=y def getx(self): return self.x def gety(self): return self.y #定义直线函数
Java基础:爬虫
最新发布
weixin_65752158的博客
07-04 842
Pattern:表示正则表达式Matcher:文本匹配器，作用按照正则表达式的规则去读取字符串，从头开始读取。在大串中去找符合匹配规则的子串。通过Pattern p = Pattern.compile("正则表达式");获得 通过Matcher m = p.matcher(str);获得 (m要在str中找符合p规则的小串)其中, m为Matcher对象, p为正则表达式规则, str为要验证的字符串. boolean b = m.find(); 表示拿着文本匹配器从头开始读取，寻找是否有
python中如何求两条直线端点距离
04-19
在Python中可以使用以下公式来计算两条直线的端点距离： d = math.sqrt((x2 - x1)**2 + (y2 - y1)**2) 其中，(x1, y1)和(x2, y2)分别是直线的两个端点的坐标，math.sqrt()是Python中的平方根函数。 需要注意的是，要先导入math模块才能使用sqrt()函数。具体实现代码如下： import math def line_distance(x1, y1, x2, y2): d = math.sqrt((x2 - x1)**2 + (y2 - y1)**2) return d x1 = 1 y1 = 2 x2 = 3 y2 = 4 distance = line_distance(x1, y1, x2, y2) print("两条直线的端点距离为：", distance)

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值