.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "_examples/01_ERT/plot_02_syscal_rollalong.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_01_ERT_plot_02_syscal_rollalong.py: Importing Syscal ERT data from roll-a-long scheme ================================================= .. GENERATED FROM PYTHON SOURCE LINES 8-9 .. code-block:: Python import reda .. GENERATED FROM PYTHON SOURCE LINES 10-16 create an ERT container and import first dataset The 'elecs_transform_reg_spacing_x' parameter transforms the spacing in x direction from the spacing stored in the data file (here 1 m) to the true spacing of 0.2 m. This is useful if the spacing on the measurement system is not changed between measurements and only changed in the postprocessing (this is common practice in some work groups). .. GENERATED FROM PYTHON SOURCE LINES 16-27 .. code-block:: Python ert_p1 = reda.ERT() ert_p1.import_syscal_bin( 'data_syscal_rollalong/profile_1.bin', elecs_transform_reg_spacing_x=(1, 0.2), ) ert_p1.compute_K_analytical(spacing=0.2) with reda.CreateEnterDirectory('output_02_rollalong'): ert_p1.pseudosection( filename='profile1_pseudosection.pdf', column='r', log10=True ) .. image-sg:: /_examples/01_ERT/images/sphx_glr_plot_02_syscal_rollalong_001.png :alt: plot 02 syscal rollalong :srcset: /_examples/01_ERT/images/sphx_glr_plot_02_syscal_rollalong_001.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none found duplicate coordinates! .. GENERATED FROM PYTHON SOURCE LINES 28-29 Print statistics .. GENERATED FROM PYTHON SOURCE LINES 29-31 .. code-block:: Python ert_p1.print_data_journal() .. rst-class:: sphx-glr-script-out .. code-block:: none --- Data Journal Start --- 2024-07-09 12:19:43.474520 Data was imported from file (1460 data points) --- Data Journal End --- .. GENERATED FROM PYTHON SOURCE LINES 32-33 This here is an activity list .. GENERATED FROM PYTHON SOURCE LINES 33-35 .. code-block:: Python ert_p1.print_log() .. rst-class:: sphx-glr-script-out .. code-block:: none 2024-07-09 12:19:42,375 - reda.containers.BaseContainer - INFO - IRIS Syscal Pro bin import 2024-07-09 12:19:43,203 - reda.containers.BaseContainer - INFO - Data sized changed from 0 to 1460 .. GENERATED FROM PYTHON SOURCE LINES 36-37 create an ERT container and import second dataset .. GENERATED FROM PYTHON SOURCE LINES 37-49 .. code-block:: Python ert_p2 = reda.ERT() ert_p2.import_syscal_bin( 'data_syscal_rollalong/profile_2.bin', elecs_transform_reg_spacing_x=(1, 0.2), ) ert_p2.compute_K_analytical(spacing=0.2) with reda.CreateEnterDirectory('output_02_rollalong'): ert_p2.pseudosection( filename='profile2_pseudosection.pdf', column='r', log10=True ) .. image-sg:: /_examples/01_ERT/images/sphx_glr_plot_02_syscal_rollalong_002.png :alt: plot 02 syscal rollalong :srcset: /_examples/01_ERT/images/sphx_glr_plot_02_syscal_rollalong_002.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none found duplicate coordinates! .. GENERATED FROM PYTHON SOURCE LINES 50-52 Again print the data journal to see what changed in terms of imported data points .. GENERATED FROM PYTHON SOURCE LINES 52-54 .. code-block:: Python ert_p2.print_data_journal() .. rst-class:: sphx-glr-script-out .. code-block:: none --- Data Journal Start --- 2024-07-09 12:19:44.600429 Data was imported from file (1460 data points) --- Data Journal End --- .. GENERATED FROM PYTHON SOURCE LINES 55-59 Now we start over again and fix the coordinates of the second profile by using the 'shift_by_xyz' parameter, which can be used to move the electrode positions by a fixed offset. In our case this offset amounts to the position of the 25th position (first electrode starts at 0, therefore 24 * 0.25 m). .. GENERATED FROM PYTHON SOURCE LINES 59-75 .. code-block:: Python ert = reda.ERT() # first profile ert.import_syscal_bin( 'data_syscal_rollalong/profile_1.bin', elecs_transform_reg_spacing_x=(1, 0.2), ) ert2 = reda.ERT() # second profile ert2.import_syscal_bin( 'data_syscal_rollalong/profile_2.bin', elecs_transform_reg_spacing_x=(1, 0.2), shift_by_xyz=[24 * 0.2], ) .. GENERATED FROM PYTHON SOURCE LINES 76-77 Electrode positions of first dataset .. GENERATED FROM PYTHON SOURCE LINES 77-80 .. code-block:: Python with reda.CreateEnterDirectory('output_02_rollalong'): fig, ax = ert.plot_electrode_positions_2d() fig.show() .. image-sg:: /_examples/01_ERT/images/sphx_glr_plot_02_syscal_rollalong_003.png :alt: plot 02 syscal rollalong :srcset: /_examples/01_ERT/images/sphx_glr_plot_02_syscal_rollalong_003.png :class: sphx-glr-single-img .. GENERATED FROM PYTHON SOURCE LINES 81-82 Electrode positions of second dataset .. GENERATED FROM PYTHON SOURCE LINES 82-85 .. code-block:: Python with reda.CreateEnterDirectory('output_02_rollalong'): fig, ax = ert2.plot_electrode_positions_2d() fig.show() .. image-sg:: /_examples/01_ERT/images/sphx_glr_plot_02_syscal_rollalong_004.png :alt: plot 02 syscal rollalong :srcset: /_examples/01_ERT/images/sphx_glr_plot_02_syscal_rollalong_004.png :class: sphx-glr-single-img .. GENERATED FROM PYTHON SOURCE LINES 86-87 Now we merge the datasets and plot the electrode positions .. GENERATED FROM PYTHON SOURCE LINES 87-93 .. code-block:: Python ert.merge_container(ert2) with reda.CreateEnterDirectory('output_02_rollalong'): fig, ax = ert.plot_electrode_positions_2d() fig.show() .. image-sg:: /_examples/01_ERT/images/sphx_glr_plot_02_syscal_rollalong_005.png :alt: plot 02 syscal rollalong :srcset: /_examples/01_ERT/images/sphx_glr_plot_02_syscal_rollalong_005.png :class: sphx-glr-single-img .. GENERATED FROM PYTHON SOURCE LINES 94-95 compute geometric factors .. GENERATED FROM PYTHON SOURCE LINES 95-97 .. code-block:: Python ert.compute_K_analytical(spacing=0.2) .. rst-class:: sphx-glr-script-out .. code-block:: none array([ -3.76991118, -3.76991118, -15.07964474, ..., 15.07964474, -3.76991118, 3.76991118]) .. GENERATED FROM PYTHON SOURCE LINES 98-99 Plot Pseudosection .. GENERATED FROM PYTHON SOURCE LINES 99-103 .. code-block:: Python with reda.CreateEnterDirectory('output_02_rollalong'): ert.pseudosection(filename='pseudosection_both_profiles.pdf', column='r', log10=True) .. image-sg:: /_examples/01_ERT/images/sphx_glr_plot_02_syscal_rollalong_006.png :alt: plot 02 syscal rollalong :srcset: /_examples/01_ERT/images/sphx_glr_plot_02_syscal_rollalong_006.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none found duplicate coordinates! .. GENERATED FROM PYTHON SOURCE LINES 104-107 .. code-block:: Python with reda.CreateEnterDirectory('output_02_rollalong'): ert.histogram(filename='hist_both_profiles.pdf', column=['r', 'rho_a', 'Iab', ]) .. image-sg:: /_examples/01_ERT/images/sphx_glr_plot_02_syscal_rollalong_007.png :alt: plot 02 syscal rollalong :srcset: /_examples/01_ERT/images/sphx_glr_plot_02_syscal_rollalong_007.png :class: sphx-glr-single-img .. _sphx_glr_download__examples_01_ERT_plot_02_syscal_rollalong.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: plot_02_syscal_rollalong.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: plot_02_syscal_rollalong.py ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_