Greenstone.

import os

# Importing gempy
import gempy as gp
import gempy_viewer as gpv

print(gp.__version__)

2023.2.0b1

data_path = os.path.abspath('../../data/input_data/tut_SandStone')

# Importing the data from csv

geo_model: gp.data.GeoModel = gp.create_geomodel(
    project_name='Greenstone',
    extent=[696000, 747000, 6863000, 6930000, -20000, 200],  # * Here we define the extent of the model
    resolution=[20, 20, 20],  # * Here we define the resolution of the voxels
    refinement=6,
    # * Here we define the number of octree levels. If octree levels are defined, the resolution is ignored.
    importer_helper=gp.data.ImporterHelper(
        path_to_orientations=data_path + "/SandStone_Foliations.csv",
        path_to_surface_points=data_path + "/SandStone_Points.csv",
        hash_surface_points=None,
        hash_orientations=None
    )
)

gpv.plot_2d(geo_model, direction=['z'])

<gempy_viewer.modules.plot_2d.visualization_2d.Plot2D object at 0x7f08ac509e40>

gp.map_stack_to_surfaces(
    gempy_model=geo_model,
    mapping_object={
        "EarlyGranite_Series": 'EarlyGranite',
        "BIF_Series": ('SimpleMafic2', 'SimpleBIF'),
        "SimpleMafic_Series": 'SimpleMafic1', 'Basement': 'basement'
    }
)

Could not find element 'basement' in any group.

Structural Groups:	StructuralGroup: Name: EarlyGranite_Series Structural Relation: StackRelationType.ERODE Elements: StructuralElement: Name: EarlyGranite StructuralGroup: Name: BIF_Series Structural Relation: StackRelationType.ERODE Elements: StructuralElement: Name: SimpleMafic2 StructuralElement: Name: SimpleBIF StructuralGroup: Name: SimpleMafic_Series Structural Relation: StackRelationType.ERODE Elements: StructuralElement: Name: SimpleMafic1
Fault Relations:	EarlyGrani... BIF_Series SimpleMafi... EarlyGranite_Series BIF_Series SimpleMafic_Series
	True False

gp.compute_model(geo_model)

Setting Backend To: AvailableBackends.numpy
/home/leguark/gempy/gempy/core/data/geo_model.py:164: UserWarning: You are using refinement and passing a regular grid. The resolution of the regular grid will be overwritten
  warnings.warn(

Solutions: 6 Octree Levels, 4 DualContouringMeshes

gpv.plot_2d(geo_model, cell_number=[-1], direction=['z'], show_data=False)

<gempy_viewer.modules.plot_2d.visualization_2d.Plot2D object at 0x7f08ac453580>

gpv.plot_2d(geo_model, cell_number=[25], direction='x')

<gempy_viewer.modules.plot_2d.visualization_2d.Plot2D object at 0x7f08ac450610>

sphinx_gallery_thumbnail_number = -1

gpv.plot_3d(geo_model)

<gempy_viewer.modules.plot_3d.vista.GemPyToVista object at 0x7f08ac657850>