如何使用PostGIS将多边形数据转换为线段

Question

我在 PostgreSQL/PostGIS 中有一个多边形数据表。现在我需要将此多边形数据转换为其相应的线段。谁能告诉我如何使用 PostGIS 查询进行转换。

提前致谢

Answer 1

一般来说，将多边形转换为线可能并不简单，因为没有人一对一映射 和多边形的各种元素映射到不同的线串（外环、内环等）。

考虑到这一点，您需要按照如下可能的方法单独分割每个数据：

SELECT ST_AsText( ST_MakeLine(sp,ep) )
FROM
   -- extract the endpoints for every 2-point line segment for each linestring
   (SELECT
      ST_PointN(geom, generate_series(1, ST_NPoints(geom)-1)) as sp,
      ST_PointN(geom, generate_series(2, ST_NPoints(geom)  )) as ep
    FROM
       -- extract the individual linestrings
      (SELECT (ST_Dump(ST_Boundary(geom))).geom
       FROM mypolygontable
       ) AS linestrings
    ) AS segments;

根据 mypolygontable 中存储的多边形数据，您可能不仅需要转储边界（如上所述使用 >ST_Boundary），还有其他元素。上面的代码具有更详细的概述，取自 postgis-users 列表： 将多边形分割为 N 条线

对于 将线串或多边形分解为PostGIS 中的各个向量

Answer 2

自 PostGIS 3.2.0（2021 年 12 月）开始，您可以< a href="https://postgis.net/docs/ST_DumpSegments.html" rel="nofollow noreferrer">ST_DumpSegments() 直接。它甚至可以处理从多边形中提取的所有环（不仅仅是外部/边界）：
_{db<>fiddle 演示}

select st_astext((st_dumpsegments('polygon((0 0,1 2,3 3,0 0))')).geom);

st_astext
LINESTRING(0 0,1 2)
LINESTRING (1 2,3 3)
LINESTRING(3 3,0 0)

上面的示例调用 SRF 在选择 部分而不是通常应位于的 from 列表中，并且它会忽略部分转储。这可能更清晰一些：

select dmp.ordinality
      ,dmp.path
      ,st_astext(dmp.geom)
from (values('polygon((0 0,1 2,3 3,0 0))'))as a(g)
cross join lateral st_dumpsegments(g)with ordinality as dmp;

序数	路径	st_astext
1	{1,1}	LINESTRING(0 0,1 2)
2	{1,2}	LINESTRING(1 2,3 3)
3	{1,3}	LINESTRING(3 3,0 0)

WITH ORDINALITY 对于单环多边形（单线串）来说非常方便，因为您可以直接获取线段编号，而无需到达 路径内部 数组。
对于多环，path 会加宽以保存该段来自哪个环的信息，而 ordinality 则不会显示这种区别。

Answer 3

当你搜索这个问题时，这是谷歌上的第一个点击。我不知道是否已经过去了这么长时间才创建了一个函数，但对于未来的谷歌用户来说 ST_ExteriorRings(geom) 对我来说非常有用。
http://postgis.net/docs/ST_ExteriorRing.html

SELECT ST_ExteriorRing(ST_Dump(geom)).geom)
FROM foo