Analiza krawędzi wielokąta PostGIS (orientacja, długość krawędzi)

Jestem raczej nowy w świecie GIS, a zwłaszcza PostGIS, więc przepraszam, jeśli odpowiedź wydaje się oczywista ...

Chciałbym przeprowadzić analizę szeregu budynków. Interesuje mnie ich powierzchnia elewacji wraz z odpowiednią orientacją. Jak pokazano na poniższym obrazku, chciałbym mieć długość i (normalną) orientację wszystkich krawędzi w szeregu wielokątów. W przykładzie wyróżniłem tylko jedną powierzchnię.

Tabela wyników może wyglądać następująco:

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

Nie jestem jednak pewien, czy jest to inteligentny sposób na przechowywanie wyników do dalszego przetwarzania (np. Oblicz odległość od krawędzi do następnego budynku itp.). Więc moje pytanie jest dwojakie:

1. Czy istnieje skuteczna funkcja PostGIS, która może analizować krawędzie wielokąta? W przypadku, gdy nie ma natywnej funkcji PostGIS, alternatywnie byłbym zainteresowany podejściem opartym na języku Python.
2. Jaki byłby inteligentny sposób przechowywania wyniku w tabeli PostGIS, ponieważ wielokąty mogą mieć różną liczbę krawędzi?

Wczoraj nie miałem czasu na szczegółowe tworzenie ... Zobacz moje rozwiązanie w 4 krokach:

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;

Ostatnie zapytanie daje identyfikatory budynku z łączeniem przestrzennym za pomocą st_touches. Mam nadzieję, że to pomoże. Aktualizacja - w qgis rozwiązanie wygląda następująco: