Analisi del bordo poligonale PostGIS (orientamento, lunghezza del bordo)

Sono piuttosto nuovo nel mondo dei GIS e in particolare PostGIS, quindi mi scusi se la risposta sembra evidente ...

Vorrei fare analisi su un certo numero di edifici. Una cosa che mi interessa sono le superfici delle loro facciate e il rispettivo orientamento. Come illustrato nella figura seguente, vorrei avere la lunghezza e l'orientamento (normale) di tutti i bordi in una serie di poligoni. Nell'esempio ho evidenziato solo una superficie.

Una tabella dei risultati potrebbe apparire così:

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

Tuttavia, non sono sicuro che sia un modo intelligente per memorizzare il risultato per ulteriori elaborazioni (ad es. Calcolare la distanza dal bordo al prossimo edificio, ecc.). Quindi la mia domanda è duplice:

1. Esiste un'efficace funzione PostGIS in grado di analizzare i bordi di un poligono? Nel caso in cui non esistesse una funzione PostGIS nativa, sarei in alternativa interessato ad un approccio basato su Python.
2. Quale sarebbe un modo intelligente per archiviare il risultato in una tabella PostGIS, dal momento che i poligoni possono avere un numero diverso di bordi?

Ieri non ho avuto il tempo di crearlo in dettaglio ... Vedi la mia soluzione in 4 passaggi:

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;

L'ultima query ti fornisce gli ID di costruzione con un join spaziale usando st_touches. Spero che sia d'aiuto. Aggiornamento - In qgis la soluzione è simile alla seguente: