¿Cómo puedo cambiar la resolución de un ráster usando GDAL?

Question

Estoy buscando la mejor manera de cambiar la resolución de un dataset ráster GDAL.

Por ejemplo, tengo un ráster que tiene un tamaño de píxel de (30, -30), y me gustaría cambiar el tamaño del píxel a (5, -5), interpolando todos los valores para un píxel dado en el ráster de salida

Por lo tanto, para cada píxel del ráster de entrada, me gustaría tener 36 píxeles en el ráster de salida que comparten el mismo valor.

Si ejecuto gdalwarp -tr 5 -5 inputRaster.tif outputRaster.tif, obtengo exactamente el resultado que estoy buscando, así que supongo que debería ser capaz de replicar esta funcionalidad con alguna función GDAL.

Preferiría evitar usar una llamada a la clase Subprocess de python, si es posible.

Answer 1

Necesita reproyectar el ráster. Por ejemplo, de una cáscara de Python interactivo: se proporciona

from osgeo import gdal 
help(gdal.ReprojectImage) 

ejemplo un pitón en el test suite.

Se proporciona documentación más completa para la función C++ GDALReprojectImage.