About GeoTrace

GeoTrace is a a QGIS plugin developed by Lachlan Grose and Sam Thiele containing some helpful tools for extracting and analysing the orientations of geological structures. It can be used to rapidly digitize structural traces in raster data, estimate their 3D orientations using an associated DEM, and then visualise the results on stereonets and rose diagrams. There is a complementary plugin, Compass, for CloudCompare that provides similar functionalities for 3D point clouds

The trace extraction method (Trace tab) uses a least-cost path algorithm to "follow" linear features in the raster. This relies on a single-channel cost raster in which the structures of interest are represented by low values, and the background by high values. A variety of functions for quickly calculating such a cost function have been included in the Cost Calculator tab.

Useage Instructions

Trace tab

The trace tab is used for computer-assisted digitization. Before starting, select: - An output layer (polyline .shp file) to write digitized traces to - mA point layer to store the control points in (optional) - A cost layer. This must be a one-channel raster, in which traces will follow low values. (Though the Invert Cost check will make the trace follow high values). The 'Cost Calculator* tab can be used to assist creation of the single-channel cost raster. - A DEM layer, used to estimate 3D orientations from the traces (optional)

Once the relevent information has been set, start interpreting by clicking the Start Digitizing button. Left-click adds control points to your trace and Right-Click completes a trace. Hit Backspace to undo.

Advanced Trace

The Advanced trace tab is used to generate traces from predefined control points. This uses a Cost layer and writes to an Output layer, as above, but rather than requiring manually inserted control points it takes a point feature layer (Control Points) and ID field defining which trace each point belongs to (Unique ID Field), and will then automatically generate traces on clicking Run.

Cost Calculator

This tab wraps a variety of python functions from the scikit-image package for easy generation of cost rasters. Please refer to the scikit-image website for detailed descriptions of each of these functions.

Stereonet

The stereonet tab can be used to plot stereonets of planar orientation estimates (strike/dip or dip/dip direction) created using this plugin, or otherwise. Simply select the layer and associated fields containing the orientation estimates and then use the plotting tools to draw the stereonet.

The stereonet utility is a wrapper around the python library mplstereonet

Rose Diagram

This tab works as above, but creates a rose diagram rather than a stereonet.

Licence

GeoTrace was developed by Lachlan Grose and Sam Thiele and is free software licenced under the GNU licence v2

Further Reading and Citation

If you found this tool useful, please cite Thiele et al., 2017. The publication (currently in prep.) also contains a more detailed description of the methods employed by this plugin.

Thiele, ST., Grose, L., Samsu, A., Micklethwaite, S., Vollgger, SA. & Cruden, S., 2017, 'A computer-assisted approach to structural interpretation of point clouds and rasters', Journal TBA