OSM道路を交差点で個々のセグメントに分割する方法は？

OpenStreetMapデータを使用してpgRoutingで使用する道路ネットワークを作成したいと思います。GeoFabrikからシェープファイルをPostgresテーブルにロードしました（PostGISを有効にした状態）。しかし、私が抱えていた問題の1つは、道路が常に交差点で終わらないということでした。そのため、すべての交差点または交差点で道路をすべて分割することにしました。

道路が交差または交差するすべての交差点を特定するために、私は次のものを使用しました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を試してください... pgroutingでルーティングするための完璧なSQLコードを作成します：http ://osm2po.de/

簡単な答え：しないでください。そのようにすべきではありません。

OSM道路シェープファイルからは、交差点と高架/地下道を区別することはできません。交差しているように見える道路をすべて分割すると、現実には存在しない交差点が作成されます。

osm2pgrouting（ネットワークが十分に小さい場合）やosm2poなどの既存のツールを使用しない場合は、元のOSMファイルで手を汚す必要があります。

pgRoutingを使用した一般的な問題について：@OSCなどを使用する方法を示している@ Uffer、@ GisStudentなどが正しいと思います。「ベストプラクティス」と「標準」の手がかりに従ってください...

質問について：「交差点で道路を個々のセグメントに分割する」または「元の分割されていないジオメトリをすべて削除する方法」ここに結果をステップバイステップで表示すると、私がお手伝いできます...

最初のステップ：トポロジー分析

 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を比較できるため、この方法の方が効率的です。

アルゴリズムで解決する1つの方法は、道路全体の始点と終点を「交差点」のセットに追加して、すべてのセグメントが2つの交差点の間にあることを確認することです。