PostGISを使用して道路の交差点を特定する

道路が交差する場所を特定し、交差点を形成する道路の数をリストして、この交差点でポイントを作成しようとしています。

これを達成するためにST_NumPointsを使用する方法があるかどうか疑問に思っていましたが、私は何をすべきかを理解できません。次のコードを使用して、線が交差するポイントのテーブルを作成しました。

CREATE TABLE test_points as
SELECT      
    ST_Intersection(a.geom, b.geom),
    a.gid
FROM
    roads as a,
    roads as b
WHERE
    ST_Touches(a.geom, b.geom);

道路のサンプルでこれを実行すると、次の点のグリッドが得られます（道路は説明のために示されています）。

ポイントの1つを調べると、多くのポイントが積み重なっていることがわかります。

ここのGIDは道路IDですが、なぜ多くのポイントがあるのか​​理解できません。中央の道路の交差点で4つのポイントがカウントされることは理解できますが、ここには12ポイントがリストされています。PostGISでこの計算を実行するより良い方法はありますか？

グループ化する場合、一意のポイントのみを取得する必要があります。

CREATE TABLE test_points as
SELECT      
    ST_Intersection(a.geom, b.geom),
    Count(Distinct a.gid)
FROM
    roads as a,
    roads as b
WHERE
    ST_Touches(a.geom, b.geom)
    AND a.gid != b.gid
GROUP BY
    ST_Intersection(a.geom, b.geom)
;

これは、予想より少し複雑です。ペア以上の関係を分析する良い方法がないためです。3行を関数に入れて、それらすべてが交差するかどうかを尋ねることはできません。

ただし、少なくとも1つの方法として、最初に交差点を見つけてから、各交差点で何本の道路が接触しているかを確認することができます（同じクエリですべて実行できます）。

道路が互いに完全に接続し、交差点を通過する道路がない場合、次のようなことができます（テストされていません）：
忘れられたグループ句で編集されています（まだテストされていません）：

SELECT distinct_crosspoints.geom as crossing, array_agg(roads.gid), count(*) FROM
  (SELECT DISTINCT (geom) geom FROM 
    (SELECT ST_Intersection(a.geom, b.geom) geom 
     FROM roads a, roads b 
     WHERE ST_Intersects(a.geom, b.geom)
    ) all_crosspoints
   ) distinct_crosspoints
   ,roads 
 WHERE ST_Intersects(distinct_crosspoints.geom, roads.geom)
 GROUP BY distinct_crosspoints.geom;

道路が適切に接続されていない場合や、一部の道路が交差点を通過している場合は、より複雑になります。

HTH

ニクラス

 CREATE TABLE test_points as
    SELECT      
        ST_Intersection(a.geom, b.geom),
        Count(Distinct a.gid)
    FROM
        roads as a,
        roads as b
    WHERE
        ST_Touches(a.geom, b.geom)
        AND a.gid < b.gid   /* !!! Changed "!=" for "<"  */
    GROUP BY
        ST_Intersection(a.geom, b.geom)
    ;

ラインA（id 1）がラインB（id 2）と交差する場合、これが必要な交点です。ただし、ラインBは同じポイントでラインAと交差しています。しかし、この点は二度必要ではありません。それが私がa.gid < b.gid 代わりに使用している理由ですa.gid != b.gid