Comment appeler gdal_translate à partir de code Python?

Est-il possible d'utiliser une API gdal pour appeler à gdal_translatepartir de code Python? Je ne parle pas simplement d'exécuter gdal_translate.exe à partir du système de fichiers, mais de l'appeler d'une manière ou d'une autre dans le code, de sorte que je n'ai pas besoin de connaître le répertoire exact dans lequel se trouve l'exécutable gdal_translate?

Depuis GDAL 2.1 (plus d’informations ici ), les utilitaires GDAL et OGR peuvent être utilisés comme fonctions de bibliothèque. Par exemple:

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

Voir le didacticiel de l'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

Si vous souhaitez davantage de contrôle de sortie, tel que redimensionnement, sous-ensemble, etc., utilisez un VRT comme entrée, voici comment gdal_translate le fait en interne.

Oui, vous pouvez appeler les utilitaires GDAL à partir de Python. Il y a de très petites différences dans l'approche selon que l'utilitaire est un exe à part entière ou aussi un morceau de code python. Quoi qu'il en soit, vous devez utiliser le module de sous - processus :

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)

Vous remarquerez que j'ajoute des guillemets échappés autour de mes chemins. En effet, sous Windows, j'ai eu des problèmes avec les chemins, en particulier ceux avec des espaces ou lorsqu'un des caractères "\" crée un autre caractère échappé par accident. Donc, je conserve juste le chemin approprié dans aspec pour ainsi dire.

Si vous utilisez un des utilitaires python, faites simplement la même chose, sauf que votre fichier exe au début de la chaîne de commande du sous-processus est maintenant "C: \ python32 \ python.exe" (ou quelle que soit la version que vous avez) et que votre deuxième élément est l'utilitaire Python que vous souhaitez utiliser.

Bien entendu, vous pouvez également parcourir votre système de fichiers plutôt que d'utiliser des constantes codées en dur, mais ceci n'est qu'un exemple.

EDIT - Généraliser pour les plugins
QGIS QGIS crée / modifie un certain nombre de variables d'environnement au démarrage. Vous pouvez donc créer des variables de chemin d'accès généralisées aux bibliothèques / utilitaires GDAL à l'aide de celles-ci (voir Paramètres-> Options-> Système) au lieu des chemins codés en dur de l'exemple ci-dessus.

Je le fais avec différentes commandes gdal utilisant os.system que vous pouvez utiliser pour appeler des fonctions à partir de la ligne de commande:

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

Il est également décrit dans la conférence 7 ici: http://www.gis.usu.edu/~chrisg/python/2009/

Voici un code rapide pour ceux qui souhaitent enregistrer des groupes d’un fichier TIF composite multibande composite dans des fichiers individuels à l’aide de GDAL Traduire 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

Cela pourrait être utile pour un traitement ultérieur (par exemple, en utilisant Rasterio, comme ici ).