在shp2pgsql中自动修复环的自交

Question

我们正在将一大堆ArcGIS shapefile导入到PostGIS中，并使用shp2pgsql进行动态转换。问题是，如果shapefile有任何环形自交点，导入就会阻塞：

NOTICE:  Ring Self-intersection at or near point -80.1338 25.8102
ERROR:  new row for relation "place_shapes" violates
  check constraint "shape_is_valid"

我们该如何解决这个问题呢？

Answer 1

虽然按照Marcelo的建议，对自交多边形(在shp文件中太常见)来说，通过零缓冲特征是一种已知的修复方法，但也有用于此目的的ST_MakeValid函数。还有一个相关的函数，ST_IsValidReason，它将通知问题出在哪里，而不是盲目地试图修复它。

实际上，使用ST_MakeValid(geom)或ST_Buffer(geom，0)可能会产生混合的几何图形类型，包括孤立的点和线串。因此，进一步的改进可能是检查返回的几何类型，并且只包括例如从ST_MakeValid产生的多边形。

create table valid_geoms as
with make_valid (id, geom) as 
   (select 
      row_number() over() as id, 
     (ST_Dump(ST_MakeValid(geom))).geom as geom from invalid_table
  )
select id, geom from make_valid where ST_GeometryType(geom)='ST_Polygon';

其中，invalid_table是原始shp2pgsql导入产生的表。

我在这里包含了一个生成的id，因为ST_MakeValid可能会从输入几何图形中生成多个多边形。可以重写查询以包括原始id字段，但不再保证它是唯一的。

Answer 2

这个查询通常会帮我解决这个问题：

UPDATE place_shapes
  SET geometry=ST_Buffer(geometry, 0.0);

Answer 3

事实证明，运行这个后处理步骤可以解决这个问题：

UPDATE place_shapes
  SET geometry=ST_SimplifyPreserveTopology(geometry, 0.0001)
  WHERE ST_IsValid(geometry) = false;

0.0001是以度为单位的容差，您可能需要根据自己的喜好对其进行调整，但对于街道地图数据似乎是正确的。

此外，如果形状表强制执行有效性(这是应该的)，则需要使用shp2pgsql创建一个临时的未验证的表，在那里修复多边形，然后才将它们复制到主表中。