Obliczanie profilu wysokości wzdłuż linii z DEM?

Biorąc pod uwagę DEM (cyfrowy model wysokości) i linię (x1, y1) -> (x2, y2) (sparowane współrzędne (x, y)), jak można obliczyć profil wysokości lub przekrój rzutowany na tę linię?

Szukam kodu open source, którego mogę użyć w moim projekcie lub pseudokodu, którego mogę użyć jako przewodnika.

Polecenie GRASS r.profilewykonuje to ( dokumentacja , źródło ) i powinno stanowić dobrą podstawę do implementacji przekroju i jest dostępne na licencji GPL.

Może kod źródłowy wtyczki QGIS „Profile” może być początkiem.

Możesz pójść znacznie dalej, łącząc kolory DEM i rastrowe, używając GRASS z R. Zobacz:

• Profil rastrowy wzdłuż dowolnych segmentów linii lub
• http://www.portailsig.org/content/grass-gis-et-r-superposition-des-couleurs-d-un-raster-quelconque-sur-un-profil-topographique oraz,
• with Python ( Python: wykorzystanie des Couches vectorielles et matricielles dans une perspektywa géologique, sans logiciel SIG ; w języku francuskim, ale skrypty są uniwersalne).

Zobacz kilka przykładów poniżej:

Z kolorami elewacji (GRASS GIS i R):

W kolorach mapy geologicznej (Python z rastrem GRASS):

I używając skryptu Python w konsoli Python QGIS:

Oto zrzuty ekranu dla GRASS: