Dato un DEM (Digital Elevation Model) e una linea (x1, y1) -> (x2, y2) (coordinate accoppiate (x, y)), come si può calcolare il profilo di elevazione o la sezione trasversale proiettata su quella linea?
Sto cercando un codice open source che posso usare nel mio progetto o pseudo-codice che posso usare come guida.
Risposte:
Il comando GRASS r.profileesegue questo ( documentazione , fonte ) e dovrebbe fornire una buona base per l'implementazione di una sezione trasversale ed è disponibile sotto GPL.
Forse il codice sorgente del plugin "Profilo" di QGIS può essere un inizio.
Puoi andare molto oltre combinando i colori DEM e raster, usando GRASS con R. Vedi:
Vedi alcuni esempi di seguito:
Con i colori di elevazione (GRASS GIS e R):
Con i colori di una mappa geologica (Python con raster GRASS):
E usando lo script Python nella console Python di QGIS:
Ecco schermate per GRASS: 
