datasets/sentinel-2-l2a/baseline-change.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "b7d2e416-6fae-4e13-9de3-149f333c0048",
   "metadata": {},
   "source": [
    "## Adjusting for the Sentinel-2 Baseline Change\n",
    "\n",
    "On January 25th, 2022 the European Space Agency implemented a [change to the processing baseline](https://sentinels.copernicus.eu/web/sentinel/-/copernicus-sentinel-2-major-products-upgrade-upcoming) for Sentinel-2 data that affects both L1C and the L2A data hosted by the Planetary Computer.\n",
    "\n",
    "The crux of the change is an offset to the L1C image reflectance data, which propagates to the L2A data hosted by the Planetary Computer.\n",
    "\n",
    "As explained by ESA:\n",
    "\n",
    "> This evolution allows avoiding the loss of information due to clamping of negative values in the predefined range `[1-32767]` occurring over dark surfaces.\n",
    "\n",
    "The presence of this offset in newer data makes it inappropriate to (naively) compare data from before and after the baseline change. Trying to apply a machine learning model trained on the old baseline to (unadjusted) imagery from the new baseline will yield incorrect results. And trying to train a machine learning model on old and (unadjusted) new data will likely result in much worse performance.\n",
    "\n",
    "This notebook demonstrates how the change affects the Sentinel-2 L2A data hosted on the Planetary Computer and how to harmonize the new data to be comparable to the old data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "b30543ce-24e0-4f84-a046-d137c06ba896",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pystac_client\n",
    "import planetary_computer\n",
    "import xarray as xr\n",
    "import matplotlib.pyplot as plt\n",
    "import skimage.exposure\n",
    "import seaborn as sns\n",
    "import stackstac\n",
    "import datetime"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e4b0be43-3b54-4188-b2e6-6097838eade4",
   "metadata": {},
   "source": [
    "The `harmonize_to_old` function harmonizes new (post January 25th, 2022) Sentinel 2 data to be comparable to older data. It expects an xarray DataArray similar to those returned by `stackstac`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "89595090-7ead-4f8f-976a-602a226141fe",
   "metadata": {},
   "outputs": [],
   "source": [
    "def harmonize_to_old(data):\n",
    "    \"\"\"\n",
    "    Harmonize new Sentinel-2 data to the old baseline.\n",
    "\n",
    "    Parameters\n",
    "    ----------\n",
    "    data: xarray.DataArray\n",
    "        A DataArray with four dimensions: time, band, y, x\n",
    "\n",
    "    Returns\n",
    "    -------\n",
    "    harmonized: xarray.DataArray\n",
    "        A DataArray with all values harmonized to the old\n",
    "        processing baseline.\n",
    "    \"\"\"\n",
    "    cutoff = datetime.datetime(2022, 1, 25)\n",
    "    offset = 1000\n",
    "    bands = [\n",
    "        \"B01\",\n",
    "        \"B02\",\n",
    "        \"B03\",\n",
    "        \"B04\",\n",
    "        \"B05\",\n",
    "        \"B06\",\n",
    "        \"B07\",\n",
    "        \"B08\",\n",
    "        \"B8A\",\n",
    "        \"B09\",\n",
    "        \"B10\",\n",
    "        \"B11\",\n",
    "        \"B12\",\n",
    "    ]\n",
    "\n",
    "    old = data.sel(time=slice(cutoff))\n",
    "\n",
    "    to_process = list(set(bands) & set(data.band.data.tolist()))\n",
    "    new = data.sel(time=slice(cutoff, None)).drop_sel(band=to_process)\n",
    "\n",
    "    new_harmonized = data.sel(time=slice(cutoff, None), band=to_process).clip(offset)\n",
    "    new_harmonized -= offset\n",
    "\n",
    "    new = xr.concat([new, new_harmonized], \"band\").sel(band=data.band.data.tolist())\n",
    "    return xr.concat([old, new], dim=\"time\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ed2742fb-08bb-45c0-92d1-34e329eedae2",
   "metadata": {},
   "source": [
    "To demonstrate the problem, we'll load an \"old\" scene with the old processing baseline of `03.00`, and a \"new\" scene with a new processing baseline of `04.00`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "bd7b7501-c1d9-43fb-ab44-b6c4b6628bbc",
   "metadata": {},
   "outputs": [],
   "source": [
    "s2 = pystac_client.Client.open(\n",
    "    \"https://planetarycomputer.microsoft.com/api/stac/v1/\",\n",
    "    modifier=planetary_computer.sign_inplace,\n",
    ").get_collection(\"sentinel-2-l2a\")\n",
    "\n",
    "old = s2.get_item(\"S2B_MSIL2A_20220106T021339_R060_T50KPC_20220107T033435\")\n",
    "new = s2.get_item(\"S2A_MSIL2A_20230116T021341_R060_T50KPC_20230116T112621\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e66dd073-dc01-401b-9f28-ad042b5c3d0b",
   "metadata": {},
   "source": [
    "We'll stack these together, cropping out a small section of the image to speed up loading."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "9eeff95b-c7cc-457d-bc30-1fd6d55d8830",
   "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[data-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-index-preview {\n",
       "  grid-column: 2 / 5;\n",
       "  color: var(--xr-font-color2);\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",
       ".xr-index-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",
       ".xr-index-data-in:checked ~ .xr-index-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-index-name div,\n",
       ".xr-index-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data,\n",
       ".xr-index-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",
       ".xr-no-icon {\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.DataArray &#x27;stackstac-7f8b9fd36f84dee2d0c759174b259b65&#x27; (time: 2,\n",
       "                                                                band: 6,\n",
       "                                                                y: 1000, x: 1000)&gt;\n",
       "dask.array&lt;fetch_raster_window, shape=(2, 6, 1000, 1000), dtype=float64, chunksize=(1, 1, 1000, 1000), chunktype=numpy.ndarray&gt;\n",
       "Coordinates: (12/44)\n",
       "  * time                                     (time) datetime64[ns] 2022-01-06...\n",
       "    id                                       (time) &lt;U54 &#x27;S2B_MSIL2A_20220106...\n",
       "  * band                                     (band) &lt;U3 &#x27;B02&#x27; &#x27;B03&#x27; ... &#x27;AOT&#x27;\n",
       "  * x                                        (x) float64 6e+05 6e+05 ... 6.1e+05\n",
       "  * y                                        (y) float64 7.7e+06 ... 7.69e+06\n",
       "    s2:reflectance_conversion_factor         (time) float64 1.034 1.034\n",
       "    ...                                       ...\n",
       "    proj:bbox                                object {600000.0, 7690240.0, 780...\n",
       "    gsd                                      (band) float64 10.0 10.0 ... 10.0\n",
       "    common_name                              (band) object &#x27;blue&#x27; ... None\n",
       "    center_wavelength                        (band) object 0.49 0.56 ... None\n",
       "    full_width_half_max                      (band) object 0.098 0.045 ... None\n",
       "    epsg                                     int64 32750\n",
       "Attributes:\n",
       "    spec:        RasterSpec(epsg=32750, bounds=(600000.0, 7690240.0, 610000.0...\n",
       "    crs:         epsg:32750\n",
       "    transform:   | 10.00, 0.00, 600000.00|\\n| 0.00,-10.00, 7700240.00|\\n| 0.0...\n",
       "    resolution:  10.0</pre><div class='xr-wrap' style='display:none'><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'>'stackstac-7f8b9fd36f84dee2d0c759174b259b65'</div><ul class='xr-dim-list'><li><span class='xr-has-index'>time</span>: 2</li><li><span class='xr-has-index'>band</span>: 6</li><li><span class='xr-has-index'>y</span>: 1000</li><li><span class='xr-has-index'>x</span>: 1000</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-011e302e-0eed-4ddc-b7fd-178d08cf48a7' class='xr-array-in' type='checkbox' checked><label for='section-011e302e-0eed-4ddc-b7fd-178d08cf48a7' title='Show/hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-array-preview xr-preview'><span>dask.array&lt;chunksize=(1, 1, 1000, 1000), meta=np.ndarray&gt;</span></div><div class='xr-array-data'><table>\n",
       "    <tr>\n",
       "        <td>\n",
       "            <table style=\"border-collapse: collapse;\">\n",
       "                <thead>\n",
       "                    <tr>\n",
       "                        <td> </td>\n",
       "                        <th> Array </th>\n",
       "                        <th> Chunk </th>\n",
       "                    </tr>\n",
       "                </thead>\n",
       "                <tbody>\n",
       "                    \n",
       "                    <tr>\n",
       "                        <th> Bytes </th>\n",
       "                        <td> 91.55 MiB </td>\n",
       "                        <td> 7.63 MiB </td>\n",
       "                    </tr>\n",
       "                    \n",
       "                    <tr>\n",
       "                        <th> Shape </th>\n",
       "                        <td> (2, 6, 1000, 1000) </td>\n",
       "                        <td> (1, 1, 1000, 1000) </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <th> Dask graph </th>\n",
       "                        <td colspan=\"2\"> 12 chunks in 3 graph layers </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <th> Data type </th>\n",
       "                        <td colspan=\"2\"> float64 numpy.ndarray </td>\n",
       "                    </tr>\n",
       "                </tbody>\n",
       "            </table>\n",
       "        </td>\n",
       "        <td>\n",
       "        <svg width=\"374\" height=\"184\" style=\"stroke:rgb(0,0,0);stroke-width:1\" >\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"0\" y1=\"0\" x2=\"25\" y2=\"0\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"0\" y1=\"25\" x2=\"25\" y2=\"25\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"0\" y1=\"0\" x2=\"0\" y2=\"25\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"12\" y1=\"0\" x2=\"12\" y2=\"25\" />\n",
       "  <line x1=\"25\" y1=\"0\" x2=\"25\" y2=\"25\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"0.0,0.0 25.412616514582485,0.0 25.412616514582485,25.412616514582485 0.0,25.412616514582485\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Text -->\n",
       "  <text x=\"12.706308\" y=\"45.412617\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" >2</text>\n",
       "  <text x=\"45.412617\" y=\"12.706308\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(0,45.412617,12.706308)\">1</text>\n",
       "\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"95\" y1=\"0\" x2=\"109\" y2=\"14\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"95\" y1=\"120\" x2=\"109\" y2=\"134\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"95\" y1=\"0\" x2=\"95\" y2=\"120\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"97\" y1=\"2\" x2=\"97\" y2=\"122\" />\n",
       "  <line x1=\"99\" y1=\"4\" x2=\"99\" y2=\"124\" />\n",
       "  <line x1=\"102\" y1=\"7\" x2=\"102\" y2=\"127\" />\n",
       "  <line x1=\"104\" y1=\"9\" x2=\"104\" y2=\"129\" />\n",
       "  <line x1=\"107\" y1=\"12\" x2=\"107\" y2=\"132\" />\n",
       "  <line x1=\"109\" y1=\"14\" x2=\"109\" y2=\"134\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"95.0,0.0 109.9485979497544,14.948597949754403 109.9485979497544,134.9485979497544 95.0,120.0\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"95\" y1=\"0\" x2=\"215\" y2=\"0\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"97\" y1=\"2\" x2=\"217\" y2=\"2\" />\n",
       "  <line x1=\"99\" y1=\"4\" x2=\"219\" y2=\"4\" />\n",
       "  <line x1=\"102\" y1=\"7\" x2=\"222\" y2=\"7\" />\n",
       "  <line x1=\"104\" y1=\"9\" x2=\"224\" y2=\"9\" />\n",
       "  <line x1=\"107\" y1=\"12\" x2=\"227\" y2=\"12\" />\n",
       "  <line x1=\"109\" y1=\"14\" x2=\"229\" y2=\"14\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"95\" y1=\"0\" x2=\"109\" y2=\"14\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"215\" y1=\"0\" x2=\"229\" y2=\"14\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"95.0,0.0 215.0,0.0 229.9485979497544,14.948597949754403 109.9485979497544,14.948597949754403\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"109\" y1=\"14\" x2=\"229\" y2=\"14\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"109\" y1=\"134\" x2=\"229\" y2=\"134\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"109\" y1=\"14\" x2=\"109\" y2=\"134\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"229\" y1=\"14\" x2=\"229\" y2=\"134\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"109.9485979497544,14.948597949754403 229.9485979497544,14.948597949754403 229.9485979497544,134.9485979497544 109.9485979497544,134.9485979497544\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Text -->\n",
       "  <text x=\"169.948598\" y=\"154.948598\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" >1000</text>\n",
       "  <text x=\"249.948598\" y=\"74.948598\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(-90,249.948598,74.948598)\">1000</text>\n",
       "  <text x=\"92.474299\" y=\"147.474299\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(45,92.474299,147.474299)\">6</text>\n",
       "</svg>\n",
       "        </td>\n",
       "    </tr>\n",
       "</table></div></div></li><li class='xr-section-item'><input id='section-29306d52-dd52-4f32-9584-f6e051732a56' class='xr-section-summary-in' type='checkbox'  ><label for='section-29306d52-dd52-4f32-9584-f6e051732a56' class='xr-section-summary' >Coordinates: <span>(44)</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'>2022-01-06T02:13:39.024000 2023-...</div><input id='attrs-377e6416-1680-45cb-9ca9-9db65e908877' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-377e6416-1680-45cb-9ca9-9db65e908877' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f8bc5c15-277f-45ea-bbe1-b953bbeafc78' class='xr-var-data-in' type='checkbox'><label for='data-f8bc5c15-277f-45ea-bbe1-b953bbeafc78' 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;2022-01-06T02:13:39.024000000&#x27;, &#x27;2023-01-16T02:13:41.024000000&#x27;],\n",
       "      dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>id</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>&lt;U54</div><div class='xr-var-preview xr-preview'>&#x27;S2B_MSIL2A_20220106T021339_R060...</div><input id='attrs-694c160a-f9d5-40a7-ba37-18ad40f0fa79' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-694c160a-f9d5-40a7-ba37-18ad40f0fa79' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-d2c4b52f-dbee-46a5-aef7-6c9037f7aa0f' class='xr-var-data-in' type='checkbox'><label for='data-d2c4b52f-dbee-46a5-aef7-6c9037f7aa0f' 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;S2B_MSIL2A_20220106T021339_R060_T50KPC_20220107T033435&#x27;,\n",
       "       &#x27;S2A_MSIL2A_20230116T021341_R060_T50KPC_20230116T112621&#x27;],\n",
       "      dtype=&#x27;&lt;U54&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>band</span></div><div class='xr-var-dims'>(band)</div><div class='xr-var-dtype'>&lt;U3</div><div class='xr-var-preview xr-preview'>&#x27;B02&#x27; &#x27;B03&#x27; &#x27;B04&#x27; &#x27;B05&#x27; &#x27;B06&#x27; &#x27;AOT&#x27;</div><input id='attrs-600f8e3a-beb5-4bfb-94b5-21e108f7a22b' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-600f8e3a-beb5-4bfb-94b5-21e108f7a22b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-eff7e96f-8f7e-424f-97bf-6752c913a868' class='xr-var-data-in' type='checkbox'><label for='data-eff7e96f-8f7e-424f-97bf-6752c913a868' 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;B02&#x27;, &#x27;B03&#x27;, &#x27;B04&#x27;, &#x27;B05&#x27;, &#x27;B06&#x27;, &#x27;AOT&#x27;], dtype=&#x27;&lt;U3&#x27;)</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'>6e+05 6e+05 ... 6.1e+05 6.1e+05</div><input id='attrs-ab6f4d9c-54b1-4f20-b4dc-038609e0a2c4' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-ab6f4d9c-54b1-4f20-b4dc-038609e0a2c4' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-558b62bb-75ed-4aaf-847f-57ebb0f8b67a' class='xr-var-data-in' type='checkbox'><label for='data-558b62bb-75ed-4aaf-847f-57ebb0f8b67a' 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([600000., 600010., 600020., ..., 609970., 609980., 609990.])</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'>7.7e+06 7.7e+06 ... 7.69e+06</div><input id='attrs-bc129dba-28da-4f8e-951a-b4b57b82cb72' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-bc129dba-28da-4f8e-951a-b4b57b82cb72' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-359732ab-ea17-48a6-a106-d7f67895d204' class='xr-var-data-in' type='checkbox'><label for='data-359732ab-ea17-48a6-a106-d7f67895d204' 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([7700240., 7700230., 7700220., ..., 7690270., 7690260., 7690250.])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:reflectance_conversion_factor</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>1.034 1.034</div><input id='attrs-cb76c8e5-bb47-497d-a1aa-c5f89ae7dbcd' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-cb76c8e5-bb47-497d-a1aa-c5f89ae7dbcd' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-45369d11-bf37-4ab2-94ac-0f9ef9919b56' class='xr-var-data-in' type='checkbox'><label for='data-45369d11-bf37-4ab2-94ac-0f9ef9919b56' 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([1.03431236, 1.03393585])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>sat:orbit_state</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U10</div><div class='xr-var-preview xr-preview'>&#x27;descending&#x27;</div><input id='attrs-45811fcb-63eb-4772-ba82-bda7fb2a3784' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-45811fcb-63eb-4772-ba82-bda7fb2a3784' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5df0a74a-1f77-4e3c-927e-2c56c758bc82' class='xr-var-data-in' type='checkbox'><label for='data-5df0a74a-1f77-4e3c-927e-2c56c758bc82' 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;descending&#x27;, dtype=&#x27;&lt;U10&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:dark_features_percentage</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.118 0.00577</div><input id='attrs-6b361c02-c2f3-411b-9f9e-d7c67be13604' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-6b361c02-c2f3-411b-9f9e-d7c67be13604' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4e199310-1cce-4c0d-9ac9-6483c856b5b9' class='xr-var-data-in' type='checkbox'><label for='data-4e199310-1cce-4c0d-9ac9-6483c856b5b9' 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([0.117992, 0.00577 ])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:processing_baseline</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>&lt;U5</div><div class='xr-var-preview xr-preview'>&#x27;03.00&#x27; &#x27;04.00&#x27;</div><input id='attrs-e11e7fd4-379b-4a06-95f8-c051c35c1773' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-e11e7fd4-379b-4a06-95f8-c051c35c1773' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-fecb4b3e-f9db-4365-a880-7f577047d552' class='xr-var-data-in' type='checkbox'><label for='data-fecb4b3e-f9db-4365-a880-7f577047d552' 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;03.00&#x27;, &#x27;04.00&#x27;], dtype=&#x27;&lt;U5&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:mean_solar_azimuth</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>99.76 96.38</div><input id='attrs-407d9588-90b8-43a7-8fe8-a81e5cecfea7' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-407d9588-90b8-43a7-8fe8-a81e5cecfea7' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-efe5c190-8bbe-4463-b1d7-84d834f704a0' class='xr-var-data-in' type='checkbox'><label for='data-efe5c190-8bbe-4463-b1d7-84d834f704a0' 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([99.75935412, 96.37513022])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>eo:cloud_cover</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.03192 1.856</div><input id='attrs-2ec206a3-a8d6-42e2-89ac-34891e357caa' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-2ec206a3-a8d6-42e2-89ac-34891e357caa' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-551de04a-4d6e-4e3d-952b-6992021f4a39' class='xr-var-data-in' type='checkbox'><label for='data-551de04a-4d6e-4e3d-952b-6992021f4a39' 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([0.031924, 1.855648])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:thin_cirrus_percentage</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.0 1.648</div><input id='attrs-31858587-aa60-401b-8b18-b11f523a34c8' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-31858587-aa60-401b-8b18-b11f523a34c8' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-93ceeebe-1e1c-44b2-aaf2-1a4c88d87c10' class='xr-var-data-in' type='checkbox'><label for='data-93ceeebe-1e1c-44b2-aaf2-1a4c88d87c10' 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([0.      , 1.647632])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:nodata_pixel_percentage</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.0 7e-06</div><input id='attrs-8bb2dd81-ec4c-440d-9a4b-f9b18cede0a7' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-8bb2dd81-ec4c-440d-9a4b-f9b18cede0a7' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f50ccf3d-17d6-4ba0-8537-74751fd3d065' class='xr-var-data-in' type='checkbox'><label for='data-f50ccf3d-17d6-4ba0-8537-74751fd3d065' 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([0.e+00, 7.e-06])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>instruments</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U3</div><div class='xr-var-preview xr-preview'>&#x27;msi&#x27;</div><input id='attrs-29f3a574-e6ae-4e10-b08f-124cad50e19a' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-29f3a574-e6ae-4e10-b08f-124cad50e19a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-35a66c5b-3409-4c2e-a637-3d85e72358f7' class='xr-var-data-in' type='checkbox'><label for='data-35a66c5b-3409-4c2e-a637-3d85e72358f7' 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;msi&#x27;, dtype=&#x27;&lt;U3&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:not_vegetated_percentage</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>54.74 55.23</div><input id='attrs-fa1a22a6-dac7-476d-a8c8-7452b2078cea' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-fa1a22a6-dac7-476d-a8c8-7452b2078cea' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-71b579c1-a37a-40ac-8754-9ce458008104' class='xr-var-data-in' type='checkbox'><label for='data-71b579c1-a37a-40ac-8754-9ce458008104' 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([54.737449, 55.233008])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:vegetation_percentage</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.386 0.3828</div><input id='attrs-bcf78b8d-6706-44c1-8942-f9dc14a2afd7' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-bcf78b8d-6706-44c1-8942-f9dc14a2afd7' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f59fc2d9-c9b0-45a7-a2fe-15fae038ffba' class='xr-var-data-in' type='checkbox'><label for='data-f59fc2d9-c9b0-45a7-a2fe-15fae038ffba' 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([0.385951, 0.382849])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:high_proba_clouds_percentage</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.009021 0.06926</div><input id='attrs-ef6f8450-ed03-4623-8200-8d98f732338f' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-ef6f8450-ed03-4623-8200-8d98f732338f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-eebe26a7-ca4e-446c-839f-f2a0dc4be81d' class='xr-var-data-in' type='checkbox'><label for='data-eebe26a7-ca4e-446c-839f-f2a0dc4be81d' 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([0.009021, 0.069263])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:product_uri</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>&lt;U65</div><div class='xr-var-preview xr-preview'>&#x27;S2B_MSIL2A_20220106T021339_N030...</div><input id='attrs-52ad32c3-3b79-4631-b339-538b613cc4b4' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-52ad32c3-3b79-4631-b339-538b613cc4b4' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-39a630ab-aec0-4bc0-bc34-7b612d0ec006' class='xr-var-data-in' type='checkbox'><label for='data-39a630ab-aec0-4bc0-bc34-7b612d0ec006' 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;S2B_MSIL2A_20220106T021339_N0300_R060_T50KPC_20220107T033435.SAFE&#x27;,\n",
       "       &#x27;S2A_MSIL2A_20230116T021341_N0400_R060_T50KPC_20230116T112621.SAFE&#x27;],\n",
       "      dtype=&#x27;&lt;U65&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:cloud_shadow_percentage</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.01462 0.2521</div><input id='attrs-9fc40dbe-a69c-4d71-aa1a-71f51984ed7a' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-9fc40dbe-a69c-4d71-aa1a-71f51984ed7a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c099a3b0-f984-4528-80b1-575192240779' class='xr-var-data-in' type='checkbox'><label for='data-c099a3b0-f984-4528-80b1-575192240779' 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([0.014618, 0.252099])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:mgrs_tile</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U5</div><div class='xr-var-preview xr-preview'>&#x27;50KPC&#x27;</div><input id='attrs-dd9c681e-686b-4701-8289-68c6fdc121d6' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-dd9c681e-686b-4701-8289-68c6fdc121d6' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-da9f2ddf-a046-4955-bdb1-9315e4bda3f5' class='xr-var-data-in' type='checkbox'><label for='data-da9f2ddf-a046-4955-bdb1-9315e4bda3f5' 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;50KPC&#x27;, dtype=&#x27;&lt;U5&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>constellation</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U10</div><div class='xr-var-preview xr-preview'>&#x27;Sentinel 2&#x27;</div><input id='attrs-bea99846-4412-4d2e-bbdb-e99a9338dde0' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-bea99846-4412-4d2e-bbdb-e99a9338dde0' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5f4752c5-efa4-4fa0-8abf-58dd7a6275ec' class='xr-var-data-in' type='checkbox'><label for='data-5f4752c5-efa4-4fa0-8abf-58dd7a6275ec' 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;Sentinel 2&#x27;, dtype=&#x27;&lt;U10&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:water_percentage</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>43.82 42.21</div><input id='attrs-636e7ead-578e-4c3e-b57a-005e1aa20539' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-636e7ead-578e-4c3e-b57a-005e1aa20539' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-11f13022-c552-46d7-9345-31d914b71ae0' class='xr-var-data-in' type='checkbox'><label for='data-11f13022-c552-46d7-9345-31d914b71ae0' 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([43.815303, 42.207554])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>sat:relative_orbit</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>60</div><input id='attrs-f7ee7019-682a-4d09-9a5e-9312fd9e9557' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-f7ee7019-682a-4d09-9a5e-9312fd9e9557' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b2dd500a-2f46-4b50-adad-80c777016e47' class='xr-var-data-in' type='checkbox'><label for='data-b2dd500a-2f46-4b50-adad-80c777016e47' 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(60)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:snow_ice_percentage</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.0</div><input id='attrs-a31a3157-a605-4ffe-8230-f080815fe415' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-a31a3157-a605-4ffe-8230-f080815fe415' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-bf895549-0c08-43d5-b0ee-ca51e091d1f9' class='xr-var-data-in' type='checkbox'><label for='data-bf895549-0c08-43d5-b0ee-ca51e091d1f9' 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(0.)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:datastrip_id</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>&lt;U64</div><div class='xr-var-preview xr-preview'>&#x27;S2B_OPER_MSI_L2A_DS_ESRI_202201...</div><input id='attrs-b7c4f324-5fa8-49cf-a49a-a0accb1a85cb' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-b7c4f324-5fa8-49cf-a49a-a0accb1a85cb' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-84148726-2188-4064-8cdf-4c0d2c316518' class='xr-var-data-in' type='checkbox'><label for='data-84148726-2188-4064-8cdf-4c0d2c316518' 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;S2B_OPER_MSI_L2A_DS_ESRI_20220107T033435_S20220106T022252_N03.00&#x27;,\n",
       "       &#x27;S2A_OPER_MSI_L2A_DS_ESRI_20230116T112623_S20230116T021431_N04.00&#x27;],\n",
       "      dtype=&#x27;&lt;U64&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:medium_proba_clouds_percentage</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.0229 0.1388</div><input id='attrs-d6842322-2569-4af5-8dbf-b179713b21ae' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-d6842322-2569-4af5-8dbf-b179713b21ae' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-440149c6-02fd-4519-9591-acc930a17d53' class='xr-var-data-in' type='checkbox'><label for='data-440149c6-02fd-4519-9591-acc930a17d53' 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([0.022903, 0.138752])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:product_type</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U7</div><div class='xr-var-preview xr-preview'>&#x27;S2MSI2A&#x27;</div><input id='attrs-5b718628-fef0-4f95-809d-800972f179ed' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-5b718628-fef0-4f95-809d-800972f179ed' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8e97baed-8b55-477a-8f26-aa41f8285472' class='xr-var-data-in' type='checkbox'><label for='data-8e97baed-8b55-477a-8f26-aa41f8285472' 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;S2MSI2A&#x27;, dtype=&#x27;&lt;U7&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:mean_solar_zenith</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>25.28 26.23</div><input id='attrs-448e61ed-d030-4112-9ef1-ac1952ead860' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-448e61ed-d030-4112-9ef1-ac1952ead860' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-1643d989-bb78-4885-9837-9bb27fb97242' class='xr-var-data-in' type='checkbox'><label for='data-1643d989-bb78-4885-9837-9bb27fb97242' 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([25.28209086, 26.22624432])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:saturated_defective_pixel_percentage</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.0</div><input id='attrs-051c146c-3a99-4f8b-9b73-cb6f9689454c' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-051c146c-3a99-4f8b-9b73-cb6f9689454c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-d51fc26c-cdb0-427e-ae45-a4c1f6af3f58' class='xr-var-data-in' type='checkbox'><label for='data-d51fc26c-cdb0-427e-ae45-a4c1f6af3f58' 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(0.)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>proj:epsg</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>32750</div><input id='attrs-37ce0118-3e6d-43a4-b9a0-a8cca24f4cb1' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-37ce0118-3e6d-43a4-b9a0-a8cca24f4cb1' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-a0886159-fb26-4f42-80e4-f4e504be0780' class='xr-var-data-in' type='checkbox'><label for='data-a0886159-fb26-4f42-80e4-f4e504be0780' 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(32750)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:datatake_type</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U8</div><div class='xr-var-preview xr-preview'>&#x27;INS-NOBS&#x27;</div><input id='attrs-8c84f7bd-3817-4bb7-8782-a033752b4bc1' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-8c84f7bd-3817-4bb7-8782-a033752b4bc1' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-efe9d97a-ee14-4cbc-bb9f-1cd7e86635a7' class='xr-var-data-in' type='checkbox'><label for='data-efe9d97a-ee14-4cbc-bb9f-1cd7e86635a7' 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;INS-NOBS&#x27;, dtype=&#x27;&lt;U8&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:datatake_id</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>&lt;U34</div><div class='xr-var-preview xr-preview'>&#x27;GS2B_20220106T021339_025254_N03...</div><input id='attrs-8147a88c-b9c0-4062-9f34-3cf0c99d22fc' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-8147a88c-b9c0-4062-9f34-3cf0c99d22fc' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-fd2e678f-9add-4ddc-b9e5-757d06a1a142' class='xr-var-data-in' type='checkbox'><label for='data-fd2e678f-9add-4ddc-b9e5-757d06a1a142' 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;GS2B_20220106T021339_025254_N03.00&#x27;,\n",
       "       &#x27;GS2A_20230116T021341_039525_N04.00&#x27;], dtype=&#x27;&lt;U34&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:unclassified_percentage</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.8968 0.06308</div><input id='attrs-79725a47-1486-4d10-ae73-6433ff5f9bc8' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-79725a47-1486-4d10-ae73-6433ff5f9bc8' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-36d883ec-b35a-4db0-83ce-ca6042da198d' class='xr-var-data-in' type='checkbox'><label for='data-36d883ec-b35a-4db0-83ce-ca6042da198d' 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([0.896762, 0.063075])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:degraded_msi_data_percentage</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.0</div><input id='attrs-0248fe22-7187-4e63-9d74-e70da0807a30' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-0248fe22-7187-4e63-9d74-e70da0807a30' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-d85e371d-6280-4926-96f1-fa8920028941' class='xr-var-data-in' type='checkbox'><label for='data-d85e371d-6280-4926-96f1-fa8920028941' 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(0.)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:generation_time</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>&lt;U27</div><div class='xr-var-preview xr-preview'>&#x27;2022-01-07T03:34:35.286504Z&#x27; &#x27;2...</div><input id='attrs-5194fa80-5da4-400a-a986-ef6fa44266b6' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-5194fa80-5da4-400a-a986-ef6fa44266b6' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-1354365c-db19-4de0-9e92-ca77a9bb0567' class='xr-var-data-in' type='checkbox'><label for='data-1354365c-db19-4de0-9e92-ca77a9bb0567' 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;2022-01-07T03:34:35.286504Z&#x27;, &#x27;2023-01-16T11:26:21.757846Z&#x27;],\n",
       "      dtype=&#x27;&lt;U27&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>s2:granule_id</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>&lt;U62</div><div class='xr-var-preview xr-preview'>&#x27;S2B_OPER_MSI_L2A_TL_ESRI_202201...</div><input id='attrs-f74a0119-bc59-4cfc-b6e3-89c82c28b784' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-f74a0119-bc59-4cfc-b6e3-89c82c28b784' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5f7225cf-7a57-4a7d-b1d8-c664dfe3d04b' class='xr-var-data-in' type='checkbox'><label for='data-5f7225cf-7a57-4a7d-b1d8-c664dfe3d04b' 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;S2B_OPER_MSI_L2A_TL_ESRI_20220107T033435_A025254_T50KPC_N03.00&#x27;,\n",
       "       &#x27;S2A_OPER_MSI_L2A_TL_ESRI_20230116T112623_A039525_T50KPC_N04.00&#x27;],\n",
       "      dtype=&#x27;&lt;U62&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>platform</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>&lt;U11</div><div class='xr-var-preview xr-preview'>&#x27;Sentinel-2B&#x27; &#x27;Sentinel-2A&#x27;</div><input id='attrs-e373675a-4481-469b-9e16-d93189ca962c' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-e373675a-4481-469b-9e16-d93189ca962c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-6d25d7c5-8092-4999-9234-a2fcd084507d' class='xr-var-data-in' type='checkbox'><label for='data-6d25d7c5-8092-4999-9234-a2fcd084507d' 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;Sentinel-2B&#x27;, &#x27;Sentinel-2A&#x27;], dtype=&#x27;&lt;U11&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>title</span></div><div class='xr-var-dims'>(band)</div><div class='xr-var-dtype'>&lt;U36</div><div class='xr-var-preview xr-preview'>&#x27;Band 2 - Blue - 10m&#x27; ... &#x27;Aeros...</div><input id='attrs-de6437fc-ab95-4b60-b95b-76e7e46f5101' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-de6437fc-ab95-4b60-b95b-76e7e46f5101' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-0e927e73-ca70-4afa-bd1b-fcad1cdef34b' class='xr-var-data-in' type='checkbox'><label for='data-0e927e73-ca70-4afa-bd1b-fcad1cdef34b' 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;Band 2 - Blue - 10m&#x27;, &#x27;Band 3 - Green - 10m&#x27;,\n",
       "       &#x27;Band 4 - Red - 10m&#x27;, &#x27;Band 5 - Vegetation red edge 1 - 20m&#x27;,\n",
       "       &#x27;Band 6 - Vegetation red edge 2 - 20m&#x27;,\n",
       "       &#x27;Aerosol optical thickness (AOT)&#x27;], dtype=&#x27;&lt;U36&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>proj:bbox</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>{600000.0, 7690240.0, 7800040.0,...</div><input id='attrs-d49979ec-d128-43fd-a810-a533c7efa9bb' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-d49979ec-d128-43fd-a810-a533c7efa9bb' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-06a7c337-43f2-48b9-9984-64204a3e0242' class='xr-var-data-in' type='checkbox'><label for='data-06a7c337-43f2-48b9-9984-64204a3e0242' 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({600000.0, 7690240.0, 7800040.0, 709800.0}, dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>gsd</span></div><div class='xr-var-dims'>(band)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>10.0 10.0 10.0 20.0 20.0 10.0</div><input id='attrs-a26d156a-cf20-4f0e-9798-c6de67a9cf18' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-a26d156a-cf20-4f0e-9798-c6de67a9cf18' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e400a70f-7f5b-43f2-b9e8-b48872747fec' class='xr-var-data-in' type='checkbox'><label for='data-e400a70f-7f5b-43f2-b9e8-b48872747fec' 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([10., 10., 10., 20., 20., 10.])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>common_name</span></div><div class='xr-var-dims'>(band)</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>&#x27;blue&#x27; &#x27;green&#x27; ... &#x27;rededge&#x27; None</div><input id='attrs-6172d3d6-7121-4b7b-a9ba-0a95d5833709' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-6172d3d6-7121-4b7b-a9ba-0a95d5833709' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5e2fcbcb-3eab-498c-b5e2-2e25601bea6e' class='xr-var-data-in' type='checkbox'><label for='data-5e2fcbcb-3eab-498c-b5e2-2e25601bea6e' 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;blue&#x27;, &#x27;green&#x27;, &#x27;red&#x27;, &#x27;rededge&#x27;, &#x27;rededge&#x27;, None], dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>center_wavelength</span></div><div class='xr-var-dims'>(band)</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>0.49 0.56 0.665 0.704 0.74 None</div><input id='attrs-1d30b895-654c-43d2-8fd8-fd09b8662e80' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-1d30b895-654c-43d2-8fd8-fd09b8662e80' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-37abddc0-2e24-43e2-a553-ad5154cc3223' class='xr-var-data-in' type='checkbox'><label for='data-37abddc0-2e24-43e2-a553-ad5154cc3223' 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([0.49, 0.56, 0.665, 0.704, 0.74, None], dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>full_width_half_max</span></div><div class='xr-var-dims'>(band)</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>0.098 0.045 0.038 0.019 0.018 None</div><input id='attrs-90976822-7070-4289-8a3e-e18e8191fabe' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-90976822-7070-4289-8a3e-e18e8191fabe' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-13fe2a13-0d19-4d22-9dd0-6829dad0a9b6' class='xr-var-data-in' type='checkbox'><label for='data-13fe2a13-0d19-4d22-9dd0-6829dad0a9b6' 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([0.098, 0.045, 0.038, 0.019, 0.018, None], dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>epsg</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>32750</div><input id='attrs-701bd12a-e402-4bf7-81c0-aada641cea2f' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-701bd12a-e402-4bf7-81c0-aada641cea2f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b7b696fe-d161-4087-a8da-61cba971082e' class='xr-var-data-in' type='checkbox'><label for='data-b7b696fe-d161-4087-a8da-61cba971082e' 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(32750)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-646dc75d-bb5d-4262-b890-e535e0ce0860' class='xr-section-summary-in' type='checkbox'  ><label for='section-646dc75d-bb5d-4262-b890-e535e0ce0860' class='xr-section-summary' >Indexes: <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-index-name'><div>time</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-07a9ab6d-afef-40b0-a7e7-a90d9352eb3d' class='xr-index-data-in' type='checkbox'/><label for='index-07a9ab6d-afef-40b0-a7e7-a90d9352eb3d' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(DatetimeIndex([&#x27;2022-01-06 02:13:39.024000&#x27;, &#x27;2023-01-16 02:13:41.024000&#x27;], dtype=&#x27;datetime64[ns]&#x27;, name=&#x27;time&#x27;, freq=None))</pre></div></li><li class='xr-var-item'><div class='xr-index-name'><div>band</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-6aea56cd-0bdd-45a5-9e30-bd0e7a96658b' class='xr-index-data-in' type='checkbox'/><label for='index-6aea56cd-0bdd-45a5-9e30-bd0e7a96658b' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Index([&#x27;B02&#x27;, &#x27;B03&#x27;, &#x27;B04&#x27;, &#x27;B05&#x27;, &#x27;B06&#x27;, &#x27;AOT&#x27;], dtype=&#x27;object&#x27;, name=&#x27;band&#x27;))</pre></div></li><li class='xr-var-item'><div class='xr-index-name'><div>x</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-82bdbca6-228d-481c-b745-2fbe82f9d9cf' class='xr-index-data-in' type='checkbox'/><label for='index-82bdbca6-228d-481c-b745-2fbe82f9d9cf' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Float64Index([600000.0, 600010.0, 600020.0, 600030.0, 600040.0, 600050.0,\n",
       "              600060.0, 600070.0, 600080.0, 600090.0,\n",
       "              ...\n",
       "              609900.0, 609910.0, 609920.0, 609930.0, 609940.0, 609950.0,\n",
       "              609960.0, 609970.0, 609980.0, 609990.0],\n",
       "             dtype=&#x27;float64&#x27;, name=&#x27;x&#x27;, length=1000))</pre></div></li><li class='xr-var-item'><div class='xr-index-name'><div>y</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-c009919d-3fa6-4fcb-bf78-93b5da3f10dd' class='xr-index-data-in' type='checkbox'/><label for='index-c009919d-3fa6-4fcb-bf78-93b5da3f10dd' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Float64Index([7700240.0, 7700230.0, 7700220.0, 7700210.0, 7700200.0, 7700190.0,\n",
       "              7700180.0, 7700170.0, 7700160.0, 7700150.0,\n",
       "              ...\n",
       "              7690340.0, 7690330.0, 7690320.0, 7690310.0, 7690300.0, 7690290.0,\n",
       "              7690280.0, 7690270.0, 7690260.0, 7690250.0],\n",
       "             dtype=&#x27;float64&#x27;, name=&#x27;y&#x27;, length=1000))</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-0893303d-aa94-4d15-941b-07e9856cb12f' class='xr-section-summary-in' type='checkbox'  checked><label for='section-0893303d-aa94-4d15-941b-07e9856cb12f' class='xr-section-summary' >Attributes: <span>(4)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>spec :</span></dt><dd>RasterSpec(epsg=32750, bounds=(600000.0, 7690240.0, 610000.0, 7700240.0), resolutions_xy=(10.0, 10.0))</dd><dt><span>crs :</span></dt><dd>epsg:32750</dd><dt><span>transform :</span></dt><dd>| 10.00, 0.00, 600000.00|\n",
       "| 0.00,-10.00, 7700240.00|\n",
       "| 0.00, 0.00, 1.00|</dd><dt><span>resolution :</span></dt><dd>10.0</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.DataArray 'stackstac-7f8b9fd36f84dee2d0c759174b259b65' (time: 2,\n",
       "                                                                band: 6,\n",
       "                                                                y: 1000, x: 1000)>\n",
       "dask.array<fetch_raster_window, shape=(2, 6, 1000, 1000), dtype=float64, chunksize=(1, 1, 1000, 1000), chunktype=numpy.ndarray>\n",
       "Coordinates: (12/44)\n",
       "  * time                                     (time) datetime64[ns] 2022-01-06...\n",
       "    id                                       (time) <U54 'S2B_MSIL2A_20220106...\n",
       "  * band                                     (band) <U3 'B02' 'B03' ... 'AOT'\n",
       "  * x                                        (x) float64 6e+05 6e+05 ... 6.1e+05\n",
       "  * y                                        (y) float64 7.7e+06 ... 7.69e+06\n",
       "    s2:reflectance_conversion_factor         (time) float64 1.034 1.034\n",
       "    ...                                       ...\n",
       "    proj:bbox                                object {600000.0, 7690240.0, 780...\n",
       "    gsd                                      (band) float64 10.0 10.0 ... 10.0\n",
       "    common_name                              (band) object 'blue' ... None\n",
       "    center_wavelength                        (band) object 0.49 0.56 ... None\n",
       "    full_width_half_max                      (band) object 0.098 0.045 ... None\n",
       "    epsg                                     int64 32750\n",
       "Attributes:\n",
       "    spec:        RasterSpec(epsg=32750, bounds=(600000.0, 7690240.0, 610000.0...\n",
       "    crs:         epsg:32750\n",
       "    transform:   | 10.00, 0.00, 600000.00|\\n| 0.00,-10.00, 7700240.00|\\n| 0.0...\n",
       "    resolution:  10.0"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "bounds = (600000, 7690240, 610000, 7700240)\n",
    "assets = [\"B02\", \"B03\", \"B04\", \"B05\", \"B06\", \"AOT\"]\n",
    "ds = stackstac.stack([old, new], assets=assets, bounds=bounds)\n",
    "ds"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ef0f882f-99e3-422a-8973-78d2ec999486",
   "metadata": {},
   "source": [
    "If we plot the red, green, and blue bands of the two datetimes, we'll see that the \"old\" data looks relatively dark compared to the \"new\" data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "49eebc5c-5d71-401c-836b-8689624dcc27",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<img src=\"https://ai4edatasetspublicassets.blob.core.windows.net/assets/notebook-output/datasets-sentinel-2-l2a-baseline-change.ipynb/5.png\"/>"
      ],
      "text/plain": [
       "<Figure size 1400x600 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, axes = plt.subplots(ncols=2, figsize=(14, 6))\n",
    "\n",
    "img = ds.sel(band=[\"B04\", \"B03\", \"B02\"]).compute()\n",
    "img = xr.apply_ufunc(skimage.exposure.rescale_intensity, img)\n",
    "\n",
    "img.isel(time=0).plot.imshow(rgb=\"band\", ax=axes[0])\n",
    "axes[0].set(title=\"Old\")\n",
    "\n",
    "img.isel(time=1).plot.imshow(rgb=\"band\", ax=axes[1])\n",
    "axes[1].set(title=\"New\")\n",
    "for ax in axes:\n",
    "    ax.set_axis_off()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "db22b0b0-b421-46c8-87ba-318bc084bfc7",
   "metadata": {},
   "source": [
    "By taking the histogram of, say, the green band, we'll see that the new data are shifted to the right."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "5e8b8eac-54bb-45a8-8677-21184958e5f0",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<img src=\"https://ai4edatasetspublicassets.blob.core.windows.net/assets/notebook-output/datasets-sentinel-2-l2a-baseline-change.ipynb/6.png\"/>"
      ],
      "text/plain": [
       "<Figure size 1200x800 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots(figsize=(12, 8))\n",
    "\n",
    "hist, bins = skimage.exposure.histogram(ds.isel(time=0, band=1).data.compute())\n",
    "ax.plot(bins, hist, label=\"Old\")\n",
    "\n",
    "hist, bins = skimage.exposure.histogram(ds.isel(time=1, band=1).data.compute())\n",
    "ax.plot(bins, hist, label=\"New\")\n",
    "plt.legend()\n",
    "sns.despine()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3d0058aa-9459-48c2-a9bb-ab817a5947fc",
   "metadata": {},
   "source": [
    "Now, let's harmonize data, shifting the new data to make it comparable to the old."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "d6d71030-e747-4c08-83f2-2809b0388031",
   "metadata": {},
   "outputs": [],
   "source": [
    "harmonized = harmonize_to_old(ds)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d966f3e0-d86a-4f1f-ab06-2518b8213e60",
   "metadata": {},
   "source": [
    "Now the RGB images look much more similar:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "f374c97b-9ca7-4437-bfef-1c0a03428a2c",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<img src=\"https://ai4edatasetspublicassets.blob.core.windows.net/assets/notebook-output/datasets-sentinel-2-l2a-baseline-change.ipynb/8.png\"/>"
      ],
      "text/plain": [
       "<Figure size 1400x600 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, axes = plt.subplots(ncols=2, figsize=(14, 6))\n",
    "\n",
    "img = harmonized.sel(band=[\"B04\", \"B03\", \"B02\"]).compute()\n",
    "img = xr.apply_ufunc(skimage.exposure.rescale_intensity, img)\n",
    "\n",
    "img.isel(time=0).plot.imshow(rgb=\"band\", ax=axes[0])\n",
    "axes[0].set(title=\"Old\")\n",
    "\n",
    "img.isel(time=1).plot.imshow(rgb=\"band\", ax=axes[1])\n",
    "axes[1].set(title=\"New\")\n",
    "for ax in axes:\n",
    "    ax.set_axis_off()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "0cae2347-76c9-4822-a481-c5ba5bb108ce",
   "metadata": {},
   "source": [
    "And the histograms now mostly overlap. Any remaining differences are a result of the two scenes being captured a year apart."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "f306d91f-626f-4814-b6a5-25b379aafeaf",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<img src=\"https://ai4edatasetspublicassets.blob.core.windows.net/assets/notebook-output/datasets-sentinel-2-l2a-baseline-change.ipynb/9.png\"/>"
      ],
      "text/plain": [
       "<Figure size 1200x800 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots(figsize=(12, 8))\n",
    "\n",
    "hist, bins = skimage.exposure.histogram(harmonized.isel(time=0, band=1).data.compute())\n",
    "ax.plot(bins, hist, label=\"Old\")\n",
    "\n",
    "hist, bins = skimage.exposure.histogram(harmonized.isel(time=1, band=1).data.compute())\n",
    "ax.plot(bins, hist, label=\"New\")\n",
    "plt.legend()\n",
    "sns.despine()"
   ]
  }
 ],
 "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.10.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}