在PostGIS中计算两个多边形共享边的长度

Question

我有一个包含几个多角体的数据库。每个多边形都有很多定义它的节点。此外，每个多聚体至少有一个邻居。为了进行计算，我需要确定两个相邻多边形的共享边的长度。在示例图片中，显示了这样的情况。这两个多边形共用红色的边。如何在PostGIS的帮助下计算红边的长度？我没找到那个功能。

​

​

Answer 1

你可以得到两个多边形的交叉点，然后得到任意LineStrings的长度。例如，以两个几何图形为例：

SELECT ST_Length(ST_CollectionExtract(ST_Intersection(a_geom, b_geom), 2))
FROM (
  SELECT
    'POLYGON((70 170,160 250,240 220,270 130,185 62,90 80,70 170))'::geometry AS a_geom,
    'POLYGON((160 250,236 314,380 290,390 150,270 130,240 220,160 250))'::geometry AS b_geom
) f;

或对表mypoly中的所有多边形执行此操作，即触摸

SELECT a.gid AS gid_a, b.gid AS gid_b,
  ST_Length(ST_CollectionExtract(ST_Intersection(a.geom, b.geom), 2))
FROM mypoly a, mypoly b
WHERE a.gid < b.gid AND ST_Touches(a.geom, b.geom);

您还可以探索拓扑扩展，它是PostGIS 2.x的一部分。