Xác định các nút giao thông đường bộ bằng PostGIS

Tôi đang cố gắng xác định nơi các con đường giao nhau và để tạo một điểm tại giao lộ này, với số lượng đường tạo thành giao lộ được liệt kê.

Tôi đã tự hỏi nếu có một số cách sử dụng ST_NumPoints để đạt được điều này nhưng tôi không thể hiểu được những gì tôi nên làm. Những gì tôi đã làm là tạo một bảng các điểm trong đó các đường giao nhau bằng cách sử dụng mã sau đây:

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

Nếu tôi chạy cái này trên một mẫu đường, tôi sẽ nhận được lưới các điểm sau (các đường được hiển thị để minh họa):

Nếu tôi kiểm tra một trong những điểm, tôi thấy rằng có nhiều điểm xếp chồng lên nhau:

GID ở đây là ID đường, nhưng tôi không hiểu tại sao có nhiều điểm. Tôi có thể hiểu 4 điểm được tính cho một ngã tư đường trung tâm, nhưng có 12 điểm được liệt kê ở đây. Có cách nào tốt hơn để thực hiện tính toán này trong PostGIS không?

Nếu bạn nhóm, bạn sẽ chỉ nhận được điểm duy nhất.

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

Đây là một chút khó khăn hơn bạn có thể mong đợi. Đó là bởi vì không có cách nào tốt trong việc phân tích các mối quan hệ cho nhiều hơn các cặp. Bạn không thể đặt ba dòng vào một hàm và hỏi xem tất cả chúng có giao nhau không.

Nhưng, ít nhất một cách tiếp cận có thể là trước tiên tìm đường giao nhau, sau đó kiểm tra xem có bao nhiêu con đường chạm vào mỗi đường giao nhau (tất cả đều có thể được thực hiện trong cùng một truy vấn).

Nếu các đường của bạn kết nối hoàn hảo với nhau và không có đường nào đi qua ngã tư, thì bạn có thể làm một cái gì đó như thế này (không được kiểm tra): được
chỉnh sửa với mệnh đề nhóm bị lãng quên (vẫn chưa được kiểm tra):

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;

Nếu đường không được kết nối đúng cách và / hoặc một số đường đi qua ngã tư thì sẽ phức tạp hơn.

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

Nếu dòng A (id 1) vượt qua dòng B (id 2) thì đó là điểm giao nhau mà chúng ta cần. Nhưng dòng B cũng vượt qua dòng A trong cùng một điểm. Nhưng chúng ta không cần điểm này hai lần. Đó là lý do tại sao tôi đang sử dụng a.gid < b.gid thay vìa.gid != b.gid