Individuazione degli incroci stradali tramite PostGIS

Sto cercando di identificare dove si intersecano le strade e di stabilire un punto a questo incrocio, con il numero di strade che formano l'incrocio elencato.

Mi chiedevo se ci fosse un modo per usare ST_NumPoints per raggiungere questo obiettivo, ma non riesco proprio a capire cosa dovrei fare. Quello che ho fatto è creare una tabella di punti in cui le linee si intersecano usando il seguente codice:

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

Se eseguo questo su un campione di strade ottengo la seguente griglia di punti (le strade sono mostrate a scopo illustrativo):

Se ispeziono uno dei punti, vedo che ci sono molti punti uno sopra l'altro:

Il GID qui è l'ID della strada, ma non capisco perché ci siano alcuni punti. Riesco a capire che 4 punti vengono contati per un incrocio stradale centrale, ma ci sono 12 punti elencati qui. Esiste un modo migliore per eseguire questo calcolo in PostGIS?

Se raggruppi, dovresti ottenere solo punti unici.

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

Questo è un po 'più complicato di quanto ci si possa aspettare. Questo perché non esiste un buon modo per analizzare le relazioni per più di coppie. Non è possibile inserire tre linee in una funzione e chiedere se si intersecano tutte.

Ma almeno un approccio potrebbe essere quello di trovare prima gli incroci, quindi controllare quante strade toccano ad ogni incrocio (tutto può essere fatto nella stessa query).

Se le tue strade si collegano perfettamente tra loro e non ci sono strade che passano da un incrocio, allora potresti fare qualcosa del genere (non testato):
modificato con clausola di gruppo dimenticata (ancora non testata):

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;

Se le strade non sono collegate correttamente e / o alcune strade attraversano un incrocio, è più complicato.

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

Se la linea A (id 1) attraversa la linea B (id 2) è un punto di intersezione di cui abbiamo bisogno. Ma la linea B incrocia anche la linea A nello stesso punto. Ma non abbiamo bisogno di questo punto due volte. Ecco perché sto usando a.gid < b.gid invece dia.gid != b.gid