とらりもん - 国土地理院数値地図50mメッシュ(標高) Diff
- Added parts are displayed like this.
- Deleted parts are displayed
like this.
履歴: 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
以下、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