Analyse de bord de polygone PostGIS (orientation, longueur de bord)

Je suis plutôt nouveau dans le monde des SIG et surtout PostGIS, alors veuillez m'excuser si la réponse semble évidente ...

Je voudrais faire une analyse sur un certain nombre de bâtiments. Une chose qui m'intéresse est leurs surfaces de façade ainsi que l'orientation respective. Comme illustré dans l'image ci-dessous, je voudrais avoir la longueur et l'orientation (normale) de toutes les arêtes dans une série de polygones. Dans l'exemple, j'ai mis en évidence une seule surface.

Un tableau de résultats pourrait ressembler à ceci:

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

Cependant, je ne sais pas s'il s'agit d'un moyen intelligent de stocker le résultat pour un traitement ultérieur (par exemple, calculer la distance entre le bord et le bâtiment suivant, etc.). Ma question est donc double:

1. Existe-t-il une fonction PostGIS efficace qui peut analyser les bords d'un polygone? Dans le cas où il n'y a pas de fonction PostGIS native, je serais également intéressé par une approche basée sur Python.
2. Quelle serait une façon intelligente de stocker le résultat dans une table PostGIS, car les polygones peuvent avoir différents nombres d'arêtes?

Hier je n'ai pas eu le temps de la créer en détails ... Voir ma solution en 4 étapes:

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;

La dernière requête vous donne les identifiants de construction avec une jointure spatiale à l'aide de st_touches. J'espère que cela aide. Mise à jour - Dans qgis, la solution ressemble à ceci: