判断线段是否相交
admin
2023-07-03 22:24:56
0
# -*- coding: utf-8 -*-
# @Time         : 2019-09-18 16:55
# @Author       : Jayce Wong
# @ProjectName  : job
# @FileName     : segment_cross.py
# @Blog         : https://blog.51cto.com/jayce1111
# @Github       : https://github.com/SysuJayce

"""
Q:给定两个线段的坐标(也就是四个点的直角坐标系坐标),判断这两个线段是否相交

### 判断线段是否相交可以利用向量的叉乘 ###
假定输入为P1、P2、Q1、Q2四个点的坐标,P1P2为一条线段,Q1Q2为另一条线段

两条线段相交只有两种情况
1. 其中一条线段的某一端点在另一条线段上;
2. 两条线段形成X形。

首先判断这四个点是否在另一条线段上,也就是说,判断P1是否在线段Q1Q2上,P2是否在线段Q1Q2上...
如果上述判断为真,那么这两条线段相交。【解决了第一种情况】

如果没有点在另一条线段上,那么进行叉乘判断。
先固定线段Q1Q2,然后以Q1为轴,计算Q1P1和Q1Q2、Q1P2和Q1Q2的叉乘是否异号;
然后固定线段P1P2,然后以P1为轴,计算P1Q1和P1P2、P1Q2和P1P2的叉乘是否异号。
当上述的叉乘都异号的时候,两条线段相交。
【解决了第二种情况】
"""

class Point:
    def __init__(self, x, y):
        self.x = x
        self.y = y

    def __sub__(self, other):
        return Point(self.x - other.x, self.y - other.y)

class Segment:
    def __init__(self, point1, point2):
        self.point1 = point1
        self.point2 = point2
        self.x = point1.x - point2.x
        self.y = point1.y - point2.y

def crossProduct(v1, v2):
    return v1.x * v2.y - v2.x * v1.y

def onSegment(p, seg):
    """
    判断点在不在一条线段上,关键在于:
    1. 三点是否共线
    2. 点p是否在线段的延长线上。

    只要满足了三点共线,且点p不在延长线上,那么点p就在线段上。

    判断三点共线可以用向量的叉乘,三点共线即两个向量平行,也就是叉乘结果为零向量(对应到二维就是零)
    当点p的横纵坐标都在线段端点之间的时候,点p不在延长线上。
    :param p:
    :param seg:
    :return:
    """
    # 先确保点p不在延长线上
    if min(seg.point1.x, seg.point2.x) <= p.x <= max(seg.point1.x, seg.point2.x)\
            and min(seg.point1.y, seg.point2.y) <= p.y <= max(seg.point1.y, seg.point2.y):
        # 然后确保这三个点形成的向量两两平行,这里只要这三个向量中任意两个平行,第三个一定也平行
        if crossProduct(p - seg.point1, p - seg.point2) == 0:
            return True
        else:
            return False
    else:
        return False

def isCross(p1, p2, q1, q2):
    p1p2 = Segment(p1, p2)
    q1q2 = Segment(q1, q2)

    p1q1 = Segment(p1, q1)
    p1q2 = Segment(p1, q2)

    q1p1 = Segment(q1, p1)
    q1p2 = Segment(q1, p2)

    # 判断是否存在端点位于另一条线段上,是的话则两条线段相交
    if any([onSegment(p1, q1q2), onSegment(p2, q1q2),
            onSegment(q1, p1p2), onSegment(q2, p1p2)]):
        return True

    # 否则固定线段P1P2,判断Q1和Q2是否在P1P2的两侧(计算叉乘)
    # 然后固定线段Q1Q2,判断P1和P2是否在Q1Q2的两侧
    # 如果上面的判断均为真,那么这两条线段形成一个X
    return (crossProduct(p1p2, p1q1) * crossProduct(p1p2, p1q2) < 0)\
           and (crossProduct(q1q2, q1p1) * crossProduct(q1q2, q1p2) < 0)

def main():
    p1 = Point(0, 0)
    p2 = Point(2, 2)
    q1 = Point(1, 1)
    q2 = Point(0, 2)
    if isCross(p1, p2, q1, q2):
        print('Yes')
    else:
        print('No')

if __name__ == '__main__':
    main()

相关内容

热门资讯

内蒙古包头:首批投运1000辆... 新华社5月31日消息,近日,我国自主研发的氢能两轮车在内蒙古自治区包头市面向公众投入运营,首批将在公...
缅甸边境发生爆炸,云南群众目击... 5月31日,缅甸掸邦北部南坎镇发生一起爆炸事故,已导致多人伤亡,多处民宅、房屋遭到严重损毁。经初步调...
大湾区打出智造新名片,高域首台... 近日,在广州黄埔区的智能制造产业园内,一架白色多旋翼飞行器缓缓驶出生产线,标志着高域(GOVY)这家...
VR眼镜与元宇宙:沉浸式体验会... VR眼镜与元宇宙:沉浸式体验会不会导致“魂不守舍”?虚拟现实对精神能量的消耗 杨明德/文 朋友们好,...
张凌赫粉丝挤爆玻璃门、活动临时... 5月31日,张凌赫原定在广西南宁万象城出席品牌活动,大批粉丝早早前往商场排队等候。现场人潮涌动,人群...
5月31日“蓝月亮”上线,还是... 据新华社:5月31日,农历四月十五,一轮满月将现身夜空。这轮满月有些特别,它是本月第二次满月,同时它...
中国的刀产量,究竟有多恐怖? 中年男人在短视频平台最爱看啥?小姐姐跳舞?那你可就想错了——锻刀大赛才是中年男人真正的“减速带”。在...
花费半生积蓄,农村自建房热背后 今年以来,在安徽安庆望江县做外墙漆生意的徐仙琴越来越忙碌了,她和员工们每天奔波在县城的各个村庄里,找...
黎以冲突再升级对中东地区影响几... 从当前的局势来看,黎以双方博弈已陷入谈判停滞、战火升级的恶性循环,地区冲突风险持续走高。根据以军发布...
易事达取得载带冷却定型装置专利... 国家知识产权局信息显示,浙江易事达电子材料有限公司取得一项名为“一种载带冷却定型装置”的专利,授权公...