QGISによる洪水想定浸水図の作成について

[黒木修身]

QGISの初心者です。

QGISを使って洪水想定氾濫図を作成する方法が分かりません。

数値標高データ（5ｍメッシュ）より、各メッシュの浸水深を別途求めます。

浸水深をレベル毎に色分けして地図に張り付けて浸水想定図を作成します。

この時に、以下の操作方法を御教授いただけたら幸いです。

　・ＱＧＩＳで5ｍメッシュを作成する方法

　・各メッシュに浸水深を設定し、レベル毎に色分けする方法

　・地図に張り付ける方法

よろしくお願いいたします。

[ありた]

こんにちは

グリッド（メッシュ）の生成はベクタメニューまたはプロセシングツールの

「グリッドを作成」でできます。

ただ、その前に、「数値標高データより求めた浸水深データ」をお持ちだと思います。

そのデータの形式等についてご確認いただき、提示してください。

すでにグリッド状ポリゴンの形式であったり、ラスタ形式であると思います。

そのデータをそのまま用いるか、あるいは加工するほうが楽です。

仮に、別途生成したメッシュ（グリッド）に浸水深データを設定（付与）するに

しても、やはりお持ちのデータの形式等の詳細情報が必要になります。

色分けについては、 GIS は同一の情報から目的に応じて色等の表現方法を変化させる

ことを得意としたツールです。レイヤをダブルクリックあるいは右クリックより

「プロパティ」を開き、「シンボロジ」で見た目のルールを設定することができます。

また GIS はパソコン上に仮想の広大な地図が広がっているイメージですが、それを

印刷等を見据えた出力物を作成するにはプロジェクトメニューの「新規印刷レイアウト」を

選択し、レイアウトの設定を行ってください。投影法（メルカトルとかモルワイデとか）や

縮尺とかを含めて調整できます。

縮尺等は気にせず、単に画像形式でほしいというのであれば、プロジェクトメニューの

「インポートとエクスポート」より「地図を画像にエクスポート」を実行してください。

現在表示されている範囲がそのまま画像として出力されます。

[ありた]

> 地図に張り付ける方法

他の地図と重ね合わせるという意味でしたら、他の地図データを用意し QGIS で

読み込み、 GIS 上で重ね合わせてください。

（既に適切な位置情報をもっているデータであれば）正しく描画をしてくれます。

ベクタデータだと

・国土地理院　基盤地図情報

　https://www.gsi.go.jp/kiban/

・NaturalEarth

　https://www.naturalearthdata.com

背景地図データだと XYZ Tiles 形式の

・国土地理院　タイル地図

　https://maps.gsi.go.jp/development/ichiran.html

・OpenStreetMap

　https://wiki.openstreetmap.org/wiki/JA:%E3%82%BF%E3%82%A4%E3%83%AB

などがあります。

[東 雲]

黒木様

氾濫解析自体はシミュレーションプログラムでやって、結果を浸水想定区域図にするか

あるいは既にある氾濫解析結果を図化するということでいいでしょうか。

まず浸水想定区域の5mメッシュと地理院の数値標高モデルの5mメッシュは基本的に別物なので

その辺を把握する必要があるんじゃないかと思います。

「浸水想定区域図データ電子化ガイドライン」で検索すると

浸水想定区域のメッシュの仕様が出てくるので、まずそこからかと。

第3版のp23以降（pdfのでいうとp32以降）にメッシュの仕様が出てるのでご一読を。

ちなみに数値標高モデルは0.2秒の正方メッシュです。

「split polygon」というQGISのプラグインがあるので、それを使って

適当な3次メッシュを分割して作るのが、手間はかかりますが技術的に簡単です。

2021年6月25日金曜日 16:02:33 UTC+9 黒木修身:

[黒木修身]

ありた様、東雲様

ご教授ありがとうございます。

プロセシングツールからグリッド作成まではできました。

浸水深データは、CSVで（No,X,Y,H)のデータです。

このグリッドに浸水深データを設定、浸水深の色分け方法が分かりません。

また、東雲さんのアドバイスにあった[Split Polygon]の使い方を書いたものはありますでしょうか。

当方はQGISの技量がほとんどなく、お手数をおかけいたします。

よろしくお願いいたします。

2021年6月26日土曜日 21:05:54 UTC+9 東 雲:

[ありた]

こんにちは

CSV でレイヤとして追加します。

その際、ジオメトリの CRS （座標系）を適切に設定してください。

（経緯度座標であれば EPSG:6668 、平面直角座標系、 UTM 等であればそれぞれ）

この座標がきっちり 5m 間隔であれば、プロセシングツールの「長方形・楕円・

ダイヤモンド」で、点から幅 5m の四角形を生成する方が楽です。

（経緯度座標等であれば、事前に「レイヤの再投影」で座標系の変換が必要）

別途作成したグリッドに割り当てたい場合は、「属性の空間結合」を行います。

[黒木修身]

ありた様

丁寧な御指導ありがとうございます。

・CSVデータのレイヤ追加

・プロセシングツールで5ｍメッシュを作成しました。

ここから、各メッシュの浸水深ごとの色分けはどうするのでしょうか。

よろしくお願いいたします。

2021年6月27日日曜日 20:58:37 UTC+9 ありた:

[sakuya_izayoi]

izayoiです。

リンクを貼るだけで恐縮ですが

https://misaki-blog.com/190620-qgis-color-grad/

を参考に調整ができると思います。

2021年6月28日月曜日 11:29:58 UTC+9 黒木修身:

[福岡]

黒木 修身 様

取り越し苦労であれば良いのですが、遣り取りをみていて少し不安に思いましたので投稿いたします。

まずCSVデータのレイヤ追加ですが、ポイント座標としてCRSを正しく設定して取り込むことができたでしょうか？

正しく取り込めたかどうかは、WEB地図を背景に敷くことで簡単に確認できます。

そしてCRSを正しく理解して設定しないと、大きく手戻りすることがあるかもしれません。

これらのことについて記事にして、昨年Qiitaに投稿したものがあります。

よろしければ、ご覧ください。

参照座標系のざっくりとしたお話し

https://qiita.com/fukuoka1960/items/2318898e47ec1e749e1c

次に5mメッシュについてです。

CSVに記載されているX,Y座標は、投影座標系で5m間隔のデータでしょうか？

東雲さんの投稿で示されていた「浸水想定区域図データ電子化ガイドライン」に従って作成されたものなら、位置は緯度経度でまとめられているように思いました。

もし緯度経度のデータを取り込んだとするのなら、X方向は0.15秒なので0.225/3600=0.0000625、Y方向は0.15秒なので0.15/3600=0.00004166･･･という数値でプロセッシングツールを使う必要があります。

以上が心配し過ぎで、順調にCSVからメッシュデータを作成できたのであれば、メッシュデータの属性テーブルにCSVの浸水深が取り込まれているはずです。

色分けは、この浸水深のフィールドを利用し、izayoiさんの示されたリンク先の手順を参考すると良いように思います。

「Split Polygon」ですが、プラグインの説明に出てくる「ホームページ」が次にリンクされています。

https://github.com/uilic/SplitPolygon

使い方は細かく説明されていないのですが、試しに使ってみての勘所は次のとおりでした。

１ 一気にメッシュは作ることができないので、0度と90度で2回分割する必要がある。

２ 分割は「Equal parts」を選び、分割数を入力する。

３ 1回0度で分割したら一時レイヤをshapeファイルなどで保存し、保存したものを使って再度分割する。

　 ※一時レイヤから再分割はできないようでした。

４ 角度の入力では分、秒に「0」を入力する。

2021年6月28日月曜日 13:26:34 UTC+9 sakuya_izayoi:

[東 雲]

福岡様

私が不安に思っているところを丁寧に説明していただいてありがとうございます。

黒木様

福岡様の説明に付け足しとして、split polygon で3次メッシュを分割して5mメッシュを作る場合は

3次メッシュ5つくらいずつに分けて作業した方がいいと思います。

作業するPCのスペックにもよるのですが、エラーが起こる原因の一つがこれです。

福岡様の説明通り、縦横それぞれに処理をする必要がありますが

2回目の処理でエラーが起こる場合は、処理数が多すぎる可能性もあります。

また、csvの内容について

・xy座標の整数部分の桁数

・都道府県

をお伝えいただければ、もう少し選択の幅を狭めたアドバイスが可能になるかと思います。

2021年6月28日月曜日 15:31:53 UTC+9 福岡:

[黒木修身]

親切な御助言ありがとうございます。

一応、浸水想定図らしきものは作成できました。

ここで、以下の点をお伺いします。

①浸水深CSVをレイヤ追加し、「長方形・楕円・ダイヤモンド」で5ｍメッシュを作成し、浸水深毎に色分けしました。

　この結果をファイル保存して、再度開いたときポリゴンが消えます。

　浸水深の属性データも入ってません。何か保存方法があるのでしょうか。

②ポリゴンで浸水想定図を描いたとき、ポリゴンの枠線を消す方法はありますか。

③浸水深ごとに色分けの指定があります。色分けをＲＧＢで設定できますか。

お手数をおかけいたします。

よろしくお願いいたします。

2021年6月29日火曜日 6:51:19 UTC+9 東 雲:

[福岡]

黒木 修身 様

こんにちは

まだ不安感が消えませんので、最初にこのことについて書かせてください。

今回行おうとしていることは、多分中級レベルの力が必要なものだと思っています。

このことを踏まえて、質問の内容が入門者が最初に学ぶべきものからそれなりの知識が必要なものまであり、アンバランスさと危うさを感じている次第です。

少なくとも入門者レベルの知識にアクセスし、確認できるような準備は必要だと思います。

また、基礎となるお手持ちのデータの仕様と最終的にどのような成果をお求めかを残念ながら投稿からは読み取ることができませんでした。

例えばQGISで浸水想定図を眺めるだけなら良いのですが、成果として印刷物やPDFでの出力、あるいは電子ファイルの提出が求められるものなら、助言が違ってくるように思います。

さてご質問に関してです。

下図をご覧ください。ポリゴンところに「一時レイヤの作成」と表示されていませんか？

この設定でポリゴンを作成した場合、表示のとおり一時的なものですので消えてしまうのは当然です。

たぶんQGISを閉じるときに、このデータが失われるという警告も出たはずです。

これを回避するためには「一時レイヤ」を保存すると良いです。

レイヤパネルの「一時レイヤ」の上で右クリックして「エクスポート － 地物の保存」でファイルを保存します。

下の「保存されたファイルを地図に追加」のチェックを入れておいて地図に追加後、一時レイヤを削除すると良いでしょう。

属性データに関しては、CSVに整理されているものなら属性テーブルに反映されているはずだと思います。

ポリゴンの枠線ですが、レイヤプロパティの「シンボロジ」タブで設定できます。

ストロークスタイルを「実線」から「ペンなし」にすると良いでしょう。

「シンボロジ」タブの「塗りつぶした色」で着色されているところの上をクリック、あるいはその右端にある三角の部分をクリックして「色の選択」を選ぶと「塗りつぶし色を選択」ダイアログが立ち上がります。

入門者レベルの知識にアクセスし、確認できるような準備として書籍の購入を強くお勧めします。

体系的に情報が整理されていて、索引から必要な事柄を探すことができるので「業務で使うＱＧＩＳ Ver.3　完全使いこなしガイド」を私個人としてはお勧めします。

AmazonでQGISをキーワードにすると最初に出てきます。

大型書店であれば店舗に在庫があるかもしれませんし、Kindle版を購入してPCやタブレットなどで利用しても良いと思います。

最後ですが、東雲さんの問いかけた「csvの内容」についてお答えになってはいかがでしょうか？

問いかけのみの一方通行だと、適切な助言を得られにくいと思います。

できれば助言者からの問いかけについて答えて頂き、会話のキャッチボールがうまくいくようになれば良いと思っています

2021年6月29日火曜日 17:12:15 UTC+9 黒木修身:

[黒木修身]

東雲様

5ｍメッシュ標高は、「基盤地図情報ビュアー」からダウンロードした、下記形式のデータです。

ちなみに、宮崎県を対象にしてます。

浸水深は、（H)を浸水深に置き替えました。

No、ｘ、Y,　H

0,-79824.543,33075.900,747.06

1,-79824.525,33081.133,750.05

2,-79824.508,33086.367,754.76

3,-79824.491,33091.600,759.48

4,-79824.474,33096.833,764.33

5,-79824.457,33102.066,768.78

福岡様

「業務で使うＱＧＩＳ Ver.3　完全使いこなしガイド」を購入しました。

　QGIS初心者である私に親切な御助言ありがとうございます。

　QGISの素晴らしさが少し分かりました。

　もっと勉強して使いこなせるようになりたいと思います。

　今後ともよろしくお願いいたします。

　感謝いたします。

2021年6月30日水曜日 9:49:38 UTC+9 福岡:

[福岡]

黒木 修身 様

ご投稿ありがとうございます。

書籍も準備頂き、なりよりと思っております。

それで気がついた点について、ご連絡いたします。

東雲さん宛てに示したデータをQGISに落としてみたのですが、宮崎県に落ちませんでした。

そこで、もしやと思いX座標とY座標を逆にしたところ、宮崎県に落ちました。

これは黒木さんが測量に通じていたからこそ、おきたことだと思いました。

海外製のGIS全てに言えるのですが、座標は数学座標が採用されています。最も普及しているArcGISもそのような仕様になっています。

つまり北方向がY軸、東方向がX軸になるわけです。

QGISに読み込む際は、１行目にあるフィールド名を書き換えるか、手作業で「X属性」と「Y属性」のフィールドを指定してください。

このようなことは背景地図敷けば容易に確認できますので、Web地図を読み込む環境を整えることをお勧めします。

今書籍を貸し出し中のため具体的なページをお示しできないのですが、「XYZTiles」に国土地理院のWeb地図を追加なさると良いと思います。

念のため「GIS実習オープン教材」での掲載箇所へのリンクを張っておきます。

https://gis-oer.github.io/gitbook/book/materials/QGIS/QGIS.html#%E5%9C%B0%E5%9B%B3%E3%82%BF%E3%82%A4%E3%83%AB%E3%81%AE%E8%AA%AD%E3%81%BF%E8%BE%BC%E3%81%BF

2021年6月30日水曜日 15:13:50 UTC+9 黒木修身: