Menghitung profil ketinggian sepanjang garis dari DEM?

Diberi DEM (Digital Elevation Model) dan sebuah garis (x1, y1) -> (x2, y2) (koordinat berpasangan (x, y)), bagaimana seseorang dapat menghitung profil elevasi atau penampang yang diproyeksikan pada garis itu?

Saya mencari kode sumber terbuka yang dapat saya gunakan dalam proyek saya atau pseudo-code yang dapat saya gunakan sebagai panduan.

Perintah GRASS r.profilemelakukan ini ( dokumentasi , sumber ) dan harus memberikan dasar yang baik untuk mengimplementasikan penampang, dan tersedia di bawah GPL.

Mungkin kode sumber plugin "Profil" QGIS bisa menjadi awal.

Anda dapat melangkah lebih jauh dengan menggabungkan DEM dan warna raster, menggunakan GRASS dengan R. Lihat:

• Profil raster di sepanjang segmen garis arbitrer , atau
• http://www.portailsig.org/content/grass-gis-et-r-superposition-des-couleurs-d-un-raster-quelconque-sur-un-profil-topographique dan,
• dengan Python ( Python: utilise des couches vectorielles et matricielles dans une géologique perspektif, sans logiciel SIG ; dalam bahasa Prancis, tetapi skripnya universal).

Lihat beberapa contoh di bawah ini:

Dengan warna elevasi (GRASS GIS dan R):

Dengan warna Peta Geologis (Python dengan GRASS raster):

Dan menggunakan skrip Python di konsol Python QGIS:

Berikut screenshot untuk GRASS: