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:

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.

If a DEM layer has been included (see above), best-fit-planes for each trace will be computed using the trace eigenvectors. It is important to note that this will often produce poor results, especially in flat topography or where the traces have variable orientations. To help identify traces with poor orientation estimates, trace eigenvalues are recorded in the output, along with a "planarity" metric that approaches 1 as as plane orientation becomes well constrained (cf. Thiele et al., 2015). For convenience, orientation estimates are also classified as "Good" (planarity <0.75), "Average" (0.5<planarity<0.75), and "Poor" (planarity<0.5). Also be aware that well-constrained ("Good") orientation estimates sub-parallel to the DEM surface can be produced if the structures being digitised have variable orientations.

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 also contains a more detailed description of the methods employed by this plugin.

Thiele, S. T., Grose, L., Samsu, A., Micklethwaite, S., Vollgger, S. A., and Cruden, A. R.: Rapid, semi-automatic fracture and contact mapping for point clouds, images and geophysical data, Solid Earth, 8, 1241-1253,10.5194/se-8-1241-2017, 2017

For further information on the plane-fitting approach and planarity metric please refer to:

Thiele, S. T., Micklethwaite, S., Bourke, P., Verrall, M., and Kovesi, P.: Insights into the mechanics of en-echelon sigmoidal vein formation using ultra-high resolution photogrammetry and computed tomography, Journal of Structural Geology, 77, 27-44, 10.1016/j.jsg.2015.05.006, 2015.

Where the derived orientation estimates are of critical importance, the following will also be a useful reference:

Seers, T. D. and Hodgetts, D.: Probabilistic constraints on structural lineament best fit plane precision obtained through numerical analysis, Journal of Structural Geology, 82, 37-47, 10.1016/j.jsg.2015.11.004, 2016.