概要
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