PostGIS多边形边缘分析（方向，边缘长度）

我对GIS尤其是PostGIS不太熟悉，因此，如果答案似乎显而易见，请原谅...

我想对许多建筑物进行分析。我感兴趣的一件事是它们的外墙表面以及各自的方向。如下图所示，我希望具有一系列多边形中所有边的长度和（法线）方向。在示例中，我仅突出显示一个表面。

结果表可能是这样的：

building_id | edge_id | orientation | edge_length
-------------------------------------------------
      1     |    1    |     315     |    10.0
      1     |    2    |      45     |     7.0
      1     |   ...   |     ...     |     ...

但是，我不确定这是否是存储结果以进行进一步处理的明智方法（例如，计算从边缘到下一栋建筑物的距离等）。所以我的问题是双重的：

1. 是否有一个有效的PostGIS功能可以分析多边形的边缘？如果没有本地PostGIS功能，我可能会对基于Python的方法感兴趣。
2. 由于多边形的边数可能不同，因此将结果存储在PostGIS表中的聪明方法是什么？

昨天，我没有时间详细创建它...分4步查看解决方案：

CREATE OR REPLACE VIEW bd_segment AS
SELECT
      ST_PointN(geom, generate_series(1, ST_NPoints(geom)-1)) AS sp,
      ST_PointN(geom, generate_series(2, ST_NPoints(geom)  )) AS ep
    FROM
       -- extract the individual linestrings
       (SELECT (ST_Dump(ST_Boundary(the_geom))).geom
       FROM bd) AS linestrings;

CREATE OR REPLACE VIEW bd_segment_geom AS
SELECT sp, ep, st_makeline(sp, ep) 
FROM bd_segment;

CREATE OR REPLACE view bd_segment_id AS 
SELECT bd.gid, row_number() 
    OVER (order by bd.gid), degrees(st_azimuth(ff.sp, ff.ep)-1.57079633) AS az_deg,
    ST_LENGTH(ff.st_makeline) , ff.st_makeline FROM bd_segment_geom ff
JOIN bd ON st_touches(ff.st_makeline, bd.the_geom)
GROUP BY bd.gid, ff.sp, ff.ep, ff.st_makeline;

UPDATE bd_segment_id
SET az_deg = az_deg + 360
WHERE az_deg < 0;

最后一个查询使用st_touches为您提供带有空间连接的建筑物ID。希望能帮助到你。更新-在qgis中，解决方案如下所示：