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

.. only:: html

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

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

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

.. _sphx_glr_auto_examples_image_plot_single_rp.py:


======================
Single recurrence plot
======================

A recurrence plot is an image obtained from a time series, representing the
pairwise Euclidean distances for each value (and more generally for each
trajectory) in the time series.
The image can be binarized using a threshold.
It is implemented as :class:`pyts.image.RecurrencePlot`.

In this example, the considered time series is the sequence of the sine
function values for 1000 equally-spaced points in the interval
:math:`[0, 4\pi]` and the threshold used is 0.1.
One can see on the recurrence plot that the sine function is periodic with
period :math:`2\pi` and that its derivative take lower absolute values around
:math:`\frac{\pi}{2} + k\pi` for any integer :math:`k` (because of the higher
density of black pixels).

Since the API is designed for machine learning, the
:meth:`~pyts.image.RecurrencePlot.transform` method of the
:class:`pyts.image.RecurrencePlot` class expects a data set of time series
as input, so the time series is transformed into a data set with a single time
series (``X = np.array([x])``) and the first element of the data set of
recurrence plots is retrieved (``ax_rp.imshow(X_rp[0], ...``).

.. GENERATED FROM PYTHON SOURCE LINES 27-86



.. image-sg:: /auto_examples/image/images/sphx_glr_plot_single_rp_001.png
   :alt: plot single rp
   :srcset: /auto_examples/image/images/sphx_glr_plot_single_rp_001.png
   :class: sphx-glr-single-img





.. code-block:: default


    # Author: Johann Faouzi <johann.faouzi@gmail.com>
    # License: BSD-3-Clause

    import numpy as np
    import matplotlib.pyplot as plt
    from pyts.image import RecurrencePlot


    # Create a toy time series using the sine function
    time_points = np.linspace(0, 4 * np.pi, 1000)
    x = np.sin(time_points)
    X = np.array([x])

    # Recurrence plot transformation
    rp = RecurrencePlot(threshold=np.pi/18)
    X_rp = rp.transform(X)

    # Plot the time series and its recurrence plot
    fig = plt.figure(figsize=(6, 6))

    gs = fig.add_gridspec(2, 2,  width_ratios=(2, 7), height_ratios=(2, 7),
                          left=0.1, right=0.9, bottom=0.1, top=0.9,
                          wspace=0.05, hspace=0.05)

    # Define the ticks and their labels for both axes
    time_ticks = np.linspace(0, 4 * np.pi, 9)
    time_ticklabels = [r'$0$', r'$\frac{\pi}{2}$', r'$\pi$',
                       r'$\frac{3\pi}{2}$', r'$2\pi$', r'$\frac{5\pi}{2}$',
                       r'$3\pi$', r'$\frac{7\pi}{2}$', r'$4\pi$']
    value_ticks = [-1, 0, 1]
    reversed_value_ticks = value_ticks[::-1]

    # Plot the time series on the left with inverted axes
    ax_left = fig.add_subplot(gs[1, 0])
    ax_left.plot(x, time_points)
    ax_left.set_xticks(reversed_value_ticks)
    ax_left.set_xticklabels(reversed_value_ticks, rotation=90)
    ax_left.set_yticks(time_ticks)
    ax_left.set_yticklabels(time_ticklabels, rotation=90)
    ax_left.invert_xaxis()

    # Plot the time series on the top
    ax_top = fig.add_subplot(gs[0, 1])
    ax_top.plot(time_points, x)
    ax_top.set_xticks(time_ticks)
    ax_top.set_xticklabels(time_ticklabels)
    ax_top.set_yticks(value_ticks)
    ax_top.set_yticklabels(value_ticks)
    ax_top.xaxis.tick_top()

    # Plot the recurrence plot on the bottom right
    ax_rp = fig.add_subplot(gs[1, 1])
    ax_rp.imshow(X_rp[0], cmap='binary', origin='lower',
                 extent=[0, 4 * np.pi, 0, 4 * np.pi])
    ax_rp.set_xticks([])
    ax_rp.set_yticks([])

    plt.show()


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

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


.. _sphx_glr_download_auto_examples_image_plot_single_rp.py:

.. only:: html

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


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

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

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

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


.. only:: html

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

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