Análisis del borde del polígono PostGIS (orientación, longitud del borde)

Soy bastante nuevo en el mundo de los SIG y especialmente en PostGIS, así que discúlpeme si la respuesta parece evidente ...

Me gustaría hacer un análisis de varios edificios. Una cosa que me interesa es su superficie de fachada junto con la orientación respectiva. Como se ilustra en la imagen a continuación, me gustaría tener la longitud y la orientación (normal) de todos los bordes en una serie de polígonos. En el ejemplo, destaqué solo una superficie.

Una tabla de resultados podría verse así:

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

Sin embargo, no estoy seguro de si es una forma inteligente de almacenar el resultado para su posterior procesamiento (por ejemplo, calcular la distancia desde el borde hasta el próximo edificio, etc.). Entonces mi pregunta es doble:

1. ¿Existe una función eficiente de PostGIS que pueda analizar los bordes de un polígono? En caso de que no haya una función nativa de PostGIS, alternativamente me interesaría un enfoque basado en Python.
2. ¿Cuál sería una forma inteligente de almacenar el resultado en una tabla PostGIS, ya que los polígonos pueden tener diferentes números de bordes?

Ayer no tuve tiempo para crearlo en detalles ... Vea mi solución en 4 pasos:

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 última consulta le proporciona los identificadores de construcción con una unión espacial usando st_touches. Espero eso ayude. Actualización: en qgis la solución se ve así: