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, 2020 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, 2020 at 10:06

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

5

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, 2020 at 9:45

Zoltan

6,72215 silver badges23 bronze badges

any drawbacks of that solution? possibilities that it won't work? this should be an answer.
– Flash Thunder
Apr 20 at 17:43
I does work using Python 3.8.10 and gdal 3.1.4.
– Zoltan
Apr 20 at 19:54

Add a comment |

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, 2020 at 9:52

BERA

51.2k13 gold badges47 silver badges106 bronze badges

Add a comment |

2 Answers 2