Examples/01_ADC_Feature_Space_Generation.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "fa3f0820",
   "metadata": {},
   "source": [
    "<center><h2>Feature Space Generation</h2></center>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "72689598",
   "metadata": {},
   "source": [
    "<img src='https://www.projectrhea.org/rhea/images/thumb/2/22/Hyperplane.png/700px-Hyperplane.png'/>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "18d43282",
   "metadata": {},
   "source": [
    "This notebook showcases how we can utilize the ADC in order to create a feature space based on parameters such as time range, bands and buffer zone for each parcel. In addition, a daily interpolation has been applied to the data. "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3a138fa0",
   "metadata": {},
   "source": [
    "### Import Libraries "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "e038cce0",
   "metadata": {},
   "outputs": [],
   "source": [
    "import sys,logging\n",
    "import xarray as xr\n",
    "import pandas as pd\n",
    "from pathlib import Path\n",
    "import h5py\n",
    "sys.path.append(str(Path('.').absolute().parent))\n",
    "sys.path.insert(0,str(Path('.').absolute().parent))\n",
    "from adc_utils import * # Please change the path for the ADC configuration (if you use the Cube)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "50e89092",
   "metadata": {},
   "source": [
    "### Load Sentinel Data\n",
    "Open Data Cube (on which ADC has been based) can be installed via the instructions given in <a href=\"https://github.com/Agri-Hub/datacap\"> DataCAP</a> repository. Nevertheless, we have uploaded <b>netCDF</b> files, extracted from the ADC, so users can easily load the required data and run these notebooks as a proof of concept. These data can be found here **https://zenodo.org/record/6458275#.YlcRz1vP2E1**\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "5cce1473",
   "metadata": {},
   "outputs": [],
   "source": [
    "import xarray as xr\n",
    "ws_datasets = '../datasets/Cyprus'\n",
    "netcdf = []\n",
    "for filename in sorted(os.listdir(ws_datasets)):\n",
    "    if not 'optical' in filename:\n",
    "        continue\n",
    "    netcdf.append(xr.open_dataset(os.path.join(ws_datasets,filename)))\n",
    "data_nc = xr.merge(netcdf)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b7ef5193",
   "metadata": {},
   "source": [
    "### Overview of Optical Data loaded from netCDF"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "8c191400",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block !important;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2 {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.Dataset&gt;\n",
       "Dimensions:      (time: 224, y: 1354, x: 1609)\n",
       "Coordinates:\n",
       "  * time         (time) datetime64[ns] 2018-10-01 ... 2020-09-30T08:27:39\n",
       "  * y            (y) float64 4.156e+06 4.156e+06 ... 4.142e+06 4.142e+06\n",
       "  * x            (x) float64 3.611e+06 3.611e+06 ... 3.627e+06 3.627e+06\n",
       "    spatial_ref  int32 3857\n",
       "Data variables:\n",
       "    B11          (time, y, x) uint16 ...\n",
       "    B12          (time, y, x) uint16 ...\n",
       "    B02          (time, y, x) uint16 ...\n",
       "    B03          (time, y, x) uint16 ...\n",
       "    B04          (time, y, x) uint16 ...\n",
       "    B05          (time, y, x) uint16 ...\n",
       "    B06          (time, y, x) uint16 ...\n",
       "    B07          (time, y, x) uint16 ...\n",
       "    B08          (time, y, x) uint16 ...\n",
       "    B8A          (time, y, x) uint16 ...\n",
       "    SCL          (time, y, x) uint8 ...\n",
       "Attributes:\n",
       "    crs:           EPSG:3857\n",
       "    grid_mapping:  spatial_ref</pre><div class='xr-wrap' style='display:none'><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-6df59d33-e484-42e0-a2c5-8c727c447b8d' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-6df59d33-e484-42e0-a2c5-8c727c447b8d' class='xr-section-summary'  title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'><ul class='xr-dim-list'><li><span class='xr-has-index'>time</span>: 224</li><li><span class='xr-has-index'>y</span>: 1354</li><li><span class='xr-has-index'>x</span>: 1609</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-b6f3d5e1-805a-4182-8659-a6a9991dbcb1' class='xr-section-summary-in' type='checkbox'  checked><label for='section-b6f3d5e1-805a-4182-8659-a6a9991dbcb1' class='xr-section-summary' >Coordinates: <span>(4)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>time</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>datetime64[ns]</div><div class='xr-var-preview xr-preview'>2018-10-01 ... 2020-09-30T08:27:39</div><input id='attrs-2f304504-d5d3-464b-9a19-5bea3adcdda8' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-2f304504-d5d3-464b-9a19-5bea3adcdda8' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-02abf189-8597-490b-b915-57f144bff765' class='xr-var-data-in' type='checkbox'><label for='data-02abf189-8597-490b-b915-57f144bff765' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([&#x27;2018-10-01T00:00:00.000000000&#x27;, &#x27;2018-10-04T00:00:00.000000000&#x27;,\n",
       "       &#x27;2018-10-06T00:00:00.000000000&#x27;, ..., &#x27;2020-09-20T08:26:29.000000000&#x27;,\n",
       "       &#x27;2020-09-25T08:27:11.000000000&#x27;, &#x27;2020-09-30T08:27:39.000000000&#x27;],\n",
       "      dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>y</span></div><div class='xr-var-dims'>(y)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>4.156e+06 4.156e+06 ... 4.142e+06</div><input id='attrs-cd696a08-101c-445c-80b8-f93e275ceb30' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-cd696a08-101c-445c-80b8-f93e275ceb30' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-2afd2ef3-2fb0-4f86-9d8b-8dd6cef5764d' class='xr-var-data-in' type='checkbox'><label for='data-2afd2ef3-2fb0-4f86-9d8b-8dd6cef5764d' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>metre</dd><dt><span>resolution :</span></dt><dd>-10.0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd></dl></div><div class='xr-var-data'><pre>array([4155685., 4155675., 4155665., ..., 4142175., 4142165., 4142155.])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>x</span></div><div class='xr-var-dims'>(x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>3.611e+06 3.611e+06 ... 3.627e+06</div><input id='attrs-0ffa8abf-2d8c-4e5e-8e7a-8d52b3256d3f' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-0ffa8abf-2d8c-4e5e-8e7a-8d52b3256d3f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-06a81c09-7ef2-41ed-aa5a-55eaf113ab0c' class='xr-var-data-in' type='checkbox'><label for='data-06a81c09-7ef2-41ed-aa5a-55eaf113ab0c' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>metre</dd><dt><span>resolution :</span></dt><dd>10.0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd></dl></div><div class='xr-var-data'><pre>array([3610515., 3610525., 3610535., ..., 3626575., 3626585., 3626595.])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>spatial_ref</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>3857</div><input id='attrs-0edc687d-82aa-469e-ad96-b265e71f3695' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-0edc687d-82aa-469e-ad96-b265e71f3695' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3f69852b-b212-4429-8a90-cf9ea13ad4c4' class='xr-var-data-in' type='checkbox'><label for='data-3f69852b-b212-4429-8a90-cf9ea13ad4c4' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>spatial_ref :</span></dt><dd>PROJCS[&quot;WGS 84 / Pseudo-Mercator&quot;,GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],PRIMEM[&quot;Greenwich&quot;,0,AUTHORITY[&quot;EPSG&quot;,&quot;8901&quot;]],UNIT[&quot;degree&quot;,0.0174532925199433,AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]],PROJECTION[&quot;Mercator_1SP&quot;],PARAMETER[&quot;central_meridian&quot;,0],PARAMETER[&quot;scale_factor&quot;,1],PARAMETER[&quot;false_easting&quot;,0],PARAMETER[&quot;false_northing&quot;,0],UNIT[&quot;metre&quot;,1,AUTHORITY[&quot;EPSG&quot;,&quot;9001&quot;]],AXIS[&quot;Easting&quot;,EAST],AXIS[&quot;Northing&quot;,NORTH],EXTENSION[&quot;PROJ4&quot;,&quot;+proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext +no_defs&quot;],AUTHORITY[&quot;EPSG&quot;,&quot;3857&quot;]]</dd><dt><span>grid_mapping_name :</span></dt><dd>??</dd></dl></div><div class='xr-var-data'><pre>array(3857, dtype=int32)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-4e45d63f-6732-4cf9-a29d-c69e9e38cea4' class='xr-section-summary-in' type='checkbox'  checked><label for='section-4e45d63f-6732-4cf9-a29d-c69e9e38cea4' class='xr-section-summary' >Data variables: <span>(11)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span>B11</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>uint16</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-7bbe1566-5576-4273-a5c6-537141d46f1a' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-7bbe1566-5576-4273-a5c6-537141d46f1a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-18edb00e-6ee9-407f-a88f-5073e92d06b8' class='xr-var-data-in' type='checkbox'><label for='data-18edb00e-6ee9-407f-a88f-5073e92d06b8' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>[488003264 values with dtype=uint16]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B12</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>uint16</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-4b4697a6-fa9d-484b-81de-ab30470ebb95' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-4b4697a6-fa9d-484b-81de-ab30470ebb95' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-06b26480-4924-42bc-bdbc-226120439b80' class='xr-var-data-in' type='checkbox'><label for='data-06b26480-4924-42bc-bdbc-226120439b80' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>[488003264 values with dtype=uint16]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B02</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>uint16</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-6c105ddf-715b-4212-9b2b-920589cb6edf' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-6c105ddf-715b-4212-9b2b-920589cb6edf' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4d5a65ba-db35-4005-8400-b3caf7504be8' class='xr-var-data-in' type='checkbox'><label for='data-4d5a65ba-db35-4005-8400-b3caf7504be8' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>[488003264 values with dtype=uint16]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B03</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>uint16</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-dec8278b-53e8-433b-ac55-197b2035e9bf' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-dec8278b-53e8-433b-ac55-197b2035e9bf' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-bf16dcdb-bab0-4c1b-b786-eff803865c4f' class='xr-var-data-in' type='checkbox'><label for='data-bf16dcdb-bab0-4c1b-b786-eff803865c4f' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>[488003264 values with dtype=uint16]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B04</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>uint16</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-f2ce3895-a2ac-4a00-8721-3e1d5e72fad1' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-f2ce3895-a2ac-4a00-8721-3e1d5e72fad1' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-0c912690-fd19-4c24-b6c4-a82ba0c4915a' class='xr-var-data-in' type='checkbox'><label for='data-0c912690-fd19-4c24-b6c4-a82ba0c4915a' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>[488003264 values with dtype=uint16]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B05</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>uint16</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-21fd0092-4abe-40fb-8d60-7f7650e70ad0' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-21fd0092-4abe-40fb-8d60-7f7650e70ad0' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4b8d9013-26f7-4885-89bb-eed8ea02d3c4' class='xr-var-data-in' type='checkbox'><label for='data-4b8d9013-26f7-4885-89bb-eed8ea02d3c4' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>[488003264 values with dtype=uint16]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B06</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>uint16</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-40358cd3-aa21-4188-b4d4-69aab2c81c76' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-40358cd3-aa21-4188-b4d4-69aab2c81c76' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5c651d1c-84c1-45ea-acd8-a74e9661322a' class='xr-var-data-in' type='checkbox'><label for='data-5c651d1c-84c1-45ea-acd8-a74e9661322a' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>[488003264 values with dtype=uint16]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B07</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>uint16</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-501becbc-70dd-41ff-a36d-2dffc0ea6732' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-501becbc-70dd-41ff-a36d-2dffc0ea6732' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-cf34eee1-3c55-4544-98d1-fe4835be4a7e' class='xr-var-data-in' type='checkbox'><label for='data-cf34eee1-3c55-4544-98d1-fe4835be4a7e' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>[488003264 values with dtype=uint16]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B08</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>uint16</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-5ec9e6ce-f306-4499-ac47-f1175704e0ec' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-5ec9e6ce-f306-4499-ac47-f1175704e0ec' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-eb109c35-0d44-4cd1-b264-3612e4d356f9' class='xr-var-data-in' type='checkbox'><label for='data-eb109c35-0d44-4cd1-b264-3612e4d356f9' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>[488003264 values with dtype=uint16]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B8A</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>uint16</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-034124c4-8672-4202-9d14-9f6d6bf2a551' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-034124c4-8672-4202-9d14-9f6d6bf2a551' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3a9d522c-2294-46d5-9c58-d944c5d16655' class='xr-var-data-in' type='checkbox'><label for='data-3a9d522c-2294-46d5-9c58-d944c5d16655' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>[488003264 values with dtype=uint16]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>SCL</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>uint8</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-7dced800-efa3-44a0-ae1e-3da6d7e34f22' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-7dced800-efa3-44a0-ae1e-3da6d7e34f22' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-2afce14e-c7f4-40d8-9fd5-c2be8cf802c8' class='xr-var-data-in' type='checkbox'><label for='data-2afce14e-c7f4-40d8-9fd5-c2be8cf802c8' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>[488003264 values with dtype=uint8]</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-f9649567-55fe-4378-ad87-86f7767d452b' class='xr-section-summary-in' type='checkbox'  checked><label for='section-f9649567-55fe-4378-ad87-86f7767d452b' class='xr-section-summary' >Attributes: <span>(2)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset>\n",
       "Dimensions:      (time: 224, y: 1354, x: 1609)\n",
       "Coordinates:\n",
       "  * time         (time) datetime64[ns] 2018-10-01 ... 2020-09-30T08:27:39\n",
       "  * y            (y) float64 4.156e+06 4.156e+06 ... 4.142e+06 4.142e+06\n",
       "  * x            (x) float64 3.611e+06 3.611e+06 ... 3.627e+06 3.627e+06\n",
       "    spatial_ref  int32 3857\n",
       "Data variables:\n",
       "    B11          (time, y, x) uint16 ...\n",
       "    B12          (time, y, x) uint16 ...\n",
       "    B02          (time, y, x) uint16 ...\n",
       "    B03          (time, y, x) uint16 ...\n",
       "    B04          (time, y, x) uint16 ...\n",
       "    B05          (time, y, x) uint16 ...\n",
       "    B06          (time, y, x) uint16 ...\n",
       "    B07          (time, y, x) uint16 ...\n",
       "    B08          (time, y, x) uint16 ...\n",
       "    B8A          (time, y, x) uint16 ...\n",
       "    SCL          (time, y, x) uint8 ...\n",
       "Attributes:\n",
       "    crs:           EPSG:3857\n",
       "    grid_mapping:  spatial_ref"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data_nc"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f07bb23b",
   "metadata": {},
   "source": [
    "### Loading Sentinel Data (if ADC is installed and configured): Setting the required parameters\n",
    "The following parameters are used in order to query the ADC. As the data have been already exported in netCDF format, there is no need for using them. "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "51c22ca3",
   "metadata": {},
   "source": [
    "<code>\n",
    "year = '2019'\n",
    "pilot = 'cyprus'\n",
    "region = 'new_area_1'\n",
    "buffer = 'buffer5'\n",
    "shapefile_bbox = 'full_path_of_desired_shapefile'\n",
    "driver = ogr.GetDriverByName(\"ESRI Shapefile\")\n",
    "ds = driver.Open(shapefile_bbox,0)\n",
    "layer = ds.GetLayer()\n",
    "xmin,xmax,ymin,ymax = layer.GetExtent()\n",
    "bbox = [xmin,xmax,ymin,ymax]\n",
    "product_ids = 'ids_{}_{}'.format(pilot,year)\n",
    "ids = getIDs(product_ids,xmin,xmax,ymin,ymax)\n",
    "</code>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f85d2714",
   "metadata": {},
   "source": [
    "###  Set the parameters "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5dbbde33",
   "metadata": {},
   "source": [
    "Except from the parameters above, both netCDF files and direct ADC query demand the following parameters:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "7761ac58",
   "metadata": {},
   "outputs": [],
   "source": [
    "year = '2019'\n",
    "\n",
    "pilot = 'cyprus'\n",
    "\n",
    "buffer = 'buffer5'\n",
    "\n",
    "\n",
    "optical_bands = ['B02','B03','B04','B05','B06','B07','B08','B8A','B11','B12','ndvi','ndwi','psri']\n",
    "\n",
    "sar_bands = ['vv','vh']\n",
    "\n",
    "# two month, as a proof of concept. Change it!\n",
    "start_date,end_date = '2019-01-01','2019-03-01' \n",
    "\n",
    "# number of breaks between the dates interval. Data loading takes place on each break reducing the usage of RAM\n",
    "d_break = 8 "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "940eec7e",
   "metadata": {},
   "source": [
    "### Load IDs of declarations and correspondive crop types\n",
    "\n",
    "Each declared parcel in LPIS comes with a unique id and it's crop type. We have converted unique ids to new ids named noa_id"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "e8e502ef",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Unnamed: 0</th>\n",
       "      <th>Decl_Code</th>\n",
       "      <th>noa_id</th>\n",
       "      <th>Description</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>Durum Wheat</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>2</td>\n",
       "      <td>Durum Wheat</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>3</td>\n",
       "      <td>Durum Wheat</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>4</td>\n",
       "      <td>Durum Wheat</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>4</td>\n",
       "      <td>3</td>\n",
       "      <td>5</td>\n",
       "      <td>Barley</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Unnamed: 0  Decl_Code  noa_id  Description\n",
       "0           0          1       1  Durum Wheat\n",
       "1           1          1       2  Durum Wheat\n",
       "2           2          1       3  Durum Wheat\n",
       "3           3          1       4  Durum Wheat\n",
       "4           4          3       5       Barley"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dbf = pd.read_csv(os.path.join(ws_datasets,'cyprus_ids_cropTypes_2019.csv'))\n",
    "dbf.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "eea88069",
   "metadata": {},
   "outputs": [],
   "source": [
    "noa_to_id = {k:v for k,v in dbf[['noa_id','Decl_Code']].values}"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "16d03aaf",
   "metadata": {},
   "source": [
    "### Load the ids for each parcel based on a layer that has been indexed on the cube"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "bb6eef23",
   "metadata": {},
   "source": [
    "**Before the process starts, the LPIS shapefile including farmers' declaration has been rasterized and indexed into the ADC. Thus, it enables fast processing by allowing xarray operations.**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "1391275a",
   "metadata": {},
   "outputs": [],
   "source": [
    "ids = xr.open_dataset(os.path.join(ws_datasets,'ids_2019.nc'))"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ff4384c5",
   "metadata": {},
   "source": [
    "### An overview of the xarray of ids\n",
    "Data variables contain instances of the same LPIS in different buffer zones. By using different buffers we can exlude mixels from data analysis "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "1c452dc0",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block !important;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2 {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.Dataset&gt;\n",
       "Dimensions:      (time: 1, y: 1354, x: 1609)\n",
       "Coordinates:\n",
       "  * time         (time) datetime64[ns] 2019-01-01\n",
       "  * y            (y) float64 4.156e+06 4.156e+06 ... 4.142e+06 4.142e+06\n",
       "  * x            (x) float64 3.611e+06 3.611e+06 ... 3.627e+06 3.627e+06\n",
       "    spatial_ref  int32 3857\n",
       "Data variables:\n",
       "    buffer0      (time, y, x) int32 ...\n",
       "    buffer5      (time, y, x) int32 ...\n",
       "Attributes:\n",
       "    crs:           EPSG:3857\n",
       "    grid_mapping:  spatial_ref</pre><div class='xr-wrap' style='display:none'><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-a13bc7fa-4d97-4062-922a-306dcc1cf28d' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-a13bc7fa-4d97-4062-922a-306dcc1cf28d' class='xr-section-summary'  title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'><ul class='xr-dim-list'><li><span class='xr-has-index'>time</span>: 1</li><li><span class='xr-has-index'>y</span>: 1354</li><li><span class='xr-has-index'>x</span>: 1609</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-a0d178ea-7270-4b7a-8499-57eb38981f7c' class='xr-section-summary-in' type='checkbox'  checked><label for='section-a0d178ea-7270-4b7a-8499-57eb38981f7c' class='xr-section-summary' >Coordinates: <span>(4)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>time</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>datetime64[ns]</div><div class='xr-var-preview xr-preview'>2019-01-01</div><input id='attrs-02ce466e-c5f9-4de1-8d7a-0906644ab186' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-02ce466e-c5f9-4de1-8d7a-0906644ab186' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8b5581c2-ed44-48b3-b127-68e755418f1b' class='xr-var-data-in' type='checkbox'><label for='data-8b5581c2-ed44-48b3-b127-68e755418f1b' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([&#x27;2019-01-01T00:00:00.000000000&#x27;], dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>y</span></div><div class='xr-var-dims'>(y)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>4.156e+06 4.156e+06 ... 4.142e+06</div><input id='attrs-b959b379-6de2-472f-b220-fff5ca7a8d57' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-b959b379-6de2-472f-b220-fff5ca7a8d57' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8dc0b914-b00d-4d34-8afd-d19e9408b792' class='xr-var-data-in' type='checkbox'><label for='data-8dc0b914-b00d-4d34-8afd-d19e9408b792' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>metre</dd><dt><span>resolution :</span></dt><dd>-10.0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd></dl></div><div class='xr-var-data'><pre>array([4155685., 4155675., 4155665., ..., 4142175., 4142165., 4142155.])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>x</span></div><div class='xr-var-dims'>(x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>3.611e+06 3.611e+06 ... 3.627e+06</div><input id='attrs-db49d6dd-57af-445d-b1a7-01c851a2da1c' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-db49d6dd-57af-445d-b1a7-01c851a2da1c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-88a0dfb8-50e5-4a68-aa42-8ca252fd8f1d' class='xr-var-data-in' type='checkbox'><label for='data-88a0dfb8-50e5-4a68-aa42-8ca252fd8f1d' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>metre</dd><dt><span>resolution :</span></dt><dd>10.0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd></dl></div><div class='xr-var-data'><pre>array([3610515., 3610525., 3610535., ..., 3626575., 3626585., 3626595.])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>spatial_ref</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-5f3c7197-99ff-4a64-a017-a3a6779a5898' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-5f3c7197-99ff-4a64-a017-a3a6779a5898' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5d5a018a-43a1-4a95-b2f9-f8992eb96347' class='xr-var-data-in' type='checkbox'><label for='data-5d5a018a-43a1-4a95-b2f9-f8992eb96347' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>spatial_ref :</span></dt><dd>PROJCS[&quot;WGS 84 / Pseudo-Mercator&quot;,GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],PRIMEM[&quot;Greenwich&quot;,0,AUTHORITY[&quot;EPSG&quot;,&quot;8901&quot;]],UNIT[&quot;degree&quot;,0.0174532925199433,AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]],PROJECTION[&quot;Mercator_1SP&quot;],PARAMETER[&quot;central_meridian&quot;,0],PARAMETER[&quot;scale_factor&quot;,1],PARAMETER[&quot;false_easting&quot;,0],PARAMETER[&quot;false_northing&quot;,0],UNIT[&quot;metre&quot;,1,AUTHORITY[&quot;EPSG&quot;,&quot;9001&quot;]],AXIS[&quot;Easting&quot;,EAST],AXIS[&quot;Northing&quot;,NORTH],EXTENSION[&quot;PROJ4&quot;,&quot;+proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext +no_defs&quot;],AUTHORITY[&quot;EPSG&quot;,&quot;3857&quot;]]</dd><dt><span>grid_mapping_name :</span></dt><dd>??</dd></dl></div><div class='xr-var-data'><pre>array(3857, dtype=int32)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-6f67c084-79fb-41e2-9e50-c6ff4b164b3a' class='xr-section-summary-in' type='checkbox'  checked><label for='section-6f67c084-79fb-41e2-9e50-c6ff4b164b3a' class='xr-section-summary' >Data variables: <span>(2)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span>buffer0</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-4241d84e-977b-479a-8d69-7f3cfd9ab5c3' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-4241d84e-977b-479a-8d69-7f3cfd9ab5c3' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-1d879dd2-b921-4e43-ae4c-4c95a2aaf150' class='xr-var-data-in' type='checkbox'><label for='data-1d879dd2-b921-4e43-ae4c-4c95a2aaf150' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>-1</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>[2178586 values with dtype=int32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>buffer5</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-ad9252d3-52d4-4f93-a609-c0a5f972eb4e' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-ad9252d3-52d4-4f93-a609-c0a5f972eb4e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c4a5bfb1-7bd0-4ec5-ac85-b26dce1af899' class='xr-var-data-in' type='checkbox'><label for='data-c4a5bfb1-7bd0-4ec5-ac85-b26dce1af899' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>-1</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>[2178586 values with dtype=int32]</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-803906c7-2667-42dc-bfaa-cbdaf388c5c6' class='xr-section-summary-in' type='checkbox'  checked><label for='section-803906c7-2667-42dc-bfaa-cbdaf388c5c6' class='xr-section-summary' >Attributes: <span>(2)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset>\n",
       "Dimensions:      (time: 1, y: 1354, x: 1609)\n",
       "Coordinates:\n",
       "  * time         (time) datetime64[ns] 2019-01-01\n",
       "  * y            (y) float64 4.156e+06 4.156e+06 ... 4.142e+06 4.142e+06\n",
       "  * x            (x) float64 3.611e+06 3.611e+06 ... 3.627e+06 3.627e+06\n",
       "    spatial_ref  int32 ...\n",
       "Data variables:\n",
       "    buffer0      (time, y, x) int32 ...\n",
       "    buffer5      (time, y, x) int32 ...\n",
       "Attributes:\n",
       "    crs:           EPSG:3857\n",
       "    grid_mapping:  spatial_ref"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ids"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f61beec5",
   "metadata": {},
   "source": [
    "### Filter Data for a specific time period"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "90e291c1",
   "metadata": {},
   "outputs": [],
   "source": [
    "start = pd.to_datetime(start_date).to_numpy()\n",
    "end = pd.to_datetime(end_date).to_numpy()\n",
    "data_nc = data_nc.where((data_nc.time>=start) & (data_nc.time<end),drop=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "870e90bc",
   "metadata": {},
   "source": [
    "### Filtered Xarray"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "8c423f8a",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block !important;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2 {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.Dataset&gt;\n",
       "Dimensions:      (time: 21, y: 1354, x: 1609)\n",
       "Coordinates:\n",
       "  * time         (time) datetime64[ns] 2019-01-04 2019-01-07 ... 2019-02-28\n",
       "  * y            (y) float64 4.156e+06 4.156e+06 ... 4.142e+06 4.142e+06\n",
       "  * x            (x) float64 3.611e+06 3.611e+06 ... 3.627e+06 3.627e+06\n",
       "    spatial_ref  int32 3857\n",
       "Data variables:\n",
       "    B11          (time, y, x) float64 1.894e+03 1.86e+03 ... 1.305e+03 1.018e+03\n",
       "    B12          (time, y, x) float64 1.189e+03 1.242e+03 ... 1.265e+03 826.0\n",
       "    B02          (time, y, x) float64 1.642e+03 1.59e+03 ... 219.0 243.0\n",
       "    B03          (time, y, x) float64 1.824e+03 1.704e+03 ... 284.0 332.0\n",
       "    B04          (time, y, x) float64 1.918e+03 1.774e+03 ... 223.0 252.0\n",
       "    B05          (time, y, x) float64 2.327e+03 2.105e+03 ... 748.0 599.0\n",
       "    B06          (time, y, x) float64 2.706e+03 2.631e+03 ... 1.262e+03\n",
       "    B07          (time, y, x) float64 2.751e+03 2.743e+03 ... 1.641e+03\n",
       "    B08          (time, y, x) float64 2.822e+03 2.798e+03 ... 1.516e+03\n",
       "    B8A          (time, y, x) float64 2.765e+03 2.828e+03 ... 1.827e+03\n",
       "    SCL          (time, y, x) float64 8.0 8.0 8.0 8.0 8.0 ... 10.0 4.0 4.0 4.0\n",
       "Attributes:\n",
       "    crs:           EPSG:3857\n",
       "    grid_mapping:  spatial_ref</pre><div class='xr-wrap' style='display:none'><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-02e45e97-d309-4232-a308-62bb76b3d5b8' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-02e45e97-d309-4232-a308-62bb76b3d5b8' class='xr-section-summary'  title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'><ul class='xr-dim-list'><li><span class='xr-has-index'>time</span>: 21</li><li><span class='xr-has-index'>y</span>: 1354</li><li><span class='xr-has-index'>x</span>: 1609</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-ac096ad5-fc93-4501-8969-dffd1bcec302' class='xr-section-summary-in' type='checkbox'  checked><label for='section-ac096ad5-fc93-4501-8969-dffd1bcec302' class='xr-section-summary' >Coordinates: <span>(4)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>time</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>datetime64[ns]</div><div class='xr-var-preview xr-preview'>2019-01-04 ... 2019-02-28</div><input id='attrs-a8012363-4a60-457d-b50d-7723e7a5c556' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-a8012363-4a60-457d-b50d-7723e7a5c556' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f8532864-271c-44a3-b0ba-0be16f94f159' class='xr-var-data-in' type='checkbox'><label for='data-f8532864-271c-44a3-b0ba-0be16f94f159' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([&#x27;2019-01-04T00:00:00.000000000&#x27;, &#x27;2019-01-07T00:00:00.000000000&#x27;,\n",
       "       &#x27;2019-01-09T00:00:00.000000000&#x27;, &#x27;2019-01-12T00:00:00.000000000&#x27;,\n",
       "       &#x27;2019-01-14T00:00:00.000000000&#x27;, &#x27;2019-01-17T00:00:00.000000000&#x27;,\n",
       "       &#x27;2019-01-19T00:00:00.000000000&#x27;, &#x27;2019-01-22T00:00:00.000000000&#x27;,\n",
       "       &#x27;2019-01-27T00:00:00.000000000&#x27;, &#x27;2019-02-01T00:00:00.000000000&#x27;,\n",
       "       &#x27;2019-02-03T00:00:00.000000000&#x27;, &#x27;2019-02-06T00:00:00.000000000&#x27;,\n",
       "       &#x27;2019-02-08T00:00:00.000000000&#x27;, &#x27;2019-02-11T00:00:00.000000000&#x27;,\n",
       "       &#x27;2019-02-13T00:00:00.000000000&#x27;, &#x27;2019-02-16T00:00:00.000000000&#x27;,\n",
       "       &#x27;2019-02-18T00:00:00.000000000&#x27;, &#x27;2019-02-21T00:00:00.000000000&#x27;,\n",
       "       &#x27;2019-02-23T00:00:00.000000000&#x27;, &#x27;2019-02-26T00:00:00.000000000&#x27;,\n",
       "       &#x27;2019-02-28T00:00:00.000000000&#x27;], dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>y</span></div><div class='xr-var-dims'>(y)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>4.156e+06 4.156e+06 ... 4.142e+06</div><input id='attrs-78dc2cd7-ac42-42a9-a305-99844ce75bd6' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-78dc2cd7-ac42-42a9-a305-99844ce75bd6' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-fc2a40f7-feb4-4710-be89-995969a9cbdf' class='xr-var-data-in' type='checkbox'><label for='data-fc2a40f7-feb4-4710-be89-995969a9cbdf' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>metre</dd><dt><span>resolution :</span></dt><dd>-10.0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd></dl></div><div class='xr-var-data'><pre>array([4155685., 4155675., 4155665., ..., 4142175., 4142165., 4142155.])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>x</span></div><div class='xr-var-dims'>(x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>3.611e+06 3.611e+06 ... 3.627e+06</div><input id='attrs-8df9ac6e-1d1c-41a0-b67d-4744c172c0ca' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-8df9ac6e-1d1c-41a0-b67d-4744c172c0ca' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-fc0584e9-a98c-4502-bd3b-978e35349578' class='xr-var-data-in' type='checkbox'><label for='data-fc0584e9-a98c-4502-bd3b-978e35349578' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>metre</dd><dt><span>resolution :</span></dt><dd>10.0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd></dl></div><div class='xr-var-data'><pre>array([3610515., 3610525., 3610535., ..., 3626575., 3626585., 3626595.])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>spatial_ref</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>3857</div><input id='attrs-80fb55c2-6460-4051-a306-cda0427b1162' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-80fb55c2-6460-4051-a306-cda0427b1162' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e733bf9b-030b-4c0a-98f2-31165baed2ab' class='xr-var-data-in' type='checkbox'><label for='data-e733bf9b-030b-4c0a-98f2-31165baed2ab' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>spatial_ref :</span></dt><dd>PROJCS[&quot;WGS 84 / Pseudo-Mercator&quot;,GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],PRIMEM[&quot;Greenwich&quot;,0,AUTHORITY[&quot;EPSG&quot;,&quot;8901&quot;]],UNIT[&quot;degree&quot;,0.0174532925199433,AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]],PROJECTION[&quot;Mercator_1SP&quot;],PARAMETER[&quot;central_meridian&quot;,0],PARAMETER[&quot;scale_factor&quot;,1],PARAMETER[&quot;false_easting&quot;,0],PARAMETER[&quot;false_northing&quot;,0],UNIT[&quot;metre&quot;,1,AUTHORITY[&quot;EPSG&quot;,&quot;9001&quot;]],AXIS[&quot;Easting&quot;,EAST],AXIS[&quot;Northing&quot;,NORTH],EXTENSION[&quot;PROJ4&quot;,&quot;+proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext +no_defs&quot;],AUTHORITY[&quot;EPSG&quot;,&quot;3857&quot;]]</dd><dt><span>grid_mapping_name :</span></dt><dd>??</dd></dl></div><div class='xr-var-data'><pre>array(3857, dtype=int32)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-335625a9-eb1d-43e9-993b-fb14314d2d02' class='xr-section-summary-in' type='checkbox'  checked><label for='section-335625a9-eb1d-43e9-993b-fb14314d2d02' class='xr-section-summary' >Data variables: <span>(11)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span>B11</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>1.894e+03 1.86e+03 ... 1.018e+03</div><input id='attrs-a799fb6a-37be-4446-881f-dbdd5224c950' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-a799fb6a-37be-4446-881f-dbdd5224c950' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-29ec243f-2963-4e3f-81e9-a3552c0bccc0' class='xr-var-data-in' type='checkbox'><label for='data-29ec243f-2963-4e3f-81e9-a3552c0bccc0' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>array([[[1894., 1860., 1860., ..., 3105., 3105., 3006.],\n",
       "        [1774., 1819., 1819., ..., 3105., 3105., 3006.],\n",
       "        [1774., 1819., 1819., ..., 3105., 3105., 3006.],\n",
       "        ...,\n",
       "        [2059., 1923., 1923., ..., 2956., 2956., 2908.],\n",
       "        [2221., 2111., 2111., ..., 2956., 2956., 2908.],\n",
       "        [2221., 2111., 2111., ..., 2956., 2956., 2908.]],\n",
       "\n",
       "       [[2540., 2583., 2583., ...,    0.,    0.,    0.],\n",
       "        [2604., 2649., 2649., ...,    0.,    0.,    0.],\n",
       "        [2604., 2649., 2649., ...,    0.,    0.,    0.],\n",
       "        ...,\n",
       "        [2774., 2734., 2734., ...,    0.,    0.,    0.],\n",
       "        [2791., 2769., 2769., ...,    0.,    0.,    0.],\n",
       "        [2791., 2769., 2769., ...,    0.,    0.,    0.]],\n",
       "\n",
       "       [[5483., 7149., 7149., ...,  557.,  557.,  541.],\n",
       "        [6647., 8665., 8665., ...,  557.,  557.,  541.],\n",
       "        [6647., 8665., 8665., ...,  557.,  557.,  541.],\n",
       "        ...,\n",
       "...\n",
       "        ...,\n",
       "        [1458., 1455., 1455., ...,  586.,  586.,  626.],\n",
       "        [1482., 1493., 1493., ...,  586.,  586.,  626.],\n",
       "        [1482., 1493., 1493., ...,  586.,  586.,  626.]],\n",
       "\n",
       "       [[5537., 5665., 5665., ...,    0.,    0.,    0.],\n",
       "        [5616., 5676., 5676., ...,    0.,    0.,    0.],\n",
       "        [5616., 5676., 5676., ...,    0.,    0.,    0.],\n",
       "        ...,\n",
       "        [5672., 5939., 5939., ...,    0.,    0.,    0.],\n",
       "        [5413., 5403., 5403., ...,    0.,    0.,    0.],\n",
       "        [5413., 5403., 5403., ...,    0.,    0.,    0.]],\n",
       "\n",
       "       [[3109., 2804., 2804., ..., 1355., 1355., 1262.],\n",
       "        [3122., 2565., 2565., ..., 1355., 1355., 1262.],\n",
       "        [3122., 2565., 2565., ..., 1355., 1355., 1262.],\n",
       "        ...,\n",
       "        [1793., 2224., 2224., ..., 1305., 1305., 1018.],\n",
       "        [2129., 2144., 2144., ..., 1305., 1305., 1018.],\n",
       "        [2129., 2144., 2144., ..., 1305., 1305., 1018.]]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B12</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>1.189e+03 1.242e+03 ... 826.0</div><input id='attrs-618b465e-a1b0-4640-af6e-144b7eddbf31' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-618b465e-a1b0-4640-af6e-144b7eddbf31' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b4651028-8c7b-40aa-b577-63f27bb0177f' class='xr-var-data-in' type='checkbox'><label for='data-b4651028-8c7b-40aa-b577-63f27bb0177f' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>array([[[1189., 1242., 1242., ..., 2335., 2335., 2444.],\n",
       "        [1063., 1192., 1192., ..., 2335., 2335., 2444.],\n",
       "        [1063., 1192., 1192., ..., 2335., 2335., 2444.],\n",
       "        ...,\n",
       "        [1664., 1554., 1554., ..., 2379., 2379., 2402.],\n",
       "        [1814., 1720., 1720., ..., 2379., 2379., 2402.],\n",
       "        [1814., 1720., 1720., ..., 2379., 2379., 2402.]],\n",
       "\n",
       "       [[1809., 1847., 1847., ...,    0.,    0.,    0.],\n",
       "        [1848., 1886., 1886., ...,    0.,    0.,    0.],\n",
       "        [1848., 1886., 1886., ...,    0.,    0.,    0.],\n",
       "        ...,\n",
       "        [2059., 2045., 2045., ...,    0.,    0.,    0.],\n",
       "        [2130., 2110., 2110., ...,    0.,    0.,    0.],\n",
       "        [2130., 2110., 2110., ...,    0.,    0.,    0.]],\n",
       "\n",
       "       [[5117., 6934., 6934., ...,  310.,  310.,  313.],\n",
       "        [6139., 8301., 8301., ...,  310.,  310.,  313.],\n",
       "        [6139., 8301., 8301., ...,  310.,  310.,  313.],\n",
       "        ...,\n",
       "...\n",
       "        ...,\n",
       "        [1595., 1590., 1590., ...,  374.,  374.,  389.],\n",
       "        [1613., 1628., 1628., ...,  374.,  374.,  389.],\n",
       "        [1613., 1628., 1628., ...,  374.,  374.,  389.]],\n",
       "\n",
       "       [[3623., 3997., 3997., ...,    0.,    0.,    0.],\n",
       "        [3747., 3927., 3927., ...,    0.,    0.,    0.],\n",
       "        [3747., 3927., 3927., ...,    0.,    0.,    0.],\n",
       "        ...,\n",
       "        [5292., 5543., 5543., ...,    0.,    0.,    0.],\n",
       "        [5031., 5155., 5155., ...,    0.,    0.,    0.],\n",
       "        [5031., 5155., 5155., ...,    0.,    0.,    0.]],\n",
       "\n",
       "       [[2889., 2366., 2366., ...,  860.,  860.,  811.],\n",
       "        [2890., 2230., 2230., ...,  860.,  860.,  811.],\n",
       "        [2890., 2230., 2230., ...,  860.,  860.,  811.],\n",
       "        ...,\n",
       "        [1132., 1689., 1689., ..., 1265., 1265.,  826.],\n",
       "        [1522., 1711., 1711., ..., 1265., 1265.,  826.],\n",
       "        [1522., 1711., 1711., ..., 1265., 1265.,  826.]]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B02</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>1.642e+03 1.59e+03 ... 219.0 243.0</div><input id='attrs-59a7eafb-d0c6-412d-8acb-6ef22c0cea33' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-59a7eafb-d0c6-412d-8acb-6ef22c0cea33' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-2cf075bc-1378-482c-a0f7-78d13c672501' class='xr-var-data-in' type='checkbox'><label for='data-2cf075bc-1378-482c-a0f7-78d13c672501' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>array([[[ 1642.,  1590.,  1520., ...,  8600.,  8600.,  8504.],\n",
       "        [ 1554.,  1520.,  1460., ...,  8552.,  8552.,  8440.],\n",
       "        [ 1442.,  1386.,  1272., ...,  8552.,  8552.,  8440.],\n",
       "        ...,\n",
       "        [ 1394.,  1670.,  1670., ...,  4792.,  4896.,  4956.],\n",
       "        [ 1430.,  1664.,  1664., ...,  4928.,  4996.,  4984.],\n",
       "        [ 1530.,  1722.,  1722., ...,  4928.,  4996.,  4984.]],\n",
       "\n",
       "       [[ 7192.,  7220.,  7264., ...,     0.,     0.,     0.],\n",
       "        [ 7272.,  7272.,  7320., ...,     0.,     0.,     0.],\n",
       "        [ 7384.,  7452.,  7416., ...,     0.,     0.,     0.],\n",
       "        ...,\n",
       "        [ 7208.,  7140.,  7140., ...,     0.,     0.,     0.],\n",
       "        [ 7132.,  7096.,  7096., ...,     0.,     0.,     0.],\n",
       "        [ 7208.,  7232.,  7232., ...,     0.,     0.,     0.]],\n",
       "\n",
       "       [[ 4908.,  5980.,  6552., ...,   401.,   401.,   379.],\n",
       "        [ 5464.,  6392.,  7276., ...,   510.,   510.,   422.],\n",
       "        [ 6032.,  6992.,  7772., ...,   510.,   510.,   422.],\n",
       "        ...,\n",
       "...\n",
       "        ...,\n",
       "        [ 8352.,  8424.,  8424., ...,   856.,   835.,   822.],\n",
       "        [ 8480.,  8560.,  8560., ...,   798.,   758.,   773.],\n",
       "        [ 8728.,  8792.,  8792., ...,   798.,   758.,   773.]],\n",
       "\n",
       "       [[12488., 11976., 11720., ...,     0.,     0.,     0.],\n",
       "        [12656., 12160., 11832., ...,     0.,     0.,     0.],\n",
       "        [12968., 12464., 12096., ...,     0.,     0.,     0.],\n",
       "        ...,\n",
       "        [ 6684.,  6260.,  6260., ...,     0.,     0.,     0.],\n",
       "        [ 6324.,  6152.,  6152., ...,     0.,     0.,     0.],\n",
       "        [ 5892.,  6012.,  6012., ...,     0.,     0.,     0.]],\n",
       "\n",
       "       [[ 1842.,  1488.,  1086., ...,   376.,   376.,   362.],\n",
       "        [ 1618.,  1294.,   984., ...,   350.,   350.,   379.],\n",
       "        [ 1268.,   977.,   902., ...,   350.,   350.,   379.],\n",
       "        ...,\n",
       "        [  425.,   402.,   402., ...,   298.,   241.,   244.],\n",
       "        [  655.,   933.,   933., ...,   268.,   219.,   243.],\n",
       "        [ 1164.,  1144.,  1144., ...,   268.,   219.,   243.]]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B03</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>1.824e+03 1.704e+03 ... 284.0 332.0</div><input id='attrs-7833cf21-0e04-47f8-a29b-194faafaa0cb' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-7833cf21-0e04-47f8-a29b-194faafaa0cb' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e5e898c5-4ee1-409a-8f9f-06846f0bc58a' class='xr-var-data-in' type='checkbox'><label for='data-e5e898c5-4ee1-409a-8f9f-06846f0bc58a' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>array([[[ 1824.,  1704.,  1644., ...,  7720.,  7720.,  7672.],\n",
       "        [ 1824.,  1668.,  1516., ...,  7628.,  7628.,  7540.],\n",
       "        [ 1894.,  1704.,  1554., ...,  7628.,  7628.,  7540.],\n",
       "        ...,\n",
       "        [ 1738.,  1752.,  1752., ...,  4764.,  4860.,  4988.],\n",
       "        [ 1844.,  1832.,  1832., ...,  4676.,  4788.,  4852.],\n",
       "        [ 1892.,  1868.,  1868., ...,  4676.,  4788.,  4852.]],\n",
       "\n",
       "       [[ 6156.,  6156.,  6192., ...,     0.,     0.,     0.],\n",
       "        [ 6220.,  6300.,  6368., ...,     0.,     0.,     0.],\n",
       "        [ 6400.,  6388.,  6444., ...,     0.,     0.,     0.],\n",
       "        ...,\n",
       "        [ 6432.,  6424.,  6424., ...,     0.,     0.,     0.],\n",
       "        [ 6440.,  6388.,  6388., ...,     0.,     0.,     0.],\n",
       "        [ 6464.,  6408.,  6408., ...,     0.,     0.,     0.]],\n",
       "\n",
       "       [[ 4696.,  5716.,  6572., ...,   466.,   466.,   398.],\n",
       "        [ 5340.,  6220.,  7348., ...,   475.,   475.,   395.],\n",
       "        [ 5968.,  6720.,  7880., ...,   475.,   475.,   395.],\n",
       "        ...,\n",
       "...\n",
       "        ...,\n",
       "        [ 8180.,  8224.,  8224., ...,  1008.,   969.,   993.],\n",
       "        [ 8280.,  8368.,  8368., ...,   929.,   890.,   962.],\n",
       "        [ 8528.,  8568.,  8568., ...,   929.,   890.,   962.]],\n",
       "\n",
       "       [[11328., 11048., 11048., ...,     0.,     0.,     0.],\n",
       "        [11368., 11016., 11024., ...,     0.,     0.,     0.],\n",
       "        [11512., 11056., 10968., ...,     0.,     0.,     0.],\n",
       "        ...,\n",
       "        [ 6276.,  6208.,  6208., ...,     0.,     0.,     0.],\n",
       "        [ 6208.,  6016.,  6016., ...,     0.,     0.,     0.],\n",
       "        [ 5936.,  5888.,  5888., ...,     0.,     0.,     0.]],\n",
       "\n",
       "       [[ 1942.,  1812.,  1500., ...,   467.,   467.,   461.],\n",
       "        [ 1730.,  1500.,  1184., ...,   465.,   465.,   523.],\n",
       "        [ 1570.,  1306.,  1016., ...,   465.,   465.,   523.],\n",
       "        ...,\n",
       "        [  582.,   744.,   744., ...,   397.,   326.,   300.],\n",
       "        [  941.,  1100.,  1100., ...,   362.,   284.,   332.],\n",
       "        [ 1368.,  1390.,  1390., ...,   362.,   284.,   332.]]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B04</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>1.918e+03 1.774e+03 ... 223.0 252.0</div><input id='attrs-740112fa-57c7-4dee-88c8-8a053dca0568' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-740112fa-57c7-4dee-88c8-8a053dca0568' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-6dd94130-64f9-4a52-8287-93f07da0169a' class='xr-var-data-in' type='checkbox'><label for='data-6dd94130-64f9-4a52-8287-93f07da0169a' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>array([[[ 1918.,  1774.,  1598., ...,  7104.,  7104.,  7080.],\n",
       "        [ 1990.,  1844.,  1700., ...,  7060.,  7060.,  7064.],\n",
       "        [ 1902.,  1798.,  1736., ...,  7060.,  7060.,  7064.],\n",
       "        ...,\n",
       "        [ 1868.,  1888.,  1888., ...,  4748.,  4800.,  4780.],\n",
       "        [ 1990.,  1956.,  1956., ...,  4696.,  4784.,  4848.],\n",
       "        [ 2122.,  2126.,  2126., ...,  4696.,  4784.,  4848.]],\n",
       "\n",
       "       [[ 5468.,  5504.,  5552., ...,     0.,     0.,     0.],\n",
       "        [ 5516.,  5592.,  5620., ...,     0.,     0.,     0.],\n",
       "        [ 5600.,  5620.,  5640., ...,     0.,     0.,     0.],\n",
       "        ...,\n",
       "        [ 5800.,  5788.,  5788., ...,     0.,     0.,     0.],\n",
       "        [ 5904.,  5776.,  5776., ...,     0.,     0.,     0.],\n",
       "        [ 5960.,  5844.,  5844., ...,     0.,     0.,     0.]],\n",
       "\n",
       "       [[ 4840.,  5440.,  6476., ...,   402.,   402.,   367.],\n",
       "        [ 5260.,  6064.,  7192., ...,   358.,   358.,   311.],\n",
       "        [ 5700.,  6680.,  7848., ...,   358.,   358.,   311.],\n",
       "        ...,\n",
       "...\n",
       "        ...,\n",
       "        [ 8240.,  8280.,  8280., ...,   961.,   959.,   940.],\n",
       "        [ 8400.,  8424.,  8424., ...,   922.,   914.,   912.],\n",
       "        [ 8632.,  8616.,  8616., ...,   922.,   914.,   912.]],\n",
       "\n",
       "       [[10568., 10408., 10360., ...,     0.,     0.,     0.],\n",
       "        [10696., 10536., 10592., ...,     0.,     0.,     0.],\n",
       "        [10664., 10496., 10744., ...,     0.,     0.,     0.],\n",
       "        ...,\n",
       "        [ 6112.,  6372.,  6372., ...,     0.,     0.,     0.],\n",
       "        [ 5976.,  6092.,  6092., ...,     0.,     0.,     0.],\n",
       "        [ 5824.,  5888.,  5888., ...,     0.,     0.,     0.]],\n",
       "\n",
       "       [[ 1876.,  1690.,  1560., ...,   414.,   414.,   397.],\n",
       "        [ 1704.,  1500.,  1256., ...,   379.,   379.,   415.],\n",
       "        [ 1622.,  1436.,  1198., ...,   379.,   379.,   415.],\n",
       "        ...,\n",
       "        [  474.,   529.,   529., ...,   564.,   250.,   263.],\n",
       "        [  809.,   989.,   989., ...,   427.,   223.,   252.],\n",
       "        [ 1352.,  1382.,  1382., ...,   427.,   223.,   252.]]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B05</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>2.327e+03 2.105e+03 ... 748.0 599.0</div><input id='attrs-712cee17-c0b0-4c97-82a9-f53ff549bfc9' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-712cee17-c0b0-4c97-82a9-f53ff549bfc9' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-006684dd-0a02-4927-b546-de1888dc16a2' class='xr-var-data-in' type='checkbox'><label for='data-006684dd-0a02-4927-b546-de1888dc16a2' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>array([[[ 2327.,  2105.,  2105., ...,  7117.,  7117.,  7053.],\n",
       "        [ 2141.,  2104.,  2104., ...,  7117.,  7117.,  7053.],\n",
       "        [ 2141.,  2104.,  2104., ...,  7117.,  7117.,  7053.],\n",
       "        ...,\n",
       "        [ 2312.,  2095.,  2095., ...,  4901.,  4901.,  4947.],\n",
       "        [ 2569.,  2362.,  2362., ...,  4901.,  4901.,  4947.],\n",
       "        [ 2569.,  2362.,  2362., ...,  4901.,  4901.,  4947.]],\n",
       "\n",
       "       [[ 5548.,  5560.,  5560., ...,     0.,     0.,     0.],\n",
       "        [ 5637.,  5659.,  5659., ...,     0.,     0.,     0.],\n",
       "        [ 5637.,  5659.,  5659., ...,     0.,     0.,     0.],\n",
       "        ...,\n",
       "        [ 5847.,  5766.,  5766., ...,     0.,     0.,     0.],\n",
       "        [ 5948.,  5783.,  5783., ...,     0.,     0.,     0.],\n",
       "        [ 5948.,  5783.,  5783., ...,     0.,     0.,     0.]],\n",
       "\n",
       "       [[ 4616.,  6137.,  6137., ...,   459.,   459.,   405.],\n",
       "        [ 5544.,  7407.,  7407., ...,   459.,   459.,   405.],\n",
       "        [ 5544.,  7407.,  7407., ...,   459.,   459.,   405.],\n",
       "        ...,\n",
       "...\n",
       "        ...,\n",
       "        [ 8327.,  8334.,  8334., ...,  1207.,  1207.,  1205.],\n",
       "        [ 8598.,  8620.,  8620., ...,  1207.,  1207.,  1205.],\n",
       "        [ 8598.,  8620.,  8620., ...,  1207.,  1207.,  1205.]],\n",
       "\n",
       "       [[10974., 10917., 10917., ...,     0.,     0.,     0.],\n",
       "        [11196., 10955., 10955., ...,     0.,     0.,     0.],\n",
       "        [11196., 10955., 10955., ...,     0.,     0.,     0.],\n",
       "        ...,\n",
       "        [ 6898.,  7418.,  7418., ...,     0.,     0.,     0.],\n",
       "        [ 6553.,  6561.,  6561., ...,     0.,     0.,     0.],\n",
       "        [ 6553.,  6561.,  6561., ...,     0.,     0.,     0.]],\n",
       "\n",
       "       [[ 2499.,  2233.,  2233., ...,   755.,   755.,   718.],\n",
       "        [ 2512.,  1912.,  1912., ...,   755.,   755.,   718.],\n",
       "        [ 2512.,  1912.,  1912., ...,   755.,   755.,   718.],\n",
       "        ...,\n",
       "        [  969.,  1505.,  1505., ...,   748.,   748.,   599.],\n",
       "        [ 1448.,  1425.,  1425., ...,   748.,   748.,   599.],\n",
       "        [ 1448.,  1425.,  1425., ...,   748.,   748.,   599.]]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B06</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>2.706e+03 2.631e+03 ... 1.262e+03</div><input id='attrs-898b0c64-668e-403a-a826-f65e2ee63af9' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-898b0c64-668e-403a-a826-f65e2ee63af9' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f1843a34-430d-4040-9f44-2f4a8081e464' class='xr-var-data-in' type='checkbox'><label for='data-f1843a34-430d-4040-9f44-2f4a8081e464' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>array([[[ 2706.,  2631.,  2631., ...,  6947.,  6947.,  6885.],\n",
       "        [ 2441.,  2540.,  2540., ...,  6947.,  6947.,  6885.],\n",
       "        [ 2441.,  2540.,  2540., ...,  6947.,  6947.,  6885.],\n",
       "        ...,\n",
       "        [ 2748.,  2475.,  2475., ...,  4789.,  4789.,  4929.],\n",
       "        [ 2923.,  2750.,  2750., ...,  4789.,  4789.,  4929.],\n",
       "        [ 2923.,  2750.,  2750., ...,  4789.,  4789.,  4929.]],\n",
       "\n",
       "       [[ 5402.,  5433.,  5433., ...,     0.,     0.,     0.],\n",
       "        [ 5424.,  5461.,  5461., ...,     0.,     0.,     0.],\n",
       "        [ 5424.,  5461.,  5461., ...,     0.,     0.,     0.],\n",
       "        ...,\n",
       "        [ 5666.,  5631.,  5631., ...,     0.,     0.,     0.],\n",
       "        [ 5774.,  5727.,  5727., ...,     0.,     0.,     0.],\n",
       "        [ 5774.,  5727.,  5727., ...,     0.,     0.,     0.]],\n",
       "\n",
       "       [[ 5329.,  6744.,  6744., ...,   659.,   659.,   634.],\n",
       "        [ 6244.,  8200.,  8200., ...,   659.,   659.,   634.],\n",
       "        [ 6244.,  8200.,  8200., ...,   659.,   659.,   634.],\n",
       "        ...,\n",
       "...\n",
       "        ...,\n",
       "        [ 8208.,  8179.,  8179., ...,  1732.,  1732.,  1933.],\n",
       "        [ 8414.,  8480.,  8480., ...,  1732.,  1732.,  1933.],\n",
       "        [ 8414.,  8480.,  8480., ...,  1732.,  1732.,  1933.]],\n",
       "\n",
       "       [[10356., 10418., 10418., ...,     0.,     0.,     0.],\n",
       "        [10608., 10452., 10452., ...,     0.,     0.,     0.],\n",
       "        [10608., 10452., 10452., ...,     0.,     0.,     0.],\n",
       "        ...,\n",
       "        [ 7140.,  7645.,  7645., ...,     0.,     0.,     0.],\n",
       "        [ 6658.,  6828.,  6828., ...,     0.,     0.,     0.],\n",
       "        [ 6658.,  6828.,  6828., ...,     0.,     0.,     0.]],\n",
       "\n",
       "       [[ 3486.,  3134.,  3134., ...,  1403.,  1403.,  1372.],\n",
       "        [ 3505.,  2876.,  2876., ...,  1403.,  1403.,  1372.],\n",
       "        [ 3505.,  2876.,  2876., ...,  1403.,  1403.,  1372.],\n",
       "        ...,\n",
       "        [ 2665.,  2554.,  2554., ...,  1594.,  1594.,  1262.],\n",
       "        [ 2736.,  2364.,  2364., ...,  1594.,  1594.,  1262.],\n",
       "        [ 2736.,  2364.,  2364., ...,  1594.,  1594.,  1262.]]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B07</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>2.751e+03 2.743e+03 ... 1.641e+03</div><input id='attrs-680d69f8-6a7c-4eb8-9cf3-08a5f3b811fa' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-680d69f8-6a7c-4eb8-9cf3-08a5f3b811fa' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8ab0c92f-9a4c-4185-b7c3-cc5ce6179853' class='xr-var-data-in' type='checkbox'><label for='data-8ab0c92f-9a4c-4185-b7c3-cc5ce6179853' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>array([[[ 2751.,  2743.,  2743., ...,  7003.,  7003.,  6904.],\n",
       "        [ 2425.,  2624.,  2624., ...,  7003.,  7003.,  6904.],\n",
       "        [ 2425.,  2624.,  2624., ...,  7003.,  7003.,  6904.],\n",
       "        ...,\n",
       "        [ 2934.,  2716.,  2716., ...,  4738.,  4738.,  4855.],\n",
       "        [ 3103.,  2996.,  2996., ...,  4738.,  4738.,  4855.],\n",
       "        [ 3103.,  2996.,  2996., ...,  4738.,  4738.,  4855.]],\n",
       "\n",
       "       [[ 5182.,  5216.,  5216., ...,     0.,     0.,     0.],\n",
       "        [ 5252.,  5274.,  5274., ...,     0.,     0.,     0.],\n",
       "        [ 5252.,  5274.,  5274., ...,     0.,     0.,     0.],\n",
       "        ...,\n",
       "        [ 5534.,  5465.,  5465., ...,     0.,     0.,     0.],\n",
       "        [ 5657.,  5567.,  5567., ...,     0.,     0.,     0.],\n",
       "        [ 5657.,  5567.,  5567., ...,     0.,     0.,     0.]],\n",
       "\n",
       "       [[ 5258.,  6923.,  6923., ...,   646.,   646.,   697.],\n",
       "        [ 6195.,  8294.,  8294., ...,   646.,   646.,   697.],\n",
       "        [ 6195.,  8294.,  8294., ...,   646.,   646.,   697.],\n",
       "        ...,\n",
       "...\n",
       "        ...,\n",
       "        [ 8071.,  8100.,  8100., ...,  1840.,  1840.,  2105.],\n",
       "        [ 8318.,  8396.,  8396., ...,  1840.,  1840.,  2105.],\n",
       "        [ 8318.,  8396.,  8396., ...,  1840.,  1840.,  2105.]],\n",
       "\n",
       "       [[ 9936., 10090., 10090., ...,     0.,     0.,     0.],\n",
       "        [10111., 10060., 10060., ...,     0.,     0.,     0.],\n",
       "        [10111., 10060., 10060., ...,     0.,     0.,     0.],\n",
       "        ...,\n",
       "        [ 7130.,  7703.,  7703., ...,     0.,     0.,     0.],\n",
       "        [ 6670.,  6992.,  6992., ...,     0.,     0.,     0.],\n",
       "        [ 6670.,  6992.,  6992., ...,     0.,     0.,     0.]],\n",
       "\n",
       "       [[ 3826.,  3356.,  3356., ...,  1545.,  1545.,  1572.],\n",
       "        [ 3812.,  3186.,  3186., ...,  1545.,  1545.,  1572.],\n",
       "        [ 3812.,  3186.,  3186., ...,  1545.,  1545.,  1572.],\n",
       "        ...,\n",
       "        [ 2947.,  2855.,  2855., ...,  1966.,  1966.,  1641.],\n",
       "        [ 3068.,  2618.,  2618., ...,  1966.,  1966.,  1641.],\n",
       "        [ 3068.,  2618.,  2618., ...,  1966.,  1966.,  1641.]]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B08</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>2.822e+03 2.798e+03 ... 1.516e+03</div><input id='attrs-d08b0d35-6ffb-418c-b858-c73f9157d80a' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-d08b0d35-6ffb-418c-b858-c73f9157d80a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-fb88ac7d-3ef8-4987-834e-d20c2ef58273' class='xr-var-data-in' type='checkbox'><label for='data-fb88ac7d-3ef8-4987-834e-d20c2ef58273' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>array([[[ 2822.,  2798.,  2696., ...,  7044.,  7044.,  7044.],\n",
       "        [ 2764.,  2686.,  2600., ...,  6960.,  6960.,  6916.],\n",
       "        [ 2790.,  2644.,  2536., ...,  6960.,  6960.,  6916.],\n",
       "        ...,\n",
       "        [ 2404.,  2570.,  2570., ...,  4604.,  4748.,  4792.],\n",
       "        [ 2464.,  2576.,  2576., ...,  4672.,  4712.,  4740.],\n",
       "        [ 2576.,  2574.,  2574., ...,  4672.,  4712.,  4740.]],\n",
       "\n",
       "       [[ 5704.,  5712.,  5804., ...,     0.,     0.,     0.],\n",
       "        [ 5760.,  5864.,  5856., ...,     0.,     0.,     0.],\n",
       "        [ 5820.,  5904.,  5912., ...,     0.,     0.,     0.],\n",
       "        ...,\n",
       "        [ 6036.,  6068.,  6068., ...,     0.,     0.,     0.],\n",
       "        [ 6128.,  6100.,  6100., ...,     0.,     0.,     0.],\n",
       "        [ 6132.,  6144.,  6144., ...,     0.,     0.,     0.]],\n",
       "\n",
       "       [[ 6456.,  7612.,  8664., ...,   703.,   703.,   653.],\n",
       "        [ 7280.,  8472.,  9720., ...,   827.,   827.,   688.],\n",
       "        [ 8160.,  9280., 10672., ...,   827.,   827.,   688.],\n",
       "        ...,\n",
       "...\n",
       "        ...,\n",
       "        [ 8488.,  8584.,  8584., ...,  1548.,  1592.,  1794.],\n",
       "        [ 8656.,  8712.,  8712., ...,  1492.,  1622.,  1980.],\n",
       "        [ 8840.,  8904.,  8904., ...,  1492.,  1622.,  1980.]],\n",
       "\n",
       "       [[11048., 10568., 10448., ...,     0.,     0.,     0.],\n",
       "        [11184., 10664., 10456., ...,     0.,     0.,     0.],\n",
       "        [11464., 10880., 10664., ...,     0.,     0.,     0.],\n",
       "        ...,\n",
       "        [ 6920.,  6668.,  6668., ...,     0.,     0.,     0.],\n",
       "        [ 6764.,  6612.,  6612., ...,     0.,     0.,     0.],\n",
       "        [ 6420.,  6440.,  6440., ...,     0.,     0.,     0.]],\n",
       "\n",
       "       [[ 3648.,  3378.,  2908., ...,  1606.,  1606.,  1668.],\n",
       "        [ 3356.,  3032.,  2670., ...,  1544.,  1544.,  1778.],\n",
       "        [ 3112.,  2688.,  2580., ...,  1544.,  1544.,  1778.],\n",
       "        ...,\n",
       "        [ 3360.,  3308.,  3308., ...,  1220.,   994.,  1056.],\n",
       "        [ 3480.,  3038.,  3038., ...,  1178.,  1034.,  1516.],\n",
       "        [ 2684.,  2574.,  2574., ...,  1178.,  1034.,  1516.]]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B8A</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>2.765e+03 2.828e+03 ... 1.827e+03</div><input id='attrs-24141796-4bb1-4cc1-8fb5-d94a3806be93' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-24141796-4bb1-4cc1-8fb5-d94a3806be93' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8e033780-2290-4e47-b931-0b076484ab0c' class='xr-var-data-in' type='checkbox'><label for='data-8e033780-2290-4e47-b931-0b076484ab0c' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>array([[[2765., 2828., 2828., ..., 6705., 6705., 6748.],\n",
       "        [2447., 2627., 2627., ..., 6705., 6705., 6748.],\n",
       "        [2447., 2627., 2627., ..., 6705., 6705., 6748.],\n",
       "        ...,\n",
       "        [3116., 2930., 2930., ..., 4707., 4707., 4787.],\n",
       "        [3306., 3197., 3197., ..., 4707., 4707., 4787.],\n",
       "        [3306., 3197., 3197., ..., 4707., 4707., 4787.]],\n",
       "\n",
       "       [[4911., 5007., 5007., ...,    0.,    0.,    0.],\n",
       "        [5040., 5133., 5133., ...,    0.,    0.,    0.],\n",
       "        [5040., 5133., 5133., ...,    0.,    0.,    0.],\n",
       "        ...,\n",
       "        [5342., 5377., 5377., ...,    0.,    0.,    0.],\n",
       "        [5449., 5451., 5451., ...,    0.,    0.,    0.],\n",
       "        [5449., 5451., 5451., ...,    0.,    0.,    0.]],\n",
       "\n",
       "       [[5391., 7110., 7110., ...,  635.,  635.,  682.],\n",
       "        [6374., 8302., 8302., ...,  635.,  635.,  682.],\n",
       "        [6374., 8302., 8302., ...,  635.,  635.,  682.],\n",
       "        ...,\n",
       "...\n",
       "        ...,\n",
       "        [7963., 7946., 7946., ..., 1922., 1922., 2219.],\n",
       "        [8195., 8243., 8243., ..., 1922., 1922., 2219.],\n",
       "        [8195., 8243., 8243., ..., 1922., 1922., 2219.]],\n",
       "\n",
       "       [[9483., 9684., 9684., ...,    0.,    0.,    0.],\n",
       "        [9732., 9671., 9671., ...,    0.,    0.,    0.],\n",
       "        [9732., 9671., 9671., ...,    0.,    0.,    0.],\n",
       "        ...,\n",
       "        [7214., 7775., 7775., ...,    0.,    0.,    0.],\n",
       "        [6777., 7035., 7035., ...,    0.,    0.,    0.],\n",
       "        [6777., 7035., 7035., ...,    0.,    0.,    0.]],\n",
       "\n",
       "       [[4174., 3629., 3629., ..., 1726., 1726., 1734.],\n",
       "        [4037., 3478., 3478., ..., 1726., 1726., 1734.],\n",
       "        [4037., 3478., 3478., ..., 1726., 1726., 1734.],\n",
       "        ...,\n",
       "        [3317., 2988., 2988., ..., 2027., 2027., 1827.],\n",
       "        [3142., 2688., 2688., ..., 2027., 2027., 1827.],\n",
       "        [3142., 2688., 2688., ..., 2027., 2027., 1827.]]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>SCL</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>8.0 8.0 8.0 8.0 ... 4.0 4.0 4.0</div><input id='attrs-efdcf265-aaea-406a-99da-74e146b57caa' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-efdcf265-aaea-406a-99da-74e146b57caa' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-9d76c4d2-5a44-4704-a4fe-c63f53610204' class='xr-var-data-in' type='checkbox'><label for='data-9d76c4d2-5a44-4704-a4fe-c63f53610204' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>1</dd><dt><span>nodata :</span></dt><dd>0</dd><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div><div class='xr-var-data'><pre>array([[[ 8.,  8.,  8., ...,  9.,  9.,  9.],\n",
       "        [ 8.,  8.,  8., ...,  9.,  9.,  9.],\n",
       "        [ 8.,  8.,  8., ...,  9.,  9.,  9.],\n",
       "        ...,\n",
       "        [ 8.,  8.,  8., ...,  9.,  9.,  9.],\n",
       "        [ 8.,  8.,  8., ...,  9.,  9.,  9.],\n",
       "        [ 8.,  8.,  8., ...,  9.,  9.,  9.]],\n",
       "\n",
       "       [[ 9.,  9.,  9., ...,  0.,  0.,  0.],\n",
       "        [ 9.,  9.,  9., ...,  0.,  0.,  0.],\n",
       "        [ 9.,  9.,  9., ...,  0.,  0.,  0.],\n",
       "        ...,\n",
       "        [ 9.,  9.,  9., ...,  0.,  0.,  0.],\n",
       "        [ 9.,  9.,  9., ...,  0.,  0.,  0.],\n",
       "        [ 9.,  9.,  9., ...,  0.,  0.,  0.]],\n",
       "\n",
       "       [[ 9.,  9.,  9., ...,  3.,  3.,  3.],\n",
       "        [ 9.,  9.,  9., ...,  3.,  3.,  3.],\n",
       "        [ 9.,  9.,  9., ...,  3.,  3.,  3.],\n",
       "        ...,\n",
       "...\n",
       "        ...,\n",
       "        [11., 11., 11., ...,  7.,  7.,  7.],\n",
       "        [11., 11., 11., ...,  7.,  7.,  7.],\n",
       "        [11., 11., 11., ...,  7.,  7.,  7.]],\n",
       "\n",
       "       [[ 9.,  9.,  9., ...,  0.,  0.,  0.],\n",
       "        [ 9.,  9.,  9., ...,  0.,  0.,  0.],\n",
       "        [ 9.,  9.,  9., ...,  0.,  0.,  0.],\n",
       "        ...,\n",
       "        [ 9.,  9.,  9., ...,  0.,  0.,  0.],\n",
       "        [ 9.,  9.,  9., ...,  0.,  0.,  0.],\n",
       "        [ 9.,  9.,  9., ...,  0.,  0.,  0.]],\n",
       "\n",
       "       [[10., 10., 10., ...,  4.,  4.,  4.],\n",
       "        [10., 10., 10., ...,  4.,  4.,  4.],\n",
       "        [10., 10., 10., ...,  4.,  4.,  4.],\n",
       "        ...,\n",
       "        [ 4.,  4.,  4., ...,  4.,  4.,  4.],\n",
       "        [ 4.,  5.,  5., ...,  4.,  4.,  4.],\n",
       "        [ 4.,  5.,  5., ...,  4.,  4.,  4.]]])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-0516f6f6-81e5-4b24-9a1b-ae7e3b2255d1' class='xr-section-summary-in' type='checkbox'  checked><label for='section-0516f6f6-81e5-4b24-9a1b-ae7e3b2255d1' class='xr-section-summary' >Attributes: <span>(2)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>crs :</span></dt><dd>EPSG:3857</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset>\n",
       "Dimensions:      (time: 21, y: 1354, x: 1609)\n",
       "Coordinates:\n",
       "  * time         (time) datetime64[ns] 2019-01-04 2019-01-07 ... 2019-02-28\n",
       "  * y            (y) float64 4.156e+06 4.156e+06 ... 4.142e+06 4.142e+06\n",
       "  * x            (x) float64 3.611e+06 3.611e+06 ... 3.627e+06 3.627e+06\n",
       "    spatial_ref  int32 3857\n",
       "Data variables:\n",
       "    B11          (time, y, x) float64 1.894e+03 1.86e+03 ... 1.305e+03 1.018e+03\n",
       "    B12          (time, y, x) float64 1.189e+03 1.242e+03 ... 1.265e+03 826.0\n",
       "    B02          (time, y, x) float64 1.642e+03 1.59e+03 ... 219.0 243.0\n",
       "    B03          (time, y, x) float64 1.824e+03 1.704e+03 ... 284.0 332.0\n",
       "    B04          (time, y, x) float64 1.918e+03 1.774e+03 ... 223.0 252.0\n",
       "    B05          (time, y, x) float64 2.327e+03 2.105e+03 ... 748.0 599.0\n",
       "    B06          (time, y, x) float64 2.706e+03 2.631e+03 ... 1.262e+03\n",
       "    B07          (time, y, x) float64 2.751e+03 2.743e+03 ... 1.641e+03\n",
       "    B08          (time, y, x) float64 2.822e+03 2.798e+03 ... 1.516e+03\n",
       "    B8A          (time, y, x) float64 2.765e+03 2.828e+03 ... 1.827e+03\n",
       "    SCL          (time, y, x) float64 8.0 8.0 8.0 8.0 8.0 ... 10.0 4.0 4.0 4.0\n",
       "Attributes:\n",
       "    crs:           EPSG:3857\n",
       "    grid_mapping:  spatial_ref"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data_nc"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "674da1c7",
   "metadata": {},
   "source": [
    "### Get the dates related to the data and data intervals"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "cd2f679c",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[('2019-01-01', '2019-01-07'),\n",
       " ('2019-01-08', '2019-01-14'),\n",
       " ('2019-01-15', '2019-01-22'),\n",
       " ('2019-01-23', '2019-01-29'),\n",
       " ('2019-01-30', '2019-02-05'),\n",
       " ('2019-02-06', '2019-02-13'),\n",
       " ('2019-02-14', '2019-02-20'),\n",
       " ('2019-02-21', '2019-02-28')]"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from datetime import timedelta\n",
    "dates_l = list(date_range(start_date,end_date,d_break))\n",
    "dates_list = []\n",
    "for i in range(len(dates_l)-1):\n",
    "    start,end = dates_l[i],dates_l[i+1]\n",
    "    start = datetime.strptime(start,\"%Y-%m-%d\")\n",
    "    end = datetime.strptime(end,\"%Y-%m-%d\")\n",
    "    end = end - timedelta(1)\n",
    "    dates_list.append((start.strftime(\"%Y-%m-%d\"),end.strftime(\"%Y-%m-%d\")))\n",
    "dates_list"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f8d784e9",
   "metadata": {},
   "source": [
    "### Creation of a hdf5 file to write in"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "3b2a365b",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<HDF5 dataset \"sar_bands\": shape (2,), type \"|S10\">"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ws_output = '/home/eouser/Desktop' \n",
    "hf = h5py.File(os.path.join(ws_output,'fs_{0}_{1}_{2}_raw5.h5'.format(pilot,year,buffer)), 'a')\n",
    "data = hf.create_group('data')\n",
    "sar_data = hf.create_group('sar_data')\n",
    "coords = hf.create_group('coords')\n",
    "meta = hf.create_group('metadata')\n",
    "meta.create_dataset('bands',(len(optical_bands)),'S10',[b.encode('ascii','ignore') for b in optical_bands])\n",
    "meta.create_dataset('sar_bands',(len(sar_bands)),'S10',[b.encode('ascii','ignore') for b in sar_bands])"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "99e22fff",
   "metadata": {},
   "source": [
    "### Iterate over time and get pixel data for each parcel"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1ccef831",
   "metadata": {},
   "source": [
    "The iteration over the times periods (set before) include the data retrieval and the zonal statistics generations. Instead of iterating over every parcel, we exploit the xarray of ids. Specifically, the command <code> grouped_data = data_cube.groupby(ids[buffer][0]) </code> generates data groups per id decreasing the time complexity and allowing aggregations functions on each group. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "ba997091",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      " 12%|█████▌                                      | 1/8 [04:10<29:12, 250.42s/it]"
     ]
    }
   ],
   "source": [
    "flag = True\n",
    "\n",
    "for d_n in tqdm(range(len(dates_list))):  \n",
    "    try:\n",
    "        d_start,d_end = dates_list[d_n][0],dates_list[d_n][1]\n",
    "        \n",
    "        # ---------------------------------------------------#\n",
    "        # ----Start of two approaches for loading data ------\n",
    "        # ---------------------------------------------------#\n",
    "        \n",
    "        ################## By using ADC ###########################\n",
    "        \n",
    "        # data_cube = getData_optical(bbox,d_start,d_end,optical_bands)\n",
    "        # data_cube = data_cube.load()\n",
    "\n",
    "         ########### By using datasets provided in git ###########\n",
    "\n",
    "        d_start = pd.to_datetime(d_start).to_numpy()\n",
    "        d_end = pd.to_datetime(d_end).to_numpy()\n",
    "        data_cube = data_nc.where((data_nc.time>=d_start)&(data_nc.time<=d_end),drop=True)\n",
    "        \n",
    "        # No need to change the two following lines:\n",
    "        all_bands = ['B02','B03','B04','B05','B06','B07','B08','B8A','B11','B12','SCL']\n",
    "        all_indices = ['ndvi','ndwi','ndmi','psri','savi','evi','dvi','rdvi','rvi','tvi','tcari','gi','vigreen','varigreen','gari','gdvi','sipi','wdrvi','gvmi','gcvi']\n",
    "        \n",
    "        bands = [b for b in optical_bands if b in all_bands]\n",
    "        indices = [b for b in optical_bands if b in all_indices]\n",
    "  \n",
    "        # calculating new indices and drop cloudy images based on a user-defined threshold\n",
    "        \n",
    "        for i,index in enumerate(optical_bands):\n",
    "            if index in indices:\n",
    "                data_cube[index] = calculate_index(data_cube,index)\n",
    "            data_cube[index] = cloud_data(data_cube,index)\n",
    "            if i == 0:\n",
    "                to_keep = data_cube[index].dropna(dim='time',thresh=0.25).time\n",
    "                data_cube = data_cube.sel(time=to_keep)\n",
    "        for b in all_bands:\n",
    "            if b not in bands:\n",
    "                data_cube = data_cube.drop(b)        \n",
    "        # ---------------------------------------------------#\n",
    "        # ------ End of two approaches for loading data ------\n",
    "        # ---------------------------------------------------#\n",
    "        ### important step: grouping of pixels per parcel based on id raster\n",
    "        grouped_data = data_cube.groupby(ids[buffer][0])\n",
    "    except Exception as e:\n",
    "        print(e)\n",
    "        continue\n",
    "    \n",
    "    dates = data_cube.time.values\n",
    "    dates = np.array([str(t).split('T')[0] for t in dates])\n",
    "    dates_unique = sorted(set(dates))  \n",
    "    \n",
    "    all_keys = []\n",
    "    for f in grouped_data:\n",
    "        key, parcel_data = f[0],f[1]\n",
    "        if key!=-1:\n",
    "            coords_all = list(zip(parcel_data.x.values,parcel_data.y.values))\n",
    "            parcel_data = np.array([parcel_data[b].values for b in optical_bands])\n",
    "\n",
    "            if len(dates)!=dates_unique:\n",
    "                vals = []\n",
    "                for b in range(parcel_data.shape[0]):\n",
    "                    df = pd.DataFrame(data=parcel_data[b],index=dates)\n",
    "                    vals.append(df.groupby(df.index).mean().values)\n",
    "                vals = np.array(vals)\n",
    "            else:\n",
    "                vals = parcel_data.copy()\n",
    "            if flag:\n",
    "                coords.create_dataset(str(key),data=np.array(coords_all), compression='gzip',compression_opts=9)\n",
    "                data.create_dataset(str(key),data=vals.astype('float64'),compression='gzip',\n",
    "                                    compression_opts=9,maxshape=(vals.shape[0],None,vals.shape[2]))\n",
    "            else:\n",
    "                data[str(key)].resize((data[str(key)].shape[1] + vals.shape[1]), axis=1)\n",
    "                data[str(key)][:,-vals.shape[1]:,:] = vals\n",
    "            all_keys.append(key)\n",
    "\n",
    "    if flag:\n",
    "        meta.create_dataset('unique_ids',(len(all_keys)),'S40',\n",
    "                            [noa_to_id[i] for i in np.array(sorted(all_keys))])\n",
    "        meta.create_dataset('dates',(len(dates_unique)),'S10',\n",
    "                            [d.encode('ascii','ignore') for d in dates_unique],maxshape=(None,))\n",
    "        flag = False\n",
    "    else:\n",
    "        meta['dates'].resize((meta['dates'].shape[0] + len(dates_unique)), axis=0)\n",
    "        meta['dates'][-len(dates_unique):] = dates_unique\n",
    "    "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "35b5c7b7",
   "metadata": {},
   "source": [
    "### The same for SAR data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "fd41df41",
   "metadata": {},
   "outputs": [],
   "source": [
    "# def getData_sar(bbox,timeStart,timeEnd,sar_bands,resolution = 10):\n",
    "    \n",
    "#     product_sar= 'sentinel1_sar'\n",
    "#     all_sar_bands = ['vv','vh']\n",
    "    \n",
    "#     bands = [b for b in sar_bands if b in all_sar_bands]\n",
    "#     if bbox is not None:\n",
    "#         xmin,xmax,ymin,ymax = bbox[0],bbox[1],bbox[2],bbox[3]\n",
    "#     query = {\n",
    "#         'time': (timeStart,timeEnd),\n",
    "#         'product': product_sar,\n",
    "#         'x':(xmin,xmax),\n",
    "#         'y':(ymin,ymax),\n",
    "#         'crs':'EPSG:3857'\n",
    "#     }\n",
    "#     dc = datacube.Datacube(app=\"test\", config=config)\n",
    "#     data = dc.load(**query,measurements=all_sar_bands,dask_chunks={})\n",
    "#     for b in all_sar_bands:\n",
    "#         if b not in bands:\n",
    "#             data = data.drop(b)\n",
    "#     return data\n",
    "\n",
    "# flag = True\n",
    "\n",
    "# if sar_bands:\n",
    "    \n",
    "#     for d_n in tqdm(range(len(dates_list))):\n",
    "        \n",
    "#         try:        \n",
    "#             d_start,d_end = dates_list[d_n][0],dates_list[d_n][1]\n",
    "#             data_cube = getData_sar(bbox,d_start,d_end,sar_bands)\n",
    "#             data_cube = data_cube.load()\n",
    "#             grouped_data = data_cube.groupby(ids[buffer][0])\n",
    "#         except:\n",
    "#             continue\n",
    "#         sar_dates = data_cube.time.values\n",
    "#         sar_dates = np.array([str(t).split('T')[0] for t in sar_dates])\n",
    "#         dates_unique = sorted(set(sar_dates))  \n",
    "\n",
    "\n",
    "#         for f in grouped_data:\n",
    "#             key, parcel_data = f[0],f[1]\n",
    "# #             if key in noa_ids_to_keep:\n",
    "#             if key!=-1:\n",
    "#                 parcel_data = np.array([parcel_data[b].values for b in sar_bands])\n",
    "\n",
    "#                 if len(sar_dates)!=dates_unique:\n",
    "#                     vals = []\n",
    "#                     for b in range(parcel_data.shape[0]):\n",
    "#                         df = pd.DataFrame(data=parcel_data[b],index=sar_dates)\n",
    "#                         vals.append(df.groupby(df.index).mean().values)\n",
    "#                     vals = np.array(vals)\n",
    "#                 else:\n",
    "#                     vals = parcel_data.copy()\n",
    "                    \n",
    "#                 if flag:\n",
    "#                     sar_data.create_dataset(str(key),data=vals.astype('float64'),compression='gzip',\n",
    "#                                         compression_opts=9,maxshape=(vals.shape[0],None,vals.shape[2]))\n",
    "#                 else:\n",
    "#                     sar_data[str(key)].resize((sar_data[str(key)].shape[1] + vals.shape[1]), axis=1)\n",
    "#                     sar_data[str(key)][:,-vals.shape[1]:,:] = vals\n",
    "\n",
    "\n",
    "#         dates_unique = sorted(set(sar_dates))   \n",
    "#         if flag:\n",
    "#             meta.create_dataset('sar_dates',(len(dates_unique)),'S10',\n",
    "#                                 [d.encode('ascii','ignore') for d in dates_unique],maxshape=(None,))\n",
    "#             flag = False\n",
    "#         else:\n",
    "#             meta['sar_dates'].resize((meta['sar_dates'].shape[0] + len(dates_unique)), axis=0)\n",
    "#             meta['sar_dates'][-len(dates_unique):] = dates_unique"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "8a13d717",
   "metadata": {},
   "outputs": [],
   "source": [
    "hf.close()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "0d12321d",
   "metadata": {},
   "source": [
    "## Filtering and Daily Linear Interpolation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "5c764c92",
   "metadata": {},
   "outputs": [],
   "source": [
    "# ### interpolation timestamps\n",
    "\n",
    "# def date_range(start, end, intv):\n",
    "\n",
    "#     start = datetime.strptime(start,\"%Y-%m-%d\")#+timedelta(days=31)\n",
    "#     end = datetime.strptime(end,\"%Y-%m-%d\")-timedelta(days=1)\n",
    "#     diff = (end  - start ) / intv\n",
    "#     for i in range(intv):\n",
    "#         yield (start + diff * i)\n",
    "#     yield end\n",
    "\n",
    "# s2_temporal_resolution = 5\n",
    "# s1_temporal_resolution = 6\n",
    "\n",
    "# start = datetime.strptime(start_date,\"%Y-%m-%d\")#+timedelta(days=31)\n",
    "# end = datetime.strptime(end_date,\"%Y-%m-%d\")-timedelta(days=1)\n",
    "# d_break = (end-start).days//s2_temporal_resolution\n",
    "# d_break_sar = (end-start).days//s1_temporal_resolution\n",
    "    \n",
    "# dates_interp = list(date_range(start_date,end_date,d_break))\n",
    "# dates_interp = [d.date() for d in dates_interp]\n",
    "# dates_interp_sar = list(date_range(start_date,end_date,d_break_sar))\n",
    "# dates_interp_sar = [d.date() for d in dates_interp_sar]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "cd44760b",
   "metadata": {},
   "outputs": [],
   "source": [
    "# def daily_interpolation(vals,dates_timestamp,sar=False,smoothing=False):\n",
    "\n",
    "#     df_band = pd.DataFrame(data=vals.T.copy(),columns=dates_timestamp)\n",
    "#     #     df_band[df_band==-9999] = np.nan\n",
    "#     start = datetime.strptime(start_date,\"%Y-%m-%d\")-timedelta(days=1)\n",
    "#     end = datetime.strptime(end_date,\"%Y-%m-%d\")+timedelta(days=1)\n",
    "#     df_band[start] = np.nan\n",
    "#     df_band[end] = np.nan\n",
    "#     df_band = df_band[sorted(df_band.columns)]\n",
    "#     df_band = df_band.resample('D',axis=1).mean().interpolate('linear',axis=1).bfill(axis=1).ffill(axis=1)\n",
    "#     if sar:\n",
    "#         df_band = df_band[dates_interp_sar]\n",
    "#     else:\n",
    "#         df_band = df_band[dates_interp]\n",
    "#     if smoothing:\n",
    "#         df_band = df_band.rolling(window=3,center=True,axis=1).median().bfill(axis=1).ffill(axis=1)\n",
    "    \n",
    "#     return df_band.values.T"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "ce981a26",
   "metadata": {},
   "outputs": [],
   "source": [
    "# hf = h5py.File('/home/eouser/Desktop/jason_notebooks/fs/{0}/fs_{0}_{1}_{2}_{3}_raw.h5'.format(pilot,region,year,buffer), \"r\", libver='latest', swmr=True)\n",
    "# unique_ids = np.array([str(i.decode('UTF-8')) for i in hf['metadata']['unique_ids'][:]])\n",
    "# dates = np.array([str(d.decode('UTF-8')) for d in hf['metadata']['dates'][:]])\n",
    "# dates_timestamp = np.array([datetime.strptime(x,\"%Y-%m-%d\") for x in dates])\n",
    "# sar_dates = np.array([str(d.decode('UTF-8')) for d in hf['metadata']['sar_dates'][:]])\n",
    "# sar_dates_timestamp = np.array([datetime.strptime(x,\"%Y-%m-%d\") for x in sar_dates])\n",
    "# bands = np.array([str(b.decode('UTF-8')) for b in hf['metadata']['bands'][:]])\n",
    "# sar_bands = np.array([str(b.decode('UTF-8')) for b in hf['metadata']['sar_bands'][:]])\n",
    "# ids = np.array(list(hf['data'].keys()))\n",
    "\n",
    "# # ids = sorted(set(ids.astype(int)).intersection(set(noa_ids_to_keep)))\n",
    "# ids = np.array(ids).astype(str)\n",
    "# unique_ids = np.array([noa_to_id[int(i)] for i in ids])\n",
    "\n",
    "# ndvi_i = np.where(bands=='ndvi')[0][0]\n",
    "# lower_ndvi_thresh = 0 \n",
    "# sample = True # put yes if you want to extract a sample of random pixels inside the parcel\n",
    "# # sample_size = 0.2 # the number the sample pixels > 1 or the portion <= 1\n",
    "# sample_size = 10\n",
    "\n",
    "\n",
    "# dates_new = [d.strftime(\"%Y-%m-%d\") for d in dates_interp]\n",
    "# sar_dates_new = [d.strftime(\"%Y-%m-%d\") for d in dates_interp_sar]\n",
    "\n",
    "\n",
    "# hf_interp = h5py.File('/home/eouser/Desktop/jason_notebooks/fs/{0}/fs_{0}_{1}_{2}_{3}_interp.h5'.format(pilot,region,year,buffer), 'w')\n",
    "# data_interp = hf_interp.create_group('data')\n",
    "# sar_data_interp = hf_interp.create_group('sar_data')\n",
    "# coords_interp = hf_interp.create_group('coords')\n",
    "# meta_interp = hf_interp.create_group('metadata')\n",
    "# meta_interp.create_dataset('bands',(len(bands)),'S10',[b.encode('ascii','ignore') for b in bands])\n",
    "# meta_interp.create_dataset('sar_bands',(len(sar_bands)),'S10',[b.encode('ascii','ignore') for b in sar_bands])\n",
    "# meta_interp.create_dataset('dates',(len(dates_new)),'S10',[d.encode('ascii','ignore') for d in dates_new])\n",
    "# meta_interp.create_dataset('sar_dates',(len(sar_dates_new)),'S10',[d.encode('ascii','ignore') for d in sar_dates_new])\n",
    "# meta_interp.create_dataset('unique_ids',(len(unique_ids)),'S40',[i.encode('ascii','ignore') for i in unique_ids])\n",
    "\n",
    "\n",
    "# for i in tqdm(ids):\n",
    "    \n",
    "#     vals = hf['data'][i][:]\n",
    "#     sar_vals = hf['sar_data'][i][:]\n",
    "#     coords = hf['coords'][i][:]\n",
    "    \n",
    "#     if sample:\n",
    "#         vals_size = vals.shape[-1]\n",
    "#         if sample_size <= 1:\n",
    "#             sample_size = int(vals_size*sample_size)\n",
    "#             s = np.sort(np.random.choice(np.arange(vals_size),sample_size,replace=False))\n",
    "#             vals = vals[:,:,s]\n",
    "#             sar_vals = sar_vals[:,:,s]\n",
    "#             coords = coords[s,:]\n",
    "#         else:\n",
    "#             vals_size = vals.shape[-1]\n",
    "#             if sample_size<=vals_size:\n",
    "#                 s = np.sort(np.random.choice(np.arange(vals_size),sample_size,replace=False))\n",
    "#             else:\n",
    "#                 s = np.sort(np.random.choice(np.arange(vals_size),sample_size,replace=True))\n",
    "#             vals = vals[:,:,s]\n",
    "#             sar_vals = sar_vals[:,:,s]\n",
    "#             coords = coords[s,:]\n",
    "    \n",
    "    \n",
    "#     vals[vals==-9999.] = np.nan\n",
    "#     sar_vals[(sar_vals==0.)|(sar_vals<=-30.)] = np.nan\n",
    "    \n",
    "#     ii = np.where(vals[ndvi_i,:,:]<=lower_ndvi_thresh) ## put nan for every band if ndvi<=ndvi_threshold\n",
    "    \n",
    "#     vals_new = np.zeros((vals.shape[0],len(dates_new),vals.shape[-1]))\n",
    "#     sar_vals_new = np.zeros((sar_vals.shape[0],len(sar_dates_new),sar_vals.shape[-1]))\n",
    "#     for b in range(len(bands)):\n",
    "#         vals[b,ii[0],ii[1]] = np.nan\n",
    "#         vals_new[b,:,:] = daily_interpolation(vals[b,:,:],dates_timestamp,sar=False)\n",
    "#     for b in range(len(sar_bands)):\n",
    "#         sar_vals_new[b,:,:] = daily_interpolation(sar_vals[b,:,:],sar_dates_timestamp,sar=True)\n",
    "\n",
    "#     coords_interp.create_dataset(str(i),data=np.array(coords), compression='gzip',compression_opts=9)\n",
    "#     data_interp.create_dataset(str(i),data=vals_new.astype('float16'),compression='gzip',compression_opts=9)\n",
    "#     sar_data_interp.create_dataset(str(i),data=sar_vals_new.astype('float16'),compression='gzip',compression_opts=9)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "8a2e0210",
   "metadata": {},
   "outputs": [],
   "source": [
    "# hf.close()\n",
    "# hf_interp.close()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "fc1b0507",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}