使用postgis数据库进行多边形裁切线
背景：有一份polyline的基础数据，有一个多边形，求出多边形内的所有polylinePostGIS参考手册： http://postgis.net/docs/reference.html
1、polyline数据表、qgis可视化
2、polygon数据表、qgis可视化
3、使用ST_Covers，求出所有完全包含在面内的线，有一部分线是被多边形穿过了，此部分线是没有在covers里的
SELECT f.geom FROM "freeway" f LEFT JOIN area r ON 1 = 1 WHERE ST_Covers (r.geom, f.geom) = 't'
4、获取多边形边界（ST_Boundary）和polyline的相交的线（ST_Intersects），并用边界分割（ST_Split）这些线，此时结果为GEOMETRYCOLLECTION，需要提取(ST_Dump)几何图形组件的集合（GEOMETRYCOLLECTION）中的图形，涉及到的函数如下：
ST_Boundary(geometry geomA)
ST_Intersects(geometry geomA , geometry geomB )
ST_Split(geometry input , geometry blade)
ST_Dump(geometry geom)
SELECT(ST_Dump (ST_Split (l.geom, ST_Boundary (rr.geom)))).geom FROM( SELECT f.geom FROM "freeway" f LEFT JOIN area r ON 1 = 1 WHEREST_Intersects (ST_Boundary (r.geom), f.geom) = 't') lLEFT JOIN area rr ON 1 =1
5、第4步骤分割完成后，一部分线在面内，一部分线在面外，此时：
面内的线与面还有交点，所以ST_Covers不可用
面外的线与面还有交点，所以ST_Intersects不可用
思路为：判断所有线的中心点，中心点在面内，则线段被面完全包含
SELECT f.geom FROM res_3 f LEFT JOIN area r ON 1 = 1 WHEREST_Covers (r.geom, ST_Centroid (f.geom)) = 't';
6、把1的查询结果和5的结果union一下，最终结果如下
页:
[1]