Qgis:福島原発避難指示区域の状況をポリゴン表示にする

はじめに

福島原発避難指示区域を地図に落として調査する仕事が入ったので、一般的に出回っている情報を使用してgisにおける福島原発避難指示区域をポリゴン化してみた
要するに以下のような地図を作成して、調査先が最寄りの避難指示区域よりどれだけ離れているのかをGISで調べるためにポリゴン情報を作成してみた

情報の入手

【難指示区域の状況】
経済産業省のページに公示情報（避難区域情報）と解除情報が入手出来るのでPDFだけど入手する
【丁町目ポリゴンの入手】
ポリゴン情報が必要になるのだが、最小の地区情報としては丁町目ポリゴンしか入手出来ないので独立行政法人統計センターが行っている政府統計の総合窓口から福島県の丁町目ポリゴンを入手します

丁町目ポリゴンの入手

政府統計の総合窓口より
地図で見る統計（統計GIS）を選択します

次に「データダウンロード」を選択します

統計表検索（ダウンロード用）に移動したら、調査名に「国勢調査」、「平成22年国勢調査（小地域）2010/10/01」を選択します。そして、「男女別人口総数及び世帯総数」にチェックを入れて「統計表各種データダウンロードへ」を押下します。

都道府県に「福島県」を選択し、市区町村（複数選択可）とあるので、福島県域全て選択します。実際に使うには避難区域の市区町ですが、全て選択します。検索を押下すると、統計データと境界データが表示されますので境界データ（(GIS)で利用するための境界データ及び定義書）をダウンロードします。
境界データの保存形式には以下の種類がダウンロード出来ますが、今回は世界測地系平面直角座標系・Shape形式をダウンロードします

日本測地系平面直角座標系・Shape形式
世界測地系平面直角座標系・Shape形式
世界測地系平面直角座標系・GML形式
日本測地系緯度経度・Shape形式
世界測地系緯度経度・Shape形式

情報の整形

ダウンロードしたファイルは複数に分かれていますのでqgis機能の「複数のshapeファイルを一つに統合する」機能を使って、１ファイルに統合します。
また、CRSが異なっていますので自分が主に使用しているCRSに変換して保存します。そうしてqgis上に福島県の丁町目ポリゴンが作られることになります

shapeファイルと避難区域情報の突合

ここからが結構面倒くさいです。
やり方としては、公示に書かれている避難指示と町名を福島県丁町目ポリゴンから見つけ出して、どのような避難指示になっているのかを見つけなければラなりません。
具体的には・・・
公示に帰還困難区域：富岡町桜１丁目の全ての区域（2016.08.24現在）と有った場合
shapeファイルから福島県双葉郡富岡町桜一丁目を探しだし、該当するKEYCODE1（桜一丁目の場合、543013001）を抽出します
KEYCODE1というのは今回ダウンロードしたshapeファイルの属性テーブルに書かれている丁町目を一意に抽出できる数字ですので、この数字の一覧を公示データから該当するshapeファイル内の属性テーブルから見つけ出します。
住所の数字がアラビアだったり、漢数字だったり、結構面倒くさくて大変です。

自分が突合した避難区域情報

自分がなんとか突合した避難区域情報情報は以下の通りです

通しNo.	県	市	郡	KEYCODE1	住所
2015/12/01 12時台	宮崎県宮崎市島之内付近 (1)	iPhone 6s	LTE	33Mbps	7Mbps
2015/12/01 13時台	宮崎県宮崎市島之内付近 (2)	iPhone 6s	LTE	67Mbps	12Mbps
2015/12/01 13時台	宮崎県宮崎市島之内付近 (3)	iPhone 6s	LTE	79Mbps	10Mbps
2015/12/01 13時台	宮崎県宮崎市島之内付近 (4)	iPhone 6s	LTE	48Mbps	10Mbps
2015/12/01 14時台	宮崎県宮崎市島之内付近 (5)	iPhone 6s	LTE	78Mbps	12Mbps
2015/12/02 12時台	宮崎県宮崎市新名爪付近 (1)	iPhone 6s	LTE	60Mbps	12Mbps
2015/12/02 12時台	宮崎県宮崎市新名爪付近 (2)	iPhone 6s	LTE	47Mbps	8Mbps
2015/12/02 13時台	宮崎県宮崎市新名爪付近 (3)	iPhone 6s	LTE	69Mbps	12Mbps
2015/12/02 13時台	宮崎県宮崎市新名爪付近 (4)	iPhone 6s	LTE	64Mbps	12Mbps