LoopStructural 1.6.2__py3-none-any.whl → 1.6.3__py3-none-any.whl

LoopStructural/utils/colours.py

@@ -0,0 +1,26 @@
+from LoopStructural.utils import rng
+
+
+def random_colour(n: int, cmap='tab20'):
+    """
+    Generate a list of random colours
+
+    Parameters
+    ----------
+    n : int
+        Number of colours to generate
+    cmap : str, optional
+        Name of the matplotlib colour map to use, by default 'tab20'
+
+    Returns
+    -------
+    list
+        List of colours in the form of (r,g,b,a) tuples
+    """
+    import matplotlib.cm as cm
+
+    colours = []
+    for _i in range(n):
+        colours.append(cm.get_cmap(cmap)(rng.random()))
+
+    return colours

LoopStructural/utils/features.py

@@ -0,0 +1,5 @@
+from ..modelling.features import LambdaGeologicalFeature
+
+X = LambdaGeologicalFeature(lambda pos: pos[:, 0], name="x")
+Y = LambdaGeologicalFeature(lambda pos: pos[:, 1], name="y")
+Z = LambdaGeologicalFeature(lambda pos: pos[:, 2], name="z")

LoopStructural/utils/maths.py

@@ -244,3 +244,54 @@ def get_dip_vector(strike, dip):
         ]
     )
     return v
+
+
+def regular_tetraherdron_for_points(xyz, scale_parameter):
+    regular_tetrahedron = np.array(
+        [
+            [np.sqrt(8 / 9), 0, -1 / 3],
+            [-np.sqrt(2 / 9), np.sqrt(2 / 3), -1 / 3],
+            [-np.sqrt(2 / 9), -np.sqrt(2 / 3), -1 / 3],
+            [0, 0, 1],
+        ]
+    )
+    regular_tetrahedron *= scale_parameter
+    tetrahedron = np.zeros((xyz.shape[0], 4, 3))
+    tetrahedron[:] = xyz[:, None, :]
+    tetrahedron[:, :, :] += regular_tetrahedron[None, :, :]
+
+    return tetrahedron
+
+
+def gradient_from_tetrahedron(tetrahedron, value):
+    """
+    Calculate the gradient from a tetrahedron
+    """
+    tetrahedron = tetrahedron.reshape(-1, 4, 3)
+    m = np.array(
+        [
+            [
+                (tetrahedron[:, 1, 0] - tetrahedron[:, 0, 0]),
+                (tetrahedron[:, 1, 1] - tetrahedron[:, 0, 1]),
+                (tetrahedron[:, 1, 2] - tetrahedron[:, 0, 2]),
+            ],
+            [
+                (tetrahedron[:, 2, 0] - tetrahedron[:, 0, 0]),
+                (tetrahedron[:, 2, 1] - tetrahedron[:, 0, 1]),
+                (tetrahedron[:, 2, 2] - tetrahedron[:, 0, 2]),
+            ],
+            [
+                (tetrahedron[:, 3, 0] - tetrahedron[:, 0, 0]),
+                (tetrahedron[:, 3, 1] - tetrahedron[:, 0, 1]),
+                (tetrahedron[:, 3, 2] - tetrahedron[:, 0, 2]),
+            ],
+        ]
+    )
+    I = np.array([[-1.0, 1.0, 0.0, 0.0], [-1.0, 0.0, 1.0, 0.0], [-1.0, 0.0, 0.0, 1.0]])
+    m = np.swapaxes(m, 0, 2)
+    element_gradients = np.linalg.inv(m)
+
+    element_gradients = element_gradients.swapaxes(1, 2)
+    element_gradients = element_gradients @ I
+    v = np.sum(element_gradients * value[:, None, :], axis=2)
+    return v

LoopStructural/version.py

@@ -1 +1 @@
-__version__ = "1.6.2"
+__version__ = "1.6.3"

LoopStructural-1.6.3.dist-info/METADATA

@@ -0,0 +1,146 @@
+Metadata-Version: 2.1
+Name: LoopStructural
+Version: 1.6.3
+Summary: 3D geological modelling
+Author-email: Lachlan Grose <lachlan.grose@monash.edu>
+License: MIT
+Project-URL: Documentation, https://Loop3d.org/LoopStructural/
+Project-URL: Bug Tracker, https://github.com/loop3d/loopstructural/issues
+Project-URL: Source Code, https://github.com/loop3d/loopstructural
+Keywords: earth sciences,geology,3-D modelling,structural geology,uncertainty
+Classifier: Development Status :: 5 - Production/Stable
+Classifier: Intended Audience :: Science/Research
+Classifier: Topic :: Scientific/Engineering :: Information Analysis
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Operating System :: Microsoft :: Windows
+Classifier: Operating System :: POSIX
+Classifier: Operating System :: MacOS
+Classifier: Programming Language :: Python :: 3.8
+Classifier: Programming Language :: Python :: 3.9
+Classifier: Programming Language :: Python :: 3.10
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
+Requires-Python: >=3.8
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: numpy >=1.18
+Requires-Dist: pandas
+Requires-Dist: scipy
+Requires-Dist: scikit-image
+Requires-Dist: scikit-learn
+Requires-Dist: tqdm
+Provides-Extra: all
+Requires-Dist: loopstructural[export,inequalities,visualisation] ; extra == 'all'
+Provides-Extra: export
+Requires-Dist: geoh5py ; extra == 'export'
+Requires-Dist: pyevtk ; extra == 'export'
+Requires-Dist: dill ; extra == 'export'
+Provides-Extra: inequalities
+Requires-Dist: loopsolver ; extra == 'inequalities'
+Provides-Extra: jupyter
+Requires-Dist: pyvista[all] ; extra == 'jupyter'
+Requires-Dist: tqdm ; extra == 'jupyter'
+Provides-Extra: visualisation
+Requires-Dist: matplotlib ; extra == 'visualisation'
+Requires-Dist: pyvista ; extra == 'visualisation'
+Requires-Dist: loopstructuralviusualisation >0.1.4 ; extra == 'visualisation'
+
+# LoopStructural: Loop3D Implicit Geological Modelling
+
+![3D model of Hamersley created using loopstructural](docs/source/images/image823.png)
+[![PyPI version](https://badge.fury.io/py/LoopStructural.svg)](https://badge.fury.io/py/LoopStructural)
+[![GitHub license](https://img.shields.io/github/license/Naereen/StrapDown.js.svg)](https://github.com/Loop3D/LoopStructural/blob/master/LICENSE)
+[![Documentation loop3d.github.io/LoopStructural/](https://img.shields.io/badge/docs-githubio-brightgreen)](https://loop3d.github.io/LoopStructural)
+
+## Overview
+
+LoopStructural is an opensource Python library for 3D geological modelling. The library has been built in the scope of the Loop project ([Loop3d.org](https://loop3d.org/)). LoopStructural can:
+
+- Model fault networks integrating fault kinematics and overprinting relationships
+- Model folds using structural elements (fold axis, fold axial surface) for multiple fold events
+- Use different implicit interpolation algorithms (Finite Difference, Piecewiese Linear, RBF)
+- Export models to vtk, geoh5, omf, gocad, csv, obj formats
+- Visualise models in an interactive python environment.
+
+## Installation
+
+### Google colab
+
+LoopStructural can be used inside a google colab notebook, however interactive visualisation using kitware trame does not work. Only static plots can be displayed.
+
+### Pip
+
+`pip install loopstructural`
+
+Optional dependencies are:
+
+- loopstructuralvisualisation (wrapper for pyvista to interface with loopstructural objects)
+- pyvista (3D visualisation)
+- geoh5py (export to .geoh5 format)
+- mira-omf (export to .omf format to load in Leapfrog)
+- pyevtk
+- dill (serialize loop models to a pickle file)
+- tqdm (progress bars)
+- loopsolver (experimental solver for admm inequality constraints)
+
+You can install all of the optional dependencies using:
+`pip install loopstructural[all]`
+
+### conda
+
+`conda install -c conda-forge -c loop3d loopstructural`
+
+to install the working 3D visualisation environment
+`conda install -c conda-forge -c loop3d loopstructural loopstructuralvisualisation pyvista trame trame-vtk trame-vuetify`
+
+## Quickstart
+
+- A basic geological model
+
+```Python
+from LoopStructural import GeologicalModel
+from LoopStructural.datatypes import BoundingBox
+from LoopStructural.visualisation import Loop3DView
+from LoopStructural.datasets import load_claudius
+
+import numpy as np
+data, bb = load_claudius()
+
+#bb constaints origin and maximum of axis aligned bounding box
+#data is a pandas dataframe with X,Y,Z,val,nx,ny,nz, feature_name
+
+model = GeologicalModel(bb[0,:],bb[1,:])
+model.data = data
+# nelements specifies the number of discrete interpolation elements
+# 'stratí' is the feature name in the data dataframe
+model.create_and_add_foliation('strati',nelements=1e5)
+model.update()
+# get the value of the interpolator at some random locations
+locations = np.array(
+    [
+        np.random.uniform(bb[0, 0], bb[1, 0],5),
+        np.random.uniform(bb[0, 1], bb[1, 1],5),
+        np.random.uniform(bb[0, 2], bb[1, 2],5),
+    ]
+).T
+val = model.evaluate_feature_value('strati', locations)
+# get the gradient of the interpolator
+gradient = model.evaluate_feature_gradient('strati',locations)
+
+#Plot the scalar field of the model
+model['strati'].scalar_field().plot()
+
+
+```
+
+## Documentation
+
+The LoopStructural documentation can be found [here](https://loop3d.github.io/LoopStructural)
+
+## Problems
+
+Any bugs/feature requests/comments please create a new [issue](https://github.com/Loop3D/LoopStructural/issues).
+
+## Acknowledgements
+
+_The Loop platform is an open source 3D probabilistic geological and geophysical modelling platform, initiated by Geoscience Australia and the OneGeology consortium. The project is funded by Australian territory, State and Federal Geological Surveys, the Australian Research Council and the MinEx Collaborative Research Centre._

{LoopStructural-1.6.2.dist-info → LoopStructural-1.6.3.dist-info}/RECORD

@@ -1,5 +1,5 @@
 LoopStructural/__init__.py,sha256=nlJ0csTZ_xlTZrWvTKiTp3ER_XRQN6iqQ5sR94RD5QA,1331
-LoopStructural/version.py,sha256=iaZ7XZ-SdVNT4ZACUYJL_5bHKP1lOYNnZyPeapZHkxk,22
+LoopStructural/version.py,sha256=X-erfDySXG8dR_oX8H_L2S8f9Y_VDgdltPb-jH2ruKg,22
 LoopStructural/datasets/__init__.py,sha256=ylb7fzJU_DyQ73LlwQos7VamqkDSGITbbnoKg7KAOmE,677
 LoopStructural/datasets/_base.py,sha256=FB_D5ybBYHoaNbycdkpZcRffzjrrL1xp9X0k-pyob9Y,7618
 LoopStructural/datasets/_example_models.py,sha256=Zg33IeUyh4C-lC0DRMLqCDP2IrX8L-gNV1WxJwBGjzM,113
@@ -29,77 +29,80 @@ LoopStructural/datasets/data/geological_map_data/stratigraphic_order.csv,sha256=
 LoopStructural/datasets/data/geological_map_data/stratigraphic_orientations.csv,sha256=RysyqUAIjY6iIDUfTh11n9QUQWXB_qxKnZeN_DqNzlY,26745
 LoopStructural/datasets/data/geological_map_data/stratigraphic_thickness.csv,sha256=pnSmG-wL8-kxuoHo_pgpJrfTmsZOzc8L0vxpBRh3r8A,355
 LoopStructural/datatypes/__init__.py,sha256=lVg64DnynMm58qvYTjLrcyWH7vk2ngr9JGMo5FaiALI,160
-LoopStructural/datatypes/_bounding_box.py,sha256=kK2hbI9NCBcR9K6S7FNZs9cbaUPVcT41Okg_MZqObcA,14101
-LoopStructural/datatypes/_point.py,sha256=DH-s8GJeomOwRmalqooNqeV5AOaEHOxVeW7a8yrb14A,5251
-LoopStructural/datatypes/_structured_grid.py,sha256=OFE7dog3TSBSUDGbSEV7ZPUekG0ToNAkFROw8pcm014,3184
-LoopStructural/datatypes/_surface.py,sha256=J40crDWkbGAbaDNYP0J8koO_eWdZucmbXRsGxLPVHDU,5964
+LoopStructural/datatypes/_bounding_box.py,sha256=nTtvNj0UiQWX-61tBH6iWS8YDzHBIqTKWHUxStkKBfU,14732
+LoopStructural/datatypes/_point.py,sha256=rYolPcd134DCj5kvGkTG9_T2S68XOPaBuLfq6-Ud7MU,6140
+LoopStructural/datatypes/_structured_grid.py,sha256=AwSzkoAySU8ZHGbbGTyUoU3fCXuBVQAOD_6th6v0J1o,3665
+LoopStructural/datatypes/_surface.py,sha256=3RLA6-zkWE1L3IvXdrq92Dd5asnI9Lf6ebcNmrMnn98,6445
 LoopStructural/export/exporters.py,sha256=BniZu-PqQvHqCU6GIuJQ5FPzI9Dx_T6rI8EW1pykois,17209
 LoopStructural/export/file_formats.py,sha256=0xKyYSW4Jv_4jsXwusg-WO6PNUhZKd6HdWSqGSaPve8,232
 LoopStructural/export/geoh5.py,sha256=TByfnHul1Rg4oB5KPnD-yCcZ-uNr4cxuYnRB_6mclnA,4253
 LoopStructural/export/gocad.py,sha256=cQ6v7ZD0CVubt3c2f9EwAYrziu5bEFSWBtx0uade5mg,3370
-LoopStructural/export/omf_wrapper.py,sha256=faswHT4k1ZMaJLdZIyxGJp-IlXBdwGREnlcOSooJh3o,2375
-LoopStructural/interpolators/__init__.py,sha256=pBxqCCIJ3WS8W179EnFK9QVmIzgzqqt0T3ehu5absWA,2944
-LoopStructural/interpolators/_api.py,sha256=m1o5oza7kOXDPgasDqKY1AQE-Tz0pCzuTJXtqgA5S88,4858
+LoopStructural/export/omf_wrapper.py,sha256=4vcF4WOQIVEYsWrfATgKCDh8nUybLTbrlXnCxZ_3fkU,3392
+LoopStructural/interpolators/__init__.py,sha256=HcQ-wOnTKFtpS02Pxd_aP9s6bgjCb7iECYYFzX1nd44,3396
+LoopStructural/interpolators/_api.py,sha256=EC4ogG2uPq-z_pgNGd_eTieTl92eaZ-rjyoFwXiHL_s,7839
 LoopStructural/interpolators/_builders.py,sha256=6M8__gGDO-dyu3LPIwcPia1knLR5-EH0GBMLSF0Kuwc,6476
-LoopStructural/interpolators/_discrete_fold_interpolator.py,sha256=4yfXIKiFTllomWF0nSLDtOk0S6MflkG7uZbvVcGW6KE,6247
-LoopStructural/interpolators/_discrete_interpolator.py,sha256=NOIlIosMPozSS63coGM1zwJTioezLuDzbZZ5sbX8Dzk,23671
-LoopStructural/interpolators/_finite_difference_interpolator.py,sha256=kwu5L9kTEelB0ibiKTZBtWSZxClJTP9PRMmmlfODbO8,16505
-LoopStructural/interpolators/_geological_interpolator.py,sha256=FGyc8AcuTh70kedbMK4QO9KMTdaA81bELh587n8bJHQ,10408
+LoopStructural/interpolators/_discrete_fold_interpolator.py,sha256=eDe0R1lcQ0AuMcv7zlpu5c-soCv7AybIqQAuN2vFE3M,6542
+LoopStructural/interpolators/_discrete_interpolator.py,sha256=EPiEWtx8RuE75xp7umy7Py4vTiFB0xgN3EFnl44ygKM,25577
+LoopStructural/interpolators/_finite_difference_interpolator.py,sha256=oFchMqCSWjRquJmI0rNrtjFhre4O7PLyOLxZFE64Se4,15819
+LoopStructural/interpolators/_geological_interpolator.py,sha256=4qPkJKWW1Kqe-U78CLR8ayLdFXq9aJZ2FB5_ntITxHk,11028
 LoopStructural/interpolators/_interpolator_factory.py,sha256=LSxrf7tEWek_0tNxioGHeWWwu3Ywf8Ar79M6gC75tBY,3468
-LoopStructural/interpolators/_non_linear_discrete_interpolator.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 LoopStructural/interpolators/_operator.py,sha256=PZOUzq9OMaJdG151dSLIo7AxRuhTj6-zEAzFZo-EOJU,1114
-LoopStructural/interpolators/_p1interpolator.py,sha256=et-8sRX442XVNZK14XHx7GYtkCaW7OSkawn9zzftg3g,8753
+LoopStructural/interpolators/_p1interpolator.py,sha256=4rjj4iaw2c8hOfBS9u8ycxzijYdmvpeijvhYRwUwZg0,8736
 LoopStructural/interpolators/_p2interpolator.py,sha256=UT-As5RNsmOwHOzO_6FiRcAwlNHfi4ILbJw2LGpwKAw,10274
-LoopStructural/interpolators/_surfe_wrapper.py,sha256=OOf8KevdBGRnrYtkw3EPTWR9Iy43gpH4Vj4UrbEiJmA,5909
+LoopStructural/interpolators/_surfe_wrapper.py,sha256=uwqABnixan9bcyU08hbxO9ATO4DawItuGSGj7iAKa9U,6772
 LoopStructural/interpolators/_cython/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-LoopStructural/interpolators/supports/_2d_base_unstructured.py,sha256=IPDdXGn8QVPugqVGbepj2rKQjbQinjZtduus240wvb0,11197
+LoopStructural/interpolators/supports/_2d_base_unstructured.py,sha256=maHzpHnPRoo5IOLcp4Si7kj0oaxFox8XYrknwErbyoo,11796
 LoopStructural/interpolators/supports/_2d_p1_unstructured.py,sha256=okcy56nyjuedmknQn_95V2tm0kdMA-oJcD3U2jU8u0w,2637
 LoopStructural/interpolators/supports/_2d_p2_unstructured.py,sha256=TeBVtT1PMV7CKzmnFZke37acMoFxouer20cskS7pVoE,10422
-LoopStructural/interpolators/supports/_2d_structured_grid.py,sha256=8jG0ASTHtCiEohvsc8TzwBVwJ0iyJJID4SVJQwgk2nk,14978
+LoopStructural/interpolators/supports/_2d_structured_grid.py,sha256=TbAV1iw2gmZ4ovSMj04GLKaDZjbxK7YAceJ4wBYn8m8,16209
 LoopStructural/interpolators/supports/_2d_structured_tetra.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-LoopStructural/interpolators/supports/_3d_base_structured.py,sha256=IDORlhwyCYP3paBEd1Dtwc7u6Oszs6CxcaNgHRH6wCI,14539
+LoopStructural/interpolators/supports/_3d_base_structured.py,sha256=keZnVSK0_wmXCLkWicyIlIV7Wsy2p67OC8B2AOUIMDM,15531
 LoopStructural/interpolators/supports/_3d_p2_tetra.py,sha256=CqGVJRUMxbPQZDhhopNt_s9gVhMqh4YbjQyDZonoyxc,11574
-LoopStructural/interpolators/supports/_3d_structured_grid.py,sha256=S6w0z7GpxS6EaJ-PW18_onnrTUf3U25QxkMw56MI1aU,16361
-LoopStructural/interpolators/supports/_3d_structured_tetra.py,sha256=x032Bipt5r3st5JWska8dpq-FU32rnaDQhKXJQqyYCg,26279
-LoopStructural/interpolators/supports/_3d_unstructured_tetra.py,sha256=rtMSdy4cg0bSmNtvYa8-Sue8WVhsK57F5YKBiBJr58M,23265
-LoopStructural/interpolators/supports/__init__.py,sha256=PU68529KecFJS0QiH6GzWj2q4SaCzkBEZhryO_RNZ3Q,1485
+LoopStructural/interpolators/supports/_3d_structured_grid.py,sha256=iHrZIXQk80vAxkkTzmBRLz3yA6MJEZFKm1RaRrkNlog,17446
+LoopStructural/interpolators/supports/_3d_structured_tetra.py,sha256=5zUNtvEXDvbCHZCu6Fz9WjGbnrMaq-sYJqNUufyLcq8,26505
+LoopStructural/interpolators/supports/_3d_unstructured_tetra.py,sha256=vc9ZffCYk7m1Ae5-H_m1P5OCHgBa9jztrklugMPPiHc,23500
+LoopStructural/interpolators/supports/__init__.py,sha256=V0JjixoBIUZVAo5MmqARR67xDOoQwnb4G3SXeOMRSyQ,1603
 LoopStructural/interpolators/supports/_aabb.py,sha256=Z-kH_u6c6izak0aHG3Uo14PEKQeZmYlevLDC32Q06xk,3208
-LoopStructural/interpolators/supports/_base_support.py,sha256=KgSBW3YhRSX1Tm8h6a-2dQiuJxeKIhqVCmSYgASgMX8,2275
+LoopStructural/interpolators/supports/_base_support.py,sha256=5Cwq16EhA2rx4-DFUmoM5z6DRnDfRLkachBpGyzShdM,2420
 LoopStructural/interpolators/supports/_face_table.py,sha256=Hyj4Io63NkPRN8ab9uDHyec-2Kb8BLY_xBF6STNlvBw,3095
 LoopStructural/interpolators/supports/_support_factory.py,sha256=c86NnM4azWhS2ajPApePap0sFI82mZC8siVAU1fCOn4,1175
-LoopStructural/modelling/__init__.py,sha256=pPJSZMlyE_xliv0c1C4M5QlqYuBOUTRLQQtmXKQ1xpU,754
+LoopStructural/modelling/__init__.py,sha256=IaIRrSvTezC1_RHZpGiVuxZ9tv8XZ4yRLHC2AgXi0kE,787
 LoopStructural/modelling/core/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-LoopStructural/modelling/core/geological_model.py,sha256=izhq1plPWLjJb8qwJbpqf6PPN4ndC1_2-5ujv6x4dQg,66173
+LoopStructural/modelling/core/geological_model.py,sha256=8TFQho3QolAMwb21uWsr_Ys5rSWHOyrXSVBgQPOnoS4,66092
 LoopStructural/modelling/features/__init__.py,sha256=qTl-jTau-nnSZ4G4Np6JkBcXB01ViXhouUo_-ImoDvk,877
-LoopStructural/modelling/features/_analytical_feature.py,sha256=L0GQvOPAk27KNejaqvbLXOJo_jbJJS_V3m4c9Dymheg,2701
+LoopStructural/modelling/features/_analytical_feature.py,sha256=7-wFD_h6XDE89jn_trMkLeL98nHAj_5mj802mWTIp7E,2825
 LoopStructural/modelling/features/_base_geological_feature.py,sha256=SP-lWDJDdXZY0UAht33sJkOpW-lGN9NZqP2RuZDgyxg,11102
 LoopStructural/modelling/features/_cross_product_geological_feature.py,sha256=IOcJBgQeGBOgwymcyyEdsQHHbVs23sYPpl8xavBaPZE,2781
-LoopStructural/modelling/features/_geological_feature.py,sha256=pdLum5tTn4aImUhe-47F6jyqbSb73-AIWhegLNG7FP4,9699
+LoopStructural/modelling/features/_geological_feature.py,sha256=Migi29Qr_hjIhNYuhclSK4I_MyETwa27lHeClhDKyu0,11221
 LoopStructural/modelling/features/_lambda_geological_feature.py,sha256=pEfLN9HcYGluxRB7LJpcKdYqXP_Mkyo1pLZ3kZZ4mvA,2599
 LoopStructural/modelling/features/_region.py,sha256=TB4qnoTDQM2VgRjgyODN839fKe3kuRYLllJj0xnDKXo,478
 LoopStructural/modelling/features/_structural_frame.py,sha256=I6W-YkQiRADX7xcKMFzDefYGDlXWVAp5nAVRFVcqNow,5062
 LoopStructural/modelling/features/_unconformity_feature.py,sha256=-0shSjYdT8IuNDNDYnqJ2niFjFFiBt-Izc-oMBYa0MQ,2405
 LoopStructural/modelling/features/builders/__init__.py,sha256=Gqld1C-PcaXfJ8vpkWMDCmehmd3hZNYQk1knPtl59Bk,266
-LoopStructural/modelling/features/builders/_base_builder.py,sha256=cafVl-KImwSZmHxTuzlpnkh1lW9ARJXJLE0NQ4T_ook,3221
+LoopStructural/modelling/features/builders/_base_builder.py,sha256=ImlKZs4sFSrtAAp1NEpQD35AKF_mLEHr7jj2DzSyAd4,3561
 LoopStructural/modelling/features/builders/_fault_builder.py,sha256=EyDLyovrnxqSk11zWFnViok7lXMt_E-pG2QqgRgvOLk,25653
-LoopStructural/modelling/features/builders/_folded_feature_builder.py,sha256=oeJjCnY_L4FXRvH-bPBySWlyups68H0cM1eXjXef3aY,5242
-LoopStructural/modelling/features/builders/_geological_feature_builder.py,sha256=uBPcwW87SQVH7nYQIZCXFqCpWrKhI1QJ70WC0rTTdlk,20960
-LoopStructural/modelling/features/builders/_structural_frame_builder.py,sha256=jVIV1wUo1MxcpFx6RxBtl8IgXJt2USq7XX1m-NAcMNU,7635
-LoopStructural/modelling/features/fault/__init__.py,sha256=tkKhGJDyuzzV4WoXsoknwUBaLY_WmpZYGMCY9NK7dfM,170
-LoopStructural/modelling/features/fault/_fault_function.py,sha256=eIZCsGp2-zNgFWjtlmgq3JGRvQ8EeUyU6iVOuVAeLVU,12083
+LoopStructural/modelling/features/builders/_folded_feature_builder.py,sha256=ZIhzQV0ZRFfZF_3zJGMyMPk-xVmCro5k2HJCzB790gA,6535
+LoopStructural/modelling/features/builders/_geological_feature_builder.py,sha256=Js3NgfbFipZrTmKrIMKsNpRE5EeBxy0_TmGz3ytO5YA,21628
+LoopStructural/modelling/features/builders/_structural_frame_builder.py,sha256=1SS_wrnkmaP7NTvTyHGTkuYaaU2Vc-_XmpmxS-CQrwg,7762
+LoopStructural/modelling/features/fault/__init__.py,sha256=4u0KfYzmoO-ddFGo9qd9ov0gBoLqBiPAUsaw5zhEOAQ,189
+LoopStructural/modelling/features/fault/_fault_function.py,sha256=5IzVyvv1tS5Z6geMB8FdTZpGrBYbyu38cOH03NACm0o,12797
 LoopStructural/modelling/features/fault/_fault_function_feature.py,sha256=1F_xGGRtLZ7mBnXutEGkwWoJImSEPUuxstf0Nu-7k5w,2452
-LoopStructural/modelling/features/fault/_fault_segment.py,sha256=3JXRGq2cz3rqfl3OhDE3dWzZE2-g_nihHtmS0muxRgk,18221
-LoopStructural/modelling/features/fold/__init__.py,sha256=taSRFnKibB0KdvFlvpJKrL_ZtmFaiQzC_U4qL7n4K9o,244
+LoopStructural/modelling/features/fault/_fault_segment.py,sha256=cAnf8cKYJd-3Mb8Q6F1ZVAuez2jkHZhRN4yAz_aMNfI,16947
+LoopStructural/modelling/features/fold/__init__.py,sha256=mOq_K5IHjE3FoK7dhYpljUMGj0gvEX2Fyqtzjd4EYfQ,176
 LoopStructural/modelling/features/fold/_fold.py,sha256=bPnnLUSiF4uoMRg8aHoOSTPRgaM0JyLoRQPu5_A-J3w,5448
-LoopStructural/modelling/features/fold/_fold_rotation_angle.py,sha256=9kKMQsOmJJGAs5vGDp7FSYgKL71yH04FslILRG3Beb8,4456
-LoopStructural/modelling/features/fold/_fold_rotation_angle_feature.py,sha256=SQ7HxlHf_QrnlpK6EJ6lpn4yLcw5t820hjAxRs6EQpk,1222
-LoopStructural/modelling/features/fold/_foldframe.py,sha256=sX_Akae4JLlB4rC3oMQ6KEtEX57euqInNLF8hc8-AoE,7721
-LoopStructural/modelling/features/fold/_svariogram.py,sha256=uatBG12QUxenNUEiayUlRQMoLREFAc0mxyaZVPIVz5k,6326
+LoopStructural/modelling/features/fold/_fold_rotation_angle_feature.py,sha256=CXLbFRQ3CrTMAcHmfdbKcmSvvLs9_6TLe0Wqi1pK2tg,892
+LoopStructural/modelling/features/fold/_foldframe.py,sha256=Rgf5aofN0OVDTZ2pzqLzAGlJUO2rnNm3aFvLSnH77yo,7669
+LoopStructural/modelling/features/fold/_svariogram.py,sha256=uzGaKZ5HGh8xZcsGGg68GUKVjkd5udLy7-4lh0NQc2Y,7765
+LoopStructural/modelling/features/fold/fold_function/__init__.py,sha256=VqMjabsBd5GnPnDMXeKwXqtd0te2iXnvHxpf6jCC9YU,830
+LoopStructural/modelling/features/fold/fold_function/_base_fold_rotation_angle.py,sha256=EC9q-zKX86AoDruXg_6MOD7lo1kVshfFcLlZ25EfFp0,8048
+LoopStructural/modelling/features/fold/fold_function/_fourier_series_fold_rotation_angle.py,sha256=Cjb6Pt6cdRoH3WGqFyJ2rHxnMe6SKvzRayA6hTuwZA8,4069
+LoopStructural/modelling/features/fold/fold_function/_lambda_fold_rotation_angle.py,sha256=U-IL7DETQhKh7a-cHlEKtsIg2kd8JZqO2gWyb_PKB8k,1446
+LoopStructural/modelling/features/fold/fold_function/_trigo_fold_rotation_angle.py,sha256=v-3YkBbsqPdbPZD2ykPB3vLKz3W_BDSF2m0uxtpP2vo,4904
 LoopStructural/modelling/input/__init__.py,sha256=HhJM3V5b-8_64LiRbF3Bd1pjWhJlcknxMSMPRrqZ0-I,153
 LoopStructural/modelling/input/fault_network.py,sha256=0uxl7lOySdhMhNXoiOkuiHIXqAz1Ls0j-W65cmdQoP8,2348
 LoopStructural/modelling/input/map2loop_processor.py,sha256=T7Fgqd7FNJWylLKvfIniRZBMRMeAoP8iU330-WYU8Fg,7031
-LoopStructural/modelling/input/process_data.py,sha256=M02AzLbHd7X6SIvIQG-zXBf71FIzMIjI1ITqwIW9AEk,25160
-LoopStructural/modelling/input/project_file.py,sha256=mnK6Ts2z_XquYO_JviM-cXTH8T3YA_ZpqrC2qyLaU8M,3305
+LoopStructural/modelling/input/process_data.py,sha256=Y3wwUigqc9tqaX4EOzAOAigepzbtO6VWFhl_YmF00nA,25441
+LoopStructural/modelling/input/project_file.py,sha256=BYV9KBkAFHubP3QOWbnPjOPpaoyhh4BuZU1I194wJUs,4221
 LoopStructural/modelling/intrusions/__init__.py,sha256=EpZK3cHJwGQhPUYIwKCKu8vkNdt_nOgWF0zfhiqDYDA,712
 LoopStructural/modelling/intrusions/geom_conceptual_models.py,sha256=jwTlhYySUj7z4DEnJoi4AINZB_N3-SW6ONRFL66OsW0,3665
 LoopStructural/modelling/intrusions/geometric_scaling_functions.py,sha256=PK3qf0TiK-WYIBGG7fYhTD1hwlUN0s75BK8d53SLYuQ,3209
@@ -108,22 +111,24 @@ LoopStructural/modelling/intrusions/intrusion_feature.py,sha256=5DFQ1uccX7Rp3PMc
 LoopStructural/modelling/intrusions/intrusion_frame_builder.py,sha256=Q1TPHxREcrO7Rw71nUfACZHfYnISLjqlgkUNTPT324k,40143
 LoopStructural/modelling/intrusions/intrusion_support_functions.py,sha256=wodakheMD62WJyoKnyX8UO-C1pje0I-5kHQEoDqShzo,13951
 LoopStructural/utils/__init__.py,sha256=K6OwdnC8RE-sDhG4GIbVyoJNUD0WnM1VGG6os1l_Woc,870
-LoopStructural/utils/_surface.py,sha256=U2cM9nPqAjdgRlFyeejmekJmRyG0mjm7HcGntk0Uudc,5364
+LoopStructural/utils/_surface.py,sha256=8pUn11Pe3ZlYoHkjy06v3QWJ7UJORQ7yWsgwFQptkxs,5530
 LoopStructural/utils/_transformation.py,sha256=7iDPMIBoZ73ZGa_slbC9XzdMQ-ONOxvF7UkY4dJaw54,2443
+LoopStructural/utils/colours.py,sha256=NOYTN-fh-XTfQdyqRVzO5bdt0kJQrnxH25ZxEnqp8ZQ,548
 LoopStructural/utils/config.py,sha256=ITGOtZTo2_QBwXkG_0AFANfE90J9siCXLzxypVmg9QA,414
 LoopStructural/utils/dtm_creator.py,sha256=-yqGG0wyEJfTCCDghz058wull1q3zGFASjeu8oDgYnk,535
 LoopStructural/utils/exceptions.py,sha256=SJboJ7ncMqVX-ib7MMizClwMrFZRHQhjZr2eCnVwnQE,500
+LoopStructural/utils/features.py,sha256=WCatS4lYBrURNvWvWwhOsDVUod9KIPNq3x0OHPbWctU,241
 LoopStructural/utils/helper.py,sha256=An9NuRH16cASUWq2ZakHc1tZt_AvUpgx8tv4cyWZEQk,6581
 LoopStructural/utils/json_encoder.py,sha256=5YNouf1TlhjEqOYgthd07MRXc0JLgxern-nyKSZ__ws,403
 LoopStructural/utils/linalg.py,sha256=tBXyu6NXcG2AcPuzUMnkVI4ncZWtE_MPHGj2PLXRwfY,123
 LoopStructural/utils/logging.py,sha256=dIUWEsS2lT4G1dsf4ZYXknTR7eQkrgvGA4b_E0vMIRU,2402
-LoopStructural/utils/maths.py,sha256=q8KF2lUoyXIavmlH7s0jeRwkoGG99X40A6co5-tgqfg,6693
+LoopStructural/utils/maths.py,sha256=8iqdQdB2-bf14SzIzfFxvjWbzmPknqK9DI7CWEcW6XU,8402
 LoopStructural/utils/regions.py,sha256=LvcOCPudF4u95-GKBOZqXVxOEcR3cOFgFpcs5x43sMk,3914
 LoopStructural/utils/typing.py,sha256=29uVSTZdzXXH-jdlaYyBWZ1gQ2-nlZ2-XoVgG_PXNFY,157
 LoopStructural/utils/utils.py,sha256=2Z4zVE6G752-SPmM29zebk82bROJxEwi_YiiJjcVED4,2438
 LoopStructural/visualisation/__init__.py,sha256=5BDgKor8-ae6DrS7IZybJ3Wq_pTnCchxuY4EgzA7v1M,318
-LoopStructural-1.6.2.dist-info/LICENSE,sha256=ZqGeNFOgmYevj7Ld7Q-kR4lAxWXuBRUdUmPC6XM_py8,1071
-LoopStructural-1.6.2.dist-info/METADATA,sha256=CxK-vbkaxg1TIa1Uwim0KInk_9vy-AHw_Ufv6n5W9ow,4958
-LoopStructural-1.6.2.dist-info/WHEEL,sha256=R0nc6qTxuoLk7ShA2_Y-UWkN8ZdfDBG2B6Eqpz2WXbs,91
-LoopStructural-1.6.2.dist-info/top_level.txt,sha256=QtQErKzYHfg6ddxTQ1NyaTxXBVM6qAqrM_vxEPyXZLg,15
-LoopStructural-1.6.2.dist-info/RECORD,,
+LoopStructural-1.6.3.dist-info/LICENSE,sha256=ZqGeNFOgmYevj7Ld7Q-kR4lAxWXuBRUdUmPC6XM_py8,1071
+LoopStructural-1.6.3.dist-info/METADATA,sha256=AQ55C8hbLHvryNgxtgqAnp--xB9ZjKWmrBhg3nGGfCQ,5793
+LoopStructural-1.6.3.dist-info/WHEEL,sha256=OVMc5UfuAQiSplgO0_WdW7vXVGAt9Hdd6qtN4HotdyA,91
+LoopStructural-1.6.3.dist-info/top_level.txt,sha256=QtQErKzYHfg6ddxTQ1NyaTxXBVM6qAqrM_vxEPyXZLg,15
+LoopStructural-1.6.3.dist-info/RECORD,,

{LoopStructural-1.6.2.dist-info → LoopStructural-1.6.3.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (72.1.0)
+Generator: setuptools (75.2.0)
 Root-Is-Purelib: true
 Tag: py3-none-any
 

LoopStructural/modelling/features/fold/_fold_rotation_angle.py

@@ -1,149 +0,0 @@
-import numpy as np
-from scipy.optimize import curve_fit
-
-from ....modelling.features.fold._fold_rotation_angle_feature import (
-    fourier_series,
-)
-from ....modelling.features.fold import SVariogram
-
-from ....utils import getLogger
-
-logger = getLogger(__name__)
-
-
-class FoldRotationAngle:
-    def __init__(self, rotation_angle, fold_frame_coordinate, svario=False):
-        """
-
-        Parameters
-        ----------
-        rotation_angle
-        fold_frame_coordinate
-        svario
-        """
-        self.rotation_angle = rotation_angle
-        self.fold_frame_coordinate = fold_frame_coordinate
-        self.fold_rotation_function = None
-        self.svario = None
-        if svario:
-            self.svario = SVariogram(self.fold_frame_coordinate, self.rotation_angle)
-        self.fitted_params = None
-
-    def fit_fourier_series(
-        self, wl=None, lags=None, nlag=None, lag=None, skip_variogram=False, **kwargs
-    ):
-        """
-
-        Parameters
-        ----------
-        wl
-        lags
-        nlag
-        lag
-
-        Returns
-        -------
-
-        """
-        if self.svario is None:
-            self.svario = SVariogram(self.fold_frame_coordinate, self.rotation_angle)
-        if not skip_variogram:
-            self.svario.calc_semivariogram(lags=lags, nlag=nlag, lag=lag)
-        if wl is None:
-            wl = self.svario.find_wavelengths(lags=lags, nlag=nlag, lag=lag)
-            # for now only consider single fold wavelength
-            wl = wl[0]
-        guess = np.zeros(4)
-        guess[3] = wl  # np.max(limb_wl)
-        logger.info(f"Guess: {guess[0]} {guess[1]} {guess[2]} {guess[3]}")
-        # mask nans
-        mask = np.logical_or(~np.isnan(self.fold_frame_coordinate), ~np.isnan(self.rotation_angle))
-        logger.info(
-            f"There are {np.sum(~mask)} nans for the fold limb rotation angle and { np.sum(mask)} observations"
-        )
-
-        if np.sum(mask) < len(guess):
-            logger.error(
-                "Not enough data points to fit Fourier series setting " "fold rotation angle" "to 0"
-            )
-            self.fold_rotation_function = lambda x: np.zeros(x.shape)
-        else:
-            try:
-                # try fitting using wavelength guess
-                popt, pcov = curve_fit(
-                    fourier_series,
-                    self.fold_frame_coordinate[mask],
-                    np.tan(np.deg2rad(self.rotation_angle[mask])),
-                    guess,
-                )
-            except RuntimeError:
-                try:
-                    # if fitting failed, try with just 0s
-                    logger.info("Running curve fit without initial guess")
-                    popt, pcov = curve_fit(
-                        fourier_series,
-                        self.fold_frame_coordinate[mask],
-                        np.tan(np.deg2rad(self.rotation_angle[mask])),
-                    )
-                except RuntimeError:
-                    # otherwise set the fourier series parameters to 0
-                    popt = guess
-                    logger.error("Could not fit curve to S-Plot, check the wavelength")
-            logger.info(f"Fitted: {popt[0]} {popt[1]} {popt[2]} {popt[3]}")
-            self.fold_rotation_function = lambda x: np.rad2deg(
-                np.arctan(fourier_series(x, popt[0], popt[1], popt[2], popt[3]))
-            )
-            self.fitted_params = popt
-
-    def __call__(self, fold_frame_coordinate):
-        """
-
-        Parameters
-        ----------
-        fold_frame_coordinate
-
-        Returns
-        -------
-
-        """
-        return self.fold_rotation_function(fold_frame_coordinate)
-
-    def calculate_misfit(self):
-        """
-
-        Returns
-        -------
-
-        """
-        return np.tan(np.deg2rad(self.rotation_angle)) - np.tan(
-            np.deg2rad(self.__call__(self.fold_frame_coordinate))
-        )
-
-    def set_function(self, function):
-        """
-
-        Parameters
-        ----------
-        function
-
-        Returns
-        -------
-
-        """
-        self.fold_rotation_function = function
-
-    def find_hinges(self, range, step):
-
-        import scipy.optimize as optimize
-
-        def fra(x):
-            x = np.array([x])
-            return self.__call__(x)
-
-        roots = []
-        x = range[0]
-        while x < range[1]:
-            result = optimize.root_scalar(fra, bracket=[x, x + step])
-            roots.append(result.root)
-            x += step
-        return roots

LoopStructural-1.6.2.dist-info/METADATA

@@ -1,81 +0,0 @@
-Metadata-Version: 2.1
-Name: LoopStructural
-Version: 1.6.2
-Summary: 3D geological modelling
-Author-email: Lachlan Grose <lachlan.grose@monash.edu>
-License: MIT
-Project-URL: Documentation, https://Loop3d.org/LoopStructural/
-Project-URL: Bug Tracker, https://github.com/loop3d/loopstructural/issues
-Project-URL: Source Code, https://github.com/loop3d/loopstructural
-Keywords: earth sciences,geology,3-D modelling,structural geology,uncertainty
-Classifier: Development Status :: 5 - Production/Stable
-Classifier: Intended Audience :: Science/Research
-Classifier: Topic :: Scientific/Engineering :: Information Analysis
-Classifier: License :: OSI Approved :: MIT License
-Classifier: Operating System :: Microsoft :: Windows
-Classifier: Operating System :: POSIX
-Classifier: Operating System :: MacOS
-Classifier: Programming Language :: Python :: 3.8
-Classifier: Programming Language :: Python :: 3.9
-Classifier: Programming Language :: Python :: 3.10
-Classifier: Programming Language :: Python :: 3.11
-Classifier: Programming Language :: Python :: 3.12
-Requires-Python: >=3.8
-Description-Content-Type: text/markdown
-License-File: LICENSE
-Requires-Dist: numpy >=1.18
-Requires-Dist: pandas
-Requires-Dist: scipy
-Requires-Dist: scikit-image
-Requires-Dist: scikit-learn
-Requires-Dist: tqdm
-Provides-Extra: all
-Requires-Dist: loopstructural[export,inequalities,visualisation] ; extra == 'all'
-Provides-Extra: export
-Requires-Dist: geoh5py ; extra == 'export'
-Requires-Dist: pyevtk ; extra == 'export'
-Requires-Dist: dill ; extra == 'export'
-Provides-Extra: inequalities
-Requires-Dist: loopsolver ; extra == 'inequalities'
-Provides-Extra: jupyter
-Requires-Dist: pyvista[all] ; extra == 'jupyter'
-Requires-Dist: tqdm ; extra == 'jupyter'
-Provides-Extra: visualisation
-Requires-Dist: matplotlib ; extra == 'visualisation'
-Requires-Dist: pyvista ; extra == 'visualisation'
-Requires-Dist: loopstructuralviusualisation >0.1.4 ; extra == 'visualisation'
-
-# LoopStructural: Loop3D Geological Forward Modeling Engine.
-![3D model of Hamersley created using loopstructural](docs/source/images/image823.png)
-![Continuous integration and deployment](https://github.com/Loop3D/LoopStructural/workflows/Continuous%20integration%20and%20deployment/badge.svg)
-![Publish Docker Hub](https://github.com/Loop3D/LoopStructural/workflows/Publish%20Docker%20Hub/badge.svg)
-[![PyPI version](https://badge.fury.io/py/LoopStructural.svg)](https://badge.fury.io/py/LoopStructural)
-[![GitHub license](https://img.shields.io/github/license/Naereen/StrapDown.js.svg)](https://github.com/Loop3D/LoopStructural/blob/master/LICENSE)
-[![Documentation loop3d.github.io/LoopStructural/](https://img.shields.io/badge/docs-githubio-brightgreen)](https://loop3d.github.io/LoopStructural)
-
- 
-LoopStructural is the 3D geological modelling library for Loop ([Loop3d.org](https://loop3d.org/)). The development of LoopStructural is lead by **Lachlan Grose** as an ARC (LP170100985) post-doc at Monash University. **Laurent Ailleres** and **Gautier Laurent** have made significant contributions to the conceptual design and integration of geological concepts into the geological modelling workflow. **Roy Thomson** and **Yohan de Rose** have contributed to the implementation and integration of LoopStructural into the Loop workflow. 
-
-Loop is led by Laurent Ailleres (Monash University) with a team of Work Package leaders from:
-* Monash University: Roy Thomson, Lachlan Grose and Robin Armit
-* University of Western Australia: Mark Jessell, Jeremie Giraud, Mark Lindsay and Guillaume Pirot
-* Geological Survey of Canada: Boyan Brodaric and Eric de Kemp
-
-The project benefits from in-kind contributions from the Geological Survey of Canada, the British Geological Survey, the French Bureau de Recherches Geologiques et Minieres, the RING group at the Universite de Lorraine, the RWTH in Aachen, Germany and AUSCOPE
-
-* Python/cython implementation of a Discrete interpolatiors  
-* Fold interpolation using constraints outlined in Laurent 2016 with fold geostatistical tools shown in Grose et al., 2017
-* Fault interpolation 
-
-If you want to use LoopStructural the easiest way to get started is to use a docker container and a jupyter notebook environment
-
-1. Pull the loopstructural docker image `docker pull lachlangrose/loopstructural`
-2. Start a docker container `docker run -it -p 8888:8888 lachlangrose/loopstructural` 
-
-## Documentation
-The LoopStructural documentation can be found [here](https://loop3d.github.io/LoopStructural)
-## Problems
-Any bugs/feature requests/comments please create a new [issue](https://github.com/Loop3D/LoopStructural/issues). 
-
-## Acknowledgements
-*The Loop platform is an open source 3D probabilistic geological and geophysical modelling platform, initiated by Geoscience Australia and the OneGeology consortium. The project is funded by Australian territory, State and Federal Geological Surveys, the Australian Research Council and the MinEx Collaborative Research Centre.*