Como chamar gdal_translate a partir do código Python?

É possível usar alguma API gdal para chamar gdal_translatedo código Python? Não quero dizer simplesmente executar o gdal_translate.exe no sistema de arquivos, mas chamá-lo de alguma forma no código, para que eu não precise saber o diretório exato em que o executável do gdal_translate está?

Desde o GDAL 2.1 (mais informações aqui ), os utilitários GDAL e OGR podem ser usados ​​como funções de biblioteca. Por exemplo:

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

Veja o Tutorial da 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 você deseja mais controle de saída, como redimensionamento, subconjunto, etc ... use uma VRT como entrada, é assim que o gdal_translate faz isso internamente.

Sim, você pode chamar os Utilitários GDAL de dentro do Python. Existem diferenças muito pequenas na abordagem, dependendo se o utilitário é um exe por si só ou também um trecho de código python. De qualquer forma, você precisa usar o módulo de subprocesso :

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)

Você notará que adiciono aspas aos meus caminhos. Isso ocorre porque, no Windows, tive problemas com caminhos, especialmente com espaços ou onde um dos caracteres '\' cria outro caractere de escape acidental. Portanto, apenas preservo o caminho correto no aspec por assim dizer.

Se você estiver usando um dos utilitários python, faça o mesmo, exceto que o exe no início da cadeia de comando do subprocesso agora é "C: \ python32 \ python.exe" (ou qualquer versão que você tenha) e seu segundo elemento é o utilitário python que você deseja usar.

Obviamente, você também pode interagir com o sistema de arquivos em vez de usar constantes codificadas, mas este é apenas um exemplo.

EDIT - Generalizando para plugins
QGIS O QGIS cria / modifica várias variáveis ​​de ambiente na inicialização. Portanto, você pode criar variáveis ​​de caminho generalizadas para as bibliotecas / utilitários GDAL usando estas (consulte Configurações-> Opções-> Sistema) em vez dos caminhos codificados no exemplo acima.

Eu faço isso com vários comandos gdal usando os.system, que você pode usar para chamar funções, assim como na linha de comando:

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

Também está descrito na aula 7 aqui: http://www.gis.usu.edu/~chrisg/python/2009/

Aqui está um código rápido para quem deseja salvar bandas de um TIF composto de várias bandas em arquivos individuais usando o GDAL Translate em 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

Isso pode ser útil para processamento adicional (por exemplo, usando o Rasterio, como aqui ).