计算 PostGIS 中两个多边形共享边的长度
Calculate length of a shared edge of two polygones in PostGIS
我有一个包含多个多边形的数据库。每个多边形都有很多定义它的节点。此外,每个多边形至少有一个邻居。对于计算,我需要确定两个相邻多边形的共享边的长度。在示例图片中显示了这种情况。两个多边形共享红色边缘。如何借助 PostGIS 计算红边的长度?我没有找到相应的功能。
您可以获取两个多边形的 intersection,然后获取任何 LineString 的长度。例如,取两个几何图形:
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;
或者对 table mypoly
中的所有多边形对执行此操作 touch:
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);
您还可以探索 Topology extension,它是 PostGIS 2.x 的一部分。
我有一个包含多个多边形的数据库。每个多边形都有很多定义它的节点。此外,每个多边形至少有一个邻居。对于计算,我需要确定两个相邻多边形的共享边的长度。在示例图片中显示了这种情况。两个多边形共享红色边缘。如何借助 PostGIS 计算红边的长度?我没有找到相应的功能。
您可以获取两个多边形的 intersection,然后获取任何 LineString 的长度。例如,取两个几何图形:
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;
或者对 table mypoly
中的所有多边形对执行此操作 touch:
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);
您还可以探索 Topology extension,它是 PostGIS 2.x 的一部分。