DEMからのラインに沿った標高プロファイルの計算？

DEM（デジタル標高モデル）とライン（x1、y1）->（x2、y2）（ペアの座標（x、y））が与えられた場合、そのラインに投影された標高プロファイルまたは断面をどのように計算できますか？

私は自分のプロジェクトで使用できるオープンソースコード、またはガイドとして使用できる擬似コードを探しています。

GRASSコマンドr.profileはこれを実行し（documentation、source）、断面を実装するための適切な基盤を提供する必要があり、GPLの下で利用可能です。

たぶん、QGIS「プロファイル」プラグインの ソースコードがスタートかもしれません。

RでGRASSを使用して、DEMとラスターカラーを組み合わせることで、さらに先へ進むことができます。以下を参照してください。

• 任意の線分に沿ったラスタープロファイル、または
• http://www.portailsig.org/content/grass-gis-et-r-superposition-des-couleurs-d-un-raster-quelconque-sur-un-profil-topographiqueおよび、
• Pythonを使用（Python：couches vectorielles et matricielles dans une perspectivegéologique、sans logiciel SIG ;フランス語ですが、スクリプトは普遍的です）。

以下の例を参照してください。

標高色（GRASS GISおよびR）：

地質図の色（GRASSラスターを含むPython）：

QGISのPythonコンソールでPythonスクリプトを使用します。

GRASSのスクリーンショット：