概述
多段線切割多邊形，直觀上來說首先就是求多段線與多邊形之間的交點，然后利用這些交點與多邊形頂點構成多邊形，這樣一個多邊形就被切割成多個多邊形。如下圖所示：
多段線與多邊形
切割后的多邊形
主要過程
以下圖所示的分割為例，做簡短說明：
示例
難點
- 如何連接交點與頂點？如上圖所示，切割的多邊形為AB21，D34，E56，12C3456F；當我們求出交點之后，如何組織交點與頂點的連接很重要，一個解決方法是將多邊形的邊界與多段線組成有向邊，依據(jù)他們的方向來進行組織連接；
- 剔除不合理的多邊形?很明顯，C23這個多變形不符合要求，因為它位于多邊形的外側，這個要去除掉；
過程
對于上述難點，本文主要是利用geos中構建多邊形的方法來解決，尤其是難點1，具體的原理可參考直線切割凹多邊形
1.獲取帶分割的幾何圖形，注意該幾何圖形是線形的，主要包括待切割的閉合環(huán)(如多邊形的邊界)以及分割線；
geos::geom::Geometry *tempG = _geos_geometry->getBoundary();
geos::geom::Geometry *tempG1 = tempG->Union(_asGeos(splitline));
2.構建多邊形，在geos主要是Polygonizer類，獲取所有分割的多邊形，當然也包括類似C23這樣不合理的多邊形，下一步是剔除這些多邊形；
geos::operation::polygonize::Polygonizer pz;
pz.add(tempG1);
std::vector< geos::geom::Polygon*> *result_geom = pz.getPolygons();
3.剔除不合理多邊形，理論上來說，切割出來的多變?nèi)绻辉级噙呅伟涂梢员Ａ?，調(diào)用geos中的Cover或Contains等拓撲計算便可以，但是在實際操作過程中有兩個問題：
(1)拓撲計算不正確，每次返回值都是false,猜測可能是精度問題，因為如果坐標都是整型，可以返回正確結果;
(2)效率不高，相對于contain和within,建議在geos中優(yōu)先使用cover或coverby。具體可參考九交模型
所以選取一個折中方案，如下所示：
//過濾掉不在多邊形內(nèi)部的切割面;
//不能直接用包含和覆蓋直接判斷，存在精度問題;
//利用切割面與原始多邊形相交后的多邊形面積與切割面的面積之比判斷;
//比值在0.99~1.01之間;
4.最后得到的就是一個vector集合，數(shù)據(jù)類型是geos::geom::Polygon*
