.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "gallery/3Dvispy/05-surfaces.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_05-surfaces.py>`
        to download the full example code

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

.. _sphx_glr_gallery_3Dvispy_05-surfaces.py:


Display 3D seismic data and multiple horizons
================================================

层位可以是一个 (n1, n2) 大小的 Z 值,
也可以是 (N, 3) 大小的 N 个点

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

.. GENERATED FROM PYTHON SOURCE LINES 17-47

.. code-block:: default


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

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

    sxp = root / 'data/co2/sx.dat'
    sfp1 = root / 'data/co2/mh21.dat'
    sfp2 = root / 'data/co2/mh22.dat'
    ni, nx, nt = 192, 192, 240

    sx = np.fromfile(sxp, np.float32).reshape(ni, nx, nt)
    sf1 = np.fromfile(sfp1, np.float32).reshape(ni, nx)
    sf2 = np.fromfile(sfp2, np.float32).reshape(ni, nx)

    nodes = cigvis.create_slices(sx)

    # show amplitude
    nodes += cigvis.create_surfaces([sf1, sf2],
                                    volume=sx,
                                    value_type='amp',
                                    cmap='Petrel',
                                    clim=[sx.min(), sx.max()])

    # add two points
    nodes += cigvis.create_points(np.array([[70, 50, 158], [20, 100, 80]]), r=3)

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


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

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


.. _sphx_glr_download_gallery_3Dvispy_05-surfaces.py:

.. only:: html

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




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

      :download:`Download Python source code: 05-surfaces.py <05-surfaces.py>`

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

      :download:`Download Jupyter notebook: 05-surfaces.ipynb <05-surfaces.ipynb>`


.. only:: html

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

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