How can I change the resolution of a raster using GDAL?

python gdal
James picture James · Oct 11, 2011 · Viewed 12.3k times · Source

I am looking for the best way to change the resolution of a GDAL raster dataset.

For example, I have a raster that has a pixel size of (30, -30), and I would like to change the pixel size to (5, -5), interpolating all values for a given pixel into the output raster.

So for each pixel of the input raster, I would like to have 36 pixels in the output raster that all share the same value.

If I run gdalwarp -tr 5 -5 inputRaster.tif outputRaster.tif, I get exactly the result that I'm looking for, and so I would assume that I should be able to replicate this functionality with some GDAL function.

I would prefer to avoid using a call to python's Subprocess class, if possible.

Answer

Mike T picture Mike T · Oct 16, 2011

You need to reproject the raster. For example, from an interactive Python shell:

from osgeo import gdal
help(gdal.ReprojectImage)

A Python example is provided in the test suite.

More complete documentation is provided for the C++ function GDALReprojectImage.

Related questions

Numpy - Replace a number with NaN

I am looking to replace a number with NaN in numpy and am looking for a function like numpy.nan_to_num, except in reverse. The number is likely to change as different arrays are processed because each can have …

Read More
python GDAL 2.1 installation on Ubuntu 16.04

Here will be my sequence of command lines while trying to install gdal2.1 in a UBUNTU virtual machine. My virtual machine is a UBUNTU 16.04 LTS(64bit) I would need gdal2.1 and especially the Python bindings to work with it in …

Read More
Using pip to install modules in python failing

I'm having trouble installing python modules using pip. Below is the output from the command window: Note that I installed pip immediately before trying to install GDAL module. I am on a w7 64bit machine running python 2.7 Microsoft Windows [Version 6.1.7601] …

Read More