qgisのgdal_translate

# 概要
qgisでgdalを利用する方法

# qgis
いわずとしれた,opensourceの地図解析ソフト
企業は高いお金をだして、ArcGISで地図を解析するが
安くあげるにはQgisでしょうな

# gdal
昔からあるラスター解析用のツール。qgisをインストールするとなんとbinのディレクトリにコピーされている。まあ、unix系であれば、なんなくapt-getなんかで使えるのだが、windowsはなかなかそうはいかない。(windows10はapt-getができると聞いたがまだ使ってない）

qgisのラスタのツールのほとんどがこれであり、実際はコマンドで使ってしまった方が使いやすい。

# gdal_translate
: gdal_translate [--help-general] [--long-usage]
  [-ot {Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/
        CInt16/CInt32/CFloat32/CFloat64}] [-strict]
  [-of format] [-b band] [-mask band] [-expand {gray|rgb|rgba}]
  [-outsize xsize[%] ysize[%]]
  [-unscale] [-scale[_bn] [src_min src_max [dst_min dst_max]]]* [-exponent[_bn] exp_val]*
  [-srcwin xoff yoff xsize ysize] [-projwin ulx uly lrx lry] [-epo] [-eco]
  [-a_srs srs_def] [-a_ullr ulx uly lrx lry] [-a_nodata value]
  [-gcp pixel line easting northing [elevation]]*
  [-mo "META-TAG=VALUE"]* [-q] [-sds]
  [-co "NAME=VALUE"]* [-stats] [-norat]
  src_dataset dst_dataset

文法はこんな感じである。
最後のsrc_datasetとdst_datasetの２つのファイルとオプションを指定する。

-of formatはファイル形式変換でこの種類には下に示すようになんでもある。
例えば、srtmのtifからdtedに変換するような場合は

gdal_translate -of dted -srcwin 0 1200 1201 1201 srtm_64_03.tif e130\n30.dt1

とかすればいい。(e130\n30はちょっとうろ覚えなんで計算してね）

上は1201なのでdt1だったが、srtmを圧縮して121ドットづつのdt0を作成するには
gdal_translate -of dted -srcwin 0 1200 1201 1201 -outsize 121 121 srtm_64_03.tif e130\n30.dt1

とすればいいわけである。

便利である。まあ、同じことでArcGisはお金をとっていたりする。

# gdalinfo
別のプログラムとしてgdalinfo.exeがあり、ファイル名を指定するとそのファイルの情報をダンプしてくれる。

# 対応しているformat
The following format drivers are configured and support output:
  VRT: Virtual Raster
  GTiff: GeoTIFF
  NITF: National Imagery Transmission Format
  HFA: Erdas Imagine Images (.img)
  ELAS: ELAS
  AAIGrid: Arc/Info ASCII Grid
  DTED: DTED Elevation Raster
  PNG: Portable Network Graphics
  JPEG: JPEG JFIF
  MEM: In Memory Raster
  GIF: Graphics Interchange Format (.gif)
  XPM: X11 PixMap Format
  BMP: MS Windows Device Independent Bitmap
  PCIDSK: PCIDSK Database File
  PCRaster: PCRaster Raster File
  ILWIS: ILWIS Raster Map
  SGI: SGI Image File Format 1.0
  SRTMHGT: SRTMHGT File Format
  Leveller: Leveller heightfield
  Terragen: Terragen heightfield
  ISIS2: USGS Astrogeology ISIS cube (Version 2)
  ERS: ERMapper .ers Labelled
  FIT: FIT Image
  RMF: Raster Matrix Format
  RST: Idrisi Raster A.1
  INGR: Intergraph Raster
  GSAG: Golden Software ASCII Grid (.grd)
  GSBG: Golden Software Binary Grid (.grd)
  GS7BG: Golden Software 7 Binary Grid (.grd)
  R: R Object Data Store
  PNM: Portable Pixmap Format (netpbm)
  ENVI: ENVI .hdr Labelled
  EHdr: ESRI .hdr Labelled
  PAux: PCI .aux Labelled
  MFF: Vexcel MFF Raster
  MFF2: Vexcel MFF2 (HKV) Raster
  BT: VTP .bt (Binary Terrain) 1.3 Format
  LAN: Erdas .LAN/.GIS
  IDA: Image Data and Analysis
  LCP: FARSITE v.4 Landscape File (.lcp)
  GTX: NOAA Vertical Datum .GTX
  NTv2: NTv2 Datum Grid Shift
  CTable2: CTable2 Datum Grid Shift
  KRO: KOLOR Raw
  ARG: Azavea Raster Grid format
  USGSDEM: USGS Optional ASCII DEM (and CDED)
  ADRG: ARC Digitized Raster Graphics
  BLX: Magellan topo (.blx)
  SAGA: SAGA GIS Binary Grid (.sdat)
  KMLSUPEROVERLAY: Kml Super Overlay
  XYZ: ASCII Gridded XYZ
  HF2: HF2/HFZ heightfield raster
  PDF: Geospatial PDF
  ZMap: ZMap Plus Grid