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

.. only:: html

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

        Click :ref:`here <sphx_glr_download_examples_Data_conversion_plot_tifs_to_netcdf.py>`
        to download the full example code

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

.. _sphx_glr_examples_Data_conversion_plot_tifs_to_netcdf.py:


Multiple TIFs to NetCDF conversion
----------------------------------
Dataset References:
Minsley, B.J., James, S.R., Bedrosian, P.A., Pace, M.D., Hoogenboom, B.E., and Burton, B.L., 2021, Airborne electromagnetic, magnetic, and radiometric survey of the Mississippi Alluvial Plain, November 2019 - March 2020: U.S. Geological Survey data release, https://doi.org/10.5066/P9E44CTQ.
James, S.R., and Minsley, B.J., 2021, Combined results and derivative products of hydrogeologic structure and properties from airborne electromagnetic surveys in the Mississippi Alluvial Plain: U.S. Geological Survey data release, https://doi.org/10.5066/P9382RCI.

.. GENERATED FROM PYTHON SOURCE LINES 10-14

.. code-block:: default

    import matplotlib.pyplot as plt
    from os.path import join
    from gspy import Survey








.. GENERATED FROM PYTHON SOURCE LINES 15-17

Convert the TIFs data to netcdf
+++++++++++++++++++++++++++++++

.. GENERATED FROM PYTHON SOURCE LINES 17-45

.. code-block:: default


    # Path to example files
    data_path = '..//..//supplemental//'

    # Define ..//supplemental information file
    ..//supplemental = data_path + "region//MAP//data//Tempest_survey_information.json"

    # Add ..//supplemental information to the survey
    survey = Survey(..//supplemental)

    # Define input ASEG-format data file and associated variable mapping file
    d_data = data_path + 'region//MAP//data//Tempest.dat'
    d_supp = data_path + 'region//MAP//data//Tempest_data_information.json'

    # Read data and format as tabular class object
    # survey.add_tabular(type='aseg', data_filename=d_data, metadata_file=d_supp)

    # Define input TIF-format data file and associated variable mapping file
    d_grid_path = data_path + 'region//MAP//data//'
    d_grid_supp = data_path + 'region//MAP//data//Tempest_rasters_Attributes.json'

    # Read data and format as Griddata class object
    survey.add_raster(metadata_file=d_grid_supp)

    # Write NetCDF
    d_out = data_path + 'region//MAP//data//tifs.nc'
    survey.write_netcdf(d_out)








.. GENERATED FROM PYTHON SOURCE LINES 46-47

Read in the netcdf files

.. GENERATED FROM PYTHON SOURCE LINES 47-49

.. code-block:: default

    new_survey = Survey.read_netcdf(d_out)








.. GENERATED FROM PYTHON SOURCE LINES 50-51

Plotting

.. GENERATED FROM PYTHON SOURCE LINES 51-53

.. code-block:: default

    plt.figure()
    new_survey.raster.pcolor('resistivity', stack=0, vmin=0, vmax=3, cmap='jet')
    plt.show()


.. image-sg:: /examples/Data_conversion/images/sphx_glr_plot_tifs_to_netcdf_001.png
   :alt: stack = 0 [not_defined], spatial_ref = 0
   :srcset: /examples/Data_conversion/images/sphx_glr_plot_tifs_to_netcdf_001.png
   :class: sphx-glr-single-img






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

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


.. _sphx_glr_download_examples_Data_conversion_plot_tifs_to_netcdf.py:


.. only :: html

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



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

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



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

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


.. only:: html

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

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