How to get raster resolution using GDAL

Question

I'm using gdal lib for Python in a Jupyter Notebook environment.

I'd like to retrieve the resolution of a raster to use it in a gdal.warp call.

At the moment I do the following :

src = gdal.Open(mask)
ulx, xres, xskew, uly, yskew, yres  = src.GetGeoTransform()

but I get too much information. Is there a way to only retreive xres and yres?

because I really don't like to instantiate useless variables. I try to clarify and organize a working code ;-) — Pierrick Rambaud, Jul 17 '20 at 9:35
Not adding as an answer as it's ugly, but you could use slice notation: xres, yres = src.GetGeoTransform()[1:6:4] — user2856, Jul 17 '20 at 10:06

BERA · Accepted Answer · 2020-07-17 10:11:47Z

3

Google says no. But a "workaround" could be to use operator.itemgetter:

import gdal, operator
src = gdal.Open(mask)
xres, yres = operator.itemgetter(1,5)(src.GetGeoTransform())

answered Jul 17 '20 at 9:52

BERA

35.6k10 gold badges37 silver badges77 bronze badges

Add a comment |

Zoltan · Accepted Answer · 2020-07-17 09:45:08Z

4

If you only want to avoid the unused variables use this code:

import gdal
...
src = gdal.Open(mask)
_, xres, _, _, _, yres  = src.GetGeoTransform()

answered Jul 17 '20 at 9:45

Zoltan

6,37713 silver badges22 bronze badges

Add a comment |

2 Answers 2