Análise de arestas de polígono PostGIS (orientação, comprimento da aresta)

Sou um pouco novo no mundo dos SIG e, principalmente, do PostGIS, então, desculpe-me se a resposta parecer evidente ...

Eu gostaria de fazer uma análise em vários edifícios. Uma coisa que me interessa são as superfícies das fachadas, juntamente com a respectiva orientação. Como ilustrado na figura abaixo, eu gostaria de ter o comprimento e a orientação (normal) de todas as arestas em uma série de polígonos. No exemplo, destaquei apenas uma superfície.

Uma tabela de resultados pode ficar assim:

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

No entanto, não tenho certeza se é uma maneira inteligente de armazenar o resultado para processamento adicional (por exemplo, calcular a distância da borda até o próximo edifício, etc.). Então, minha pergunta é dupla:

1. Existe uma função PostGIS eficiente que pode analisar as arestas de um polígono? Caso não exista uma função PostGIS nativa, eu estaria interessado em uma abordagem baseada em Python.
2. Qual seria uma maneira inteligente de armazenar o resultado em uma tabela PostGIS, já que os polígonos podem ter diferentes números de arestas?

Ontem não tive tempo de criá-lo em detalhes ... Veja minha solução em 4 etapas:

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;

A última consulta fornece os IDs de construção com uma junção espacial usando st_touches. Espero que ajude. Atualização - No qgis, a solução fica assim: