Вычисление профиля высот вдоль линии от матрицы высот?

Учитывая ЦМР (цифровую модель рельефа) и линию (x1, y1) -> (x2, y2) (парные координаты (x, y)), как можно рассчитать профиль высот или сечение, спроецированное на эту линию?

Я ищу открытый исходный код, который я могу использовать в своем проекте, или псевдокод, который я могу использовать в качестве руководства.

Команда GRASS r.profileвыполняет это ( документация , источник ) и должна обеспечить хорошую основу для реализации поперечного сечения и доступна под лицензией GPL.

Может быть, исходный код плагина QGIS «Профиль» может быть началом.

Вы можете пойти гораздо дальше, комбинируя DEM и растровые цвета, используя GRASS с R. См .:

• Растровый профиль вдоль произвольных отрезков или
• http://www.portailsig.org/content/grass-gis-et-r-superposition-des-couleurs-d-un-raster-quelconque-sur-un-profil-topographique и,
• с Python ( Python: использование заданных векторов и матриц) с точки зрения геологии, без логики SIG ; на французском языке, но сценарии универсальны).

Смотрите некоторые примеры ниже:

С цветами возвышения (GRASS GIS и R):

С цветами геологической карты (Python с растром GRASS):

И используя скрипт Python в консоли Python QGIS:

Вот скриншоты для GRASS: