การระบุทางแยกถนนโดยใช้ 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

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 (id 1) ข้ามเส้น B (id 2) มันเป็นจุดตัดที่เราต้องการ แต่บรรทัด B ก็ตัดกันบรรทัด A ในจุดเดียวกัน แต่เราไม่ต้องการจุดนี้สองครั้ง นั่นเป็นเหตุผลที่ฉันใช้a.gid < b.gid แทนa.gid != b.gid