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.
Você pode ir muito além combinando DEM e cores raster, usando GRASS com R. Veja:
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:
fonte
Aqui capturas de tela para GRASS:
fonte