# -*- 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：給定兩個(gè)線(xiàn)段的坐標(biāo)(也就是四個(gè)點(diǎn)的直角坐標(biāo)系坐標(biāo))，判斷這兩個(gè)線(xiàn)段是否相交

### 判斷線(xiàn)段是否相交可以利用向量的叉乘 ###
假定輸入為P1、P2、Q1、Q2四個(gè)點(diǎn)的坐標(biāo)，P1P2為一條線(xiàn)段，Q1Q2為另一條線(xiàn)段

兩條線(xiàn)段相交只有兩種情況
1. 其中一條線(xiàn)段的某一端點(diǎn)在另一條線(xiàn)段上；
2. 兩條線(xiàn)段形成X形。

首先判斷這四個(gè)點(diǎn)是否在另一條線(xiàn)段上，也就是說(shuō)，判斷P1是否在線(xiàn)段Q1Q2上，P2是否在線(xiàn)段Q1Q2上...
如果上述判斷為真，那么這兩條線(xiàn)段相交?！窘鉀Q了第一種情況】

如果沒(méi)有點(diǎn)在另一條線(xiàn)段上，那么進(jìn)行叉乘判斷。
先固定線(xiàn)段Q1Q2，然后以Q1為軸，計(jì)算Q1P1和Q1Q2、Q1P2和Q1Q2的叉乘是否異號(hào)；
然后固定線(xiàn)段P1P2，然后以P1為軸，計(jì)算P1Q1和P1P2、P1Q2和P1P2的叉乘是否異號(hào)。
當(dāng)上述的叉乘都異號(hào)的時(shí)候，兩條線(xiàn)段相交。
【解決了第二種情況】
"""

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):
    """
    判斷點(diǎn)在不在一條線(xiàn)段上，關(guān)鍵在于：
    1. 三點(diǎn)是否共線(xiàn)
    2. 點(diǎn)p是否在線(xiàn)段的延長(zhǎng)線(xiàn)上。

    只要滿(mǎn)足了三點(diǎn)共線(xiàn)，且點(diǎn)p不在延長(zhǎng)線(xiàn)上，那么點(diǎn)p就在線(xiàn)段上。

    判斷三點(diǎn)共線(xiàn)可以用向量的叉乘，三點(diǎn)共線(xiàn)即兩個(gè)向量平行，也就是叉乘結(jié)果為零向量(對(duì)應(yīng)到二維就是零)
    當(dāng)點(diǎn)p的橫縱坐標(biāo)都在線(xiàn)段端點(diǎn)之間的時(shí)候，點(diǎn)p不在延長(zhǎng)線(xiàn)上。
    :param p:
    :param seg:
    :return:
    """
    # 先確保點(diǎn)p不在延長(zhǎng)線(xiàn)上
    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):
        # 然后確保這三個(gè)點(diǎn)形成的向量?jī)蓛善叫校@里只要這三個(gè)向量中任意兩個(gè)平行，第三個(gè)一定也平行
        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)

    # 判斷是否存在端點(diǎn)位于另一條線(xiàn)段上，是的話(huà)則兩條線(xiàn)段相交
    if any([onSegment(p1, q1q2), onSegment(p2, q1q2),
            onSegment(q1, p1p2), onSegment(q2, p1p2)]):
        return True

    # 否則固定線(xiàn)段P1P2，判斷Q1和Q2是否在P1P2的兩側(cè)(計(jì)算叉乘)
    # 然后固定線(xiàn)段Q1Q2，判斷P1和P2是否在Q1Q2的兩側(cè)
    # 如果上面的判斷均為真，那么這兩條線(xiàn)段形成一個(gè)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()

分享名稱(chēng)：判斷線(xiàn)段是否相交
轉(zhuǎn)載來(lái)源：http://aaarwkj.com/article30/iipiso.html

成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián)，為您提供Google、網(wǎng)站收錄、移動(dòng)網(wǎng)站建設(shè)、域名注冊(cè)、電子商務(wù)、App設(shè)計(jì)

廣告

聲明：本網(wǎng)站發(fā)布的內(nèi)容（圖片、視頻和文字）以用戶(hù)投稿、用戶(hù)轉(zhuǎn)載內(nèi)容為主，如果涉及侵權(quán)請(qǐng)盡快告知，我們將會(huì)在第一時(shí)間刪除。文章觀點(diǎn)不代表本網(wǎng)站立場(chǎng)，如需處理請(qǐng)聯(lián)系客服。電話(huà)：028-86922220；郵箱：631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載，或轉(zhuǎn)載時(shí)需注明來(lái)源： 創(chuàng)新互聯(lián)