Определение пересечения дорог с помощью 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);

Если я запускаю это на примере дорог, я получаю следующую сетку точек (дороги показаны для иллюстрации):

Если я осматриваю одну из точек, я вижу, что есть много точек, наложенных друг на друга:

GID здесь - это идентификатор дороги, но я не понимаю, почему есть много точек. Я могу понять, что 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)
;

Это немного сложнее, чем вы могли ожидать. Это потому, что не существует хорошего способа анализа отношений для более чем пары. Вы не можете поместить три строки в функцию и спросить, пересекаются ли все они.

Но, по крайней мере, один из подходов может заключаться в том, чтобы сначала найти перекрестки, а затем проверить, сколько дорог касается каждого перекрестка (все это можно сделать в одном запросе).

Если ваши дороги идеально соединены друг с другом, и нет никаких дорог, проходящих через перекресток, то вы можете сделать что-то вроде этого (не проверено):
отредактировано с забытым предложением группы (все еще не проверено):

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;

Если дороги не соединены надлежащим образом и / или некоторые дороги проходят через перекресток, это более сложно.

НТН

Никлас

 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