数値地図の合成について

[Ryuiki]

はじめましてQGIS初心者です。

geotiffのDEM5ｍ（5ピクセル）が手元にあります。このデータはある範囲の全体をカバーしておらず、不足箇所が多数あります。この不足部分をDEM10ｍ（11ピクセル）で補ってやりたいのですが、ピクセルが異なるデーター同士の合成手法が分かりません。

現在、QGIS Desktop 3.16.6 with GRASS 7.8.5をWindows10　64bitで使っております。

ご存じの方がおりましたら、是非ともご教授をお願い致します。

[福岡]

Ryuiki 様

こんにちは

4年ほど前に、おそらく類似事例と思われることをこの掲示板で議論したことがあります。

データ無しの水面の標高を基盤地図情報数値標高モデルのGeotiffより推察して埋めたい

https://groups.google.com/g/qgisshitumon01/c/mTnnGJGk0-U

最後から2番目のラスタ計算機を使う方法が、お役に立つように思いました。

2021年4月26日月曜日 14:20:10 UTC+9 Ryuiki:

[Ryuiki]

福岡様

早速のご教授に感謝いたします。

＞最後から2番目のラスタ計算機を使う方法が、お役に立つように思いました。

同じように5ｍメッシュの穴埋めに悩んだ方がおられたのですね。

有難うございました！

2021年4月26日月曜日 14:39:45 UTC+9 福岡:

[Ryuiki]

福岡様

申し訳ありません。肝心な疑問が残っております。

地理院からDLLできるのは、DEM5ｍは3次メッシュのデータ、DEM10mは2次メッシュのデータです。

ラスターですので、それぞれの大きさが異なると思うのですが、計算により合成して出来上がったデータは、高さ・距離等に問題は生じないのでしょうか？

合成についてまだよく理解しておらず、反省しております・・。

2021年4月27日火曜日 10:32:15 UTC+9 Ryuiki:

[福岡]

Ryuiki 様

こんにちは

ラスタ計算機の右側に「X最小値」「X最大値」「Y最小値」「Y最大値」「カラム」「行」があると思います。

「X最小値」「X最大値」「Y最小値」「Y最大値」でラスタデータの範囲、「カラム」「行」で縦横のピクセル数が決まります。

そして、これらのすぐ上のところに「選択レイヤの領域」というボタンがあります。

バンドという窓内にあるレイヤを選んでアクティブにして、このボタンを押すと対象レイヤの範囲とピクセル数に設定されます。

これらを踏まえて5mDEMの範囲内で穴埋めをするのなら、5mDEMのレイヤを選んで「選択レイヤの領域」ボタンをクリックすると良いと思います。

この操作で、新たに生成されるラスタは5mDEMと同サイズで同様の細かさになります。

そして例示した式を適用すると、5mDEMの値が「-9999」となっているセルに、セルの中心が重なる10mDEMの値を持ってくることになります。

2021年4月27日火曜日 11:15:42 UTC+9 Ryuiki:

[Ryuiki]

福岡様

丁寧なご説明、恐縮しております。

ついでということでご勘弁願いたいのですが、DEM5、DEM10ともにオリジナルは地理座標のJGD2011ですが、最終的には投影座標10系に統一する予定です。

ですので座標の再投影（warp）については、合成前に一致させておいた方が良いと思いますが、入力する際のパラメータに注意点があればお聞かせ願えればと思います。

度々申し訳ありません。

2021年4月27日火曜日 11:49:49 UTC+9 福岡:

[福岡]

Ryuiki 様

こんにちは

自信を持って推奨はできませんが、私ならオリジナルで合成して、平面直角座標への変換はその後に行うと思います。

5mメッシュの正体は0.2秒（約5m）メッシュ、10mメッシュも0.4秒（約10m）メッシュです。

私の感覚的なものですが、今回のケースは生（オリジナル）から調理（合成）して、盛り付け（再投影）はその後というイメージです。

また、失敗をするのなら合成の方だと重いますので、合成前の手間がない方が試行錯誤しやすいというのも理由の1つです。

他の人はどうなんでしょう？

ご意見がありましたら、投稿をよろしくお願いいたします。m(__)m

2021年4月27日火曜日 13:19:26 UTC+9 Ryuiki:

[Ryuiki]

福岡様

早速のご返答有難うございます！

まずはオリジナルからそのままラスタ計算で、4年前の式を参考にさせて頂いて合成したところ、成功したようです。

というのは、陰影図で見た限りでは5ｍと10ｍの差異が良くわからないのです。

試しにシンボロジを等高線にして描画したところ、詳細な間隔に変化させても1/2.5万地形図にほぼ整合しました。

これは成功と判断しました。

取り急ぎ御礼を申し上げます。

2021年4月27日火曜日 14:27:44 UTC+9 福岡:

[福岡]

Ryuiki 様

うまく行ったようでよかったです。

ラスタを扱っての結果を確認するには「Value Tool」というプラグインを使うのがお勧めです。
これを使うとマウスオーバしたところの各レイヤの値を確認できます。

今後もラスタを扱う機会が多いのでしたら、インストールして損はないと思います。

※間違ってメールの方へ投稿しましたので、掲示板への再投稿です。

　ダブった方がおられましたら、ごめんなさい。m(__)m

2021年4月27日火曜日 14:38:18 UTC+9 Ryuiki:

[Ryuiki]

福岡様

Value Tool 、是非使ってみたいです。

重ね重ね有難うございました。

ずっとArcGISを使ってきたのですが、このところ業務でQGISを使うようになり、その機能の凄さに驚いております。

しかもフリー。

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

2021年4月27日火曜日 14:52:44 UTC+9 福岡: