1 votos

Python GDAL cambiando geotransformación sin abrir como array numpy

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)

1voto

Joe Puntos 16

Esto debería ser suficiente para cambiar la geotransformación (capturada de los autotests de GDAL https://github.com/OSGeo/gdal/blob/master/autotest/gcore/tiff_write.py ):

ds = gdal.Open('tmp/tiff57.tif', gdal.GA_Update)
ds.SetGeoTransform([100, 1, 3, 200, 3, 1])
ds = None

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X