Plugin QGIS - qgSurf, calcul des intersections entre des surfaces géologiques (planes) et un MNT

dim 20-01-2013 Martin Laloux

SIG OpenSource

En cartographie géologique, en analyse structurale ou en analyse géomorphologique, il est nécessaire de déterminer la trace topographique d'une surface géologique, c'est-à-dire son intersection avec la topographie,

intersection entre des couches ou un volume géologique et la topographie, réalisé avec GRASS GIS

tracé de la limite géologique sur une carte topographique (intersection)

en fonction de son orientation et de son pendage (angle d'inclinaison de la surface considérée localement comme un plan), mesurés en un ou divers points:

Si l'on dispose d'un MNT, est-il possible de calculer automatiquement cette ligne d'intersection ?

C'est le travail realisé par Mauro Alberti qui a créé le plugin qgSurf pour Quantum GIS, basé sur son application en pur Python, gSurf. Les 2 implémentations utilisent les modules Python PyQT4, osgeo/gdal-ogr et matplotlib. Les algorithmes utilisés ont été expliqués et détaillés pour la version en Python seul à Intersezioni tra DEM e superfici planari, un tema di interesse in geologia et gSurf: una applicazione Python per calcolare interattivamente l'intersezione fra piani e DEM.

application Python gSurf

Préambules

Supposons que l'on a mesuré un plan structural sur un affleurement géologique (stratification d'une couche mesurée en pendage orienté, 20° vers 164, en bleu).

fonds topographique OpenStreetMap

On dispose aussi d'un MNT (au format ESRI .asc ou GeoTIFF):

DEM SRTM

Extension qgSurf

L'extension se lance depuis le menu Extension et son interface est la suivante:

Après avoir complété tous les éléments, le résultat est obtenu par le bouton Calculate -> ligne en blanc (ici):

Les points et la ligne résultants peuvent ensuite être exportés sous forme de fichiers shapefiles et récupérés dans QGIS:

L'option Lineaments, permet d'analyser des structures linéaires dont on voudrait connaître la direction et le pendage.

Ajustement à une des fractures:

Petits problèmes et solutions

Le principal petit problème provient du fait que les lignes et les points résultants sont extrapolés sur l'emprise du MNT. Le fait de pouvoir exporter les résultats permet facilement de régler ce problème dans le cas de plusieurs mesures.