geogst

Structural geology tools with Python.

geogst is a Python 3 module library for structural geology geodata processing.

It allows to perform geological operations, such as determining the intersections between a geological plane and a topographic surface, creating multiple parallel geological profiles with attitudes, geological outcrop and line intersections.

Plane-DEM intersections

Geological profiles

Moreover, it allows to perform geometric operations on orientation data (planes, axes) and to plot the attitudes of geological plane and axes in stereonets.

Installation

You should have python-setuptools already installed, otherwise you can install it via apt or pip3, e.g.:

Linux:

sudo apt-get install python-setuptools

To install geogst, from the shell (you should be in the geogst directory) you can run:

pip install .

You can then check if the module was successfully installed in Python3 by importing the module via:

import geogst

Usage example

Orientations

Orientations can refer to linear or planar features. When referring to linear orientations, we subdivide them into directional, with a defined direction, or axial, without direction. We import all classes/methods from the orientations sub-module:

from geogst.core.geology.orientations import *

A Direct is equivalent to a unit vector in the 3D space, with orientation expressed by polar coordinates: a direction defined by a trend (from the North, 0°-360°) and a plunge (-90° to 90°, where positive values are downward-directed while negative ones are upward-directed):

or1 = Direct(130, 10)  # orientation defined by its trend and plunge. Pointing downward, since positive plunge.
or2 = Direct(312, -45)  # or2 points upwards, since it has a negative plunge
or3 = Direct(300, -20)  # as or2

We can plot using the plot function. The data to plot have to be inserted into a list.

plot([or1, or2, or3]) # we provide the function arguments in a list

P-T focal mechanisms

Creation of a P-T axes instance from a FaultSlick instance:

ptbx_from_fltsl = PTBAxes.from_fault_slickenline(Fault(90, 45, slickenlines=[Slickenline(90, 45)]))
print(ptbx_from_fltsl)

PTBAxes(P Axis(az: 180.00°, pl: 90.00°), T Axis(az: 90.00°, pl: 0.00°))

For tutorials, please refer to the Jupyter Notebooks in the doc/tutorials folder:

Geometries

• Basic spatial data

Intersections

• DEM-plane intersections

Orientations

• Geological data
• Vectors and orientations

Profiles

• Geological profiles - Valnerina
• Geological profiles with animation- Timpa San Lorenzo
• Profile with attitudes - Monte Alpi
• Seismological profiles - Colfiorito

Meta

Mauro Alberti – alberti.m65@gmail.com

Distributed under the MIT license. See LICENSE for more information.

https://gitlab.com/mauroalberti/geogst