等高線からDEMデータ作成について

[kenke...@gmail.com]

皆様のお知恵を拝借させてください。使用バージョン　QGIS ver.3.4.15 、GRASS GIS ver.7.8.2

一般的等高線データからDEMを作成したいと考えており、検索の結果２通りの方法があったので試してみたのですが、いずれも途中でつっかえてしまいます。

①「GRASS GIS v.surf.rstコマンドで等高線からDEMを作成する」https://qiita.com/ishiijunpei/items/bf404865c113c0be62e8

この方法では「GRASSツール→ファイル管理→GRASSへベクタをインポートします→ロードされたベクタをインポートする」の段階でインポート完了となりますが、いずれも一部エラーとでてしまい、結局GRASSで表示されていません。処理が正しく行われていない？

②https://groups.google.com/forum/#!searchin/qgisshitumon01/%E7%AD%89%E9%AB%98%E7%B7%9A%E3%82%92DEM|sort:date/qgisshitumon01/oIALa6eao5c/zOPa2F8hCgAJ

こちらでは4「プロセッシングツールのGRASS7から「v.to.rast.attribute」を使って、等高線を標高を値に持つラスタに変換します。」の段階がうまくいきません。そもそも「v.to.rast.attribute」ではなく、「v.to.rast」しかありませんが同じもの？　各種パラメータの設定がよく分かりません。5の「r.surf.contour」も然りですが・・・

もっと簡略な方法で作成できれば良いのですが・・・　宜しくお願い致します。

[福岡]

おはようございます。

これから出かけますので、簡単な記述になりますがご容赦ください。m(__)m

ざっとですが、先ほど動作を確認しました。

使ったのは基盤地図情報からダウンロードしたDEMから作成した等高線で、それを素材にDEMを作成できるかを確認してみました。

手順としては次のとおりです。

１．「v.to.rast」でベクタの等高線をラスタにする。

　念のため「Name of column for 'attr' parameter」に標高が記録されているテーブルを指定しました。

２。１で作成したラスタから「r.surf.contour」で、DEMを作成する。

　デフォルトのまま実行しました。

結果としてはセルサイズを指定しないデフォルトのままでは、セルサイズが粗すぎると感じましたが、出力はできました。

詳しく検証せず、情報提供程度の内容となりましたが、よろしくお願いいたします。

[福岡]

おはようございます。

帰ってきてから、確認してみたことについて投稿いたします。

セルサイズの数値は直感で設定しています、できればご自分でもお試しになることをお勧めします。

まず基盤地図情報からダウンロードしたDEM（10mメッシュ）をラスタに変換したものを掲載します。

これと等高線から作成したDEMと比較することによって、できあがりの評価ができると思います。

今回の作業に使う等高線を作成しました。間隔は10mにしてあります。

これを材料として、DEM作成を試してみます。

最初に「v.to.rast」をデフォルト設定のまま実行します。

前の投稿にも記載しましたが「Name of column for 'attr' parameter」だけ標高が記録されているテーブルを指定しています。

これを元に「r.surf.contour」で、DEMを作成します。「GRASS GIS 7 領域のセルサイズ」を「10」としました。

できあがりをみるとダウンロードしたDEMと乖離しています。これでは不満です。

そこで「v.to.rast」で「GRASS GIS 7 領域のセルサイズ」を「10」として出力してみます。

これを元に「r.surf.contour」で、DEMを作成します。先程と同様に「GRASS GIS 7 領域のセルサイズ」を「10」としています。

このできあがりでしたら、ある程度満足できるのではないかと思います。

[kenke...@gmail.com]

福岡様

ご検討いただきありがとうございます。

自分で試してみたのですが、以下の点で躓いております。

試しに用いたデータは基盤地図情報からダウンロードした標高点データ(xml)をshpに変換したものを用いています。

１.「v.to.rast」でラスタ変換する際、「Name of column for 'attr' parameter」に標高が記録されているテーブルを指定とありますが、

　　ドロップダウンで選べる選択肢が整備データと整備デーAしか選べません・・・

２．レイヤプロパティを見るとソースフィールドではID abc9に標高と記されています。

　　これを整備データと同じIDに変更できれば、先ほどの1で選択可能となるのでしょうか？　またIDの変更はどのように行えば良いのでしょうか・・・。

　

　　　　

[福岡]

おはようございます。

私の投稿したものは、等高線のラインデータからDEMを作成する方法でした。

基盤地図情報のxmlからDEMを作成するのであれば、もっと簡単な方法があります。

