I'm looking to generate elevation profiles from a number of GPX traces (which I also have GeoJSON conversions of). I have 20m resolution DEMs of the area, in GeoTIFF format. Ultimately I want to generate a static image that combines a number of these elevation profiles, but I'll be happy if I can get something like a CSV of elevation values sampled at regular intervals.
Summarised, how can I turn:
- 1+ vector paths (.gpx or .geojson)
- 1 DEM (.tif)
- 1+ .svg (raw plot, no axes or labels); or
- 1+ .csv
using only command line tools. I have managed to achieve this manually using QGIS (and the Profile Tool plugin), but need to automate it. I'd prefer to avoid writing code - but Python is much preferable to R. A workflow that added elevations to the GeoJSON file would be better than nothing. I'm really not looking to implement a sampling algorithm from scratch.
OS X or Ubuntu are ok.
(Note to moderators: there seem to be lots of similar questions, but many of them are just "how do I get an elevation profile from this GPX file". The key differences are that I have the suitable terrain file already, and I need to automate this for dozens of files.)