end0tknr's kipple - 新web写経開発

http://d.hatena.ne.jp/end0tknr/ から移転しました

国土交通省にあるGISデータをPostGISへインポート(2017年版)

先程のエントリの続きです。 postgisの動作確認の為、gisデータをimportし、画像ファイルとして表示します。

mac osx 10.6.3にpostgis, mapserverをinstall - end0tknr's kipple - 新web写経開発

国土交通省にあるGISデータをPostGISへインポート(改訂版) - end0tknr's kipple - 新web写経開発

GISデータのダウンロード & 解凍

http://nlftp.mlit.go.jp/ksj/index.html

上記urlより、行政区域面 (東京 N03-170101_13_GML.zip)をダウンロード。

$ unzip N03-170101_13_GML.zip
$ ls -l
-rw-rw-r-- 1    16232 Mar  4  2017 KS-META-N03-17_13_170101.xml
-rw-r--r-- 1  7497800 Sep  9  2017 N03-170101_13_GML.zip
-rw-rw-r-- 1   472154 Feb 17  2017 N03-17_13_170101.dbf
-rw-rw-r-- 1      145 Feb 17  2017 N03-17_13_170101.prj
-rw-rw-r-- 1  7486656 Feb 17  2017 N03-17_13_170101.shp
-rw-rw-r-- 1    49780 Feb 17  2017 N03-17_13_170101.shx
-rw-rw-r-- 1 18624682 Feb 17  2017 N03-17_13_170101.xml

前回行った際は、 shape形式への変換する必要がありましたが、 既に変換済のデータも提供されるようになったようです。

国土交通省にあるGISデータをPostGISへインポート(改訂版) - end0tknr's kipple - 新web写経開発

shapeファイルから、create table文 , insert文作成、更にimport

$ /usr/local/pgsql/bin/shp2pgsql -p \
  ~/tmp/gis/N03-17_13_170101 gyosei_kuiki > ~/tmp/gis/create_gyosei_kuiki.sql
$ /usr/local/pgsql/bin/shp2pgsql -W cp932 -a \
  ~/tmp/gis/N03-17_13_170101 gyosei_kuiki > ~/tmp/gis/insert_gyosei_kuiki.sql

$ /usr/local/pgsql/bin/psql -U postgres gis_test < ~/tmp/gis/create_gyosei_kuiki.sql
$ /usr/local/pgsql/bin/psql -U postgres gis_test < ~/tmp/gis/insert_gyosei_kuiki.sql

mapファイルの作成と、画像への出力

$ vi gyosei_kuiki.map

MAP
    SIZE 800 800           #画像size
    EXTENT 138.8 35.5 140 36   #出力範囲の座標
    STATUS ON              #地図を表示するか
    UNITS DD               #地図の単位(DD は緯度経度)
    IMAGECOLOR 255 255 255 #背景色R G B
    IMAGETYPE PNG          #地図画像を保存する形式
    LAYER
         NAME "gyosei"
         CONNECTIONTYPE POSTGIS
         CONNECTION "user=postgres password='' dbname=gis_test host=localhost"
         DATA "geom FROM gyosei_kuiki" #select文
         TYPE LINE
         STATUS ON
         CLASS
             COLOR 0 0 0
         END
    END
END

$ /usr/local/bin/shp2img -m gyosei_kuiki.map -o gyosei_kuiki.png

上記のmapファイル作成し、shp2img を実行すると、次のように出力されます。 f:id:end0tknr:20170909131918p:plain