使用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是道路ID，但我不明白为什么会有很多点。我可以理解中央道路交叉路口被计算出的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）交叉，则这是我们需要的交点。但是线B在同一点也越过线A。但是我们不需要两次。这就是为什么我用a.gid < b.gid 而不是a.gid != b.gid