.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "gallery/3Dvispy/04-overlay_rgt_fault.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

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

        :ref:`Go to the end <sphx_glr_download_gallery_3Dvispy_04-overlay_rgt_fault.py>`
        to download the full example code

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

.. _sphx_glr_gallery_3Dvispy_04-overlay_rgt_fault.py:


Overlay multiple 3D data bodies (RGT, fault) on slices of 3D seismic data bodies
====================================================================================

``create_overlay`` 的第一个参数是背景(background), 
第二个参数是叠加的前景(foreground), 其可以包含多个3D数据体

.. Note::
    foreground 需要合理设置透明度和mask

.. image:: ../../_static/cigvis/3Dvispy/04.png
    :alt: image
    :align: center

.. GENERATED FROM PYTHON SOURCE LINES 20-51

.. code-block:: default


    # sphinx_gallery_thumbnail_path = '_static/cigvis/3Dvispy/04.png'

    import numpy as np
    import cigvis
    from cigvis import colormap
    from pathlib import Path
    root = Path(__file__).resolve().parent.parent.parent

    sxp = root / 'data/rgt/sx.dat'
    uxp = root / 'data/rgt/ux.dat'
    fxp = root / 'data/rgt/fx.dat'
    ni, nx, nt = 128, 128, 128

    sx = np.fromfile(sxp, np.float32).reshape(ni, nx, nt)
    rgt = np.fromfile(uxp, np.float32).reshape(ni, nx, nt)
    fx = np.fromfile(fxp, np.float32).reshape(ni, nx, nt)

    rgt_cmap = colormap.set_alpha('jet', 0.4)
    # mask min value (0), 0 means no fault
    fx_cmap = colormap.set_alpha_except_min('jet', alpha=1)
    fg_cmap = [rgt_cmap, fx_cmap]

    # fx is discrete data, set interpolation as 'nearest'
    nodes = cigvis.create_slices(sx, pos=[[36], [28], [84]], cmap='gray')
    nodes = cigvis.add_mask(nodes, [rgt, fx],
                            cmaps=[rgt_cmap, fx_cmap],
                            interpolation=['cubic', 'nearest'])
    nodes += cigvis.create_colorbar_from_nodes(nodes, 'RGT', select='mask', idx=0) # idx = 0 means the first mask

    cigvis.plot3D(nodes, size=(750, 600), savename='example.png')


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

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


.. _sphx_glr_download_gallery_3Dvispy_04-overlay_rgt_fault.py:

.. only:: html

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




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

      :download:`Download Python source code: 04-overlay_rgt_fault.py <04-overlay_rgt_fault.py>`

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

      :download:`Download Jupyter notebook: 04-overlay_rgt_fault.ipynb <04-overlay_rgt_fault.ipynb>`


.. only:: html

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

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