.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "tutorials/ch1_fundamentals/ch1_3b_cross_sections.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_tutorials_ch1_fundamentals_ch1_3b_cross_sections.py>`
        to download the full example code.

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_tutorials_ch1_fundamentals_ch1_3b_cross_sections.py:


1.3b: 2-D sections
==================

.. GENERATED FROM PYTHON SOURCE LINES 8-9

Importing

.. GENERATED FROM PYTHON SOURCE LINES 9-15

.. code-block:: Python

    import gempy as gp
    import gempy_viewer as gpv
    import numpy as np

    np.random.seed(1234)








.. GENERATED FROM PYTHON SOURCE LINES 16-19

Setup the model
---------------


.. GENERATED FROM PYTHON SOURCE LINES 21-22

Importing the data from CSV-files and setting extent and resolution

.. GENERATED FROM PYTHON SOURCE LINES 22-51

.. code-block:: Python

    data_path = 'https://raw.githubusercontent.com/cgre-aachen/gempy_data/master/'

    geo_model: gp.data.GeoModel = gp.create_geomodel(
        project_name='Tutorial_ch1_1_Basics',
        extent=[0, 2000, 0, 2000, 0, 750],
        resolution=[20, 20, 20],  # * Here we define the resolution of the voxels
        refinement=4,  # * 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 + "/data/input_data/getting_started/simple_fault_model_orientations.csv",
            path_to_surface_points=data_path + "/data/input_data/getting_started/simple_fault_model_points.csv",
            hash_surface_points="4cdd54cd510cf345a583610585f2206a2936a05faaae05595b61febfc0191563",
            hash_orientations="7ba1de060fc8df668d411d0207a326bc94a6cdca9f5fe2ed511fd4db6b3f3526"
        )
    )

    gp.map_stack_to_surfaces(
        gempy_model=geo_model,
        mapping_object=  # TODO: This mapping I do not like it too much. We should be able to do it passing the data objects directly
        {
            "Fault_Series": 'Main_Fault',
            "Strat_Series": ('Sandstone_2', 'Siltstone', 'Shale', 'Sandstone_1')
        }
    )

    gp.set_is_fault(
        frame=geo_model.structural_frame,
        fault_groups=['Fault_Series']
    )





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    Surface points hash:  4cdd54cd510cf345a583610585f2206a2936a05faaae05595b61febfc0191563
    Orientations hash:  7ba1de060fc8df668d411d0207a326bc94a6cdca9f5fe2ed511fd4db6b3f3526


.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">

            <table>
              <tr><td>Structural Groups:</td><td>
        <table style="border-left:1.2px solid black;>
          <tr><th colspan="2"><b>StructuralGroup:</b></th></tr>
          <tr><td>Name:</td><td>Fault_Series</td></tr>
          <tr><td>Structural Relation:</td><td>StackRelationType.FAULT</td></tr>
          <tr><td>Elements:</td><td>
        <table width="50%" style="border-left:15px solid #527682;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>Main_Fault</td></tr>
        </table>
            </td></tr>
        </table>
            <br>
        <table style="border-left:1.2px solid black;>
          <tr><th colspan="2"><b>StructuralGroup:</b></th></tr>
          <tr><td>Name:</td><td>Strat_Series</td></tr>
          <tr><td>Structural Relation:</td><td>StackRelationType.ERODE</td></tr>
          <tr><td>Elements:</td><td>
        <table width="50%" style="border-left:15px solid #ffbe00;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>Sandstone_2</td></tr>
        </table>
            <br>
        <table width="50%" style="border-left:15px solid #443988;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>Siltstone</td></tr>
        </table>
            <br>
        <table width="50%" style="border-left:15px solid #728f02;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>Shale</td></tr>
        </table>
            <br>
        <table width="50%" style="border-left:15px solid #9f0052;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>Sandstone_1</td></tr>
        </table>
            </td></tr>
        </table>
            </td></tr>
              <tr><td>Fault Relations:</td><td><table style="border-collapse: collapse; table-layout: fixed;"><th></th><th style="transform: rotate(-35deg); height:150px; vertical-align: bottom; text-align: center;">Fault_Seri...</th><th style="transform: rotate(-35deg); height:150px; vertical-align: bottom; text-align: center;">Strat_Seri...</th><tr><th>Fault_Series</th><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #527682; width: 20px; height: 20px; border: 1px solid black;"></td></tr><tr><th>Strat_Series</th><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td></tr></table></td></tr>
              <tr><td></td><td>
            <table>
              <tr>
                <td><div style="display: inline-block; background-color: #527682; width: 20px; height: 20px; border: 1px solid black;"></div> True</td>
                <td><div style="display: inline-block; background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></div> False</td>
              </tr>
            </table>
            </td></tr>
            </table>
        
    </div>
    <br />
    <br />

.. GENERATED FROM PYTHON SOURCE LINES 52-55

Add sections
~~~~~~~~~~~~


.. GENERATED FROM PYTHON SOURCE LINES 58-61

pass section dictionary with startpoint, endpoint and resolution for
every section:


.. GENERATED FROM PYTHON SOURCE LINES 63-72

.. code-block:: Python

    gp.set_section_grid(
        grid=geo_model.grid,
        section_dict={
            'section1': ([0, 0], [2000, 2000], [100, 80]),
            'section2': ([800, 0], [800, 2000], [150, 100]),
            'section3': ([0, 200], [1500, 500], [200, 150])
        }  # p1,p2,resolution
    )





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    Active grids: GridTypes.DENSE|SECTIONS|NONE


.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">
    <table border="1" class="dataframe">
      <thead>
        <tr style="text-align: right;">
          <th></th>
          <th>start</th>
          <th>stop</th>
          <th>resolution</th>
          <th>dist</th>
        </tr>
      </thead>
      <tbody>
        <tr>
          <th>section1</th>
          <td>(0.0, 0.0)</td>
          <td>(2000.0, 2000.0)</td>
          <td>(100, 80)</td>
          <td>2828.427125</td>
        </tr>
        <tr>
          <th>section2</th>
          <td>(800.0, 0.0)</td>
          <td>(800.0, 2000.0)</td>
          <td>(150, 100)</td>
          <td>2000.000000</td>
        </tr>
        <tr>
          <th>section3</th>
          <td>(0.0, 200.0)</td>
          <td>(1500.0, 500.0)</td>
          <td>(200, 150)</td>
          <td>1529.705854</td>
        </tr>
      </tbody>
    </table>
    </div>
    <br />
    <br />

.. GENERATED FROM PYTHON SOURCE LINES 73-76

Add topography
~~~~~~~~~~~~~~


.. GENERATED FROM PYTHON SOURCE LINES 78-85

.. code-block:: Python

    gp.set_topography_from_random(
        grid=geo_model.grid,
        fractal_dimension=1.2,
        d_z=np.array([300, 750]),
        topography_resolution=np.array([50, 50])
    )





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    Active grids: GridTypes.DENSE|TOPOGRAPHY|SECTIONS|NONE

    Topography(_regular_grid=RegularGrid(resolution=array([20, 20, 20]), extent=array([   0., 2000.,    0., 2000.,    0.,  750.]), values=array([[  50.  ,   50.  ,   18.75],
           [  50.  ,   50.  ,   56.25],
           [  50.  ,   50.  ,   93.75],
           ...,
           [1950.  , 1950.  ,  656.25],
           [1950.  , 1950.  ,  693.75],
           [1950.  , 1950.  ,  731.25]], shape=(8000, 3)), mask_topo=array([], shape=(0, 3), dtype=bool), _transform=None, _base_resolution=array([2, 2, 2])), values_2d=array([[[   0.        ,    0.        ,  525.81364595],
            [   0.        ,   40.81632653,  534.4311714 ],
            [   0.        ,   81.63265306,  545.21138351],
            ...,
            [   0.        , 1918.36734694,  493.45417828],
            [   0.        , 1959.18367347,  505.51902125],
            [   0.        , 2000.        ,  516.62320964]],

           [[  40.81632653,    0.        ,  520.02419799],
            [  40.81632653,   40.81632653,  529.19849936],
            [  40.81632653,   81.63265306,  540.73889549],
            ...,
            [  40.81632653, 1918.36734694,  485.70563878],
            [  40.81632653, 1959.18367347,  498.13591443],
            [  40.81632653, 2000.        ,  509.94870633]],

           [[  81.63265306,    0.        ,  512.29034438],
            [  81.63265306,   40.81632653,  522.38242549],
            [  81.63265306,   81.63265306,  533.74749759],
            ...,
            [  81.63265306, 1918.36734694,  475.96047534],
            [  81.63265306, 1959.18367347,  488.33556325],
            [  81.63265306, 2000.        ,  500.75092545]],

           ...,

           [[1918.36734694,    0.        ,  536.8572301 ],
            [1918.36734694,   40.81632653,  549.32197157],
            [1918.36734694,   81.63265306,  561.72341993],
            ...,
            [1918.36734694, 1918.36734694,  503.24054915],
            [1918.36734694, 1959.18367347,  513.73925164],
            [1918.36734694, 2000.        ,  525.47837607]],

           [[1959.18367347,    0.        ,  534.68908826],
            [1959.18367347,   40.81632653,  545.65518713],
            [1959.18367347,   81.63265306,  558.27313536],
            ...,
            [1959.18367347, 1918.36734694,  501.0502999 ],
            [1959.18367347, 1959.18367347,  512.95973605],
            [1959.18367347, 2000.        ,  524.20361683]],

           [[2000.        ,    0.        ,  530.78028142],
            [2000.        ,   40.81632653,  540.15206274],
            [2000.        ,   81.63265306,  551.79148373],
            ...,
            [2000.        , 1918.36734694,  498.47649577],
            [2000.        , 1959.18367347,  510.05476038],
            [2000.        , 2000.        ,  521.16890329]]], shape=(50, 50, 3)), source=None, values=array([[   0.        ,    0.        ,  525.81364595],
           [   0.        ,   40.81632653,  534.4311714 ],
           [   0.        ,   81.63265306,  545.21138351],
           ...,
           [2000.        , 1918.36734694,  498.47649577],
           [2000.        , 1959.18367347,  510.05476038],
           [2000.        , 2000.        ,  521.16890329]], shape=(2500, 3)), resolution=(50, 50), raster_shape=())



.. GENERATED FROM PYTHON SOURCE LINES 86-88

Active grids:


.. GENERATED FROM PYTHON SOURCE LINES 90-92

.. code-block:: Python

    geo_model.grid.active_grids





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <GridTypes.DENSE|TOPOGRAPHY|SECTIONS|NONE: 1050>



.. GENERATED FROM PYTHON SOURCE LINES 93-95

.. code-block:: Python

    gpv.plot_section_traces(geo_model)




.. image-sg:: /tutorials/ch1_fundamentals/images/sphx_glr_ch1_3b_cross_sections_001.png
   :alt: Cell Number: -1 Direction: z
   :srcset: /tutorials/ch1_fundamentals/images/sphx_glr_ch1_3b_cross_sections_001.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <function plot_section_traces at 0x7f582f2ef8a0>



.. GENERATED FROM PYTHON SOURCE LINES 98-101

.. code-block:: Python

    geo_model.interpolation_options.mesh_extraction = False
    sol = gp.compute_model(geo_model)





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    Setting Backend To: AvailableBackends.PYTORCH
    GPU requested but unavailable; falling back to CPU (GEMPY_GPU_FALLBACK=True)
    Setting Backend To: AvailableBackends.PYTORCH
    Chunking done: 8 chunks




.. GENERATED FROM PYTHON SOURCE LINES 102-104

.. code-block:: Python

    gpv.plot_2d(geo_model, section_names=['topography'])




.. image-sg:: /tutorials/ch1_fundamentals/images/sphx_glr_ch1_3b_cross_sections_002.png
   :alt: Geological map
   :srcset: /tutorials/ch1_fundamentals/images/sphx_glr_ch1_3b_cross_sections_002.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    /opt/buildAgent/work/3a8738c25f60c3c9/venv/lib/python3.14/site-packages/gempy_viewer/API/_plot_2d_sections_api.py:112: UserWarning: Section contacts not implemented yet. We need to pass scalar field for the sections grid
      warnings.warn(

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



.. GENERATED FROM PYTHON SOURCE LINES 105-107

.. code-block:: Python

    gpv.plot_2d(geo_model, section_names=['section1'])




.. image-sg:: /tutorials/ch1_fundamentals/images/sphx_glr_ch1_3b_cross_sections_003.png
   :alt: section1
   :srcset: /tutorials/ch1_fundamentals/images/sphx_glr_ch1_3b_cross_sections_003.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    /opt/buildAgent/work/3a8738c25f60c3c9/venv/lib/python3.14/site-packages/gempy_viewer/API/_plot_2d_sections_api.py:112: UserWarning: Section contacts not implemented yet. We need to pass scalar field for the sections grid
      warnings.warn(

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



.. GENERATED FROM PYTHON SOURCE LINES 108-116

.. code-block:: Python

    gpv.plot_2d(
        model=geo_model,
        section_names=['section1', 'section2', 'section3', 'topography'],
        show_topography=True
    )


    # sphinx_gallery_thumbnail_number = 4



.. image-sg:: /tutorials/ch1_fundamentals/images/sphx_glr_ch1_3b_cross_sections_004.png
   :alt: section1, section2, section3, Geological map
   :srcset: /tutorials/ch1_fundamentals/images/sphx_glr_ch1_3b_cross_sections_004.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    /opt/buildAgent/work/3a8738c25f60c3c9/venv/lib/python3.14/site-packages/gempy_viewer/API/_plot_2d_sections_api.py:112: UserWarning: Section contacts not implemented yet. We need to pass scalar field for the sections grid
      warnings.warn(

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




.. rst-class:: sphx-glr-timing

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


.. _sphx_glr_download_tutorials_ch1_fundamentals_ch1_3b_cross_sections.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: ch1_3b_cross_sections.ipynb <ch1_3b_cross_sections.ipynb>`

    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: ch1_3b_cross_sections.py <ch1_3b_cross_sections.py>`

    .. container:: sphx-glr-download sphx-glr-download-zip

      :download:`Download zipped: ch1_3b_cross_sections.zip <ch1_3b_cross_sections.zip>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_