Визначення перехресть доріг за допомогою 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;

Якщо дороги не підключені належним чином та / або деякі дороги проходять через переправу, це складніше.

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), це нам потрібна точка перетину. Але лінія В також перетинає лінію А в тій же точці. Але нам ця точка не потрібна двічі. Ось чому я використовую a.gid < b.gid замість цьогоa.gid != b.gid