Mi raster está al revés, pero es demasiado grande para abrirlo como array y reescribirlo con otra geotransformación. Conoces alguna forma de cambiar la geotransformación sin abrir el raster como un array numpy?
mi enfoque normal
import gdal
import numpy
r = gdal.Open(fileLocation)
band = r.GetRasterBand(1) # bands start at one
b = band.ReadAsArray().astype(np.float16)
dst_ds = driver.Create(output_file, band.XSize, band.YSize, 1, gdal.GDT_Int16)
dst_ds.SetGeoTransform(geotransform)
dst_ds.GetRasterBand(1).WriteArray(a)