PostGISポリゴンエッジ分析（方向、エッジ長）

私はGISの世界、特にPostGISの世界にはかなり慣れていないので、答えが明らかな場合は失礼します...

多くの建物の分析をしたいと思います。私が興味を持っていることの1つは、それぞれの方向に沿ったファサードの表面です。下の図に示すように、一連のポリゴンのすべてのエッジの長さと（通常の）方向を設定したいと思います。この例では、1つのサーフェスのみを強調表示しています。

結果テーブルは次のようになります。

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

しかし、それをさらに処理するために結果を格納するスマートな方法であるかどうかはわかりません（たとえば、エッジから次の建物までの距離を計算するなど）。だから私の質問は2つあります：

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を使用した空間結合で建物IDを提供します。それが役に立てば幸い。更新-qgisでは、ソリューションは次のようになります。