-コンテンツへスキップこのwebサイトの先進的な機能は、ブラウザーのJavaScriptを有効にする必要があります。コメントを残す異なる解像度の数値標高モデルを組み合わせる方法についての覚書です。ここでは国土地理院が提供する高解像度の5mメッシュデータと10mメッシュデータを組み合わせて等高線を生成する手順を紹介します。はじめに国土地理院のサイトを覗いてみると5mメッシュの数値標高モデルが公開されています。以前GarminGPS向けに地形図を作成した際は日本全土をカバーしている10mメッシュデータを利用しました。これでも十分地形が分かる等高線が作れましたが、さらに高解像度の5mメッシュだとどうなるだろうかと気になります。5mメッシュだと単純に解像度が倍になるので、これを元にすれば更に精度の良い地形図が作れるはず。残念ながら現時点(2018/08/10)ではまだ日本全国をカバーしていません。そこで不足している部分を他のデータで補って等高線を生成出来るかやってみることにしました。進んでいるようで山岳地帯のカバー範囲も広がってきました。作業環境は以下の通りです。今回はいつもの仮想PCではなくWindowsとのデュアルブート環境を構築して作業しました(作業過程でメモリを猛烈に食うので)。18.04使用する数値標高モデル以下の数値標高モデルを等高線生成に用います。DEM5A(5mメッシュ,航空レーザ測量)DEM5B(5mメッシュ,写真測量)DEM10A(10mメッシュ)数値標高モデルは以下のサイトからダウンロード出来ます。前準備データをダウンロード出来たらツールを使ってGeoTIFFに変換します。この辺りの流れは以下の記事を参照ください。変換したGeoTIFFは種類毎にディレクトリに分けておきます。./DEM5A/./DEM5B/./DEM10B/解像度を合わせる(リサンプリング)データが揃ったので3種類ある数値標高モデルを組み合わせます。ここではgdalbuildvrtコマンドを使ってvrt形式ファイルにしました。以下のようにコマンド入力してください。../DEM5A/geotiff/FG-GML-*.tif引数のファイルは後ろになるほど優先順位が高くなります。すると、DEM5A->DEM5B->DEM10Bの順に高度情報が参照されるので、DEM5Aが提供されていない座標ではDEM5B、それも無ければDEM10Bの高度情報が読み込まれます。"-resolution"オプションは解像度を指定します。highestとすると、解像度をGeoTIFFファイルの中で最も高いもの(ここでは5mメッシュ)にあわせます。無指定だと"average"となり、解像度が(10+5)/2=7.5mになってしまいます。"-r"オプションはリサンプリングのアルゴリズムを指定します。ここでは"cubic"を指定しました。どれがベストなのかはっきりとはわかりませんが、少なくとも無指定(nearestが選択される)だとDEM10Aを参照する場所の等高線がカクカクになります。前項のvrtファイルを使って等高線を生成し地形図を作ります。ここではGarminGPS向けの地図ファイルを作成しました。詳細な手順は以下の記事を参照してください。合計となりDEM10Bを使ったときは5GByteくらいだったのでほぼ4倍のファイルサイズになりました。分解能が倍になったのでここは想定通りです。試しに作ってみたGarmin地図では、10mメッシュ版:2.8GByte5mメッシュ版:3.8GByteとかなりサイズが増えました。GarminGPSはFAT32の制約でファイルサイズ上限4GByteなのでギリギリです。3D表示対応のためにDEMをくっつけたら恐らく容量オーバーになりますね。指定したアルゴリズムによって等高線がどう変化するか比べてみました。以下の画像中央付近にある”旗のアイコン”を境界にして左側が5mメッシュ、右側が10mメッシュのデータです。”nearest”を使ったケースなどでは10mメッシュ側の等高線がリサンプリングの影響でガクガクになっているので分かりやすいと思います。"cubic","cubicspline","bilinear"辺りが無難かな。参考までにリサンプリング無し(DEM10B、10mメッシュのデータだけ)を掲載します。出来上がった地形図と見比べると5mメッシュの部分はより細かい地形が判別出来そうです。その一方で10mメッシュと5mメッシュデータの境目となる場所の等高線には違和感がありますね。異なる解像度の数値標高モデルを組み合わせて等高線を生成してみました。作業手順としては問題なさそうですが、もし登山で利用するとなると、うーんどうでしょうか。これから登る場所がちょうど継ぎ目に該当しなければ問題ないかもですが、富士山頂付近なんかがちょうど玉虫色になっているんですよね。次回の地形図更新では試しに10mメッシュ版と10m/5mメッシュ混合版のふたつを作ってみようかな。参考この記事は以下の内容を参考にさせて頂きました。地図sukeについて自転車で日本一周旅行に出ようとー思い立ったエンジニア。ガジェット類やセミナーに目がない。コメントをキャンセルメールアドレスが公開されることはありません。コメントデータの処理方法の詳細はこちらをご覧ください。