ASTER GDEM
入手方法
以下のウェブサイトから無料ダウンロードできる。ユーザ登録が必要。
- 日本 ERSDAC: http://www.ersdac.or.jp/GDEM/J/index.html
- 米国 USGS WIST: https://wist.echo.nasa.gov/~wist/api/imswelcome/
データの詳細
- ここに詳しく書いてある>ERSDACのサイト: http://www.ersdac.or.jp/GDEM/J/2.html
- 解像度30m
- Latlon(geographic; 等緯度経度)、WGS84/EGM96
- GeoTIFFの1度x1度のタイル形式
- 全球をカバーする無料のDEMでは、これが一番解像度が高い。
- 米国で配られているREADME(pdf): http://www.viewfinderpanoramas.org/GDEM/GDEM-README.pdf
- 30m解像度で作られているが、実際の地形表現は100-120m解像度らしい。
モザイクしてGRASSにとりこみの例
- モザイクは、GDALを使うと簡単。以下、解凍したtifファイルをディレクトリ下に集めて実行する。ただ、たくさんいっぺんにやると時間かかる。
$ gdalwarp *dem*tif out.tif
- GeoTIFFなので、GRASSとりこみもgdalで。
latlon+wgs84のlocation上で、 $ r.in.gdal input=out.tif output=ASTER_GDEM
欠測の補間
- GRASSで、
$ r.mapcalc dummy="if(ASTER_GDEM==-9999, null(), ASTER_GDEM)" $ r.fillnulls input=dummy output=ASTER_GDEM_Japan_fill $ g.remove rast=dummy
Keyword(s):
References:[とらりもんHOME]