【水域などの窪地処理について】

[akiyosea]

初めて質問させていただきます。

QGISバージョン：3.16　パソコンのOS：21H2

★経緯

国土地理院「基盤地図情報ダウンロードサービス」から、

5mメッシュDEMデータを取得したところ、

データ中に水域（標高値-9999，QGIS上で黒塗り表示）があると思われたため、

テキスト「業務で使うQGIS（全国林業改良普及協会出版）」を参考に下記作業を実施しました。

①プロセッシングツール「SAGA/Terrain Analysis-Hydrology/Fill sinks xxl(wang&liy)」でminimun Slope0.01とし窪地の除去。

②窪地除去した標高DEMと元の標高DEMを地図に重ね、ラスタ計算機で、水域（-9999）の部分のみ、窪地処理した標高値を入力した新しいDEMを作成する。

　ラスタ計算機の式：("元のDEM@1"=-9999)*"窪地除去したDEM@1"+("元のDEM@1"!=-9999)*"元のDEM@1"

★質問『①時点の窪地除去した標高DEMレイヤでは、水域の黒塗りがなくなっていますが、②で最後に作成されたレイヤ（窪地処理した標高値を入力した新しいDEM）は黒塗りされたままになっており，窪地除去できていないように思えます。

　⇒この事象や対応策等についてご教示ください。初歩的な質問で申し訳ありません。

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

[福岡]

こんにちは

おそらく「元のDEM」の黒塗り表示の場所が、「-9999」ではなく「nodata」として捉えられているのだと思います。
これを回避するためには「元のDEM」に対して「プロセッシングツール」の「ラスタツール －　nodataセルを値で埋める」を使って、「nodata」を「-9999」で埋めると良いと思います。

実は似たような経験をしていまして「5mメッシュの欠損値を10mメッシュDEMの値で補完する」という記事では、上記の手順を使っています。
参考になることもあるかもしれませんので、ご関心がありましたらご覧ください。

5mメッシュDEMの欠測値を10mメッシュDEMの値で補完する
https://qiita.com/fukuoka1960/items/bc75ac8ec2dea23fd098

2022年11月1日火曜日 9:18:09 UTC+9 akiyo...@gmail.com:

[akiyosea]

福岡様

早速のご対応ありがとうございます。

改めて5mメッシュデータをGeoTIFFファイルに変換したら、当該箇所が黒塗りではなく白塗りになりました。（変換ツールで海抜を0にしてしまったため黒になっていたか。）

それが福岡様のおっしゃるとおり「-9999」ではなく「nodata」で欠損しているように思えました。

リンクいただいた10mメッシュの値で補完する方法を試したら、それらしいものができました。

濃淡の調整等色々と試してみます。

大変助かりました。

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

2022年11月1日火曜日 9:42:02 UTC+9 福岡: