如何在交叉路口将OS​​M道路分割为单个路段？

我想使用OpenStreetMap数据创建一个用于pgRouting的道路网络。我从GeoFabrik将shapefile加载到Postgres表中（启用了PostGIS）。但是，我遇到的一个问题是道路并不总是在十字路口处终止，因此我决定在每个十字路口或十字路口将其全部拆分。

要找出所有地方的道路交叉或相交我用下面的interesections SQL（类似于以前的问题）：

CREATE TABLE split_points as
SELECT DISTINCT    
   ST_GeometryN(ST_Intersection(a.geom, b.geom),1) as geom      
FROM
   roads as a,
   roads as b
WHERE
    ST_Touches(a.geom, b.geom)
OR
    ST_Crosses(a.geom, b.geom)    
    AND a.gid != b.gid
GROUP BY
   ST_Intersection(a.geom, b.geom);

我现在想利用这些点来分路。我使用以下方法：

CREATE TABLE split_roads as
SELECT     
    ST_GeomFromEWKB((ST_Dump(ST_Split(g.geom, blade.geom))).geom) As geom,
    generate_series(1,ST_NumGeometries((ST_Split(g.geom, blade.geom)))) as gid
FROM    
    split_points as blade,
    roads as g
WHERE
    ST_Intersects(g.geom, blade.geom);

这种分割方法的问题在于，除了所有分割件之外，还保留了完整的道路长度。为了删除其中包含的这些未拆分的道路几何图形，我使用了该ST_Equals()函数来识别它们并删除它们：

DELETE FROM split_roads USING roads
WHERE ST_Equals(split_roads.geom, roads.geom)

但是，此方法不会删除所有原始的未拆分几何（尽管它确实删除了其中的一些）。是否有更好的删除方法（或整体删除方法），以便表中只有拆分的几何图形？

不是真正解决您的问题的方法，但是请尝试osm2po ...它会创建完美的SQL代码以进行pgrouting中的路由：http ://osm2po.de/

简单的答案：不要。你不应该那样做。

从OSM道路Shapefile中，无法区分交叉口和高架/地下通道。如果分割所有看似交叉的道路，则将创建实际上不存在的交叉路口。

如果您不想使用现有的工具，例如osm2pgrouting（网络足够小）或osm2po，则需要动手处理原始OSM文件。

关于使用pgRouting的一般问题：我认为@ Uffer，@ GisStudent和其他正在展示如何使用“ OSC＆etc.”的人，它们是正确的。遵循“最佳做法”和“标准”的线索...

关于您的问题：“将道路分割成交叉路口的各个路段”或“如何删除所有原始的未分割几何图形”。如果您在这里逐步显示结果，我们可以为您提供帮助。

第一步：拓扑分析

 CREATE TABLE split_points_topo as
  SELECT     
    a.gid as gid_a, b.gid  as gid_b, ST_Relation(a.geom, b.geom) as DE9IM_code
  FROM
    roads as a,
    roads as b
  WHERE a.gid != b.gid AND a.geom && b.geom;

 SELECT DISTINCT st_geometryType(geom) FROM roads;
 SELECT DISTINCT DE9IM_code FROM split_points_topo;
 -- list here the results o these two queries!  ... after we can continue.

另一个“不是您问题的真正解决方案”，但是我们的OSM转换器在从OSM转换为SHP的同时在交叉点分裂。这样可以更高效，因为它可以比较节点的ID，而不是进行几何计算。

一种算法上解决问题的方法是将每条道路的起点和终点添加到“路口”集合中，以便可以确定每个路段都在两个路口之间。