XYZ形式のDEMデータをラスタに変換する

2013.11.27追記:XYZ形式をQGISだけでラスタに変換する方法が分かりました。
こちらの記事を参照ください。

メッシュDEMデータがXYZ形式の場合,直接QGISに取り込めない(直接取り込めました。詳しくは上記リンク先の記事で。)のでラスタデータに変換する.

全メッシュのデータがある場合:
file=DEMデータ(XYZ形式)
name=保存ファイル名
として,
$ gdal_translate $file $name.tif
でいけるらしい.試してません.
gdalが必要.QGISと同時にインストールされるはず.

Z=NaNのメッシュデータが欠けている場合:
GMTのxyz2grdを使う.図郭の左下・右上のXY座標とメッシュサイズが必要.メッシュの縦横比が1:1でなければいけない.そうなっていない場合,grdsampleとかを使ってメッシュの切り直しをする.
range=x_min/x_max/y_min/y_max
mesh=メッシュサイズ
NaN=データなしの値(陸域のみのデータなら-9999とか)
として,
$ xyz2grd -V -R$range -I$mesh $file -G$name.grd -N$NaN
生成される.grdファイルは,QGISで直接読み込み可能.

※GRASS GISにはr.in.xyzというそのものズバリなコマンドが用意されている.上記の方法と比較する予定.

リンク:lat, lon points to a GeoTiff

コメントを残す