Cálculo del perfil de elevación a lo largo de la línea desde un DEM?

Dado un DEM (Modelo de elevación digital) y una línea (x1, y1) -> (x2, y2) (coordenadas emparejadas (x, y)), ¿cómo se puede calcular el perfil de elevación o la sección transversal proyectada en esa línea?

Estoy buscando un código fuente abierto que pueda usar en mi proyecto o pseudocódigo que pueda usar como guía.

El comando GRASS r.profilerealiza esto ( documentación , fuente ) y debe proporcionar una buena base para implementar una sección transversal, y está disponible bajo la GPL.

Quizás el código fuente del complemento "Perfil" de QGIS puede ser un comienzo.

Puede ir mucho más allá combinando DEM y colores de trama, utilizando GRASS con R. Consulte:

• Perfil ráster a lo largo de segmentos de línea arbitrarios , o
• http://www.portailsig.org/content/grass-gis-et-r-superposition-des-couleurs-d-un-raster-quelconque-sur-un-profil-topographique y,
• con Python ( Python: utilización de couches vectorielles et matricielles dans une perspective geéologique, sans logiciel SIG ; en francés, pero los guiones son universales).

Vea algunos ejemplos a continuación:

Con colores de elevación (GRASS GIS y R):

Con colores de un mapa geológico (Python con raster GRASS):

Y usando el script Python en la consola Python de QGIS:

Aquí capturas de pantalla para GRASS: