PostGIS-Polygonkantenanalyse (Ausrichtung, Kantenlänge)

Ich bin ziemlich neu in der Welt von GIS und insbesondere von PostGIS. Entschuldigen Sie mich bitte, wenn die Antwort offensichtlich erscheint ...

Ich möchte einige Gebäude analysieren. Eine Sache, die mich interessiert, sind ihre Fassadenoberflächen zusammen mit der jeweiligen Ausrichtung. Wie im Bild unten dargestellt, möchte ich die Länge und (normale) Ausrichtung aller Kanten in einer Reihe von Polygonen haben. Im Beispiel habe ich nur eine Oberfläche hervorgehoben.

Eine Ergebnistabelle könnte folgendermaßen aussehen:

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

Ich bin mir jedoch nicht sicher, ob es eine clevere Möglichkeit ist, das Ergebnis für die weitere Verarbeitung zu speichern (z. B. Abstand von der Kante zum nächsten Gebäude berechnen usw.). Meine Frage ist also zweifach:

1. Gibt es eine effiziente PostGIS-Funktion, mit der die Kanten eines Polygons analysiert werden können? Falls es keine native PostGIS-Funktion gibt, wäre ich alternativ an einem Python-basierten Ansatz interessiert.
2. Was wäre eine clevere Möglichkeit, das Ergebnis in einer PostGIS-Tabelle zu speichern, da die Polygone möglicherweise eine unterschiedliche Anzahl von Kanten haben?

Gestern hatte ich keine Zeit, es im Detail zu erstellen ... Sehen Sie meine Lösung in 4 Schritten:

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;

Bei der letzten Abfrage erhalten Sie die Gebäude-IDs mit einer räumlichen Verknüpfung mithilfe von st_touches. Ich hoffe es hilft. Update - In qgis sieht die Lösung folgendermaßen aus: