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

[Joaquim Manuel Freire Luís]

This is for the Lisbon Airport.
The first time it downloads the grid(s). Here something is probably failing because it shouldn't need to download all tiles that make up the global grid (8 tiles and 30 Mb for this resolution)

echo -9.1409 38.7718 | grdtrack -G at earth_relief_03m
grdblend [NOTICE]: Remote data courtesy of GMT data server oceania [http://oceania.generic-mapping-tools.org]

grdblend [NOTICE]: Earth Relief at 3x3 arc minutes from Gaussian Cartesian filtering (5.6 km fullwidth) of SRTM15+V2.1 [Tozer et al., 2019].
grdblend [NOTICE]:   -> Download 90x90 degree grid tile (earth_relief_03m_p): S90W180

-9.1409 38.7718 87.2690486581

-----Original Message-----
From: gdal-dev <gdal-dev-bounces at lists.osgeo.org> On Behalf Of Joaquim Manuel Freire Luís
Sent: Tuesday, December 28, 2021 6:00 PM
To: Luca Bertoncello <lucabert at lucabert.de>; gdal dev <gdal-dev at lists.osgeo.org>
Subject: Re: [gdal-dev] Get max elevation on "squares"

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
_______________________________________________
gdal-dev mailing list
gdal-dev at lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev