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

.. only:: html

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

        :ref:`Go to the end <sphx_glr_download_python_module_10_tracking_diagnostics_pet_geographic_stats.py>`
        to download the full example code. or to run this example in your browser via Binder

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

.. _sphx_glr_python_module_10_tracking_diagnostics_pet_geographic_stats.py:


Geographical statistics
=======================

.. GENERATED FROM PYTHON SOURCE LINES 6-22

.. code-block:: Python


    from matplotlib import pyplot as plt
    import py_eddy_tracker_sample

    from py_eddy_tracker.observations.tracking import TrackEddiesObservations


    def start_axes(title):
        fig = plt.figure(figsize=(13.5, 5))
        ax = fig.add_axes([0.03, 0.03, 0.90, 0.94])
        ax.set_xlim(-6, 36.5), ax.set_ylim(30, 46)
        ax.set_aspect("equal")
        ax.set_title(title)
        return ax









.. GENERATED FROM PYTHON SOURCE LINES 23-24

Load an experimental med atlas over a period of 26 years (1993-2019), we merge the 2 datasets

.. GENERATED FROM PYTHON SOURCE LINES 24-36

.. code-block:: Python

    a = TrackEddiesObservations.load_file(
        py_eddy_tracker_sample.get_demo_path(
            "eddies_med_adt_allsat_dt2018/Anticyclonic.zarr"
        )
    )
    c = TrackEddiesObservations.load_file(
        py_eddy_tracker_sample.get_demo_path("eddies_med_adt_allsat_dt2018/Cyclonic.zarr")
    )
    a = a.merge(c)

    step = 0.1








.. GENERATED FROM PYTHON SOURCE LINES 37-38

Mean of amplitude in each box

.. GENERATED FROM PYTHON SOURCE LINES 38-45

.. code-block:: Python

    ax = start_axes("Amplitude mean by box of %s°" % step)
    g = a.grid_stat(((-7, 37, step), (30, 46, step)), "amplitude")
    m = g.display(ax, name="amplitude", vmin=0, vmax=10, factor=100)
    ax.grid()
    cb = plt.colorbar(m, cax=ax.figure.add_axes([0.94, 0.05, 0.01, 0.9]))
    cb.set_label("Amplitude (cm)")




.. image-sg:: /python_module/10_tracking_diagnostics/images/sphx_glr_pet_geographic_stats_001.png
   :alt: Amplitude mean by box of 0.1°
   :srcset: /python_module/10_tracking_diagnostics/images/sphx_glr_pet_geographic_stats_001.png
   :class: sphx-glr-single-img





.. GENERATED FROM PYTHON SOURCE LINES 46-47

Mean of speed radius in each box

.. GENERATED FROM PYTHON SOURCE LINES 47-54

.. code-block:: Python

    ax = start_axes("Speed radius mean by box of %s°" % step)
    g = a.grid_stat(((-7, 37, step), (30, 46, step)), "radius_s")
    m = g.display(ax, name="radius_s", vmin=10, vmax=50, factor=0.001)
    ax.grid()
    cb = plt.colorbar(m, cax=ax.figure.add_axes([0.94, 0.05, 0.01, 0.9]))
    cb.set_label("Speed radius (km)")




.. image-sg:: /python_module/10_tracking_diagnostics/images/sphx_glr_pet_geographic_stats_002.png
   :alt: Speed radius mean by box of 0.1°
   :srcset: /python_module/10_tracking_diagnostics/images/sphx_glr_pet_geographic_stats_002.png
   :class: sphx-glr-single-img





.. GENERATED FROM PYTHON SOURCE LINES 55-56

Percent of virtual on the whole obs in each box

.. GENERATED FROM PYTHON SOURCE LINES 56-63

.. code-block:: Python

    ax = start_axes("Percent of virtual by box of %s°" % step)
    g = a.grid_stat(((-7, 37, step), (30, 46, step)), "virtual")
    g.vars["virtual"] *= 100
    m = g.display(ax, name="virtual", vmin=0, vmax=15)
    ax.grid()
    cb = plt.colorbar(m, cax=ax.figure.add_axes([0.94, 0.05, 0.01, 0.9]))
    cb.set_label("Percent of virtual (%)")



.. image-sg:: /python_module/10_tracking_diagnostics/images/sphx_glr_pet_geographic_stats_003.png
   :alt: Percent of virtual by box of 0.1°
   :srcset: /python_module/10_tracking_diagnostics/images/sphx_glr_pet_geographic_stats_003.png
   :class: sphx-glr-single-img






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

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


.. _sphx_glr_download_python_module_10_tracking_diagnostics_pet_geographic_stats.py:

.. only:: html

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

    .. container:: binder-badge

      .. image:: images/binder_badge_logo.svg
        :target: https://mybinder.org/v2/gh/AntSimi/py-eddy-tracker/master?urlpath=lab/tree/notebooks/python_module/10_tracking_diagnostics/pet_geographic_stats.ipynb
        :alt: Launch binder
        :width: 150 px

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

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

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

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

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

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


.. only:: html

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

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