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

[Joaquim Manuel Freire Luís]

With GMT that is a single command using the grdtrack module

https://docs.generic-mapping-tools.org/dev/grdtrack.html

-----Original Message-----
From: gdal-dev <gdal-dev-bounces at lists.osgeo.org> On Behalf Of Luca Bertoncello
Sent: Tuesday, December 28, 2021 5:48 PM
To: gdal dev <gdal-dev at lists.osgeo.org>
Subject: Re: [gdal-dev] Get max elevation on "squares"

Am 28.12.2021 um 16:45 schrieb Alan Snow:
> This resource is useful to help you get started exploring raster data 
> with Python:

Maybe there is other possibility than write a Python code? I really don't know Python and if I first need to learn it to extract the information I need, it will take a while...

I discovered this:
	https://www.eea.europa.eu/data-and-maps/data/digital-elevation-model-of-europe/

The resolution is enough for my purpose and the file size is OK.
So as I understand in this file every pixel has a size of 3x3 km, is it correct?

So my problem is now to get the "pixel coordinate" from latitude and longitude and extract the elevation from the TIF.

Am I right, until here?

If yes: can someone suggest me a way to extract the data from this file?
Best solution is a PERL or BASH script, no matter if it runs for hours or days, since I do the export just once.

Thanks
Luca Bertoncello
(lucabert at lucabert.de)
_______________________________________________
gdal-dev mailing list
gdal-dev at lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev