gdalwarpでタイル状のラスター型GISデータを接合

GDALのプログラム群のひとつである、gdalwarpを使うと、ラスター型（メッシュ）のGISデータの接合が簡単にできるので自分へのメモを兼ねて紹介したい。

国土地理院の基盤地図情報やASTER GDEMなどで公開されている地盤高データ（Digital Elevation Model; DEM）は、タイル状に分割されて公開されている。そのため、境界部を対象にする場合や、広域な利用を考える場合は、データを接合する必要がある。

基盤地図情報については、以前のエントリでラスター型GISデータへ変換できることを紹介した。先日公開がはじまったASTER GDEMはGeotiff形式のラスター型GISデータとして公開されている。これらのデータは、GDALのgdalwarpを使うことで、簡単に接合処理ができてしまう。

gdalwarp ファイル1 ファイル2 出力ファイル

ファイル1とファイル2をつなぎ合わせたものが出力ファイルに保存される。大量のファイルの場合は、アスタリスク"*"を使いワイルドカードによって大量のファイルを指定することができれば、大量のタイル状のファイルを接合することができる。つまり、

gdalwarp ワイルドカード付きのファイル 出力ファイル

とすることで、大量のタイルを接合することが可能となる。

今回、基盤地図情報の10mのDEMを以前エントリの方法でラスター型GISデータへ変換し、gdalwarpで接合した。以下に接合後のファイルを示した。