Python'da GDAL ile bir geotiff okuma, değiştirme ve yazma

Python GDAL bağları ve numpy kullanarak Uzaktan Algılama görüntü işlemenin iplerini öğrenmeye çalışıyorum. İlk denemede, Landsat8 geotiff dosyasını okuyorum, basit bir manipülasyon yapıyorum ve sonucu yeni bir dosyaya yazıyorum. Aşağıdaki kod, orijinal rasterin manipüle edilmiş raster yerine çıktı dosyasına dökülmesi dışında iyi çalışıyor gibi görünüyor.

Herhangi bir yorum veya öneri kabul edilir, ancak özellikle manipüle edilen rasterin sonuçta neden görünmediğine dair notlar.

import os
import gdal

gdal.AllRegister()

file = "c:\~\LC81980242015071LGN00.tiff"
(fileRoot, fileExt) = os.path.splitext(file)
outFileName = fileRoot + "_mod" + fileExt

ds = gdal.Open(file)
band = ds.GetRasterBand(1)
arr = band.ReadAsArray()

[cols, rows] = arr.shape
arr_min = arr.Min()
arr_max = arr.Max()
arr_mean = int(arr.mean())

arr_out = numpy.where((arr < arr_mean), 10000, arr)

driver = gdal.GetDriverByName("GTiff")
outdata = driver.Create(outFileName, rows, cols, 1, gdal.GDT_UInt16)
outband = outdata.GetRasterBand(1)
outband.WriteArray(arr_out)
outdata = None

print arr_min
> 0
print arr_max
> 65535
print arr_mean
> 4856

Windows 7 32 bit makinede Python 2.7.1 kullanıyorum.

Komut dosyanızda, değişikliklerin sonunda belleğe sahip olduklarınızı diske kaydeden ds.FlushCache yöntemi eksik. Aşağıda örneğinizin düzeltilmiş bir sürümüne bakın. Projeksiyon ve geotransform'u girdi olarak ayarlamak için iki satır eklediğime dikkat edin

import os
import gdal

file = "path+filename"
ds = gdal.Open(file)
band = ds.GetRasterBand(1)
arr = band.ReadAsArray()
[cols, rows] = arr.shape
arr_min = arr.min()
arr_max = arr.max()
arr_mean = int(arr.mean())
arr_out = numpy.where((arr < arr_mean), 10000, arr)
driver = gdal.GetDriverByName("GTiff")
outdata = driver.Create(outFileName, rows, cols, 1, gdal.GDT_UInt16)
outdata.SetGeoTransform(ds.GetGeoTransform())##sets same geotransform as input
outdata.SetProjection(ds.GetProjection())##sets same projection as input
outdata.GetRasterBand(1).WriteArray(arr_out)
outdata.GetRasterBand(1).SetNoDataValue(10000)##if you want these values transparent
outdata.FlushCache() ##saves to disk!!
outdata = None
band=None
ds=None