GDALでモザイク処理の変更点

追加された行はこの色です。
削除された行はこの色です。
GDALでモザイク処理へ行く。
GDALでモザイク処理の差分を削除

[[GIS/RS一般]]
20101202
~#　処理1：JGD2000緯度経度のままで全道のDEMをモザイク
--処理1：JGD2000緯度経度のままで全道のDEMをモザイク
~#　コマンドプロンプトを起動しカレントディレクトリをデータ格納場所に変更~
cd /d D:\DEM\GSI_50m~
~#　.img拡張子のついたディレクトリ内のすべてのファイルのフルパスをlist.txtに出力~
dir /b /s *.img>list.txt~
~#　OSGeo4wシェルからgdal_merge.pyを使用してモザイク処理を実施~
~#　--optfile　でlist.txtを読み込むことで、処理対象ファイルを指定する~
~#　-init　でNodata(海の部分)は"-9999"と指定~
python C:\OSGeo4W\bin\gdal_merge.py -o D:\DEM\GSI50m\merge.img -init -9999 --optfile D:\DEM\GSI50m\list.txt~
~#　モザイク処理の所要時間：全道分で3分程度~
~#　処理2：投影変換~
--処理2：投影変換~
~#　gdalwarpを使用してJGD2000緯度経度からJGD2000UTM54帯のデータへ変換~
~#　-tr で解像度を50m×50mと指定~
gdalwarp -s_srs EPSG:4612 -t_srs EPSG:3100 -dstnodata -9999 -of HFA -tr 50 50 D:\DEM\GSI50m\merge.img D:\DEM\GSI50m\merge_j2ku54.img~
~#　投影変換処理の所要時間2分程度~

GDALでモザイク処理 の変更点