私が使っているのはエコリスさんの「基盤地図情報 標高DEMデータ変換ツール」で、愛用者も多いと思います。

このツールを使えば、xmlからGeoTiff形式のDEMを簡単に作成することができます。

株式会社エコリスさんのホームページ

https://www.ecoris.co.jp/

基盤地図情報 標高DEMデータ変換ツールのページ

https://www.ecoris.co.jp/contents/demtool.html

[kenke...@gmail.com]

福岡様

おはようございます。

ご返信ありがとうございます。

エコリスさんの「基盤地図情報 標高DEMデータ変換ツール」については知っていたのですが、今回は等高線からDEMを作成するということを何とかできないかと試行錯誤しています。

というのも、とある海外地点のDEMを作成したいのですが、どうにも等高線のshpファイルしか無いらしく、そのデータはまだ購入していないため、試しに基盤地図情報のshpファイルから作成できないかと実験をしておりました。

[福岡]

お返事ありがとうございます。

実験ということですので、私と同じようにトライしてみてはいかがでしょうか？

１　基盤地図情報から標高データをダウンロード

　　お勧めは「10mメッシュの10B（地形図の等高線」です。5mメッシュの場合、データ抜けがある場合があります。

　　それにxmlのファイル数が5mメッシュよりも少なくなるので、作業が楽になると思います。

２　エコリスさんの「基盤地図情報 標高DEMデータ変換ツール」でDEMを作成

３　QGISのプロセッシングツールの「GDL － ラスタ抽出 － 等高線」を使って等高線を作成

　　デフォルトでは、標高は属性テーブルの「ELEV」へ書き込まれます。

４　「v.to.rast」でラスタ版の等高線を作成

　　以前の投稿に書いたように、標高データが書かれているところ（ELEV）の指定と領域のセルサイズを指定してください。

５　「r.surf.contour」でDEMを作成

　　こちらも以前の投稿に書いたように、領域のセルサイズを指定してください。

ところで、私はポイントデータから試していないのですが、「v.to.rast」で標高データを選択できない理由についてお伝えします。

「Name of column for 'attr' parameter」は数値でなくてはいけないのですが、貼って頂いた画像をみると「String」（文字列）になっています。

これを「Integer」（整数）または「Real」（実数）にする必要があると思います。

この修正は、プロセッシングツールの「ベクタテーブル － フィールドのリファクタリング」でできます。

標高点のデータで実験をお続けになる場合には、お試しください。

[kenke...@gmail.com]

福岡様

懇切丁寧な回答ありがとうございます。

福岡様のxml→DEM→等高線→DEMは滞りなく作成できました。

しかし今度用いる海外のデータはshpのみなので、xmlデータは恐らくないことが想定されます。

そのためなるべく等高線shpからDEMを作成したいと考えております。

ご教授いただいたとおり、プロセッシングツールの「ベクタテーブル － フィールドのリファクタリング」で標高データを「Integer」（整数）に修正し、

「v.to.rast」の「Name of column for 'attr' parameter」ドロップダウンに表示されましたので、

これで実行を選択すると下記のような表示でエラーが出てしまいます。

やはり難しいのでしょうか・・・

[kenke...@gmail.com]

福岡様

検討いただきありがとうございます。

「標高」を「level」に表記を変更したところうまく変換してくれました！！

ありがとうございました！！

[sakuya_izayoi]

こちらizayoiです。

以前ちょっと興味があったので実施してみました。

その時の条件ですが

・5m解像度のラスタまではたいていの場合は作成可能。

・2.5m解像度はラスタ化範囲次第では可能。

・1.0m解像度は実用的な範囲では作成しきってくれない

という事でした。

あんまり細かいラスタ化は注意が必要でしたという状況提供でした。

[takaha...@gmail.com]

izayoi様

情報提供ありがとうございます。

地形図から無理やり起こすのは、詳細なDEM作成の場合には適さないということになりますね。

[sakuya_izayoi]

izayoiです。

ちょっと実験的に細かいのをやってみました。

30m×30mでならば1m以下の解像度のラスターDEMも作れました。(0.5m解像度まではいけたけど、めっちゃくちゃ時間がかかる。)

自分の場合はこの面積では実用性が皆無なのですが、農地10a分でいいよという場合には細かいのも作れる感じです。

私の場合はどうしているかというと、有料製品ソフトなのですが、AutodeskのCivil3Dを使って詳細DEM作っております。(1km×1kmの0.5m解像度とかになるとかなり時間かかるけど)

オープンソースソフトの方がどのPCでも詳細DEMを作れるので,

できれば詳細DEMを作る方法が知りたいです。