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 09: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

score 9 · Answer 1 · answered Jul 17 '20 at 09:45

9

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 09:45

Zoltan

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

BERA · Accepted Answer · 2020-07-17T10:11:47.860

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())

edited Jul 17 '20 at 10:11

answered Jul 17 '20 at 09:52

BERA

score 0 · Answer 3 · edited Jul 12 '23 at 11:05

0

This worked for me in gdal version 3.4.1:

src = gdal.Open(mask)
xres = src.RasterXSize
yres = src.RasterYSize

edited Jul 12 '23 at 11:05

Kadir Şahbaz

answered Jul 12 '23 at 11:04

Zipzap

3 Answers3