end0tknr's kipple - web写経開発

太宰府天満宮の狛犬って、妙にカワイイ

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

「更に」先程の上記entryの続き & 前回2017年の上記entryの2022年版です。

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

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

2017年と比較すると、画面構成や zipファイルの名称が変更されていますが、上記urlより、行政区域面 (東京 N03-20210101_13_GML.zip)をダウンロード。

$ unzip N03-20210101_13_GML.zip
$ ls -lh N03-20210101_13_GML
total 47M
-rw-rw-r-- 1 end0tknr end0tknr  12K Mar  9  2021 KS-META-N03-21_13_210101.xml
-rw-rw-r-- 1 end0tknr end0tknr 459K Feb 25  2021 N03-21_13_210101.dbf
-rw-rw-r-- 1 end0tknr end0tknr  21M Mar  3  2021 N03-21_13_210101.geojson
-rw-rw-r-- 1 end0tknr end0tknr  145 Feb 25  2021 N03-21_13_210101.prj
-rw-rw-r-- 1 end0tknr end0tknr 7.3M Feb 25  2021 N03-21_13_210101.shp
-rw-rw-r-- 1 end0tknr end0tknr  49K Feb 25  2021 N03-21_13_210101.shx
-rw-rw-r-- 1 end0tknr end0tknr  18M Mar  3  2021 N03-21_13_210101.xml

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

$ pwd
/home/end0tknr/tmp/GIS/N03-20210101_13_GML
$ /usr/local/pgsql/bin/shp2pgsql -p \
  ./N03-21_13_210101 gyosei_kuiki > ./create_gyosei_kuiki.sql
$ /usr/local/pgsql/bin/shp2pgsql -W cp932 -a \
  ./N03-21_13_210101 gyosei_kuiki > ./insert_gyosei_kuiki.sql

$ /usr/local/pgsql/bin/psql -U postgres gis_test < ./create_gyosei_kuiki.sql
$ /usr/local/pgsql/bin/psql -U postgres gis_test < ./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

↑このように実行すると、以下のようなpngが作成されます。

f:id:end0tknr:20220110153731p:plain