多邊形相交判斷算法

一、規(guī)則

判斷兩多邊形是否相交顷级,排除邊邊重合凫乖、以及斷點(diǎn)與邊重合的情況

二、示意圖

示意圖.png

三弓颈、算法代碼(js)

/**
     * 判斷兩個多邊形是否相交(邊邊重合帽芽,點(diǎn)邊重合除外)
     * 核心算法:1-快速排除算法 2-矢量叉乘跨立算法 3-射線算法
     *
     * @param polygon1
     * @param polygon2
     */
    function intersectsPolygonAndPolygonNew(polygon1, polygon2) {
        // 兩多邊形是否相交標(biāo)記
        var _ifIntersect;
        // 經(jīng)緯度轉(zhuǎn)換成平面坐標(biāo)
        var polygon1Points = pathConvert(polygon1);
        var polygon2Points = pathConvert(polygon2);
        // 快速排除法判斷是否相交
        _ifIntersect = fastExclude(polygon1Points, polygon2Points);
        if (_ifIntersect == false) {
            return _ifIntersect;
        }
        // 獲取線段集合
        var _polygon1Segs = getSegs(polygon1Points);
        var _polygon2Segs = getSegs(polygon2Points);

        // 跨立法判斷相交
        _ifIntersect = judgeIntersect(_polygon1Segs, _polygon2Segs);
        if(_ifIntersect) {
            return _ifIntersect;
        }else {
            // 判斷是否包含
            var _containPointNum = 0;
            for (var i = 0; i < polygon1Points.length - 1; i++) {
                var _ifContain = judgePointInPolygon(polygon1Points[i], _polygon2Segs);
                if (_ifContain) {
                    _containPointNum++;
                }
            }
            if (_containPointNum == polygon1Points.length - 1) {
                return _ifIntersect = true;
            } else {
                return _ifIntersect = false;
            }
        }

        // 獲取多邊形線段集合
        function getSegs(points) {
            var _segs = [];
            for (var i = 0; i < points.length - 1; i++) {
                var _element1 = points[i];
                var _element2 = points[i + 1];
                _segs.push({x1: _element1.x, y1: _element1.y, x2: _element2.x, y2: _element2.y});
            }
            return _segs;
        }
    }

    /**
     * 快速排除算法:初步判斷多邊形是否相交
     *
     * @param polygon1Points
     * @param polygon2Points
     */
    function fastExclude(polygon1Points, polygon2Points) {
        // 多邊形1最大X坐標(biāo)
        var _polygon1MaxX;
        // 多邊形1最大Y坐標(biāo)
        var _polygon1MaxY;
        // 多邊形1最小X坐標(biāo)
        var _polygon1MinX;
        // 多邊形1最小Y坐標(biāo)
        var _polygon1MinY;

        // 多邊形2最大X坐標(biāo)
        var _polygon2MaxX;
        // 多邊形2最大Y坐標(biāo)
        var _polygon2MaxY;
        // 多邊形2最小X坐標(biāo)
        var _polygon2MinX;
        // 多邊形2最小Y坐標(biāo)
        var _polygon2MinY;

        for (var i = 0; i < polygon1Points.length; i++) {
            var _polygon1Point = polygon1Points[i];
            _polygon1MaxX = _polygon1MaxX > _polygon1Point.x ? _polygon1MaxX : _polygon1Point.x;
            _polygon1MinX = _polygon1MinX < _polygon1Point.x ? _polygon1MinX : _polygon1Point.x;
            _polygon1MaxY = _polygon1MaxY > _polygon1Point.y ? _polygon1MaxY : _polygon1Point.y;
            _polygon1MinY = _polygon1MinY < _polygon1Point.y ? _polygon1MinY : _polygon1Point.y;
        }

        for (var i = 0; i < polygon2Points.length; i++) {
            var _polygon2Point = polygon2Points[i];
            _polygon2MaxX = _polygon2MaxX > _polygon2Point.x ? _polygon2MaxX : _polygon2Point.x;
            _polygon2MinX = _polygon2MinX < _polygon2Point.x ? _polygon2MinX : _polygon2Point.x;
            _polygon2MaxY = _polygon2MaxY > _polygon2Point.y ? _polygon2MaxY : _polygon2Point.y;
            _polygon2MinY = _polygon2MinY < _polygon2Point.y ? _polygon2MinY : _polygon2Point.y;
        }

        if (_polygon1MaxX <= _polygon2MinX || _polygon2MaxX <= _polygon1MinX || _polygon1MinY >= _polygon2MaxY || _polygon2MinY >= _polygon1MaxY) {
            return false;
        } else {
            return true;
        }
    }

    /**
     * 跨立算法:
     * 利用矢量叉乘的物理意義判斷兩線段是否相交
     *
     * 有向量 a,b翔冀,c
     * 若 (a*b)*(b*c) = 0 則兩線段重合
     * 若 (a*b)*(b*c) < 0 則兩線段相交
     * 若 (a*b)*(b*c) > 0 則兩線段不相交
     *
     * @param polygon1Segs
     * @param polygon2Segs
     */
    function judgeIntersect(polygon1Segs, polygon2Segs) {
        var _ifIntersect;
        math.config({
            number: 'BigNumber'
        });
        for (var i = 0; i < polygon1Segs.length; i++) {
            var polygon1Seg = polygon1Segs[i];
            for (var j = 0; j < polygon2Segs.length; j++) {
                var polygon2Seg = polygon2Segs[j];
                // 向量1
                var _vector1 = new vector(math.parser().eval(polygon2Seg.x2 + "-" + polygon2Seg.x1), math.parser().eval(polygon2Seg.y2 + "-" + polygon2Seg.y1));
                // 向量2
                var _vector2 = new vector(math.parser().eval(polygon1Seg.x1 + "-" + polygon2Seg.x1), math.parser().eval(polygon1Seg.y1 + "-" + polygon2Seg.y1));
                // 向量3
                var _vector3 = new vector(math.parser().eval(polygon1Seg.x2 + "-" + polygon2Seg.x1), math.parser().eval(polygon1Seg.y2 + "-" + polygon2Seg.y1));
                // 向量_vector2與_vector1叉乘的解
                var _calculationVal1 = math.parser().eval(_vector2.x + "*" + _vector1.y + "-" + _vector1.x + "*" + _vector2.y);
                // 向量_vector3與_vector1叉乘的解
                var _calculationVal2 = math.parser().eval(_vector3.x + "*" + _vector1.y + "-" + _vector1.x + "*" + _vector3.y);

                // 向量4
                var _vector4 = new vector(math.parser().eval(polygon1Seg.x2 + "-" + polygon1Seg.x1), math.parser().eval(polygon1Seg.y2 + "-" + polygon1Seg.y1));
                // 向量5
                var _vector5 = new vector(math.parser().eval(polygon2Seg.x1 + "-" + polygon1Seg.x1), math.parser().eval(polygon2Seg.y1 + "-" + polygon1Seg.y1));
                // 向量6
                var _vector6 = new vector(math.parser().eval(polygon2Seg.x2 + "-" + polygon1Seg.x1), math.parser().eval(polygon2Seg.y2 + "-" + polygon1Seg.y1));
                // 向量_vector5與_vector4叉乘的解
                var _calculationVal3 = math.parser().eval(_vector5.x + "*" + _vector4.y + "-" + _vector4.x + "*" + _vector5.y);
                // 向量_vector6與_vector4叉乘的解
                var _calculationVal4 = math.parser().eval(_vector6.x + "*" + _vector4.y + "-" + _vector4.x + "*" + _vector6.y);

                if (_calculationVal1 * _calculationVal2 < 0 && _calculationVal3 * _calculationVal4 < 0) {
                    return _ifIntersect = true;
                }
            }
        }
        return _ifIntersect = false;

        // 向量對象导街,因?yàn)槭峭黄矫婢€段,所以省略Z坐標(biāo)
        function vector(x, y) {
            this.x = x;
            this.y = y;
        }
    }

    /**
     * 射線法判斷點(diǎn)是否在多邊形內(nèi)
     * 點(diǎn)射線(向右水平)與多邊形相交點(diǎn)的個數(shù)為奇數(shù)則認(rèn)為該點(diǎn)在多邊形內(nèi)
     * 點(diǎn)射線(向右水平)與多邊形相交點(diǎn)的個數(shù)為偶數(shù)則認(rèn)為該點(diǎn)不在多邊形內(nèi)
     *
     * @param point
     * @param polygonSegs
     */
    function judgePointInPolygon(point, polygonSegs) {
        debugger
        // 是否在多邊形內(nèi)
        var _ifContain;
        // 點(diǎn)坐標(biāo)
        var _X = point.x;
        var _Y = point.y;
        // 交點(diǎn)個數(shù)
        var _intersecNum = 0;
        // 判斷射線與多邊形的交點(diǎn)個數(shù)
        for (var i = 0; i < polygonSegs.length; i++) {
            var _seg = polygonSegs[i];
            var _maxY = _seg.y1 < _seg.y2 ? _seg.y2 : _seg.y1;
            var _minY = _seg.y1 < _seg.y2 ? _seg.y1 : _seg.y2;
            if (_Y > _minY && _Y < _maxY) {
                math.config({
                    number: 'BigNumber'
                });
                // 計(jì)算交點(diǎn)X坐標(biāo)
                var intersecPointX = math.parser().eval("(" + _seg.x1 + "*" + _seg.y2 + "-" + _seg.x2 + "*" + _seg.y1 + "-" + "(" + _seg.x1 + "-" + _seg.x2 + ")" + "*" + _Y + ")" + "/" + "(" + _seg.y2 + "-" + _seg.y1 + ")");
                _X <= intersecPointX ? _intersecNum++ : _intersecNum + 0;
            }
        }
        return _intersecNum % 2 == 0 ? _ifContain = false : _ifContain = true;
    }
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末纤子,一起剝皮案震驚了整個濱河市搬瑰,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌控硼,老刑警劉巖泽论,帶你破解...
    沈念sama閱讀 217,826評論 6 506
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場離奇詭異卡乾,居然都是意外死亡翼悴,警方通過查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,968評論 3 395
  • 文/潘曉璐 我一進(jìn)店門幔妨,熙熙樓的掌柜王于貴愁眉苦臉地迎上來鹦赎,“玉大人,你說我怎么就攤上這事误堡」呕埃” “怎么了?”我有些...
    開封第一講書人閱讀 164,234評論 0 354
  • 文/不壞的土叔 我叫張陵锁施,是天一觀的道長陪踩。 經(jīng)常有香客問我杖们,道長,這世上最難降的妖魔是什么膊毁? 我笑而不...
    開封第一講書人閱讀 58,562評論 1 293
  • 正文 為了忘掉前任胀莹,我火速辦了婚禮,結(jié)果婚禮上婚温,老公的妹妹穿的比我還像新娘描焰。我一直安慰自己,他們只是感情好栅螟,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,611評論 6 392
  • 文/花漫 我一把揭開白布荆秦。 她就那樣靜靜地躺著,像睡著了一般力图。 火紅的嫁衣襯著肌膚如雪步绸。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 51,482評論 1 302
  • 那天吃媒,我揣著相機(jī)與錄音瓤介,去河邊找鬼。 笑死赘那,一個胖子當(dāng)著我的面吹牛刑桑,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播募舟,決...
    沈念sama閱讀 40,271評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼祠斧,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了拱礁?” 一聲冷哼從身側(cè)響起琢锋,我...
    開封第一講書人閱讀 39,166評論 0 276
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎呢灶,沒想到半個月后吴超,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 45,608評論 1 314
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡鸯乃,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,814評論 3 336
  • 正文 我和宋清朗相戀三年烛芬,在試婚紗的時候發(fā)現(xiàn)自己被綠了。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片飒责。...
    茶點(diǎn)故事閱讀 39,926評論 1 348
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖仆潮,靈堂內(nèi)的尸體忽然破棺而出宏蛉,到底是詐尸還是另有隱情,我是刑警寧澤性置,帶...
    沈念sama閱讀 35,644評論 5 346
  • 正文 年R本政府宣布拾并,位于F島的核電站,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏嗅义。R本人自食惡果不足惜屏歹,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,249評論 3 329
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望之碗。 院中可真熱鬧蝙眶,春花似錦、人聲如沸褪那。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,866評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽博敬。三九已至友浸,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間偏窝,已是汗流浹背收恢。 一陣腳步聲響...
    開封第一講書人閱讀 32,991評論 1 269
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留祭往,地道東北人伦意。 一個月前我還...
    沈念sama閱讀 48,063評論 3 370
  • 正文 我出身青樓,卻偏偏與公主長得像链沼,于是被迫代替她去往敵國和親默赂。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,871評論 2 354

推薦閱讀更多精彩內(nèi)容