Berechnung des Höhenprofils entlang einer Linie aus einem DEM?

Wie kann man bei einem DEM (Digital Elevation Model) und einer Linie (x1, y1) -> (x2, y2) (gepaarte Koordinaten (x, y)) das Höhenprofil oder den auf diese Linie projizierten Querschnitt berechnen?

Ich suche einen Open Source Code, den ich in meinem Projekt verwenden kann, oder einen Pseudocode, den ich als Leitfaden verwenden kann.

Der Befehl GRASS r.profileführt dies aus ( Dokumentation , Quelle ) und sollte eine gute Grundlage für die Implementierung eines Querschnitts bieten. Er ist unter der GPL verfügbar.

Vielleicht kann der Quellcode des QGIS "Profile" Plugins ein Anfang sein.

Sie können noch viel weiter gehen, indem Sie DEM- und Rasterfarben mit GRASS und R kombinieren. Siehe:

• Rasterprofil entlang beliebiger Liniensegmente oder
• http://www.portailsig.org/content/grass-gis-et-r-superposition-des-couleurs-d-un-raster-quelconque-sur-un-profil-topographique und
• mit Python ( Python: Verwendung der Couches vectorielles und matricielles in einer Perspektive géologique, sans logiciel SIG ; auf Französisch, aber die Skripte sind universell).

Siehe einige Beispiele unten:

Mit Höhenfarben (GRASS GIS und R):

Mit Farben einer geologischen Karte (Python mit GRASS-Raster):

Und mit dem Python-Skript in der Python-Konsole von QGIS:

Hier Screenshots zu GRASS: