【GISツール】GDAL / OGR ~ GDALコマンド編 ~

■ はじめに

https://blogs.yahoo.co.jp/dk521123/38092840.html
の続き。

今回は、主に使うGDALコマンドをまとめる。

■ gdalwarp (投影法変換)

https://gdal.org/programs/gdalwarp.html
 * 地図投影法の変換および歪みに関するユーティリティ
 # Image reprojection and warping utility

 * warp(ワープ) : ゆがむ、ねじれ、ひずみ
https://eow.alc.co.jp/search?q=warp

構文

gdalwarp [オプション] [入力画像] [出力画像]
オプション
 * -r [リサンプリング方法] : リサンプリング方法の指定(ex. near/bilinear etc)
 * -of [画像フォーマット]: (デフォルト:GeoTIFF(GTiff))
 * -overwrite : 上書き
 * -multi:マルチスレッドの有無
  => !メモ! 実験してみたが効果が分からなかった...
 * -tr [X方向] [Y方向]:解像度
 * -ts [幅] [高さ]:画像のサイズ(tsを指定すると、trは指定できない)

コマンド例

gdalwarp -of GTiff -overwrite IMG-HH-ALOS2014410740-140829-UBSL1.5GUA IMG-HH-ALOS2014410740-140829-UBSL1.5tiff

■ gdal_translate (画像フォーマット変換)

https://gdal.org/programs/gdal_translate.html
 * ラスタデータを異なる画像フォーマットに変換

# Converts raster data between different formats.

■ gdalinfo (画像データ内容の表示)

https://gdal.org/programs/gdalinfo.html
 * ラスタデータセットに関する情報一覧

# Lists information about a raster dataset.

構文

gdalinfo [オプション] [入力画像]

コマンド例

gdalinfo IMG-HH-ALOS2014410740-140829-UBSL1.5GUA
出力例
Driver: SAR_CEOS/CEOS SAR Image
Files: C:\IMG-HH-ALOS2014410740-140829-UBSL1.5GUA
Size is 29125, 32803
Coordinate System is `'
GCP Projection = 
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]],
    AUTHORITY["EPSG","4326"]]
GCP[  0]: Id=1, Info=
          (0.5,0.5) -> (139.421872,36.081964,0)
・・・略・・・
GCP[ 14]: Id=15, Info=
          (29124.5,32800.5) -> (140.237509,35.35061,0)
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,32803.0)
Upper Right (29125.0,    0.0)
Lower Right (29125.0,32803.0)
Center      (14562.5,16401.5)
Band 1 Block=29125x1 Type=UInt16, ColorInterp=Undefined