先程のエントリの続きです。 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 を実行すると、次のように出力されます。