Determining the Unit of Spatial Resolution of a raster image in Python

Question

I have some raster images, and I would like to determine the GSD using Python.

I know that I can use rasterio.transform() or gdal.GetGeoTransform() to extract the spatial resolution of the image.

How do I determine the units of the spatial resolution.

I don't think I've seen this called GSD before, and for anyone else it stands for "Ground Sampling Distance". (My default for GSD is German Shepherd Dog) — Spacedman, Nov 08 '20 at 11:32
Sorry if it wasn't clear enough I updated the question. I was hoping there would be a way to use the spatial resolution to determine the GSD. — A Merii, Nov 08 '20 at 14:15
Thanks a lot, in this case I will rephrase the question so that others can benefit. — A Merii, Nov 08 '20 at 15:45
@Aaron, If you can, re-upload your answer, so that I can accept it as the accepted answer then that would be great! — A Merii, Nov 08 '20 at 15:49

score 4 · Accepted Answer · answered Nov 08 '20 at 04:36

4

You can use the linear_units property in the crs class. For example:

import rasterio
raster = rasterio.open('/path/to/your/raster.tif' )
print(raster.crs.linear_units)

answered Nov 08 '20 at 04:36

Aaron

1

What if this just sayis 'unknown' ? :( Showing this rater in QGis, for example, knows where to put each pixel, so surely it can calculate GSD, no? – Ciprian Tomoiagă Jul 30 '21 at 14:16
Is there an equivalent for GDAL? – dazonic Jan 25 '23 at 11:56
@dazonic This may be helpful: https://gis.stackexchange.com/a/368008/8104 – Aaron Jan 25 '23 at 15:55
@Aaron that only get the x and y res in numbers, I’m looking for a way to get units, ie degrees or metres – dazonic Feb 07 '23 at 11:33

1 Answers1