とらりもんHOME  Index  Search  Changes  Login

国土地理院数値地図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
Last modified:2008/12/30 12:08:39
Keyword(s):
References:[トレーニングコース] [とらりもんHOME]