等高線の表示のさせ方

[TAISUKE]

はじめて投稿します。

兵庫県のDEMデータから、１ｍ単位の等高線を表示させたいと思っています。

こちらのブログを参考にしたのですが、

https://note.com/kinari_iro/n/n25db0d04ba64

ラスタ→メニュー→等高線

までは行けたのですが、そこ先のパラメーターのところがmerge（投影座標系）が表示されず、DEMの選択肢が出てきます。

もとのデータはtxtデータなので、txtデータをxmlデータに変換しようとフリーのコンバーターでやろうとしてもエラーになります。

こちらどのようにすれば等高線が表示できるでしょうか。。。

とても困っていて教えていただけますと幸いです。

[sakuya_izayoi]

izayoiです。

まず、読み込ませた点群データをラスタに変換してからとなります。

https://groups.google.com/g/qgisshitumon01/c/aH12DlAHDao/m/ZYtT2AhrHAAJ

の記事を参考に

プロセッシングツールボックスの内挿コマンドの部分で、

TIN内挿かIDW補間を選んでください。

兵庫県が配っているデータというと、1m点群だと思われるので、ピクセルサイズは1.0で良いと思います。

また、点群の値をラスタピクセルの代表値としたい場合は、自動でセットされた「内挿する範囲」を0.5mずつオフセットしておくと良いかもです。

例えば、Xmin,Xmax,Ymin,Ymax=10 , 100 , 20 , 120だとしたら、9.5 , 100.5 , 19.5 , 120.5 のように、上下左右に0.5m分余白をもった四角形を範囲としてやるという感じです。

それで作成したラスタデータに対して等高線コマンドを使ってください。

2021年8月6日金曜日 14:29:22 UTC+9 TAISUKE:

[TAISUKE]

izayoiさん、お返事いただきありがとうございます。

・点群データをラスタに変換

・プロセッシングツールボックスの内挿コマンドの部分で、TIN内挿を選択

・ 作成したラスタデータに対して等高線コマンドを使う

→こちらのブログの抽出までを参考にしました。

https://note.com/kinari_iro/n/n25db0d04ba64

を行いましたが、QGIS上に表示されません。。

パネルは添付画像のようになっています。

2021年8月6日金曜日 14:56:19 UTC+9 sakuya_izayoi:

[sakuya_izayoi]

おそらくですが、点群からラスタに変換したのが【DSN_05NF072_1g】というラスタだと思いますが、いかがでしょうか？

（兵庫県の1m点群）

TIN内挿を実施する際の図は下記の図 のとおりにです。（添付1.png）

なお、内挿する範囲は、「…」をクリックして、レイヤの範囲を選択した後に、0.5mオフセットかけております。（添付図中では、yminのオフセット量を間違えております）

結果は下記の図です。（添付図2.png)

その図から等高線の抽出を行いましたが、問題なく出力ができました。

添付図をご参照頂き、再度トライしてみてください。

追伸

また、OUTPUTというラスタは数値がエラーを起こしており、また、座標系が不明となっていますが、それは何でしょうか？もしかして、その図より等高線を抽出しましたか？

[TAISUKE]

お返事いただきありがとうございます!

明日またそのとおりにしてみて、お返事させていただきます。

2021年8月6日金曜日 17:22:58 UTC+9 sakuya_izayoi:

[TAISUKE]

izayoiさん

先日はお返事くださりありがとうございました。

使われている元データはDEMデータではなくDSMデータですか？

ベクタレイヤがDSMになっていたのでそう思いました。

私はDEMデータを使っているのですが、それで等高線抽出できますでしょうか？

又、仰るとおりOUTPUTから等高線を抽出いたしました。

2021年8月6日金曜日 17:27:56 UTC+9 TAISUKE:

[sakuya_izayoi]

TAISUKE様

すみません、TAISEUKE様の場合、【DEM__05NF072_1g】とDEMデータを使われておりましたね。

もちろん、DEMでもDSMでも同様の方法で出せます。

やりかたは先日投稿した内容と全く同じです。

これで試してみてうまくいかないようでしたら、また投稿お願いします。

2021年8月10日火曜日 11:48:38 UTC+9 TAISUKE:

[sakuya_izayoi]

