とらりもんHOME  Index  Search  Changes  Login

とらりもん - 国土数値情報 Diff

  • Added parts are displayed like this.
  • Deleted parts are displayed like this.

! 土地利用細分メッシュ
[[https://nlftp.mlit.go.jp/ksj/gml/datalist/KsjTmplt-L03-b.html]]
* 昭和51年度、昭和62年度、平成3年度、平成9年度、平成18年度、平成21年度、平成26年度、平成28年度
* ベクター (shape型式)で提供されている(平成26年度版は別途ラスター型式でも提供あり)。
* 日本測地系(旧測地系; tokyo datum)と世界測地系(jgd2000)の両方で提供されている。
(ベクター版)
* ここでは平成28年度(2016年度)版の処理手順を紹介する。Linuxで。
* まずダウンロードと解凍:
for i in \
3036 3622 3623 3624 3631 3641 3653 3724 3725 3741 3823 3824 \
3831 3841 3926 3927 3928 3942 4027 4028 4040 4042 4128 4129 \
4142 4229 4230 4328 4329 4429 4440 4529 4530 4531 4540 4629 \
4630 4631 4728 4729 4730 4731 4739 4740 4828 4829 4830 4831 \
4839 4928 4929 4930 4931 4932 4933 4934 4939 5029 5030 5031 \
5032 5033 5034 5035 5036 5038 5039 5129 5130 5131 5132 5133 \
5134 5135 5136 5137 5138 5139 5229 5231 5232 5233 5234 5235 \
5236 5237 5238 5239 5240 5332 5333 5334 5335 5336 5337 5338 \
5339 5340 5432 5433 5435 5436 5437 5438 5439 5440 5531 5536 \
5537 5538 5539 5540 5541 5636 5637 5638 5639 5640 5641 5738 \
5739 5740 5741 5839 5840 5841 5939 5940 5941 5942 6039 6040 \
6041 6139 6140 6141 6239 6240 6241 6243 6339 6340 6341 6342 \
6343 6439 6440 6441 6442 6443 6444 6445 6540 6541 6542 6543 \
6544 6545 6546 6641 6642 6643 6644 6645 6646 6647 6741 6742 \
6747 6748 6840 6841 6842 6847 6848; do
wget https://nlftp.mlit.go.jp/ksj/gml/data/L03-b/L03-b-16/L03-b-16_${i}-jgd_GML.zip
done
for i in L03*zip; do unar $i; done
* 各zipファイルに対応するディレクトリができ, その中にベクターファイル(.shp)が入っている。
* 試しに覗いてみるとこんなかんじ:
ogrinfo -summary -al L03-b-16_6342-jgd_GML/L03-b-16_6342.shp

INFO: Open of `L03-b-16_6342-jgd_GML/L03-b-16_6342.shp'
      using driver `ESRI Shapefile' successful.

Layer name: L03-b-16_6342
Metadata:
  DBF_DATE_LAST_UPDATE=2020-02-15
Geometry: Polygon
Feature Count: 430000
Extent: (142.000000, 42.083333) - (143.000000, 42.666667)
Layer SRS WKT:
GEOGCS["GCS_JGD_2000",
    DATUM["Japanese_Geodetic_Datum_2000",
        SPHEROID["GRS_1980",6378137.0,298.257222101]],
    PRIMEM["Greenwich",0.0],
    UNIT["Degree",0.0174532925199433]]
メッシュ: String (10.0)
土地利用種: String (4.0)
撮影年月日: String (8.0)
* 多くの場合, この最後の3行が文字化けしているだろう。文字コードが合っていないのである(このデータの文字コードはSHIFT-JISであり, 最近の多くのLinuxの文字コードはUTF-8である)。
* そこで, GDAL/OGRライブラリの文字コードをこのデータに合わせてやろう。そのためには環境変数を設定する:
export SHAPE_ENCODING="SHIFT_JIS"
* この設定が後で効いてくる。shapeファイルをラスター化するときに, 属性名の指定で効いてくる。
* そのあとに上のコマンドを打てば文字化けは消えるはずだ。
* これをラスターファイルに変換するには:
gdal_rasterize -a 土地利用種 -tr 0.00041666666 0.00041666666 L03-b-16_6342-jgd_GML/L03-b-16_6342.shp L03-b-16_6342.tif

以下は, 地域ごとにタイルをまとめて(.vrtというのはGDALヴァーチャルフォーマットというやつ), ラスターに変換するコマンド
fnvrt="Yaeyama.vrt"
for i in 3622 3623 3624 3724 3725 3823 3824; do
ls *${i}*/*.shp; done > dummy.txt
ogrmerge.py -single -overwrite_ds -f VRT -o $fnvrt `cat dummy.txt`

fnvrt="Daito.vrt"
for i in 3831 3631; do
ls *${i}*/*.shp; done > dummy.txt
ogrmerge.py -single -overwrite_ds -f VRT -o $fnvrt `cat dummy.txt`

fnvrt="Ogasawara.vrt"
for i in 3641 3741 3841 3942 4040 4042 4142 4440; do
ls *${i}*/*.shp; done > dummy.txt
ogrmerge.py -single -overwrite_ds -f VRT -o $fnvrt `cat dummy.txt`

fnvrt="Okinawa.vrt"
for i in 3926 3927 3928 4027 4028 4128 4129 4229 4230 4328 4329 4429; do
ls *${i}*/*.shp; done > dummy.txt
ogrmerge.py -single -overwrite_ds -f VRT -o $fnvrt `cat dummy.txt`

fnvrt="Kyushu.vrt"
for i in 4529 4530 4531 4629 4630 4631 4728 4729 4730 4731 4828 4829 4830 4831 4928 4929 4930 4931 5029 5030 5031 5129 5229 5130 5131; do
ls *${i}*/*.shp; done > dummy.txt
ogrmerge.py -single -overwrite_ds -f VRT -o $fnvrt `cat dummy.txt`

fnvrt="Izu.vrt"
for i in 4739 4740 4540 4839 4939 5038 5039 5139; do
ls *${i}*/*.shp; done > dummy.txt
ogrmerge.py -single -overwrite_ds -f VRT -o $fnvrt `cat dummy.txt`

fnvrt="ShikokuChugoku.vrt"
for i in 4932 4933 4934 5032 5033 5034 5132 5133 5134 5231 5232 5233 5234 5332 5333 5334 5432 5433 5531 5130 5131
do ls *${i}*/*.shp; done > dummy.txt
ogrmerge.py -single -overwrite_ds -f VRT -o $fnvrt `cat dummy.txt`

fnvrt="KinkiChubu.vrt"
for i in 5035 5036 5135 5136 5137 5138 5235 5236 5237 5238 5335 5336 5337 5338 5435 5436 \
          5437 5438 5536 5537 5538 5636 5637 5638 5738 5134 5234 5334
do ls *${i}*/*.shp; done > dummy.txt
ogrmerge.py -single -overwrite_ds -f VRT -o $fnvrt `cat dummy.txt`

fnvrt="KantoHokuriku.vrt"
for i in 5239 5240 5339 5340 5439 5440 5539 5540 5541 5639 5640 5641 \
          5238 5338 5438 5538 5638 5738 5739 5740 5741
do ls *${i}*/*.shp; done > dummy.txt
ogrmerge.py -single -overwrite_ds -f VRT -o $fnvrt `cat dummy.txt`

fnvrt="Tohoku.vrt"
for i in 5839 5840 5841 5939 5940 5941 5942 6039 6040 6041 6139 6140 6141 5738 5739 5740 5741 6240 6241
do ls *${i}*/*.shp; done > dummy.txt
ogrmerge.py -single -overwrite_ds -f VRT -o $fnvrt `cat dummy.txt`

fnvrt="Hokkaido.vrt"
for i in 6239 6240 6241 6243 6339 6340 6341 6342 6343 6439 \
6440 6441 6442 6443 6444 6445 6540 6541 6542 6543 \
6544 6545 6546 6641 6642 6643 6644 6645 6646 6647 \
6741 6742 6747 6748 6840 6841 6842 6847 6848 \
do ls *${i}*/*.shp; done > dummy.txt
ogrmerge.py -single -overwrite_ds -f VRT -o $fnvrt `cat dummy.txt`

# 南鳥島 3653は除外
# 沖ノ鳥島 3036は除外

export SHAPE_ENCODING="SHIFT_JIS"
for i in *.vrt; do
  gdal_rasterize -a HANREI_C -co "COMPRESS=LZW" -tr 0.00041666666 0.00041666666 $i ${i%.vrt}_50m.tif
done


! 土地利用細分メッシュ(ラスタ版)
* [[https://nlftp.mlit.go.jp/ksj/gml/datalist/KsjTmplt-L03-b_r.html]]
* ラスタ版は扱いが楽なのでありがたい。
* しかし, 古いバージョン(平成26年度 = 2014年度)しか公開されていない。最新版(平成28年版)は別途ベクターで提供されている(2022/03/04現在)
* "全国の土地利用の状況について、3次メッシュ1/10細分区画(100mメッシュ)毎に、各利用区分(田、その他の農用地、森林、荒地、建物用地、幹線交通用地、湖沼、河川等)を整備したもの"
* 全体は東経122-154度, 北緯20-46.0度。それが175枚のタイルに分かれている。
* 以下のコマンド(Linux)で一括ダウンロードできる:
for i in \
3036 3622 3623 3624 3631 3641 3653 3724 3725 3741 3823 3824 \
3831 3841 3926 3927 3928 3942 4027 4028 4040 4042 4128 4129 \
4142 4229 4230 4328 4329 4429 4440 4529 4530 4531 4540 4629 \
4630 4631 4728 4729 4730 4731 4739 4740 4828 4829 4830 4831 \
4839 4928 4929 4930 4931 4932 4933 4934 4939 5029 5030 5031 \
5032 5033 5034 5035 5036 5038 5039 5129 5130 5131 5132 5133 \
5134 5135 5136 5137 5138 5139 5229 5231 5232 5233 5234 5235 \
5236 5237 5238 5239 5240 5332 5333 5334 5335 5336 5337 5338 \
5339 5340 5432 5433 5435 5436 5437 5438 5439 5440 5531 5536 \
5537 5538 5539 5540 5541 5636 5637 5638 5639 5640 5641 5738 \
5739 5740 5741 5839 5840 5841 5939 5940 5941 5942 6039 6040 \
6041 6139 6140 6141 6239 6240 6241 6243 6339 6340 6341 6342 \
6343 6439 6440 6441 6442 6443 6444 6445 6540 6541 6542 6543 \
6544 6545 6546 6641 6642 6643 6644 6645 6646 6647 6741 6742 \
6747 6748 6840 6841 6842 6847 6848; do
wget https://nlftp.mlit.go.jp/ksj/gml/data/L03-b_r/L03-b_r-14/L03-b-14_${i}.zip
done
* これら↑をこのコマンド↓で一気に解凍してモザイクする:
for i in L03*zip; do unar $i; done
gdal_merge.py L*/*.tif -o L03-b-14_merge.tif -ul_lr 122 46 154 20 -n 0 -a_nodata 0 -pct
# -pctはカラーテーブルを最初のファイルから流用する。
# -nはオリジナルのno data値
# -a_nodataは出来上がるファイルのno data値

* 凡例は以下の通り(zip解凍した中にhtmファイルで入ってる。全国共通)
コード 種別 定義
10 田 湿田・乾田・沼田・蓮田及び田とする。
20 その他の農用地 麦・陸稲・野菜・草地・芝地・りんご・梨・桃・ブドウ・茶・桐・はぜ・こうぞ・しゅろ等を栽培する土地とする。
50 森林 多年生植物の密生している地域とする。
60 荒地 しの地・荒地・がけ・岩・万年雪・湿地・採鉱地等で旧土地利用データが荒地であるところとする。
70 建物用地 住宅地・市街地等で建物が密集しているところとする。
91 道路 道路などで、面的に捉えられるものとする。
92 鉄道 鉄道・操車場などで、面的にとらえられるものとする。
100 その他の用地 運動競技場、空港、競馬場・野球場・学校港湾地区・人工造成地の空地等とする。
110 河川地及び湖沼 人工湖・自然湖・池・養魚場等で平水時に常に水を湛えているところ及び河川・河川区域の河川敷とする。
140 海浜 海岸に接する砂、れき、岩の区域とする。
150 海水域 隠顕岩、干潟、シーパースも海に含める。
160 ゴルフ場 ゴルフ場のゴルフコースの集まっている部分のフェアウエイ及びラフの外側と森林の境目を境界とする。
0 解析範囲外
* モザイクがおわったら表示してみよう:
qgis L03-b-14_merge.tif

!国土数値情報データ変換ツール
(2015/02/04 記事作成)

こいつをUbuntu Linux 14.04で動かす。[[参考|http://blogs.yahoo.co.jp/igproj_fusion/17430504.html]] [[参考|http://azemichi2000.blogspot.jp/2011/01/linuxxmlshape.html]]

mkdir tmp
cd tmp
wget http://nlftp.mlit.go.jp/ksj/jpgis/ksjtool_32_v1.8.zip
unzip ksjtool_32_v1.8.zip
cd JRE
wine jre-6u26-windows-i586-s.exe
cd ..
msiexec -i KsjToolInstaller.msi
cd ~/.wine/


メッセージが文字化けするときは? ... Windows用のJAVAの日本語フォントがうまく設定されていない。対処法:

cd ~/.wine/drive_c/Program\ Files/Java/jre6/lib/fonts/
mkdir fallback
cd fallback
sudo ln -s /usr/share/fonts/truetype/takao-gothic/TakaoPGothic.ttf ./

[[参考|https://cosmo0920.wordpress.com/2011/01/29/wine%E4%B8%8A%E3%81%A7java%E3%81%AEexe%E3%82%92%E5%8B%95%E3%81%8B%E3%81%99%E3%82%82%E3%81%AF%E3%82%84%E3%83%8D%E3%82%BF/]]

うまくいかない...
N03-14_08_010401.xmlを変換しています
[ Input Error ] 変換対象の地物が見つかりません