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

2021-12-29 Thread Luca Bertoncello
Am 29.12.2021 um 11:42 schrieb Andrew C Aitchison: > Is this the max or the average height of the square ? Mmmm... Now I'm unsure... > Do you need to warn of pylons, skyscrapers, bridges or radio masts ? > These are unlikely to be included in such map data. None of them. Gound elevation > I no

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

2021-12-29 Thread Andrew C Aitchison
On Tue, 28 Dec 2021, Luca Bertoncello wrote: 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 ex

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

2021-12-29 Thread Luca Bertoncello
Am 29.12.2021 um 10:30 schrieb Volker Wichmann: Hi Volker > I haven't followed the discussion in detail, but you could use > > https://gdal.org/programs/gdallocationinfo.html > > to read out the values from your resampled raster. I did that. And got a little Python script to export the data as

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

2021-12-29 Thread Volker Wichmann
Hi Luca, I haven't followed the discussion in detail, but you could use https://gdal.org/programs/gdallocationinfo.html to read out the values from your resampled raster. The elevation differences between high and coarse resolution DEMs strongly depend on the relief in the area of investigati

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

2021-12-28 Thread Joaquim Manuel Freire Luís
rid tile (earth_relief_03m_p): S90W180 -9.1409 38.7718 87.2690486581 -Original Message- From: gdal-dev On Behalf Of Joaquim Manuel Freire Luís Sent: Tuesday, December 28, 2021 6:00 PM To: Luca Bertoncello ; gdal dev Subject: Re: [gdal-dev] Get max elevation on "squares" With GM

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

2021-12-28 Thread Luca Bertoncello
Am 28.12.2021 um 04:54 schrieb Alan Snow: > After some thought, you may achieve what you want by resampling the file > to a 5km grid size with the Resampling.max option: > > https://rasterio.readthedocs.io/en/latest/topics/resampling.html > > https://rasterio.readthedocs.io/en/latest/api/rasterio

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

2021-12-28 Thread 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 On Behalf Of Luca Bertoncello Sent: Tuesday, December 28, 2021 5:48 PM To: gdal dev Subject: Re: [gdal-dev] Get max elevation on

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

2021-12-28 Thread Luca Bertoncello
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 w

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

2021-12-28 Thread Alan Snow
This resource is useful to help you get started exploring raster data with Python: https://carpentries-incubator.github.io/geospatial-python/05-raster-structure/index.html https://www.earthdatascience.org/courses/use-data-open-source-python/intro-raster-data-python/raster-data-processing/ If th

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

2021-12-27 Thread Luca Bertoncello
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

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

2021-12-27 Thread Alan Snow
After some thought, you may achieve what you want by resampling the file to a 5km grid size with the Resampling.max option: https://rasterio.readthedocs.io/en/latest/topics/resampling.html https://rasterio.readthedocs.io/en/latest/api/rasterio.enums.html#rasterio.enums.Resampling.max You can the

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

2021-12-27 Thread Luca Bertoncello
Am 27.12.2021 um 20:20 schrieb Alan Snow: Hi Alan, thank you very much for your answer! > You may be looking for clip/mask: > > https://rasterio.readthedocs.io/en/latest/topics/masking-by-shapefile.html > > https://corteva.github.io/rioxarray/stable/examples/clip_geom.html > > After clipping, yo

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

2021-12-27 Thread Alan Snow
You may be looking for clip/mask: https://rasterio.readthedocs.io/en/latest/topics/masking-by-shapefile.html https://corteva.github.io/rioxarray/stable/examples/clip_geom.html After clipping, you can use numpy to do the sum. On Mon, Dec 27, 2021, 12:58 PM Luca Bertoncello wrote: > Hi list! >

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

2021-12-27 Thread Luca Bertoncello
Hi list! I'm very new to GDAL, so please excuse me if the question is too simple... My problem: I'm writing an App to fly VFR and I'd like to display a vertical profile of the route. Of course, it must work offline, so I cannot use a service like Google Maps or similar... I tryed OpenTopoData wi