追記です。

おそらく、OUTPUTという名前のラスタファイルは、計算エラーを起こしております。（値が何も入っていないと思います。）

添付図の（添付1.png）の上から二番目の「内挿対象フィールド」にfield_3を選んでやり、内挿にZ座標を使用するにチェック、で、その下の右側にある「+」ボタンを押してやり、図のような感じに設定してやる必要があります。

2021年8月10日火曜日 15:57:52 UTC+9 sakuya_izayoi:

[ありた]

こんにちは

izayoi さんがおっしゃるように、 OUTPUT というラスタの生成の段階で、適切な

処理ができていないと思われます。

なお、兵庫県の 1m 間隔の標高（DEM）データであれば、ラスタ形式で配布

されています。位置情報はワールドファイルで用意されていますが、

座標系情報はもっていないので、記載にあるように平面直角座標系 V 系として

認識させてやる必要があります。

https://www.geospatial.jp/ckan/dataset/2010-2018-hyogo-geo-cs

[TAISUKE]

ありた様

お返事いただき、ありがとうございます。

ラスタ形式のtifデータをCRS変換しようとしたのですが、

不明なCRS と EPSG:4326 の間に可能な変換がありません

という表示で変換むずかしいです、、、

どのようにすればいいでしょうか、、、

2021年8月10日火曜日 21:03:42 UTC+9 ありた:

[TAISUKE]

izayoiさま

お返事遅くなりすみません。

大雨で山の補修作業に大わらわでした。

添付写真のようにしたのですが、表示されずでして、、

今度はシンプルに、

1、兵庫県からDSN(またはDEMのtxtデータをダウンロード)

２、TIN内挿で頂いた写真のように設定

という手順で行ったのですが表示されずでした、、

原因は何でしょうか、、

izayoiさん、お返事いただきありがとうございます。

2021年8月18日水曜日 16:37:09 UTC+9 TAISUKE:

[TAISUKE]

添付はこちらです。

2021年8月18日水曜日 16:37:40 UTC+9 TAISUKE:

[福岡]

TAISUKE 様

こんにちは

横入りになって申し訳ありませんが、気づいた点について投稿いたします。m(__)m

１　ラスタ形式のTiffについて

　ありたさんがおっしゃった様に、このデータを使うのが最も手間が掛からない方法だと思います。

　CRSは変換というよりも指定するということになります。

　レイヤのプロパティの「ソース」タブで、「設定されたCRS」を「JGD2000/Japan Plane Rectangular CS V」に指定することで解決すると思いました。

２　TIN内挿について

　画像を拝見すると、ピクセルサイズが0.1mになっているようです。

　izayoiさんの説明では1ｍとしていたと思いますので、設定が細かすぎるという印象です。

2021年8月18日水曜日 16:38:31 UTC+9 TAISUKE:

[sakuya_izayoi]

TAISUKEさま

izayoiです。

おはようございます。

福岡さま、ご助言ありがとうございます。

1.福岡さまの記載通り、ラスタファイルがある場合には、ラスタファイルから抽出した方が楽ですね。その際のCRS設定は福岡さん記載通りです。

2.TIN内挿についてですが、もともとの点群密度が1m間隔でしか存在していないので、ピクセルサイズは1mとしてください。それ以上細かいパラメータは今回のような格子型点群の場合、指定してもほぼエラーとなります。非格子型の場合はいけるかもですが、精度も何もあったもんではありません。

今回の例ではないですが、仮にレーザー測量などで0.1m間隔の点群が得られている場合には、ピクセルサイズは0.1と設定しても理論上TIN内挿が完了するはずです。

（理論上としているのは、あまりにも点数が多いとQGISの当該機能がハングアップする可能性があるから完成するとは言い切れないためです。）

3.「2.」で作ったラスタファイルから等高線を抽出する際には、等高線間隔は正直何mに設定してもQGIS（GDALやOGR2OGR）が追い付く限り計算はしてくれるため、お好みで設定すれば大丈夫です。

点群からラスタファイルを作成して等高線を抽出という作業については、結構躓きやすいところなので、また不明点があれば投稿お願いします。

2021年8月19日木曜日 9:07:53 UTC+9 福岡: