Analisis tepi poligon PostGIS (orientasi, panjang tepi)

Saya agak baru di dunia GIS dan terutama PostGIS, jadi tolong permisi jika jawabannya tampak jelas ...

Saya ingin melakukan analisis pada sejumlah bangunan. Satu hal yang saya tertarik adalah permukaan fasad mereka bersama dengan orientasi masing-masing. Seperti yang diilustrasikan dalam gambar di bawah ini, saya ingin memiliki panjang dan orientasi (normal) semua tepi dalam serangkaian poligon. Dalam contoh saya hanya menyoroti satu permukaan.

Tabel hasil dapat terlihat seperti ini:

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

Namun, saya tidak yakin apakah ini cara yang cerdas untuk menyimpan hasil untuk diproses lebih lanjut (mis. Menghitung jarak dari tepi ke gedung berikutnya, dll.). Jadi pertanyaan saya ada dua:

1. Apakah ada fungsi PostGIS yang efisien yang dapat menganalisis tepi poligon? Dalam hal tidak ada fungsi PostGIS asli saya atau akan tertarik pada pendekatan berbasis Python.
2. Apa yang akan menjadi cara cerdas untuk menyimpan hasilnya dalam tabel PostGIS, karena poligon mungkin memiliki jumlah sisi yang berbeda?

Kemarin saya tidak punya waktu untuk membuatnya secara detail ... Lihat solusi saya dalam 4 langkah:

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;

Kueri terakhir memberi Anda id bangunan dengan gabungan spasial menggunakan st_touches. Semoga ini bisa membantu. Perbarui - Dalam qgis solusinya terlihat seperti ini: