[PostGIS] 数値地図2500の道路情報から頂点情報を抽出する

国土地理院が提供する数値地図2500には、道路線の情報が含まれています。しかし、交差点や各道路のつなぎ目の情報は収録されていません。

そこで、次の SQL 文を実行すると、各道路線のつなぎ目が取得できます。ここで、道路情報は road テーブルの the_geom 列に格納されているとします。

SELECT DISTINCT ON ( ST_AsEWKB( location ) ) location
FROM (
  ( SELECT ST_StartPoint( the_geom ) AS location FROM road ) UNION ALL
  ( SELECT ST_EndPoint( the_geom ) AS location FROM road ) )
  AS vertex_source;

（；^ω^）道路線の両端点を抽出して、重複する点を省いているだけです。

参考文献

• 数値地図2500を MapServer で表示させる 数値地図2500を PostGIS に格納するには、このページを参照してください。