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

.. only:: html

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

        :ref:`Go to the end <sphx_glr_download__auto_examples_multisensors_cube.py>`
        to download the full example code.

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

.. _sphx_glr__auto_examples_multisensors_cube.py:


Create a multisensor cube
=================================================================

With Sentinel-2 and Landsat, using Sentinel-2 spatial resolution.

.. GENERATED FROM PYTHON SOURCE LINES 9-11

Import librairies
-------------------------------------------

.. GENERATED FROM PYTHON SOURCE LINES 11-16

.. code-block:: Python


    from matplotlib import pyplot as plt

    from earthdaily import EarthDataStore, datasets








.. GENERATED FROM PYTHON SOURCE LINES 17-19

Set parameters
-------------------------------------------

.. GENERATED FROM PYTHON SOURCE LINES 21-23

Init earthdatastore with environment variables or default credentials
----------------------------------------------------------------------------

.. GENERATED FROM PYTHON SOURCE LINES 23-35

.. code-block:: Python


    eds = EarthDataStore()

    collections = ["sentinel-2-l2a", "landsat-c2l2-sr"]
    datetime = ["2022-07-20", "2022-09-01"]
    intersects = datasets.load_pivot_corumba()
    assets = ["blue", "green", "red", "nir"]
    mask_with = "native"
    clear_cover = 50
    resampling = "cubic"
    cross_calibration_collection = "sentinel-2-l2a"



.. rst-class:: sphx-glr-script-out

.. code-block:: pytb

    Traceback (most recent call last):
      File "/private/var/folders/_t/czv_cfrd613_6_gd_p8qpqd80000gp/T/tmpmap_d3kq/7261d0bf6bc6b22a0278281791f6623d19ddf89a/examples/multisensors_cube.py", line 24, in <module>
        eds = EarthDataStore()
              ^^^^^^^^^^^^^^^^
      File "/private/var/folders/_t/czv_cfrd613_6_gd_p8qpqd80000gp/T/tmpmap_d3kq/7261d0bf6bc6b22a0278281791f6623d19ddf89a/earthdaily/__init__.py", line 51, in EarthDataStore
        return earthdatastore.Auth.from_credentials(
               ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
      File "/private/var/folders/_t/czv_cfrd613_6_gd_p8qpqd80000gp/T/tmpmap_d3kq/7261d0bf6bc6b22a0278281791f6623d19ddf89a/earthdaily/earthdatastore/__init__.py", line 464, in from_credentials
        config = cls.read_credentials(
                 ^^^^^^^^^^^^^^^^^^^^^
      File "/private/var/folders/_t/czv_cfrd613_6_gd_p8qpqd80000gp/T/tmpmap_d3kq/7261d0bf6bc6b22a0278281791f6623d19ddf89a/earthdaily/earthdatastore/__init__.py", line 529, in read_credentials
        raise NotImplementedError("Credentials weren't found.")
    NotImplementedError: Credentials weren't found.




.. GENERATED FROM PYTHON SOURCE LINES 36-38

Create the multisensors datacube
-------------------------------------------

.. GENERATED FROM PYTHON SOURCE LINES 38-55

.. code-block:: Python


    datacube = eds.datacube(
        collections,
        assets=assets,
        datetime=datetime,
        intersects=intersects,
        mask_with=mask_with,
        clear_cover=clear_cover,
        cross_calibration_collection=cross_calibration_collection,
    )

    # Add the NDVI
    datacube = datacube.ed.add_indices(["NDVI"])

    # Load in memory
    datacube = datacube.load()


.. GENERATED FROM PYTHON SOURCE LINES 56-58

See the evolution in RGB
-------------------------------------------

.. GENERATED FROM PYTHON SOURCE LINES 58-62

.. code-block:: Python


    datacube.ed.plot_rgb(col_wrap=3)
    plt.show()


.. GENERATED FROM PYTHON SOURCE LINES 63-65

See the NDVI evolution
-------------------------------------------

.. GENERATED FROM PYTHON SOURCE LINES 65-69

.. code-block:: Python


    datacube["NDVI"].ed.plot_band(col_wrap=3, vmin=0, vmax=0.8, cmap="Greens")
    plt.show()


.. GENERATED FROM PYTHON SOURCE LINES 70-72

See the NDVI mean evolution
-------------------------------------------

.. GENERATED FROM PYTHON SOURCE LINES 72-76

.. code-block:: Python


    datacube["NDVI"].groupby("time").mean(...).plot.line(x="time")
    plt.title("NDVI evolution")
    plt.show()


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

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


.. _sphx_glr_download__auto_examples_multisensors_cube.py:

.. only:: html

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

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

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

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

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

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

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


.. only:: html

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

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