PostGIS kullanarak yol kavşaklarını belirleme

Yolların birbiriyle nerede kesiştiğini belirlemeye ve bu kesişimde, kesişmeyi oluşturan yolların sayısı ile bir noktaya değinmeye çalışıyorum.

Bunu başarmak için ST_NumPoints kullanmanın bir yolu olup olmadığını merak ediyordum ama ne yapmam gerektiğini tam olarak anlayamıyorum . Ne yaptım hatları aşağıdaki kodu kullanarak kesiştiği noktaların bir tablo oluşturmaktır:

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);

Bunu bir yol örneğinde çalıştırırsam aşağıdaki nokta ızgarasını alırım (yollar gösterim için gösterilir):

Noktalardan birini incelersem, üst üste yığılmış birçok nokta olduğunu görürüm:

Buradaki GID yol kimliğidir, ancak neden çok fazla nokta olduğunu anlamıyorum. Merkezi bir yol kavşağı için 4 puan sayıldığını anlayabiliyorum, ancak burada 12 puan var. PostGIS'te bu hesaplamayı yapmanın daha iyi bir yolu var mı?

Gruplandırırsanız, yalnızca benzersiz puanlar almalısınız.

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)
;

Bu beklediğinizden biraz daha hileli. Çünkü çiftlerden daha fazla ilişkiyi analiz etmenin iyi bir yolu yoktur. Bir işleve üç satır koyamaz ve hepsinin kesişip kesişmediğini soramazsınız.

Ancak, en azından bir yaklaşım önce geçişleri bulmak, daha sonra her geçişte kaç yolun temas ettiğini kontrol etmek olabilir (hepsi aynı sorguda yapılabilir).

Yollarınız birbirine mükemmel bir şekilde bağlanıyorsa ve bir geçişten geçen yol yoksa, böyle bir şey yapabilirsiniz (test edilmedi):
unutulmuş grup maddesiyle düzenlenmiş (hala test edilmemiştir):

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;

Yollar doğru şekilde bağlanmazsa ve / veya bazı yollar bir geçişten geçerse daha karmaşıktır.

HTH

Nicklas

 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 çizgisi (id 1) B çizgisini (id 2) geçerse ihtiyacımız olan bir kavşak noktasıdır. Ancak B çizgisi de A çizgisini aynı noktada keser. Ancak bu noktaya iki kez ihtiyacımız yok. Bu yüzden kullanmak a.gid < b.gid yerinea.gid != b.gid