Skip to Content

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


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.

Résultat final

Pour aller plus loin

Il est ensuite possible d'utiliser la technique exposée dans QGIS, représentation 3D des couches vectorielles (shapefiles dits 3D ou shapefiles avec attributs z) avec les modules Python Matplotlib ou Visvis à partir de la console Python pour visualiser les résultats en 3D:

Tous les traitements ont été effectués sur Mac OS X  avec QGIS, version 1.9 "master branch"  de Larry Shaffer disponible à qgis.dakotacarto.com/


Site officiel : qgSurf, a Python program for the interactive determination of intersections between DEM and planes
Site officiel : QGIS
Site officiel : plugin qgSurf


Creative Commons License
licence Creative Commons Paternité-Pas d'Utilisation Commerciale-Partage des Conditions Initiales à l'Identique Commerciale 2.0 France

Commentaires

Poster un nouveau commentaire

Le contenu de ce champ sera maintenu privé et ne sera pas affiché publiquement.