Morphometric Protection Indexの数値の意味（単位）について

[ron tabata]

お世話になっております。

プロセッシングのMorphometric Protection Indexのラスタの値は何を表しているのかを教えていただきたいです。

開度を求めたいのですが、Morphometric Protection Indexの使用後、値の大きい地点を見ると凹部（谷）になっておりました。（横山ほか，1999）によると開度の値が小さくなるほどより凹部（谷）になると思うのですが、反転する理由も知りたいです。

予想としては、仰角のradianなどかなとは思うのですが、そのような記載も見つけられずにいます。

ラスタの値が何を表しているのか、ソースも含めて教えていただきますと幸いです。

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

[wata909]

SAGAのマニュアルのページを見ると、
https://saga-gis.sourceforge.io/saga_tool_doc/2.2.0/ta_morphometry_7.html

「This algorithm analyses the immediate surrounding of each cell up to an given distance and evaluates how the relief protects it.」とあるので、視界を阻害する凹部の値が高くなるような指標化のかもしれません。

一方で、横山先生の2002年の論文を引用して、「This algorithm analyses the immediate surrounding of each cell up to an given distance and evaluates how the relief protects it.」とも書いてあるので、地上開度の逆数的な指標かもしれません。

https://www.asprs.org/wp-content/uploads/pers/2002journal/march/2002_mar_257-265.pdf

一応、元のコードは以下で
https://github.com/saga-gis/saga-gis/blob/master/saga-gis/src/tools/terrain_analysis/ta_morphometry/ProtectionIndex.cpp

このコードについてchatGPTに解読を依頼したやり取りは以下になります
https://chat.openai.com/c/c6966444-e2c2-4409-8a0b-2efef320ed13

とはいえ、途中で追加の質問を行うと、真逆の回答になっているので、ちゃんと自分でコードを確認したほうがいいと思います。

ご参考までに。

2023年10月2日月曜日 12:02:47 UTC+9 ron tabata:

[futoshi kidani]

「ソースも含めて」というのがソースコードではなくどこを根拠に…という事だとして、真偽はご自分でご確認願いますが

SAGAのフォーラムで過去にあったやりとり

https://sourceforge.net/p/saga-gis/discussion/790705/thread/137680bc/

Morphometric Protection Indexを使った解析事例

https://forswearer66.rssing.com/chan-59216107/article30.html

ここを見ると計算結果の値は、そのセルから8方向に考慮距離まで検索して、最大の角度のものの8方向平均ということになっています。

よって、谷地形でも考慮距離が短ければ開度が大きくなることも考えられますし、一方が開けていれば平均が変わるので開度がそれに引っ張られることもあるかもです。

どういうデータをどういう条件で解析されたか書いて下さると、追試・確認できるかもです。

2023年10月2日月曜日 12:49:25 UTC+9 wata909:

[ron taba]

>wata909様

ありがとうございます。山頂の値が0になっているため、地上開度の逆数ではなさそうなので「視界を阻害する凹部の値が高くなるような指標」の可能性の方が高いかなと思いました。

chatGPTに解読を依頼したやり取りについてですが、リンクに飛んだところUnable to load conversationと出てしまい読み込めませんでした。

>futoshi kidani様

ありがとうございます。データは国土地理院の10mメッシュ数値標高DEMで、radius200にして解析しました。

山頂部の開度が0になってしまうところからも反転しているのかなと思った次第です。

2023年10月2日月曜日 13:04:23 UTC+9 futoshi kidani:

[futoshi kidani]

こんなもの作って調べてみました。

200ｍ*200ｍのグリッドを作る

↓

全体を高さ1にする

↓

島を作って、その周りを1高さを上げる

↓

ＭＰＩをかける

結果として、8方向の水平からの角度（ラジアン）の平均と一致しました。

どうも一般的に開度というと天頂や天底（Ｚ軸）から角度を測るように思えますが、水平から測っているようです。

昔と挙動が変わったのかしらん？

2023年10月2日月曜日 18:16:24 UTC+9 ron taba:

[ron taba]

futoshi kidani様

お世話になっております。やはり水平からの角度だったのですね。

ラスタまで作成して実際に測っていただき大変感謝しております。

この度はありがとうございました。

2023年10月2日月曜日 20:36:13 UTC+9 futoshi kidani:

[wata909]

taba さま

>chatGPTに解読を依頼したやり取りについてですが、リンクに飛んだところUnable to load conversationと出てしまい読み込めませんでした。

おっと、リンクが間違えておりました。
https://chat.openai.com/share/8a225273-7b90-4d9e-82e0-93317a88c067

でいかがでしょうか？

2023年10月2日月曜日 18:16:24 UTC+9 ron taba:

>wata909様

[ron taba]

wata909様

お世話になっております。上記のリンクで無事アクセスできました。

単位についてですが、chatGPT的にはコード内に記載がないので単位なしと判断されていますね。

自分でもコードを確認しながら見てみます。

この度は大変ありがとうございました。

2023年10月3日火曜日 9:28:55 UTC+9 wata909: