Phân tích cạnh đa giác PostGIS (định hướng, chiều dài cạnh)

Tôi còn khá mới mẻ với thế giới của GIS và đặc biệt là PostGIS, vì vậy xin vui lòng cho tôi biết nếu câu trả lời có vẻ hiển nhiên ...

Tôi muốn làm phân tích trên một số tòa nhà. Một điều tôi quan tâm là bề mặt của chúng cùng với hướng tương ứng. Như minh họa trong hình dưới đây, tôi muốn có chiều dài và hướng (bình thường) của tất cả các cạnh trong một loạt các đa giác. Trong ví dụ tôi chỉ nhấn mạnh một bề mặt.

Một bảng kết quả có thể trông như thế này:

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

Tuy nhiên, tôi không chắc liệu đó có phải là một cách thông minh để lưu trữ kết quả để xử lý thêm hay không (ví dụ: tính khoảng cách từ cạnh đến tòa nhà tiếp theo, v.v.). Vì vậy, câu hỏi của tôi là gấp đôi:

1. Có chức năng PostGIS hiệu quả nào có thể phân tích các cạnh của đa giác không? Trong trường hợp không có chức năng PostGIS riêng, tôi cũng sẽ quan tâm đến cách tiếp cận dựa trên Python.
2. Điều gì sẽ là một cách thông minh để lưu trữ kết quả trong bảng PostGIS, vì các đa giác có thể có số cạnh khác nhau?

Hôm qua tôi không có thời gian để tạo ra nó một cách chi tiết ... Xem giải pháp của tôi trong 4 bước:

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;

Truy vấn cuối cùng cung cấp cho bạn các id tòa nhà với phép nối không gian bằng cách sử dụng st_touches. Hy vọng nó giúp. Cập nhật - Trong qgis giải pháp trông như thế này: