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@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-boun...@lists.osgeo.org> On Behalf Of Joaquim Manuel Freire Luís Sent: Tuesday, December 28, 2021 6:00 PM To: Luca Bertoncello <lucab...@lucabert.de>; gdal dev <gdal-dev@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-boun...@lists.osgeo.org> On Behalf Of Luca Bertoncello Sent: Tuesday, December 28, 2021 5:48 PM To: gdal dev <gdal-dev@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 (lucab...@lucabert.de) _______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev _______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev _______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev
