GDAL のコマンドで GeoTIFF ファイルを扱う

【概要】

コマンドで GeoTIFF ファイルを扱う方法を説明する．GeoTIFF ファイルの拡大（縦横の画素数を増やす）などを行う．

【その他情報】

キーワード： GeoTIFF 情報の取得, 投影座標系への投影, GeoTIFF ファイルの拡大, GeoTIFF から png へ変換,png から GeoTIFF へ変換, EPSG コードを用いた座標系の変換

【目次】

前準備

あわせて，環境変数 GDAL_DATA を「C:\OSGeo4W\share\gdal」に設定しておくこと．

GeoTIFF について練習したい場合には， GeoTIFF サンプルデータファイルを入手する．例えば，次の Web ページから cea.tif をダウンロード

ダウンロードした .tif ファイルを，分かりやすいディレクトリ（例えばd:\)に保存

gdalinfo コマンドを使用．サイズ，緯度経度、バンド数、データ型などが確認できる．

gdalinfo cea.tif

あるGeoTIFF画像での実行結果例

別のGeoTIFF画像での実行結果例

UTM座標系に投影
http://spatialreference.org/ref/?search=JGD2000&srtext=Search （現存しない）

EPSG:3099 (JGD2000 UTM座標系ゾーン53N)に変換したいときは次のように実行する．
```
gdalwarp -t_srs EPSG:3099 ＜変換元ファイル.tif＞ ＜変換先ファイル.tif＞
```
実行結果例

次のようなエラーが出たときは，環境変数の設定が正しくない． OSGeo4W のインストールの Web ページを参考に，環境変数の設定を行うこと
平面直角座標系に投影
http://spatialreference.org/ref/?search=JGD2000&srtext=Search （現存しない）によれば次の通り．

国土地理院のWebページ https://www.gsi.go.jp/sokuchikijun/jpc.html などを参考に，平面直角座標系の番号を調べる
EPSG:2445 (JGD2000 平面直角座標系III系)に変換したいときは次のように実行する．
```
gdalwarp -t_srs EPSG:2445 ＜変換元ファイル.tif＞ ＜変換先ファイル.tif＞
```

GeoTIFFファイルの画素の値の範囲を調べる
```
gdalinfo -mm cea.tif
```
下の実行例では，画素の範囲は 0 から 255だと分かる
png に変換する前に決めておくこと
png ファイルは，1画素 8ビットなので, 0 から 255の値．
GeoTiff から png に変換するとき，次のことを決めておく
- ファイル名: ここでは cea.pngに設定
- スケール: ここでは，GeoTIffでは 0 から 255の値を，png での0 から 255の値に変換するとする．
GeoTIFF から png への変換
```
gdal_translate -scale 0 255 0 255 -ot Byte -of PNG cea.tif cea.png
```
結果は，画像ビューワで簡単に確認できる．

png ファイルは，1画素 8ビットなので, 0 から 255の値．

GeoTiff に変換するとき，次のことを決めておく

ファイル名: ここでは cea3.tif
スケール: ここでは，png での0 から 255の値を，GeoTIffでは 0 から 24700に変換するとする．
GeoTIFFの画素のデータ型: ここでは Int16に設定
データ型については「gdal_translate --long-usage」で確認できる．

png から GeoTIFF への変換例

gdal_translate -scale 0 255 0 24700 -ot Int16 -of GTiff cea.png cea3.tif

結果は，画像ビューワで簡単に確認できる．

拡大率 200% にしたいときの操作手順例

gdal_translate コマンドを用いる

gdal_translate -outsize 200% 200% d:\cea.tif d:\cea2.tif

gdalinfo コマンドを用いて確認してみる

GeoTIFF ファイルは，多くの場合，画像ビューワで簡単に表示できる．

拡大前

拡大後

gdalwarp -s_srs EPSG:＜EPSGコード＞ -t_srs EPSG:＜EPSGコード＞ ＜ラスタファイル名＞ ＜ラスタファイル名＞