DEMから山の尾根線を抽出する方法

[清掃奉仕]

-質問-

「DEMから山の尾根線を抽出する方法」

尾根を線としてshapefileで出力したいです。

-環境-

PCはWindows10でQGIS3.2を使用しています。

谷は上手く出力されるのですが尾根は上手くいきませんでした。

ちなみに谷(流路)の抽出方法は、

窪地を埋めたDEMを作成: プロセッシングツールボックス > SAGA > Fill sinks(wang & liu)

流路を埋め立てたDEMから抽出: プロセッシングツールボックス > SAGA > Terrain Analysis - Hydrology > Channel network and drainage basins

この工程で作成すると、谷は線としてshapefileで出力できたので似たような方法で作成したいと思っています。

ラスタデータで出力可能な方法でも良いので、ご存じであれば教えていただけますでしょうか。

[キタ]

こんにちは

谷を抽出できるのであれば、ラスタ計算機で「-1」を標高にかけて、谷と尾根を反転させたラスタを作成し、谷と同じようにすると抽出できるのではないでしょうか。

キタ

2022年6月4日土曜日 11:25:57 UTC+9 清掃奉仕:

[mooya]

清掃奉仕 様、キタ 様

こんにちは。

いろいろためしてみまして、、、

最初、この文献を参考に曲率で判定しようとしました
https://www.jstage.jst.go.jp/article/sabo1973/49/6/49_6_9/_article/-char/ja/

GRASS のr.slope.aspect で曲率のラスタを作成して、適当に敷地を設定して抽出しましたが、幅の広い尾根の抽出がいまひとつでした。
文献では計算範囲は25～30mが良いと書いているのですが、プロセッシングではその指定ができないのが原因かなと思います。








次に、キタ様の方法も試してみましたが、尾根の途中にピークがあると、反転させたピークが窪地となって埋められてしまいますで、思わしくありません。

次に、下記を参考にSAGA GISのConvergence Indexをためしたところ、比較的良かったので、その方法を紹介します。
https://ryotts.wordpress.com/2020/06/23/qgis%E3%81%A7%E3%83%A9%E3%82%B9%E3%82%BF%E8%A7%A3%E6%9E%90%EF%BC%9A%E5%B0%BE%E6%A0%B9%E8%B0%B7%E5%BA%A6%E3%81%A8%E8%B0%B7%E5%AF%86%E5%BA%A6/

1)Convergence index
　　　方法；[0]Aspect
　　　勾配計算；[1]3 x 3

2)ラスタ計算機　値が10以上を抽出
　　　式の例；"コンバージェンス指数@1" > 10

3）再投影（warp）
　　　出力バンドのnodata値；0　　
　　　出力データの型；int16
　次のthin（細線化）の前処理で、nodata値の指定と型変換をします。

4)grass　の r.thin （細線化）

5)grass　の r.to.vect　　（ラスタ→ライン変換）
　　　地物のタイプ；line


適切なパラメータは、標高ラスタのソースによるのではないかと思います。今回使用したので地理院の10mDEMです。

2022年6月4日土曜日 13:23:53 UTC+9 キタ:

[Bun]

こんにちは

QGISではなく、個別のアプリでG空間情報センターにて「数値地形解析アプリケーション」というものが公開されております。

そのなかの「TopoChar：地形特性図作成アプリ」では、谷部、尾根部、尾根線、鞍部深さなど解析できるようです。

利用条件等ありますので、条件合致されるようなら試してみるとよいかと思います。

2022年6月6日月曜日 0:24:47 UTC+9 mooya:

[Bun]

清掃奉仕 様 より

プライベートメッセージで返信がきておりましたので、以下貼付いたします。

ｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰｰ

ご記入いただいた、キタ様、bun様、mooya様

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

キタ様の方法でshapefileに書き出すことができました。現状、尾根が枝分かれしたり精度としては今一つですがパラメーターの設定を変えながらうまく抽出していきたいと思います。

bun様、mooya様の方法も順次試していきたいと思います。

このスレに書き込む形で報告しようと思いますのでよろしくお願いします。

もし他にも尾根線の抽出方法をご存じの方がいればお知らせください。

2022年6月6日月曜日 9:44:50 UTC+9 Bun: