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'])
Cell Number: mid 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_SeriesSimpleMafi...
EarlyGranite_Series
BIF_Series
SimpleMafic_Series
True
False


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)
Cell Number: -1 Direction: z
<gempy_viewer.modules.plot_2d.visualization_2d.Plot2D object at 0x7f08ac453580>
gpv.plot_2d(geo_model, cell_number=[25], direction='x')
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)
Greenstone
<gempy_viewer.modules.plot_3d.vista.GemPyToVista object at 0x7f08ac657850>

Total running time of the script: (0 minutes 8.624 seconds)

Gallery generated by Sphinx-Gallery