OSMデータから道路のジオメトリを抽出する最良の方法は何ですか？

OSMデータから道路のジオメトリを抽出する最良の方法は何ですか？目標は、天候のために閉鎖されている高速道路のセグメントを表示することです。気象データは、主要都市からの位置（道路に沿って測定されたKM単位）、および影響を受ける道路の長さで構成されます。主な問題は、これらの距離に基づいて道路を切り取ることです。その後のレンダリングは簡単に思えます。

osm2pgsqlやimposm（.osm抽出のみの場合）またはshp2pgsql（シェープファイルがある場合）などを使用して、確実にそのデータをPostgisに取り込みたいと思います。データがPostgisに保存されると、ここで説明されている線形参照関数の一部を使用できます。ST_Line_Interpolate_Point（）およびST_Line_Locate_Point（）：

http://postgis.net/docs/reference.html#Linear_Referencing

線形参照のエラーが大きな問題になると思います。距離がkm単位であり、OSM道路が気象サービスが使用するものとは異なるジオメトリである場合、先に行くほどエラーが拡大されます。

最も簡単な方法は、http: //download.geofabrik.de/osm/からSHAPEファイルをダウンロードすることです。「道路」レイヤーが含まれています（これらのファイルは夜間に更新されます）。

接続された属性テーブルを使用すると、街路を階層的な方法ですばやく分類できます。

パッケージの別のレイヤーは（一部の都市では）エクステントです。

したがって、後続のクリッピングは、たとえばQGISまたはGRASS GISを使用して実行できます。