Postgres / postgis - 确保所有多边形都已关闭

Question

我需要编写一个查询来查找任何未闭合的多边形并通过复制第一个点并创建一个额外的端点来关闭它们。

我能够 select 无效行：

SELECT delivery_zone_id, polygon from delivery_zone WHERE ST_IsClosed(polygon::geometry) = FALSE;

而且我能够从每个多边形中转储各个点：

SELECT delivery_zone_id, ST_AsText((dp).geom) FROM
    (SELECT delivery_zone_id, ST_DumpPoints(polygon::geometry) AS dp 
        FROM delivery_zone 
        WHERE ST_IsClosed(polygon::geometry) = FALSE
    ) AS coords;

结果如下所示：

1   POINT(-96.80037 33.09812)  ## Copy this point and add it to the set
1   POINT(-96.80427 33.0956)
1   POINT(-96.80401 33.09219)
1   POINT(-96.79603 33.09222)
1   POINT(-96.79346 33.09647)
1   POINT(-96.80037 33.09857)

4   POINT(-96.80037 33.099)    ## Copy this point and add it to the set
4   POINT(-96.80427 33.0956)
4   POINT(-96.80401 33.09219)
4   POINT(-96.79603 33.09222)
4   POINT(-96.79346 33.09647)
4   POINT(-96.80037 33.09923)

这是我 sql 技能不足的地方。我需要一些帮助来复制第一个点并使用该数据创建一个新的终点。欢迎伪查询 - 我只需要看看它可能是什么样子，我可以填补空白。

---

更新：最终解决方案

感谢下面 JGH 的回答，我能够创建以下更新查询。这将找到任何未闭合的多边形并通过复制第一个点添加一个新的端点。

Note: this will only work with simple "single" polygons. If you have complex outer and inner polygons, you will need to make some drastic changes to this query.

UPDATE delivery_zone dz
SET polygon=ST_MakePolygon(ST_AddPoint(subquery.openline, ST_PointN(subquery.openline, 1), -1))
FROM (
  SELECT delivery_zone_id, ST_ExteriorRing(polygon::geometry) AS openline 
  FROM delivery_zone WHERE ST_IsClosed(polygon::geometry) = FALSE
) AS subquery
WHERE dz.delivery_zone_id = subquery.delivery_zone_id;

Answer 1

您可以尝试使用线添加点，然后转换为多边形。

请注意，创建非闭合多边形是不可能的...不太确定您是如何获得的，希望您能首先将它们转换为一条线。

所以，我们的想法是获取直线，然后在最后一个位置 (-1) 添加一个点。这个点与这条线的第一个点（位置 1）相同。最后你可以转换为多边形

WITH src AS (
    SELECT ST_GeomFromText('LINESTRING(0 0, 0 1, 1 1, 1 0)') As openline)
SELECT st_asText(openline), 
        st_asText(ST_MakePolygon(st_addPoint(openline,st_PointN(openline,1),-1)))
FROM src;



      st_astext          |           st_astext
-----------------------------+--------------------------------
 LINESTRING(0 0,0 1,1 1,1 0) | POLYGON((0 0,0 1,1 1,1 0,0 0))
(1 row)