.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "examples/real/Alesmodel.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code .. rst-class:: sphx-glr-example-title .. _sphx_glr_examples_real_Alesmodel.py: Alesmodel: Plotting sections and maps. ====================================== # %% # .. admonition:: Explanation # # This model is quite unstable in general and requires float64 to find a solution. The lack of data on # one of the corners for the TRIAS and LIAS series makes that the model bends in an unrealistic # way and erodes CARBO that disappears on that section. The easy way to solve this is to add more data in that area but # I leave as it is since I did no constructed the model. # .. GENERATED FROM PYTHON SOURCE LINES 15-22 .. code-block:: Python import gempy as gp import gempy_viewer as gpv import os import numpy as np .. GENERATED FROM PYTHON SOURCE LINES 23-33 .. code-block:: Python cwd = os.getcwd() if 'examples' not in cwd: data_path = os.getcwd() + '/examples' else: data_path = cwd + '/../..' path_interf = data_path + "/data/input_data/AlesModel/2018_interf.csv" path_orient = data_path + "/data/input_data/AlesModel/2018_orient_clust_n_init5_0.csv" path_dem = data_path + "/data/input_data/AlesModel/_cropped_DEM_coarse.tif" .. GENERATED FROM PYTHON SOURCE LINES 34-45 .. code-block:: Python geo_model: gp.data.GeoModel = gp.create_geomodel( project_name='Claudius', extent=[729550.0, 751500.0, 1913500.0, 1923650.0, -1800.0, 800.0], resolution=None, refinement=6, importer_helper=gp.data.ImporterHelper( path_to_orientations=path_orient, path_to_surface_points=path_interf, ) ) .. GENERATED FROM PYTHON SOURCE LINES 46-53 .. code-block:: Python gp.set_section_grid( grid=geo_model.grid, section_dict={ 'section1': ([732000, 1916000], [745000, 1916000], [200, 150]) } ) .. rst-class:: sphx-glr-script-out .. code-block:: none Active grids: ['sections'] .. raw:: html
start stop resolution dist
section1 [732000, 1916000] [745000, 1916000] [200, 150] 13000.0


.. GENERATED FROM PYTHON SOURCE LINES 54-55 sorting of lithologies .. GENERATED FROM PYTHON SOURCE LINES 55-68 .. code-block:: Python gp.map_stack_to_surfaces( gempy_model=geo_model, mapping_object={ 'fault_left' : 'fault_left', 'fault_right' : 'fault_right', 'fault_lr' : 'fault_lr', 'Trias_Series' : ('TRIAS', 'LIAS'), 'Carbon_Series' : 'CARBO', 'Basement_Series': 'basement' }, remove_unused_series=True ) .. rst-class:: sphx-glr-script-out .. code-block:: none Could not find element 'basement' in any group. .. raw:: html
Structural Groups: StructuralGroup:
Name:fault_left
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:fault_left

StructuralGroup:
Name:fault_right
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:fault_right

StructuralGroup:
Name:fault_lr
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:fault_lr

StructuralGroup:
Name:Trias_Series
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:TRIAS

StructuralElement:
Name:LIAS

StructuralGroup:
Name:Carbon_Series
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:CARBO
Fault Relations:
fault_leftfault_righ...fault_lrTrias_Seri...Carbon_Ser...
fault_left
fault_right
fault_lr
Trias_Series
Carbon_Series
True
False


.. GENERATED FROM PYTHON SOURCE LINES 69-70 Change colors .. GENERATED FROM PYTHON SOURCE LINES 70-74 .. code-block:: Python geo_model.structural_frame.get_element_by_name("LIAS").color = "#015482" geo_model.structural_frame.get_element_by_name("TRIAS").color = "#9f0052" geo_model.structural_frame.get_element_by_name("CARBO").color = "#ffbe00" .. GENERATED FROM PYTHON SOURCE LINES 75-77 .. code-block:: Python a = gpv.plot_2d(geo_model, direction='y') .. image-sg:: /examples/real/images/sphx_glr_Alesmodel_001.png :alt: Cell Number: mid Direction: y :srcset: /examples/real/images/sphx_glr_Alesmodel_001.png :class: sphx-glr-single-img .. GENERATED FROM PYTHON SOURCE LINES 78-80 .. code-block:: Python gpv.plot_section_traces(geo_model) .. image-sg:: /examples/real/images/sphx_glr_Alesmodel_002.png :alt: Cell Number: -1 Direction: z :srcset: /examples/real/images/sphx_glr_Alesmodel_002.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none .. GENERATED FROM PYTHON SOURCE LINES 81-84 Faults '''''' .. GENERATED FROM PYTHON SOURCE LINES 86-96 .. code-block:: Python gp.set_is_fault( frame=geo_model.structural_frame, fault_groups=[ geo_model.structural_frame.get_group_by_name('fault_left'), geo_model.structural_frame.get_group_by_name('fault_right'), geo_model.structural_frame.get_group_by_name('fault_lr') ], change_color=True ) .. raw:: html
Structural Groups: StructuralGroup:
Name:fault_left
Structural Relation:StackRelationType.FAULT
Elements:
StructuralElement:
Name:fault_left

StructuralGroup:
Name:fault_right
Structural Relation:StackRelationType.FAULT
Elements:
StructuralElement:
Name:fault_right

StructuralGroup:
Name:fault_lr
Structural Relation:StackRelationType.FAULT
Elements:
StructuralElement:
Name:fault_lr

StructuralGroup:
Name:Trias_Series
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:TRIAS

StructuralElement:
Name:LIAS

StructuralGroup:
Name:Carbon_Series
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:CARBO
Fault Relations:
fault_leftfault_righ...fault_lrTrias_Seri...Carbon_Ser...
fault_left
fault_right
fault_lr
Trias_Series
Carbon_Series
True
False


.. GENERATED FROM PYTHON SOURCE LINES 97-100 Topography ~~~~~~~~~~ .. GENERATED FROM PYTHON SOURCE LINES 102-111 .. code-block:: Python gp.set_topography_from_file( grid=geo_model.grid, filepath=path_dem, crop_to_extent=[729550.0, 751500.0, 1_913_500.0, 1923650.0] ) gpv.plot_3d(geo_model, show_topography=True, ve=1, image=True) .. image-sg:: /examples/real/images/sphx_glr_Alesmodel_003.png :alt: Alesmodel :srcset: /examples/real/images/sphx_glr_Alesmodel_003.png :class: sphx-glr-single-img .. image-sg:: /examples/real/images/sphx_glr_Alesmodel_004.png :alt: Alesmodel :srcset: /examples/real/images/sphx_glr_Alesmodel_004.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none Active grids: ['topography' 'sections'] .. GENERATED FROM PYTHON SOURCE LINES 112-114 .. code-block:: Python carbo = geo_model.structural_frame.get_group_by_name("Carbon_Series") .. GENERATED FROM PYTHON SOURCE LINES 115-123 .. code-block:: Python geo_model.interpolation_options.number_octree_levels_surface = 4 geo_model.interpolation_options.kernel_options.range = .8 gp.modify_surface_points( geo_model=geo_model, elements_names=["CARBO", "LIAS", "TRIAS"], nugget=0.005 ) .. raw:: html
Structural Groups: StructuralGroup:
Name:fault_left
Structural Relation:StackRelationType.FAULT
Elements:
StructuralElement:
Name:fault_left

StructuralGroup:
Name:fault_right
Structural Relation:StackRelationType.FAULT
Elements:
StructuralElement:
Name:fault_right

StructuralGroup:
Name:fault_lr
Structural Relation:StackRelationType.FAULT
Elements:
StructuralElement:
Name:fault_lr

StructuralGroup:
Name:Trias_Series
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:TRIAS

StructuralElement:
Name:LIAS

StructuralGroup:
Name:Carbon_Series
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:CARBO
Fault Relations:
fault_leftfault_righ...fault_lrTrias_Seri...Carbon_Ser...
fault_left
fault_right
fault_lr
Trias_Series
Carbon_Series
True
False


.. GENERATED FROM PYTHON SOURCE LINES 124-127 .. code-block:: Python print(geo_model.structural_frame) geo_model.structural_frame .. rst-class:: sphx-glr-script-out .. code-block:: none StructuralFrame( structural_groups=[ StructuralGroup( name=fault_left, structural_relation=StackRelationType.FAULT, elements=[ Element( name=fault_left, color=#527682, is_active=True ) ] ), StructuralGroup( name=fault_right, structural_relation=StackRelationType.FAULT, elements=[ Element( name=fault_right, color=#527682, is_active=True ) ] ), StructuralGroup( name=fault_lr, structural_relation=StackRelationType.FAULT, elements=[ Element( name=fault_lr, color=#527682, is_active=True ) ] ), StructuralGroup( name=Trias_Series, structural_relation=StackRelationType.ERODE, elements=[ Element( name=TRIAS, color=#9f0052, is_active=True ), Element( name=LIAS, color=#015482, is_active=True ) ] ), StructuralGroup( name=Carbon_Series, structural_relation=StackRelationType.ERODE, elements=[ Element( name=CARBO, color=#ffbe00, is_active=True ) ] ) ], fault_relations= [[False, False, False, True, True], [False, False, False, True, True], [False, False, False, True, True], [False, False, False, False, False], [False, False, False, False, False]], .. raw:: html
Structural Groups: StructuralGroup:
Name:fault_left
Structural Relation:StackRelationType.FAULT
Elements:
StructuralElement:
Name:fault_left

StructuralGroup:
Name:fault_right
Structural Relation:StackRelationType.FAULT
Elements:
StructuralElement:
Name:fault_right

StructuralGroup:
Name:fault_lr
Structural Relation:StackRelationType.FAULT
Elements:
StructuralElement:
Name:fault_lr

StructuralGroup:
Name:Trias_Series
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:TRIAS

StructuralElement:
Name:LIAS

StructuralGroup:
Name:Carbon_Series
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:CARBO
Fault Relations:
fault_leftfault_righ...fault_lrTrias_Seri...Carbon_Ser...
fault_left
fault_right
fault_lr
Trias_Series
Carbon_Series
True
False


.. GENERATED FROM PYTHON SOURCE LINES 128-140 .. code-block:: Python geo_model.interpolation_options.mesh_extraction = False _ = gp.compute_model( geo_model, engine_config=gp.data.GemPyEngineConfig( backend=gp.data.AvailableBackends.PYTORCH, use_gpu=True, dtype="float64" )) gpv.plot_2d(geo_model, show_topography=False, section_names=['topography'], show_lith=True) .. image-sg:: /examples/real/images/sphx_glr_Alesmodel_005.png :alt: Geological map :srcset: /examples/real/images/sphx_glr_Alesmodel_005.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none Setting Backend To: AvailableBackends.PYTORCH .. GENERATED FROM PYTHON SOURCE LINES 141-144 .. code-block:: Python gpv.plot_2d(geo_model, cell_number=[4], direction=['y'], show_topography=True, show_data=True) gpv.plot_2d(geo_model, cell_number=[-4], direction=['y'], show_topography=True, show_data=True) .. rst-class:: sphx-glr-horizontal * .. image-sg:: /examples/real/images/sphx_glr_Alesmodel_006.png :alt: Cell Number: 4 Direction: y :srcset: /examples/real/images/sphx_glr_Alesmodel_006.png :class: sphx-glr-multi-img * .. image-sg:: /examples/real/images/sphx_glr_Alesmodel_007.png :alt: Cell Number: -4 Direction: y :srcset: /examples/real/images/sphx_glr_Alesmodel_007.png :class: sphx-glr-multi-img .. rst-class:: sphx-glr-script-out .. code-block:: none .. GENERATED FROM PYTHON SOURCE LINES 145-146 sphinx_gallery_thumbnail_number = -1 .. GENERATED FROM PYTHON SOURCE LINES 146-147 .. code-block:: Python gpv.plot_3d(geo_model, show_lith=True, show_topography=True, kwargs_plot_structured_grid={'opacity': 0.8}) .. image-sg:: /examples/real/images/sphx_glr_Alesmodel_008.png :alt: Alesmodel :srcset: /examples/real/images/sphx_glr_Alesmodel_008.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none .. rst-class:: sphx-glr-timing **Total running time of the script:** (0 minutes 31.237 seconds) .. _sphx_glr_download_examples_real_Alesmodel.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: Alesmodel.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: Alesmodel.py ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_