I use QGis, Python and OGR/Gdal libraries.
I have two WGS 84 points with elevation that define a 3D line, for example [8.11, 47.3, 800] and [8.08, 47.2, 400].
I have a tif file representing an WGS 84 Elevation Model.
A want the [x,y,h] coordinate of the intersection between Elevation Model and this line.
How do I programmatically achieve this with gdal/ogr/python/qgis libraries?