.. 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_birth_and_death.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_birth_and_death.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_birth_and_death.py:


Birth and death
===============

Following figures are based on https://doi.org/10.1016/j.pocean.2011.01.002

.. GENERATED FROM PYTHON SOURCE LINES 8-14

.. code-block:: Python

    from matplotlib import pyplot as plt
    import py_eddy_tracker_sample

    from py_eddy_tracker.observations.tracking import TrackEddiesObservations









.. GENERATED FROM PYTHON SOURCE LINES 15-30

.. code-block:: Python

    def start_axes(title):
        fig = plt.figure(figsize=(13, 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


    def update_axes(ax, mappable=None):
        ax.grid()
        if mappable:
            plt.colorbar(mappable, cax=ax.figure.add_axes([0.95, 0.05, 0.01, 0.9]))









.. GENERATED FROM PYTHON SOURCE LINES 31-32

Load an experimental med atlas over a period of 26 years (1993-2019)

.. GENERATED FROM PYTHON SOURCE LINES 32-52

.. code-block:: Python

    kwargs_load = dict(
        include_vars=(
            "longitude",
            "latitude",
            "observation_number",
            "track",
            "time",
            "speed_contour_longitude",
            "speed_contour_latitude",
        )
    )
    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")
    )








.. GENERATED FROM PYTHON SOURCE LINES 53-58

.. code-block:: Python

    t0, t1 = a.period
    step = 0.125
    bins = ((-10, 37, step), (30, 46, step))
    kwargs = dict(cmap="terrain_r", factor=100 / (t1 - t0), name="count", vmin=0, vmax=1)








.. GENERATED FROM PYTHON SOURCE LINES 59-61

Cyclonic
--------

.. GENERATED FROM PYTHON SOURCE LINES 61-66

.. code-block:: Python

    ax = start_axes("Birth cyclonic frenquency (%)")
    g_c_first = c.first_obs().grid_count(bins, intern=True)
    m = g_c_first.display(ax, **kwargs)
    update_axes(ax, m)




.. image-sg:: /python_module/10_tracking_diagnostics/images/sphx_glr_pet_birth_and_death_001.png
   :alt: Birth cyclonic frenquency (%)
   :srcset: /python_module/10_tracking_diagnostics/images/sphx_glr_pet_birth_and_death_001.png
   :class: sphx-glr-single-img





.. GENERATED FROM PYTHON SOURCE LINES 67-72

.. code-block:: Python

    ax = start_axes("Death cyclonic frenquency (%)")
    g_c_last = c.last_obs().grid_count(bins, intern=True)
    m = g_c_last.display(ax, **kwargs)
    update_axes(ax, m)




.. image-sg:: /python_module/10_tracking_diagnostics/images/sphx_glr_pet_birth_and_death_002.png
   :alt: Death cyclonic frenquency (%)
   :srcset: /python_module/10_tracking_diagnostics/images/sphx_glr_pet_birth_and_death_002.png
   :class: sphx-glr-single-img





.. GENERATED FROM PYTHON SOURCE LINES 73-75

Anticyclonic
------------

.. GENERATED FROM PYTHON SOURCE LINES 75-80

.. code-block:: Python

    ax = start_axes("Birth anticyclonic frequency (%)")
    g_a_first = a.first_obs().grid_count(bins, intern=True)
    m = g_a_first.display(ax, **kwargs)
    update_axes(ax, m)




.. image-sg:: /python_module/10_tracking_diagnostics/images/sphx_glr_pet_birth_and_death_003.png
   :alt: Birth anticyclonic frequency (%)
   :srcset: /python_module/10_tracking_diagnostics/images/sphx_glr_pet_birth_and_death_003.png
   :class: sphx-glr-single-img





.. GENERATED FROM PYTHON SOURCE LINES 81-85

.. code-block:: Python

    ax = start_axes("Death anticyclonic frequency (%)")
    g_a_last = a.last_obs().grid_count(bins, intern=True)
    m = g_a_last.display(ax, **kwargs)
    update_axes(ax, m)



.. image-sg:: /python_module/10_tracking_diagnostics/images/sphx_glr_pet_birth_and_death_004.png
   :alt: Death anticyclonic frequency (%)
   :srcset: /python_module/10_tracking_diagnostics/images/sphx_glr_pet_birth_and_death_004.png
   :class: sphx-glr-single-img






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

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


.. _sphx_glr_download_python_module_10_tracking_diagnostics_pet_birth_and_death.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_birth_and_death.ipynb
        :alt: Launch binder
        :width: 150 px

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

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

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

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

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

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


.. only:: html

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

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