温暖化に伴い洪水被害が増加しつつある．各自治体のハザードマップの重要性が増しているところであるが，今回は国土数値情報ダウンロードサイトから洪水浸水想定区域のデータをダウンロードし，QGISで表現してみた．
洪水浸水想定区域のデータをダウンロード
データはここにある．ブラウザで「洪水浸水想定区域」と検索する．メッシュとポリゴンが検索されるが，今回はポリゴンを使う．
スクロールしていくと茨城県のファイルがあるのでアンケートに答えてダウンロードする．
QGISに読み込む
「レイヤ」「レイヤを追加」「ベクタレイヤを追加」と進んでいく．
複数のシェープファイルをまとめて追加できる．計画規模，想定最大規模，浸水継続時間の3種類があるが，ここでは想定最大規模を使う．
この状態で「追加」をクリックする．
下図のように17個のファイルがレイヤとして追加される．属性テーブルを開いてみるとわかるが，一つのレイヤには数百から数千ものポリゴンが格納されている．浸水深ランクコードごとにまとめてしまいたい．そんな時は融合を使う．
複数のポリゴンを浸水深コードごとに融合する
「ベクタ」「空間演算ツール」「融合(dissolve)」と進む．
融合(dissolve)ダイアログの初期画面．入力レイヤがサービス付き高齢者住宅となっているが，これは今回追加したポリゴンのレイヤに変更していく．
複数のレイヤをまとめて処理するにはバッチプロセスだが…
手動でできなくもないが同じ処理はバッチプロセスを使って自動化したい．上図の「バッチプロセスで実行」をクリックすると下図のダイアログになる．
「入力レイヤ」「オートフィル」「レイヤから選択」と進む．
「すべて選択」し不要なチェックを外し，OKをクリックする．
「基準となる属性」を「式による計算」とする．
「地図のレイヤ」から浸水深ランクコードである「A31_205」を選びダブルクリックする．
工具アイコンの「詳細オプション」で「無効地物フィルタ」では「フィルタリングしない」を選択する．これをまとめて設定する方法がない．
OKをクリックするが結果はエラーとなる．シェープファイルのフォルダ名に日本語が使われているためかと考え英語に直したが結果は変わらなかった．
その後の調査で出力したファイルのフルパスを指定すればよいことが分かった．指定方法はいくつかあるが，「自動塗りつぶしモード」で「数値で塗りつぶす」を選択するととりあえずはうまく行くようだ．
その他の方法（「パラメータの値で塗りつぶす」など）は途中でエラーが発生するか，同じファイル名で繰り返し保存してしまい，最後の一個しかファイルが残らなかった．
やむなく手動で17個のレイヤを融合
仕方がないので手動で上記の処理を行う．結果は下図のようになる．それぞれエクスポートしておく．レイヤ名は任意で良い．ここでは2桁の数字の連番とした．
浸水深ランクコードにより色を塗り分ける
それぞれのレイヤを右クリックして「プロパティ」の「シンボロジ」で「カテゴリ値による定義」とし，「値」を浸水深ランクコードである「A31_205」とする．「カラーランプ」から「カラーランプを新規作成」を選ぶ．
「カラーランプ型」を「カタログ(ColorBrewer)」とする．
「色」を6色に変更する．
「分類」をクリックする．
浸水深ランクコードごとに分類される．「シンボル」から「色を選択」する．
「シンボルの色」ダイアログが開く．
- 浸水深ランクコード１はR230, G255, B102
- 浸水深ランクコード２はR153, G255, B255
- 浸水深ランクコード３はR230, G166, B255
- 浸水深ランクコード４以上はR153, G64, B255
とする．
入力の終わった状態である．この色指定は国土交通省の「洪水ハザードマップ作成の手引き（改訂版）」に準拠している．浸水深ランクコード１は１階床下浸水，２は１階床上浸水，３は２階浸水，４以上は２階水没を示している．CMYKからRGBへの変換についてはCMYK⇔RGB変換ツール – げんちゃんの色見本にて行った．浸水深ランクコードについては同じ国土交通省の「水害ハザードマップ作成の手引」(p36)が6段階の色見本を示している．こういった手引については部署間で認識を統一してもらいたい．
結果
結果を示す．サービス付き高齢者住宅の中には洪水浸水想定区域の中に建っているものもある．こういう情報は自分で調べないと分からない．ハザードマップだけでは自分の関係する建物がどこにあるのか分かりにくい．各自治体の公表している情報も本来ならシェープファイルで公開してほしいものである．
“国土数値情報の洪水浸水想定区域データをQGISで表現する” への1件の返信