Cálculo do perfil de elevação ao longo da linha de um DEM?

Dado um DEM (Digital Elevation Model) e uma linha (x1, y1) -> (x2, y2) (coordenadas emparelhadas (x, y)), como calcular o perfil de elevação ou seção projetada nessa linha?

Estou procurando um código-fonte aberto que possa ser usado no meu projeto ou pseudo-código que possa ser usado como guia.

O comando GRASS r.profileexecuta isso ( documentação , fonte ) e deve fornecer uma boa base para a implementação de uma seção transversal, e está disponível na GPL.

Talvez o código-fonte do plugin "Profile" do QGIS possa ser um começo.

Você pode ir muito além combinando DEM e cores raster, usando GRASS com R. Veja:

• Perfil de varredura ao longo de segmentos de linha arbitrários ou
• http://www.portailsig.org/content/grass-gis-et-r-superposition-des-couleurs-d-un-raster-quelconque-sur-un-profil-topographique e,
• com Python ( Python: utilização de vetores vetoriais e matrizes em uma perspectiva geográfica, sem lógica SIG ; em francês, mas os scripts são universais).

Veja alguns exemplos abaixo:

Com cores de elevação (GRASS GIS e R):

Com cores de um mapa geológico (Python com GRASS raster):

E usando o script Python no console Python do QGIS:

Aqui capturas de tela para GRASS: