国土地理院数値地図50mメッシュ(標高)
履歴: 2008/12/15 奈佐原
以下、GRASSによる処理法です。
sudo mount /dev/cdrom /mnt/1 ## 読み込み ## 古いCD-ROMの場合はTokyo datumのlatlonに入って。 ## 新しいCD-ROMの場合はWGS84のlatlonに入って。 cd /mnt/1/data/5338 for i in *.mem; do r.in.gdal input=$i output=$i; done cd ../5238 for i in *.mem; do r.in.gdal input=$i output=$i; done cd ../5337 for i in *.mem; do r.in.gdal input=$i output=$i; done cd ../5237 for i in *.mem; do r.in.gdal input=$i output=$i; done
## モザイク g.region w=137 e=139 s=34.5 n=36 r.mapcalc "dummy=null()" g.list rast | awk '{for (i=1; i<=NF; i++) print $i}' | grep mem | \ awk '{print "r.mapcalc \"dummy1=if(isnull(dummy),"$1",dummy)\" g.remove rast=dummy; g.rename rast=dummy1,dummy"}' > ~/dummy.sh sh ~/dummy.sh
### 測地系の変換 (新しいCD-ROMの場合は必要) # WGS84のlatlon のlocationに入って、 r.proj input=dummy location=TokyoDatum1 mapset=PERMANENT output=ggg
Keyword(s):
References:[トレーニングコース] [とらりもんHOME]