geogst 1.0.0__tar.gz

geogst-1.0.0/LICENSE

@@ -0,0 +1,15 @@
+MIT License - Copyright (c) 2024 Mauro Alberti (alberti.m65@gmail.com)
+
+Permission is hereby granted, free of charge, to any person obtaining a copy of this software and associated
+documentation files (the "Software"), to deal in the Software without restriction, including without limitation
+the rights to use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the Software,
+and to permit persons to whom the Software is furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all copies or substantial portions
+of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED
+TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
+THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
+DEALINGS IN THE SOFTWARE.

geogst-1.0.0/PKG-INFO

@@ -0,0 +1,163 @@
+Metadata-Version: 2.1
+Name: geogst
+Version: 1.0.0
+Summary: geogst is a structural geology module
+Home-page: https://gitlab.com/mauroalberti/geogst
+Author: Mauro Alberti
+Author-email: alberti.m65@gmail.com
+License: MIT
+Project-URL: Bug Tracker, https://gitlab.com/mauroalberti/geogst/-/issues
+Project-URL: Documentation, https://gitlab.com/mauroalberti/geogst/-/blob/main/README.md
+Project-URL: Source Code, https://gitlab.com/mauroalberti/geogst/-/tree/main
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Operating System :: OS Independent
+Classifier: Intended Audience :: Developers
+Classifier: Intended Audience :: Science/Research
+Classifier: Programming Language :: Python :: 3
+Classifier: Programming Language :: Python :: 3 :: Only
+Classifier: Topic :: Software Development
+Classifier: Topic :: Scientific/Engineering :: GIS
+Classifier: Topic :: Utilities
+Requires-Python: >=3
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: affine>=2.3.0
+Requires-Dist: bokeh>=2.3.1
+Requires-Dist: GDAL>=3.3.2
+Requires-Dist: geopandas>=0.9.0
+Requires-Dist: imageio>=2.9.0
+Requires-Dist: matplotlib>=3.4.2
+Requires-Dist: mplstereonet>=0.6.2
+Requires-Dist: numpy>=1.19.5
+Requires-Dist: numtest>=2.0.6
+Requires-Dist: ogr>=0.35.0
+Requires-Dist: pyproj>=2.5.0
+Requires-Dist: PyYAML>=5.3.1
+Requires-Dist: scipy
+Requires-Dist: setuptools
+Requires-Dist: shapely>=1.7.1
+
+# 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**
+
+![](ims/plane_dem_intersection.png)
+
+**Geological profiles**
+
+![](ims/timpa_san_lorenzo_profiles.png)
+
+Moreover, it allows to perform geometric operations on orientation data (planes, axes) and to
+plot the attitudes of geological plane and axes in stereonets.
+
+![](ims/stereonet.png)
+
+## Installation
+
+You should have *python-setuptools* already installed, otherwise you can install it via apt or pip3, e.g.:
+
+Linux:
+
+```sh
+sudo apt-get install python-setuptools
+```
+
+To install *geogst*, from the shell (you should be in the geogst directory) you can run:
+```
+python3 setup.py install
+```
+You can then check if the module was successfully installed in Python3 by importing the module via:
+```python
+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:
+
+```python
+from geogst.core.orientations.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):
+
+```python
+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.
+
+```python
+plot([or1, or2, or3]) # we provide the function arguments in a list
+```
+
+![](ims/usages_stereoplot.png)
+
+
+**P-T focal mechanisms**
+
+Creation of a P-T axes instance from a FaultSlick instance:
+
+```python
+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][tutorials] folder:**
+
+_Geometries_
+
+* [Basic spatial data][basic-spatial-data]
+
+_Intersections_
+
+* [DEM-plane intersections][dem-plane-intersections]
+
+_Orientations_
+
+* [Geological data][geological-data]
+* [Vectors and orientations][vectors-and-orientations]
+
+_Profiles_
+
+* [Geological profiles - Valnerina][geological-profiles-valnerina]
+* [Geological profiles with animation- Timpa San Lorenzo][geological-profiles-with-animation-timpa-san-lorenzo]
+* [Profile with attitudes - Monte Alpi][profile-with-attitudes-monte-alpi]
+* [Seismological profiles - Colfiorito][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](https://gitlab.com/mauroalberti/geogst/)
+
+<!-- Markdown link & img dfn's -->
+[npm-image]: https://img.shields.io/npm/v/datadog-metrics.svg?style=flat-square
+[npm-url]: https://npmjs.org/package/datadog-metrics
+[npm-downloads]: https://img.shields.io/npm/dm/datadog-metrics.svg?style=flat-square
+[travis-image]: https://img.shields.io/travis/dbader/node-datadog-metrics/master.svg?style=flat-square
+[travis-url]: https://travis-ci.org/dbader/node-datadog-metrics
+[tutorials]: https://gitlab.com/mauroalberti/geogst/-/tree/main/doc/tutorials
+[basic-spatial-data]: https://gitlab.com/mauroalberti/geogst/-/blob/main/doc/tutorials/geometries/Basic%20spatial%20data.ipynb
+[dem-plane-intersections]: https://gitlab.com/mauroalberti/geogst/-/blob/main/doc/tutorials/intersections/DEM-plane%20intersections.ipynb
+[geological-data]: https://gitlab.com/mauroalberti/geogst/-/blob/main/doc/tutorials/orientations/Geological%20data.ipynb
+[vectors-and-orientations]: https://gitlab.com/mauroalberti/geogst/-/blob/main/doc/tutorials/orientations/Vectors%20and%20orientations.ipynb
+[geological-profiles-valnerina]: https://gitlab.com/mauroalberti/geogst/-/blob/main/doc/tutorials/profiles/Geologic%20profiles%20-%20Valnerina.ipynb
+[geological-profiles-with-animation-timpa-san-lorenzo]: https://gitlab.com/mauroalberti/geogst/-/blob/main/doc/tutorials/profiles/Geologic%20profiles%20with%20animation%20-%20Timpa%20San%20Lorenzo.ipynb
+[profile-with-attitudes-monte-alpi]: https://gitlab.com/mauroalberti/geogst/-/blob/main/doc/tutorials/profiles/Profile%20with%20attitudes%20-%20Monte%20Alpi.ipynb
+[seismological-profiles-colfiorito]: https://gitlab.com/mauroalberti/geogst/-/blob/main/doc/tutorials/profiles/Seismological%20profiles%20-%20Colfiorito.ipynb

geogst-1.0.0/README.md

@@ -0,0 +1,124 @@
+# 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**
+
+![](ims/plane_dem_intersection.png)
+
+**Geological profiles**
+
+![](ims/timpa_san_lorenzo_profiles.png)
+
+Moreover, it allows to perform geometric operations on orientation data (planes, axes) and to
+plot the attitudes of geological plane and axes in stereonets.
+
+![](ims/stereonet.png)
+
+## Installation
+
+You should have *python-setuptools* already installed, otherwise you can install it via apt or pip3, e.g.:
+
+Linux:
+
+```sh
+sudo apt-get install python-setuptools
+```
+
+To install *geogst*, from the shell (you should be in the geogst directory) you can run:
+```
+python3 setup.py install
+```
+You can then check if the module was successfully installed in Python3 by importing the module via:
+```python
+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:
+
+```python
+from geogst.core.orientations.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):
+
+```python
+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.
+
+```python
+plot([or1, or2, or3]) # we provide the function arguments in a list
+```
+
+![](ims/usages_stereoplot.png)
+
+
+**P-T focal mechanisms**
+
+Creation of a P-T axes instance from a FaultSlick instance:
+
+```python
+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][tutorials] folder:**
+
+_Geometries_
+
+* [Basic spatial data][basic-spatial-data]
+
+_Intersections_
+
+* [DEM-plane intersections][dem-plane-intersections]
+
+_Orientations_
+
+* [Geological data][geological-data]
+* [Vectors and orientations][vectors-and-orientations]
+
+_Profiles_
+
+* [Geological profiles - Valnerina][geological-profiles-valnerina]
+* [Geological profiles with animation- Timpa San Lorenzo][geological-profiles-with-animation-timpa-san-lorenzo]
+* [Profile with attitudes - Monte Alpi][profile-with-attitudes-monte-alpi]
+* [Seismological profiles - Colfiorito][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](https://gitlab.com/mauroalberti/geogst/)
+
+<!-- Markdown link & img dfn's -->
+[npm-image]: https://img.shields.io/npm/v/datadog-metrics.svg?style=flat-square
+[npm-url]: https://npmjs.org/package/datadog-metrics
+[npm-downloads]: https://img.shields.io/npm/dm/datadog-metrics.svg?style=flat-square
+[travis-image]: https://img.shields.io/travis/dbader/node-datadog-metrics/master.svg?style=flat-square
+[travis-url]: https://travis-ci.org/dbader/node-datadog-metrics
+[tutorials]: https://gitlab.com/mauroalberti/geogst/-/tree/main/doc/tutorials
+[basic-spatial-data]: https://gitlab.com/mauroalberti/geogst/-/blob/main/doc/tutorials/geometries/Basic%20spatial%20data.ipynb
+[dem-plane-intersections]: https://gitlab.com/mauroalberti/geogst/-/blob/main/doc/tutorials/intersections/DEM-plane%20intersections.ipynb
+[geological-data]: https://gitlab.com/mauroalberti/geogst/-/blob/main/doc/tutorials/orientations/Geological%20data.ipynb
+[vectors-and-orientations]: https://gitlab.com/mauroalberti/geogst/-/blob/main/doc/tutorials/orientations/Vectors%20and%20orientations.ipynb
+[geological-profiles-valnerina]: https://gitlab.com/mauroalberti/geogst/-/blob/main/doc/tutorials/profiles/Geologic%20profiles%20-%20Valnerina.ipynb
+[geological-profiles-with-animation-timpa-san-lorenzo]: https://gitlab.com/mauroalberti/geogst/-/blob/main/doc/tutorials/profiles/Geologic%20profiles%20with%20animation%20-%20Timpa%20San%20Lorenzo.ipynb
+[profile-with-attitudes-monte-alpi]: https://gitlab.com/mauroalberti/geogst/-/blob/main/doc/tutorials/profiles/Profile%20with%20attitudes%20-%20Monte%20Alpi.ipynb
+[seismological-profiles-colfiorito]: https://gitlab.com/mauroalberti/geogst/-/blob/main/doc/tutorials/profiles/Seismological%20profiles%20-%20Colfiorito.ipynb

geogst-1.0.0/geogst/configs/plots.py

@@ -0,0 +1,21 @@
+
+# General
+color_def = 'grey'
+aspect_def = 1
+superposed_def = False
+
+# Attitudes
+a_color_def = "red"
+a_labels_orientions_def = True
+a_labels_ids_def = True
+
+# Line intersections
+li_color_def = "red"
+li_size_def = 3
+li_alpha_def = 0.5
+li_labels_def = True
+
+# Polygons intersections
+pi_linewidth = 2
+pi_labels_def = True
+pi_legend_def = True

geogst-1.0.0/geogst/core/deformations/space3d/deformations.py

@@ -0,0 +1,56 @@
+
+from geogst.core.deformations.space3d.shears import *
+from geogst.core.deformations.space3d.rotations import *
+
+
+def deformation_matrices(deform_params):
+
+    def_matrices = []
+
+    for deform_param in deform_params:
+        if deform_param['type'] == 'displacement':
+            displ_x = deform_param['parameters']['delta_x']
+            displ_y = deform_param['parameters']['delta_y']
+            displ_z = deform_param['parameters']['delta_z']
+            deformation = {'increment': 'additive',
+                           'matrix': np.array([displ_x, displ_y, displ_z])}
+        elif deform_param['type'] == 'rotation':
+            rot_matr = rotation_matrix_from_trend_and_plunge(deform_param['parameters']['rotation axis trend'],
+                                                             deform_param['parameters']['rotation axis plunge'],
+                                                             deform_param['parameters']['rotation angle'])
+            deformation = {'increment': 'multiplicative',
+                           'matrix': rot_matr,
+                           'shift_pt': np.array([deform_param['parameters']['center x'],
+                                                 deform_param['parameters']['center y'],
+                                                 deform_param['parameters']['center z']])}
+        elif deform_param['type'] == 'scaling':
+            scal_matr = matrScaling(deform_param['parameters']['x factor'],
+                                       deform_param['parameters']['y factor'],
+                                       deform_param['parameters']['z factor'])
+            deformation = {'increment': 'multiplicative',
+                           'matrix': scal_matr,
+                           'shift_pt': np.array([deform_param['parameters']['center x'],
+                                                 deform_param['parameters']['center y'],
+                                                 deform_param['parameters']['center z']])}
+        elif deform_param['type'] == 'simple shear - horizontal':
+            simple_shear_horiz_matr = simple_shear_horiz_matrix(deform_param['parameters']['psi angle (degr.)'],
+                                                                deform_param['parameters']['alpha angle (degr.)'])
+            deformation = {'increment': 'multiplicative',
+                           'matrix': simple_shear_horiz_matr,
+                           'shift_pt': np.array([deform_param['parameters']['center x'],
+                                                 deform_param['parameters']['center y'],
+                                                 deform_param['parameters']['center z']])}
+        elif deform_param['type'] == 'simple shear - vertical':
+            simple_shear_vert_matr = simple_shear_vert_matrix(deform_param['parameters']['psi angle (degr.)'],
+                                                              deform_param['parameters']['alpha angle (degr.)'])
+            deformation = {'increment': 'multiplicative',
+                           'matrix': simple_shear_vert_matr,
+                           'shift_pt': np.array([deform_param['parameters']['center x'],
+                                                 deform_param['parameters']['center y'],
+                                                 deform_param['parameters']['center z']])}
+        else:
+            continue
+
+        def_matrices.append(deformation)
+
+    return def_matrices