¿Cómo llamar a gdal_translate desde el código de Python?

¿Es posible usar alguna API de gdal para llamar gdal_translatedesde código Python? No me refiero simplemente a ejecutar gdal_translate.exe desde el sistema de archivos, sino a llamarlo de alguna manera en el código para que no necesite saber el directorio exacto en el que se encuentra el ejecutable gdal_translate.

Desde GDAL 2.1 (más información aquí ), las utilidades GDAL y OGR pueden usarse como funciones de biblioteca. Por ejemplo:

from osgeo import gdal

ds = gdal.Open('input.tif')
ds = gdal.Translate('output.tif', ds, projWin = [-75.3, 5.5, -73.5, 3.7])
ds = None

Consulte el tutorial de la API de GDAL .

#Import gdal
from osgeo import gdal

#Open existing dataset
src_ds = gdal.Open( src_filename )

#Open output format driver, see gdal_translate --formats for list
format = "GTiff"
driver = gdal.GetDriverByName( format )

#Output to new format
dst_ds = driver.CreateCopy( dst_filename, src_ds, 0 )

#Properly close the datasets to flush to disk
dst_ds = None
src_ds = None

Si desea más control de salida, como cambiar el tamaño, subconjunto, etc., use un VRT como entrada, así es como gdal_translate lo hace internamente.

Sí, puede llamar a las Utilidades GDAL desde Python. Existen diferencias muy pequeñas en el enfoque dependiendo de si la utilidad es un exe por derecho propio o también una pieza de código python. De cualquier manera, aunque necesita usar el módulo de subproceso :

import subprocess

# constants
gdalTranslate = r'C:\Program Files\GDAL\gdal_translate.exe'
src = r"C:\somefolder\somefile.tif"
dst = r"C:\someotherfolder\myresul.tif"
cmd = "-ot float32 -outsize 25 25"  # just for example!

# see note below
def youCanQuoteMe(item):
    return "\"" + item + "\""

fullCmd = ' '.join([gdalTranslate, cmd, youCanQuoteMe(src), youCanQuoteMe(dst)])
subprocess.popen(fullCmd)

Notarás que agrego comillas escapadas alrededor de mis caminos. Esto se debe a que, en Windows, he tenido problemas con las rutas, especialmente las que tienen espacios o donde uno de los caracteres '\' genera otro carácter de escape accidental. Entonces, solo conservo el camino correcto en aspec por así decirlo.

Si está utilizando una de las utilidades de python, simplemente haga lo mismo, excepto que su exe al comienzo de la cadena de comando del subproceso ahora es "C: \ python32 \ python.exe" (o la versión que tenga) y su segundo elemento es la utilidad python que quieres usar.

Obviamente, también puede iterar sobre su sistema de archivos en lugar de usar constantes codificadas, pero esto es solo un ejemplo.

EDITAR: generalización de los complementos
QGIS QGIS crea / modifica una serie de variables de entorno al inicio. Por lo tanto, puede crear variables de ruta generalizadas para las bibliotecas / utilidades GDAL utilizando estas (consulte Configuración-> Opciones-> Sistema) en lugar de las rutas codificadas en el ejemplo anterior.

Hago esto con varios comandos de gdal usando os.system que puedes usar para llamar funciones como desde la línea de comandos:

os.system("gdal_translate -of GTiff " + sourcefile + " " +  destinationfile)

También se describe en la clase 7 aquí: http://www.gis.usu.edu/~chrisg/python/2009/

Aquí hay un código rápido para cualquier persona que desee guardar bandas de un TIF compuesto de múltiples bandas en archivos individuales utilizando GDAL Translate en Python.

import gdal

in_path = 'C:/GIS/Sample.tif' #input composite raster
out_path = 'C:/GIS/Output/' #output directory for individual bands as files

#Open existing raster ds
src_ds = gdal.Open(in_path)

for i in range(1,src_ds.RasterCount +1): #Save bands as individual files
    out_ds = gdal.Translate(out_path + 'band' + str(i) + '.tiff', src_ds, format='GTiff', bandList=[i])
    out_ds=None

Esto podría ser útil para un procesamiento posterior (por ejemplo, usando Rasterio, como aquí ).