[gdal-dev] Get max elevation on "squares"

Luca Bertoncello lucabert at lucabert.de
Mon Dec 27 22:56:27 PST 2021


Am 28.12.2021 um 04:54 schrieb Alan Snow:

Hi Alan

> After some thought, you may achieve what you want by resampling the file
> to a 5km grid size with the Resampling.max option:

This is a good idea... I will need to write some code in Python (I don't
really know this language), but I think I'll get it...

> You can then store the centroid of the grid cell and the value in a
> database.

Yes, I already see with gdalinfo -mm (maybe there is a possibility to
check it with Python?) I can get the maximum elevation and the centroid.

Just some question, to correctly understand the code sample:

1) I think, dataset.read get the size as parameter, correct? Are there
meter/kilometer/something other?
I searched some information about DatasetReader.read but I didn't found
anything...
2) In the example the new image will not be saved, is it correct? How
can I do that? I didn't found any function DatasetReader.save...
3) Maybe is it possible to do that in BASH? As I sayd, I'm not really
fit in Python... :(

Thanks a lot
Luca Bertoncello
(lucabert at lucabert.de)


More information about the gdal-dev mailing list