メッシュ標高差分計算(切・盛データの計算)及び地形変化量図の作成を行いたいです。

小山守

unread,

Oct 24, 2021, 7:46:45 PM10/24/21

to QGIS初心者質問グループ

お世話になります。

国土地理院の５ｍメッシュ（ラスタ）と古い地図メッシュ（ラスタ）との地形変化量図を作成したいです。

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

福岡

unread,

Oct 25, 2021, 12:24:03 AM10/25/21

to QGIS初心者質問グループ

こんにちは

ラスタ同士の差分について計算したいということでよろしいでしょうか？

それで良ければ、ラスタ計算機で差分を計算するのが良いと思います。

5mメッシュの値から古い地図のメッシュの値を引けば、増加した場合にプラス、減少した場合にマイナスで差分を得ることができると思います。

下図はラスタ計算機の設定例ですが、古いDEMデータを持っていないので、代わりに10mメッシュを古いDEMデータに見立てています。

メッシュの細かさと範囲は5mDEMに合わせた方が良いと思いますので、5mメッシュを選択してから「選択レイヤの領域」ボタンを押してみてください。

範囲とカラムが5mメッシュと同じになるはずです。

2021年10月25日月曜日 8:46:45 UTC+9 koyama...@gmail.com:

小山守

unread,

Oct 25, 2021, 2:29:34 AM10/25/21

to QGIS初心者質問グループ

福岡様

早速の回答ありがとうございます。

おかげさまで、下記のようなイメージ図はできました。

追加の質問で恐縮ですが、

体積を求めることは可能でしょうか。

例えば、3000m2以上とか。

元資料↓

https://www.gsi.go.jp/common/000200130.pdf

お手数おかけしますが

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

2021年10月25日月曜日 13:24:03 UTC+9 福岡:

tomoyuki arai

unread,

Oct 25, 2021, 12:20:31 PM10/25/21

to QGIS初心者質問グループ

koyama様

横からすみません。

少し気になったのですが、求めたいものは「体積」でしょうか？

フローチャート内でもメール文章でも「3000m2」と書いてあるように、「面積」を求める必要があるのではと思います。

盛土体積3000m3か、盛土面積3000m2かでは対象の抽出が大きく変わってくるため、よく確認した方が良いと思います。

新井智之

2021年10月25日月曜日 15:29:34 UTC+9 koyama...@gmail.com:

小山守

unread,

Oct 25, 2021, 7:32:59 PM10/25/21

to QGIS初心者質問グループ

失礼しました。

面積です。その面積に対して体積も求めることは可能でしょうか。

お手数おかけしますが

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

2021年10月26日火曜日 1:20:31 UTC+9 tarai...@gmail.com:

福岡

unread,

Oct 25, 2021, 10:06:38 PM10/25/21

to QGIS初心者質問グループ

新井　智之様

ご助言ありがとうございました。m(__)m

koyama 様

ラスタでは無く、ベクタで解析するのが良いような印象を受けました。

机上のイメージですので、抜けや勘違いがある可能性が高いように思いますが、次のように考えました。

１変化量のラスタをベクタ化する

　メニューの「ラスタ－変換－ラスタのベクタ化」

　 ※ラスタのCRSに変換後のベクタも引っ張られると思います。

　　地理座標系の場合は、メニューの「ベクタ－データ管理ツール－レイヤの再投影」などを使って、平面直角座標系などにする必要があります。

２盛土、切土のカテゴリー分類

　フィールド計算機を使って、色分けしたような分類にカテゴリー分類を行います。

　※関数は「CASE文」を使って、カテゴリー分類を行うのが良いと思います。

３盛土箇所ごとに箇所分類を行う。

　高さ分類が違っても場所が隣接しているようなところは、盛土箇所として一体と考えなければならないように思います。

　そのためのフィールドを追加して、IDなり名前なりを付けてやる作業が必要になりそうです。これは手作業になるイメージです。

４箇所分類に利用したデータを利用して、地物をまとめる

　メニューの「ベクタ－空間演算ツール－融合」で、基準となる属性に先ほどの箇所分類を設定してまとめます。

５フィールド計算機を使って、面積を計算する

　「$area」を使って、上記レイヤの面積を計算します。

６体積の計算をする

　「１」～「３」で使ったベクタに対してフィールド計算機の「$area」を使って面積を求めます。

　引き続きフィールド計算機を使って、この面積に盛土高を乗じて体積を求めます。

７エクセルに属性テーブルをコピペして、箇所別の体積を計算する。

　「sumif」で集計できるように思います。

結構な作業量になりますし、ある程度の基礎知識も必要な内容かと感じました。

もし業務で行っているのであれば、参考書籍を準備することをお勧めします。

私個人としては「業務で使うＱＧＩＳ　Ver.3　完全使いこなしガイド」が体系的に整理されていて、入門者には最適かと思っています。

2021年10月26日火曜日 8:32:59 UTC+9 koyama...@gmail.com:

Reply all

Reply to author

Forward