ST_MakePolygon — 根据外壳和可选的孔列表创建多边形。
geometry ST_MakePolygon(
geometry linestring)
;
geometry ST_MakePolygon(
geometry outerlinestring, geometry[] interiorlinestrings)
;
根据给定的外壳和可选的孔数组创建多边形。输入几何必须是闭合线串(环)。
变体 1:接受一个外壳线串。
变体 2:接受一个外壳线串和一个内部(孔)线串数组。可以使用 PostgreSQL array_agg()、ARRAY[] 或 ARRAY() 构造函数构造几何数组。
此函数不接受多线串。使用 ST_LineMerge 生成线串,或使用 ST_Dump 提取线串。 |
此函数支持 3d,不会删除 z 索引。
根据 2D 线串创建多边形。
SELECT ST_MakePolygon( ST_GeomFromText('LINESTRING(75 29,77 29,77 29, 75 29)'));
根据开线串创建多边形,使用 ST_StartPoint 和 ST_AddPoint 将其闭合。
SELECT ST_MakePolygon( ST_AddPoint(foo.open_line, ST_StartPoint(foo.open_line)) ) FROM ( SELECT ST_GeomFromText('LINESTRING(75 29,77 29,77 29, 75 29)') As open_line) As foo;
根据 3D 线串创建多边形
SELECT ST_AsEWKT( ST_MakePolygon( 'LINESTRING(75.15 29.53 1,77 29 1,77.6 29.5 1, 75.15 29.53 1)')); st_asewkt ----------- POLYGON((75.15 29.53 1,77 29 1,77.6 29.5 1,75.15 29.53 1))
根据具有度量的线串创建多边形
SELECT ST_AsEWKT( ST_MakePolygon( 'LINESTRINGM(75.15 29.53 1,77 29 1,77.6 29.5 2, 75.15 29.53 2)' )); st_asewkt ---------- POLYGONM((75.15 29.53 1,77 29 1,77.6 29.5 2,75.15 29.53 2))
创建带有一个额外孔的圆环多边形
SELECT ST_MakePolygon( ST_ExteriorRing( ST_Buffer(ring.line,10)), ARRAY[ ST_Translate(ring.line, 1, 1), ST_ExteriorRing(ST_Buffer(ST_Point(20,20),1)) ] ) FROM (SELECT ST_ExteriorRing( ST_Buffer(ST_Point(10,10),10,10)) AS line ) AS ring;
创建一组省界,其中孔表示湖泊。输入是省多边形/多边形表和水线串表。通过使用 ST_IsClosed 确定形成湖泊的线。通过使用 ST_Boundary 提取省界线。根据 ST_MakePolygon
的要求,通过使用 ST_LineMerge 将边界强制为单线串。(但是,请注意,如果一个省有多个区域或有岛屿,这将生成一个无效的多边形。)使用 LEFT JOIN 确保即使没有湖泊,也包括所有省份。
使用 CASE 构造函数,因为将空数组传递到 ST_MakePolygon 会导致 NULL 返回值。 |
SELECT p.gid, p.province_name, CASE WHEN array_agg(w.geom) IS NULL THEN p.geom ELSE ST_MakePolygon( ST_LineMerge(ST_Boundary(p.geom)), array_agg(w.geom)) END FROM provinces p LEFT JOIN waterlines w ON (ST_Within(w.geom, p.geom) AND ST_IsClosed(w.geom)) GROUP BY p.gid, p.province_name, p.geom;
另一种技术是利用相关子查询和将行集转换为数组的 ARRAY() 构造函数。
SELECT p.gid, p.province_name, CASE WHEN EXISTS( SELECT w.geom FROM waterlines w WHERE ST_Within(w.geom, p.geom) AND ST_IsClosed(w.geom)) THEN ST_MakePolygon( ST_LineMerge(ST_Boundary(p.geom)), ARRAY( SELECT w.geom FROM waterlines w WHERE ST_Within(w.geom, p.geom) AND ST_IsClosed(w.geom))) ELSE p.geom END AS geom FROM provinces p;