Étant donné un DEM (modèle numérique d'élévation) et une ligne (x1, y1) -> (x2, y2) (coordonnées appariées (x, y)), comment peut-on calculer le profil d'élévation ou la coupe projetée sur cette ligne?
Je recherche un code open source que je peux utiliser dans mon projet ou un pseudo-code que je peux utiliser comme guide.
Réponses:
La commande GRASS r.profileeffectue cela ( documentation , source ) et devrait fournir une bonne base pour la mise en œuvre d'une coupe transversale, et est disponible sous la GPL.
Peut-être que le code source du plugin QGIS "Profile" peut être un début.
Vous pouvez aller beaucoup plus loin en combinant les couleurs DEM et raster, en utilisant GRASS avec R. Voir:
Voir quelques exemples ci-dessous:
Avec couleurs d'élévation (GRASS GIS et R):
Avec les couleurs d'une carte géologique (Python avec raster GRASS):
Et en utilisant le script Python dans la console Python de QGIS:
Voici des captures d'écran pour GRASS: 
