Анализ ребер полигонов PostGIS (ориентация, длина ребер)

Я довольно новичок в мире ГИС и особенно в PostGIS, поэтому извините, если ответ кажется очевидным ...

Я хотел бы сделать анализ ряда зданий. Одна вещь, которая меня интересует, это их фасадные поверхности вместе с соответствующей ориентацией. Как показано на рисунке ниже, я хотел бы иметь длину и (нормальную) ориентацию всех ребер в серии многоугольников. В примере я выделил только одну поверхность.

Таблица результатов может выглядеть так:

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

Однако я не уверен, что это разумный способ сохранить результат для дальнейшей обработки (например, рассчитать расстояние от края до следующего здания и т. Д.). Итак, мой вопрос состоит из двух частей:

1. Существует ли эффективная функция PostGIS, которая может анализировать ребра многоугольника? В случае отсутствия встроенной функции PostGIS, я бы заинтересовался подходом на основе Python.
2. Какой разумный способ сохранить результат в таблице PostGIS, поскольку полигоны могут иметь разное число ребер?

Вчера у меня не было времени, чтобы создать его в деталях ... Смотрите мое решение в 4 этапа:

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;

Последний запрос дает вам идентификаторы здания с пространственным объединением, используя st_touches. Надеюсь, поможет. Обновление - в qgis решение выглядит так: