GIS：GDAL实现对栅格文件的转换

我们常常在图像处理过程中遇到不同软件或程序要求输入的图像格式不同（有些程序或软件支持的数据格式不是常用的Tiff，Img等数据格式），因此需要对不同的数据格式相互进行转换。 我这里以GTiff(.tif)数据转换为PCRaster(.map)数据为例。首先需要安装GDAL，我这里是在Anaconda上直接安装了基于Python的GDAL，可以在下面网站自行下载，https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal 例如下面对应的就是Python3.8版本的GDAL。

安装好后可在Anaconda的Prompt中直接使用gdal_translate.exe直接实现转换，转换方法如下：

gdal_translate [--help-general]
    [-ot {
    
    Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/
            CInt16/CInt32/CFloat32/CFloat64}] [-strict]
    [-if format]* [-of format]
    [-b band]* [-mask band] [-expand {
    
    gray|rgb|rgba}]
    [-outsize xsize[%]|0 ysize[%]|0] [-tr xres yres]
    [-r {
    
    nearest,bilinear,cubic,cubicspline,lanczos,average,rms,mode}]
    [-unscale] [-scale[_bn] [src_min src_max [dst_min dst_max]]]* [-exponent[_bn] exp_val]*
    [-srcwin xoff yoff xsize ysize] [-epo] [-eco]
    [-projwin ulx uly lrx lry] [-projwin_srs srs_def]
    [-a_srs srs_def] [-a_ullr ulx uly lrx lry] [-a_nodata value]
    [-a_scale value] [-a_offset value]
    [-nogcp] [-gcp pixel line easting northing [elevation]]*
    |-colorinterp{
    
    _bn} {
    
    red|green|blue|alpha|gray|undefined}]
    |-colorinterp {
    
    red|green|blue|alpha|gray|undefined},...]
    [-mo "META-TAG=VALUE"]* [-q] [-sds]
    [-co "NAME=VALUE"]* [-stats] [-norat] [-noxmp]
    [-oo NAME=VALUE]*
    src_dataset dst_dataset

根据需求，输入相关的参数进行转换，不需要的参数可以忽略。主要参数说明（其他参数详细见GDAL官方文档 gdal_translate — GDAL 文档）：

-ot  强制输出图像带具有驱动程序支持的特定数据类型，该数据类型可以是以下类型之一： Byte ， UInt16 ， Int16 ， UInt32 ， Int32 ， Float32 ， Float64 ， CInt16 ， CInt32 ， CFloat32 或 CFloat64 .
-if  试图打开输入文件的格式/驱动程序名称。通常不需要指定它，但当它无法选择适当的驱动程序时，可以使用它跳过自动驱动程序检测。
-of  选择输出格式。从GDAL 2.3开始，如果未指定，则从扩展名猜测格式。
-b 选择输入波段 band 输出,从1开始编号。
-mask   选择输入波段 band 创建输出数据集掩码带区。
-tr 设定目标分辨率。这些值必须用地理参考单位表示。两者都必须是正值。
-r 采样方式，有nearest (default),bilinear,cubic,cubicspline,lanczos,average,rms,mode
<src_dataset>源数据集名称。它可以是文件名、数据源的URL或多数据集文件的子数据集名称。
<dst_dataset>目标文件名。

支持的数据格式有：

 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
  SRTMHGT: SRTMHGT File Format
  Leveller: Leveller heightfield
  Terragen: Terragen heightfield
  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
  IDA: Image Data and Analysis
  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)
  USGSDEM: USGS Optional ASCII DEM (and CDED)
  ADRG: ARC Digitized Raster Graphics