使用PostGIS识别“长而窄”的多边形

我有一组代表大面积的多边形，例如城市社区。我想确定它们之间的大重叠区域。

但是存在一个问题：有时这些多边形会沿其周界重叠（因为它们绘制的精度很低）。这将产生我不关心的长而窄的重叠。

但其他时候，健壮多边形会出现大量重叠，这意味着邻域的多边形会与另一个多边形重叠的较大区域。我只选择这些。

请参见下面的重叠图片。想象一下，我只想选择左下角的蓝色多边形。

我可以看一下区域，但是有时狭窄的区域是如此之长，以至于最终它们会具有与蓝色多边形一样大的区域。我试图做一个面积/周长的比率，但这也产生了混合的结果。

我什至尝试使用ST_MinimumClearance，但是有时大区域会附上一个狭窄的部分，或者两个非常接近的顶点。

对其他方法有什么想法吗？

最后，最适合我的方法是使用负缓冲区，如下面的@Cyril和@FGreg所建议。

我使用了类似的东西：

ST_Area(ST_Buffer(geom, -10)) as neg_buffer_area

在我的情况下，单位为米，因此10 m负缓冲。

对于窄多边形，该区域返回零（同样，几何图形将为空）。然后，我使用此列来过滤出狭窄的多边形。

我会尝试创建一个负缓冲区，如果它吃了薄多边形，那很好，如果它不吃多边形，那是我的... :-)

运行此脚本，预先设置线性多边形的宽度的2/3 ...

create table name_table as
SELECT ST_Buffer(
(ST_Dump(
(ST_Union(
ST_Buffer(
(geom),-0.0001))))).geom,
0.0001)) as geom from source_table

操作系统：-）...

代替面积/周长，最好使用面积除以周长的平方（或其倒数）。

这也称为“形状指数”。周长除以面积的平方的最小值为4 * Pi（）（对于磁盘，这是最紧凑的2D几何形状），因此可以通过4 * Pi（）进行归一化，解释（接近1的归一化值意味着您拥有非常紧凑的对象，并且正方形的值大约为1.27）。

编辑：该区域的阈值将对消除非常小的伪像很有用，它可能很紧凑。然后，形状指数将显示出更好的对比度。编辑：除了此答案，使用ST_Snap可以帮助您在问题发生之前解决问题。

一种选择是使用多边形面积与可以使用其末端绘制的最长线的比率。识别狭窄的长多边形。

select * from polygons where ST_Length(ST_LongestLine(geom, geom)) < ST_Area(geom) * 4

这对于条形多边形非常有效。您可以调整比例（乘以面积）以适合您的需求和投影。

听起来这可能与您的用例匹配：消除选定的多边形

通过擦除输入多边形的公共边界，将输入层的选定多边形与某些相邻多边形合并。相邻的多边形可以是面积最大或最小的多边形，也可以是与要消除的多边形共享最大公共边界的多边形。

消除通常用于消除条形多边形，即微小多边形，这是多边形交集过程的结果，在该过程中，输入的边界相似但不相同。

听起来您想尝试使用“最大公共边界”选项。

在我看来，这似乎是PostGIS拓扑扩展的理想用例。拓扑的公差参数将确定允许顶点捕捉到其他现有多边形的距离，以应对源数据的低精度并进行清理。

简而言之，策略是：

1.启用拓扑扩展

CREATE EXTENSION postgis_topology;

2.创建一个新的空拓扑

SELECT topology.CreateTopology('neighborhoods_topo', 4326, 1e-7);

第三个参数是公差，以CRS为单位；明智地选择它。理想情况下，您需要单位为米的CRS。如果CRS单位不是米，例如WGS 84 aka 4326，则用于ST_Transform重新投影多边形。

3.将“ TopoGeometry”列添加到多边形表

SELECT topology.AddTopoGeometryColumn('neighborhoods_topo', 'public', 'neighborhoods', 'topogeom', 'POLYGON');

这将返回一个新的layer_id。保存它，以后将需要它。1如果您从头开始，它将是分层的，并在每次新呼叫时增加。

4.将所有多边形添加到拓扑中

UPDATE public.neighborhoods
SET topogeom = topology.toTopoGeom(geom, 'neighborhoods_topo', 1, 1e-7);

对于大型数据集，这可能需要几个小时，请耐心等待。1是较早返回的layer_id。

5.查找出现在多个社区中的面孔

从拓扑中找到存在于2个或更多拓扑中的所有面。我将把查询保留为练习。最简单的方法可能是使用该GetTopoGeomElements函数，然后按面部ID分组，然后查看计数为2或更大的对象。或者，您可以使用拓扑表列中的清理后的几何图形创建一个新表，将其转换为标准几何图形topogeom::geometry，然后重复执行现在已有的操作，但现在使用一个干净的数据集而不会出现条重叠。