とらりもん - ASTER GDEM Diff
- Added parts are displayed like this.
- Deleted parts are displayed
like this.
!!!入手方法
以下のウェブサイトから無料ダウンロードできる。ユーザ登録が必要。
*日本 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
以下のウェブサイトから無料ダウンロードできる。ユーザ登録が必要。
*日本 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