Come chiamare gdal_translate dal codice Python?

È possibile utilizzare alcune API gdal per chiamare gdal_translatedal codice Python? Non intendo semplicemente eseguire gdal_translate.exe dal file system, ma piuttosto chiamarlo in qualche modo nel codice, quindi non ho bisogno di conoscere la directory esatta in cui si trova l'eseguibile gdal_translate?

A partire da GDAL 2.1 (maggiori informazioni qui ), le utilità GDAL e OGR possono essere utilizzate come funzioni di libreria. Per esempio:

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

Consulta il tutorial sull'API 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

Se vuoi un maggiore controllo dell'output, come ridimensionamento, subsetting, ecc ... usa un VRT come input, ecco come gdal_translate lo fa internamente.

Sì, puoi chiamare le utility GDAL da Python. Ci sono differenze minime nell'approccio a seconda che l'utilità sia un exe a sé stante o anche un pezzo di codice Python. In ogni caso, però, è necessario utilizzare il modulo di sottoprocesso :

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)

Noterai che aggiungo virgolette sfuggite attorno ai miei percorsi. Questo perché, su Windows, ho avuto problemi con i percorsi, specialmente quelli con spazi o in cui uno dei caratteri "\" crea un altro carattere di escape accidentale. Quindi, preservo il percorso corretto in Aspec per così dire.

Se stai usando una delle utility di Python, fai la stessa cosa tranne che il tuo exe all'inizio della stringa di comando del sottoprocesso ora è "C: \ python32 \ python.exe" (o qualunque versione tu abbia) e il tuo secondo elemento è l'utilità python che si desidera utilizzare.

Ovviamente puoi anche iterare sul tuo file system piuttosto che usare costanti codificate, ma questo è solo un esempio.

EDIT - Generalizzando per i plugin
QGIS QGIS crea / modifica una serie di variabili d'ambiente all'avvio. Quindi, puoi costruire variabili di percorso generalizzate alle librerie / utilità GDAL usando queste (vedi Impostazioni-> Opzioni-> Sistema) invece dei percorsi codificati nell'esempio sopra.

Faccio questo con vari comandi gdal usando os.system che puoi usare per chiamare le funzioni proprio come dalla riga di comando:

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

È anche descritto nella lezione 7 qui: http://www.gis.usu.edu/~chrisg/python/2009/

Ecco un codice rapido per chiunque desideri salvare bande da un TIF multi-banda composito in singoli file utilizzando GDAL Translate in 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

Questo potrebbe essere utile per ulteriori elaborazioni (ad es. Utilizzando Rasterio, come qui ).