Python kodundan gdal_translate nasıl çağırılır?

gdal_translatePython kodundan çağırmak için bir gdal API kullanmak mümkün müdür ? Sadece gdal_translate.exe dosyasını dosya sisteminden çalıştırmak demek istemiyorum, ancak kodunu bir şekilde kod olarak adlandırın, böylece gdal_translate yürütülebilir dosyasının bulunduğu dizini bilmem gerekmez mi?

GDAL 2.1'den beri ( burada daha fazla bilgi var ), GDAL ve OGR programları kütüphane işlevi olarak kullanılabilir. Örneğin:

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

GDAL API Eğitimine bakınız .

#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

Yeniden boyutlandırma, alt kümeleme, vb. Gibi daha fazla çıktı kontrolü istiyorsanız, girdi olarak VRT kullanın , gdal_translate bunu dahili olarak yapar.

Evet, GDAL Yardımcı Programlarını Python içinden arayabilirsiniz. Yardımcı programın kendi başına bir exe veya bir python kodu parçası olmasına bağlı olarak yaklaşımda çok küçük farklılıklar vardır. Her iki durumda da alt işlem modülünü kullanmanız gerekse de :

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)

Yollarımın üzerine kaçış tırnak işaretleri eklediğimi fark edeceksiniz. Bunun nedeni, Windows’ta, özellikle boşluk içeren veya '\' karakterlerinden birinin yanlışlıkla kaçan bir karakter çıkardığı yollar ile ilgili sorun yaşıyorum. Bu yüzden, sadece aspec'teki doğru yolu korudum.

Python yardımcı programlarından birini kullanıyorsanız, alt işlem komut dizesinin başlangıcındaki exe'niz dışında şimdi aynı şeyleri yapın "C: \ python32 \ python.exe" (veya hangi sürümüne sahipseniz) ve ikinci öğeniz kullanmak istediğiniz python yardımcı programını kullanın.

Açıkçası sabit kodlanmış sabitleri kullanmak yerine dosya sisteminizi yineleyebilirsiniz, ancak bu sadece bir örnek.

EDIT - QGIS eklentileri için genelleme QGIS
, başlangıçta bir dizi çevre değişkeni oluşturur / değiştirir. Böylece, yukarıdaki örnekte kodlanmış yollar yerine bunları (bkz. Ayarlar-> Seçenekler-> Sistem) kullanarak GDAL kitaplıklarına / yardımcı programlarına genelleştirilmiş yol değişkenleri oluşturabilirsiniz.

Bunu komut satırından itibaren fonksiyon çağırmak için kullanabileceğiniz os.system kullanarak çeşitli komutlarla yapıyorum:

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

Ayrıca burada 7. derste açıklanmıştır: http://www.gis.usu.edu/~chrisg/python/2009/

Python'da GDAL Translate'i kullanarak kompozit çok bantlı TIF'ten bireysel dosyalara bant kaydetmek isteyen herkes için hızlı bir kod.

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

Bu daha sonraki işlemler için faydalı olabilir (örneğin, burada olduğu gibi Rasterio kullanarak ).