datasets/io-lulc/io-lulc-example.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "657431e4-bf40-4708-8639-0b014997aaea",
   "metadata": {},
   "source": [
    "## Accessing Esri 10m Land Use/Land Cover data with the Planetary Computer STAC API\n",
    "\n",
    "This dataset contains global estimates of 10-class land use/land cover for 2020, derived from ESA Sentinel-2 imagery at 10m resolution. In this notebook, we'll demonstrate how to access and work with this data through the Planetary Computer."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "dc99c9cc",
   "metadata": {},
   "source": [
    "### Environment setup\n",
    "\n",
    "This notebook works with or without an API key, but you will be given more permissive access to the data with an API key.\n",
    "The [Planetary Computer Hub](https://planetarycomputer.microsoft.com/compute) is pre-configured to use your API key."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "b87fbc90-35c8-42fd-92a8-b1a12e8676dd",
   "metadata": {},
   "outputs": [],
   "source": [
    "import dask.distributed\n",
    "from matplotlib.colors import ListedColormap\n",
    "import pystac_client\n",
    "from pystac.extensions.projection import ProjectionExtension as proj\n",
    "\n",
    "import cartopy.crs as ccrs\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "import planetary_computer\n",
    "import rasterio\n",
    "import rasterio.features\n",
    "import stackstac\n",
    "\n",
    "# Set the environment variable PC_SDK_SUBSCRIPTION_KEY, or set it here.\n",
    "# The Hub sets PC_SDK_SUBSCRIPTION_KEY automatically.\n",
    "# pc.settings.set_subscription_key(<YOUR API Key>)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c2e0bed0-d1f9-4164-b6de-bb743fdc3683",
   "metadata": {},
   "source": [
    "### Data access\n",
    "\n",
    "The datasets hosted by the Planetary Computer are available from [Azure Blob Storage](https://docs.microsoft.com/en-us/azure/storage/blobs/). We'll use [pystac-client](https://pystac-client.readthedocs.io/) to search the Planetary Computer's [STAC API](https://planetarycomputer.microsoft.com/api/stac/v1/docs) for the subset of the data that we care about, and then we'll load the data directly from Azure Blob Storage. We'll specify a `modifier` so that we can access the data stored in the Planetary Computer's private Blob Storage Containers. See [Reading from the STAC API](https://planetarycomputer.microsoft.com/docs/quickstarts/reading-stac/) and [Using tokens for data access](https://planetarycomputer.microsoft.com/docs/concepts/sas/) for more."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "a668267d-abab-45f0-bca1-c8a0f5852014",
   "metadata": {},
   "outputs": [],
   "source": [
    "catalog = pystac_client.Client.open(\n",
    "    \"https://planetarycomputer.microsoft.com/api/stac/v1\",\n",
    "    modifier=planetary_computer.sign_inplace,\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "901a784a-842d-4aea-bd80-54d256e3fa4f",
   "metadata": {},
   "source": [
    "### Select a region and find data items\n",
    "\n",
    "We'll pick an area surrounding Manila, Philippines and use the STAC API to find what data items are available. We won't select a date range since this dataset contains items from a single timeframe in 2020."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "53cece1f-5c0b-469c-94b4-5c111829d295",
   "metadata": {},
   "outputs": [],
   "source": [
    "area_of_interest = {\n",
    "    \"type\": \"Polygon\",\n",
    "    \"coordinates\": [\n",
    "        [\n",
    "            [120.88256835937499, 14.466596475463248],\n",
    "            [121.34948730468749, 14.466596475463248],\n",
    "            [121.34948730468749, 14.81737062015525],\n",
    "            [120.88256835937499, 14.817370620155254],\n",
    "            [120.88256835937499, 14.466596475463248],\n",
    "        ]\n",
    "    ],\n",
    "}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "10e74012-4398-4762-8074-ae135bb66dc5",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Returned 4 Items\n"
     ]
    }
   ],
   "source": [
    "search = catalog.search(collections=[\"io-lulc\"], intersects=area_of_interest)\n",
    "\n",
    "# Check how many items were returned\n",
    "items = search.item_collection()\n",
    "print(f\"Returned {len(items)} Items\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9d800526-aebb-44cf-a0df-6e93041d7f2f",
   "metadata": {},
   "source": [
    "We found 4 items that intersect with our area of interest, which means the data we want to work with is spread out over 4 non-overlapping GeoTIFF files stored on blob storage. In order to merge them together, we could open each item, clip to the subset of our AoI, and merge them together manually with rasterio. We'd also have to reproject each item which may span multiple UTM projections. \n",
    "\n",
    "Instead, we'll use the [stackstac](https://stackstac.readthedocs.io/en/latest/) library to read, merge, and reproject in a single step - all without loading the rest of the file data we don't need."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "2ca1a7a9-0229-4b0d-bb2f-599f4bd4e3c3",
   "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-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.DataArray &#x27;stackstac-5c28b96685cb1458a25519df2b6bfa98&#x27; (time: 1,\n",
       "                                                                band: 1,\n",
       "                                                                y: 3924, x: 5063)&gt;\n",
       "dask.array&lt;fetch_raster_window, shape=(1, 1, 3924, 5063), dtype=uint8, chunksize=(1, 1, 1024, 1024), chunktype=numpy.ndarray&gt;\n",
       "Coordinates: (12/20)\n",
       "  * time               (time) datetime64[ns] 2020-06-01\n",
       "    id                 (time) &lt;U8 &#x27;51P-2020&#x27;\n",
       "  * band               (band) &lt;U4 &#x27;data&#x27;\n",
       "  * x                  (x) float64 2.718e+05 2.718e+05 ... 3.224e+05 3.224e+05\n",
       "  * y                  (y) float64 1.639e+06 1.639e+06 ... 1.6e+06 1.6e+06\n",
       "    proj:epsg          int64 32651\n",
       "    ...                 ...\n",
       "    proj:transform     object {0.0, 137987.1216733564, 10.0, 1799862.99727356...\n",
       "    proj:bbox          object {908347.1216733564, 137987.1216733564, 790452.9...\n",
       "    label:type         &lt;U6 &#x27;raster&#x27;\n",
       "    file:size          int64 103458531\n",
       "    raster:bands       object {&#x27;nodata&#x27;: 0, &#x27;spatial_resolution&#x27;: 10}\n",
       "    epsg               int64 32651\n",
       "Attributes:\n",
       "    spec:        RasterSpec(epsg=32651, bounds=(271760.0, 1599970.0, 322390.0...\n",
       "    crs:         epsg:32651\n",
       "    transform:   | 10.00, 0.00, 271760.00|\\n| 0.00,-10.00, 1639210.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-5c28b96685cb1458a25519df2b6bfa98'</div><ul class='xr-dim-list'><li><span class='xr-has-index'>time</span>: 1</li><li><span class='xr-has-index'>band</span>: 1</li><li><span class='xr-has-index'>y</span>: 3924</li><li><span class='xr-has-index'>x</span>: 5063</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-1b1663b7-f1e0-4621-9042-74486c8a320b' class='xr-array-in' type='checkbox' checked><label for='section-1b1663b7-f1e0-4621-9042-74486c8a320b' 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, 1024, 1024), meta=np.ndarray&gt;</span></div><div class='xr-array-data'><table>\n",
       "    <tr>\n",
       "        <td>\n",
       "            <table>\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> 18.95 MiB </td>\n",
       "                        <td> 1.00 MiB </td>\n",
       "                    </tr>\n",
       "                    \n",
       "                    <tr>\n",
       "                        <th> Shape </th>\n",
       "                        <td> (1, 1, 3924, 5063) </td>\n",
       "                        <td> (1, 1, 1024, 1024) </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <th> Count </th>\n",
       "                        <td> 3 Graph Layers </td>\n",
       "                        <td> 20 Chunks </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                    <th> Type </th>\n",
       "                    <td> uint8 </td>\n",
       "                    <td> numpy.ndarray </td>\n",
       "                    </tr>\n",
       "                </tbody>\n",
       "            </table>\n",
       "        </td>\n",
       "        <td>\n",
       "        <svg width=\"374\" height=\"157\" 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=\"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\" >1</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=\"24\" x2=\"109\" y2=\"39\" />\n",
       "  <line x1=\"95\" y1=\"48\" x2=\"109\" y2=\"63\" />\n",
       "  <line x1=\"95\" y1=\"72\" x2=\"109\" y2=\"87\" />\n",
       "  <line x1=\"95\" y1=\"93\" x2=\"109\" y2=\"107\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"95\" y1=\"0\" x2=\"95\" y2=\"93\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"109\" y1=\"14\" x2=\"109\" y2=\"107\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"95.0,0.0 109.9485979497544,14.948597949754403 109.9485979497544,107.95274568824935 95.0,93.00414773849495\" 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=\"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=\"119\" y1=\"0\" x2=\"134\" y2=\"14\" />\n",
       "  <line x1=\"143\" y1=\"0\" x2=\"158\" y2=\"14\" />\n",
       "  <line x1=\"167\" y1=\"0\" x2=\"182\" y2=\"14\" />\n",
       "  <line x1=\"192\" y1=\"0\" x2=\"207\" y2=\"14\" />\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=\"39\" x2=\"229\" y2=\"39\" />\n",
       "  <line x1=\"109\" y1=\"63\" x2=\"229\" y2=\"63\" />\n",
       "  <line x1=\"109\" y1=\"87\" x2=\"229\" y2=\"87\" />\n",
       "  <line x1=\"109\" y1=\"107\" x2=\"229\" y2=\"107\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"109\" y1=\"14\" x2=\"109\" y2=\"107\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"134\" y1=\"14\" x2=\"134\" y2=\"107\" />\n",
       "  <line x1=\"158\" y1=\"14\" x2=\"158\" y2=\"107\" />\n",
       "  <line x1=\"182\" y1=\"14\" x2=\"182\" y2=\"107\" />\n",
       "  <line x1=\"207\" y1=\"14\" x2=\"207\" y2=\"107\" />\n",
       "  <line x1=\"229\" y1=\"14\" x2=\"229\" y2=\"107\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"109.9485979497544,14.948597949754403 229.9485979497544,14.948597949754403 229.9485979497544,107.95274568824935 109.9485979497544,107.95274568824935\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Text -->\n",
       "  <text x=\"169.948598\" y=\"127.952746\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" >5063</text>\n",
       "  <text x=\"249.948598\" y=\"61.450672\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(-90,249.948598,61.450672)\">3924</text>\n",
       "  <text x=\"92.474299\" y=\"120.478447\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(45,92.474299,120.478447)\">1</text>\n",
       "</svg>\n",
       "        </td>\n",
       "    </tr>\n",
       "</table></div></div></li><li class='xr-section-item'><input id='section-974f0d1f-da1d-4635-ac9e-e11069f3dfd1' class='xr-section-summary-in' type='checkbox'  checked><label for='section-974f0d1f-da1d-4635-ac9e-e11069f3dfd1' class='xr-section-summary' >Coordinates: <span>(20)</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'>2020-06-01</div><input id='attrs-4091823c-18b3-4e87-9702-3bac0a50d968' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-4091823c-18b3-4e87-9702-3bac0a50d968' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3ef00825-1d51-4433-a637-844930589f54' class='xr-var-data-in' type='checkbox'><label for='data-3ef00825-1d51-4433-a637-844930589f54' 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;2020-06-01T00:00:00.000000000&#x27;], 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;U8</div><div class='xr-var-preview xr-preview'>&#x27;51P-2020&#x27;</div><input id='attrs-57cbc01b-cac8-4412-b393-cadaaa99f3a3' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-57cbc01b-cac8-4412-b393-cadaaa99f3a3' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-21e48474-20b0-4b27-a7e9-5cfae7735a0b' class='xr-var-data-in' type='checkbox'><label for='data-21e48474-20b0-4b27-a7e9-5cfae7735a0b' 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;51P-2020&#x27;], dtype=&#x27;&lt;U8&#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;U4</div><div class='xr-var-preview xr-preview'>&#x27;data&#x27;</div><input id='attrs-30e057c5-144b-4d64-87cf-cd6ae46d163e' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-30e057c5-144b-4d64-87cf-cd6ae46d163e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8affc483-41a0-4f6c-ae1d-0207208cc6c4' class='xr-var-data-in' type='checkbox'><label for='data-8affc483-41a0-4f6c-ae1d-0207208cc6c4' 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;data&#x27;], dtype=&#x27;&lt;U4&#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'>2.718e+05 2.718e+05 ... 3.224e+05</div><input id='attrs-c0c2d08a-ae58-437e-ab01-a2759dd6d648' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-c0c2d08a-ae58-437e-ab01-a2759dd6d648' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5eea985e-0a8f-46d7-99c2-177a4fb0bc65' class='xr-var-data-in' type='checkbox'><label for='data-5eea985e-0a8f-46d7-99c2-177a4fb0bc65' 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([271760., 271770., 271780., ..., 322360., 322370., 322380.])</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'>1.639e+06 1.639e+06 ... 1.6e+06</div><input id='attrs-37244081-9846-43d6-9f34-91f96304afd7' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-37244081-9846-43d6-9f34-91f96304afd7' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b7d5265e-ea15-4626-9435-273fd6e1a2ee' class='xr-var-data-in' type='checkbox'><label for='data-b7d5265e-ea15-4626-9435-273fd6e1a2ee' 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([1639210., 1639200., 1639190., ..., 1600000., 1599990., 1599980.])</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'>32651</div><input id='attrs-1814a86e-dd93-42e0-8be0-4d5c748b1c98' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-1814a86e-dd93-42e0-8be0-4d5c748b1c98' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-580aaa86-6cc9-454e-8f2b-404d1eb6b77a' class='xr-var-data-in' type='checkbox'><label for='data-580aaa86-6cc9-454e-8f2b-404d1eb6b77a' 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(32651)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>label:classes</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>{&#x27;name&#x27;: &#x27;&#x27;, &#x27;classes&#x27;: [&#x27;nodata...</div><input id='attrs-7e61c781-ab6b-47cf-ad93-554c7162eab6' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-7e61c781-ab6b-47cf-ad93-554c7162eab6' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f891af97-bb41-43ff-b280-382223daf60d' class='xr-var-data-in' type='checkbox'><label for='data-f891af97-bb41-43ff-b280-382223daf60d' 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;name&#x27;: &#x27;&#x27;, &#x27;classes&#x27;: [&#x27;nodata&#x27;, &#x27;water&#x27;, &#x27;trees&#x27;, &#x27;grass&#x27;, &#x27;flooded veg&#x27;, &#x27;crops&#x27;, &#x27;scrub&#x27;, &#x27;built area&#x27;, &#x27;bare&#x27;, &#x27;snow/ice&#x27;, &#x27;clouds&#x27;]},\n",
       "      dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>start_datetime</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U20</div><div class='xr-var-preview xr-preview'>&#x27;2020-01-01T00:00:00Z&#x27;</div><input id='attrs-29bf2a6b-c672-4f02-98cb-c1659ed5fe96' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-29bf2a6b-c672-4f02-98cb-c1659ed5fe96' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-2d62335b-fbac-44ea-84df-c4650e42422b' class='xr-var-data-in' type='checkbox'><label for='data-2d62335b-fbac-44ea-84df-c4650e42422b' 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;2020-01-01T00:00:00Z&#x27;, dtype=&#x27;&lt;U20&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>end_datetime</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U20</div><div class='xr-var-preview xr-preview'>&#x27;2021-01-01T00:00:00Z&#x27;</div><input id='attrs-307fd12d-79e3-481c-964c-c20f07b8b5ba' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-307fd12d-79e3-481c-964c-c20f07b8b5ba' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8c4a8240-ce20-4284-a837-b8eda31d54c3' class='xr-var-data-in' type='checkbox'><label for='data-8c4a8240-ce20-4284-a837-b8eda31d54c3' 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;2021-01-01T00:00:00Z&#x27;, dtype=&#x27;&lt;U20&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>label:description</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U4</div><div class='xr-var-preview xr-preview'>&#x27;lulc&#x27;</div><input id='attrs-7e942afe-426c-425a-92e1-7ea2d8c371ec' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-7e942afe-426c-425a-92e1-7ea2d8c371ec' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-299d5e16-74ae-4f58-add6-037353eaa7d9' class='xr-var-data-in' type='checkbox'><label for='data-299d5e16-74ae-4f58-add6-037353eaa7d9' 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;lulc&#x27;, dtype=&#x27;&lt;U4&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>io:tile_id</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U3</div><div class='xr-var-preview xr-preview'>&#x27;51P&#x27;</div><input id='attrs-0d2f7ea3-c017-4a92-a184-f6bfe614187a' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-0d2f7ea3-c017-4a92-a184-f6bfe614187a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-a85aa373-7ac6-4bc8-93d7-f03d5f684331' class='xr-var-data-in' type='checkbox'><label for='data-a85aa373-7ac6-4bc8-93d7-f03d5f684331' 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;51P&#x27;, dtype=&#x27;&lt;U3&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>label:properties</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>None</div><input id='attrs-baa9f27a-028e-42a9-99da-096b5715ab69' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-baa9f27a-028e-42a9-99da-096b5715ab69' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-19b7140e-e25d-45ca-9510-2d5246189f36' class='xr-var-data-in' type='checkbox'><label for='data-19b7140e-e25d-45ca-9510-2d5246189f36' 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(None, dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>proj:shape</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>{77036, 100941}</div><input id='attrs-7109ea2f-8cf7-4409-87ba-8d9a469d63dd' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-7109ea2f-8cf7-4409-87ba-8d9a469d63dd' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-257b70b6-cc18-4175-a50d-432f19f8dad0' class='xr-var-data-in' type='checkbox'><label for='data-257b70b6-cc18-4175-a50d-432f19f8dad0' 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({77036, 100941}, dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>io:supercell_id</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U3</div><div class='xr-var-preview xr-preview'>&#x27;51P&#x27;</div><input id='attrs-6fb143f6-09c5-46cc-b9bc-8fe471f300b9' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-6fb143f6-09c5-46cc-b9bc-8fe471f300b9' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-03be4b78-037d-4092-a1bc-41f1fe7c1ec9' class='xr-var-data-in' type='checkbox'><label for='data-03be4b78-037d-4092-a1bc-41f1fe7c1ec9' 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;51P&#x27;, dtype=&#x27;&lt;U3&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>proj:transform</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>{0.0, 137987.1216733564, 10.0, 1...</div><input id='attrs-645cad6b-62de-46df-a297-5687126ed9b0' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-645cad6b-62de-46df-a297-5687126ed9b0' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-2e688de7-9cc1-4bd7-a629-107688492dbc' class='xr-var-data-in' type='checkbox'><label for='data-2e688de7-9cc1-4bd7-a629-107688492dbc' 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.0, 137987.1216733564, 10.0, 1799862.9972735613, -10.0},\n",
       "      dtype=object)</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'>{908347.1216733564, 137987.12167...</div><input id='attrs-c5c3d935-dfbd-419d-aab5-7e3e44324903' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-c5c3d935-dfbd-419d-aab5-7e3e44324903' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-96e3b271-bb5d-4fbf-9152-31f1fc60224c' class='xr-var-data-in' type='checkbox'><label for='data-96e3b271-bb5d-4fbf-9152-31f1fc60224c' 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({908347.1216733564, 137987.1216733564, 790452.9972735613, 1799862.9972735613},\n",
       "      dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>label:type</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U6</div><div class='xr-var-preview xr-preview'>&#x27;raster&#x27;</div><input id='attrs-02d747ef-a788-4459-ab7c-e2da9a9eba59' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-02d747ef-a788-4459-ab7c-e2da9a9eba59' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-ddaa4b7d-5a4c-4d2a-9694-75c1d925a815' class='xr-var-data-in' type='checkbox'><label for='data-ddaa4b7d-5a4c-4d2a-9694-75c1d925a815' 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;raster&#x27;, dtype=&#x27;&lt;U6&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>file:size</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>103458531</div><input id='attrs-8ae6adfd-5d6c-4ea5-a287-4a67e7954e48' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-8ae6adfd-5d6c-4ea5-a287-4a67e7954e48' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3bf400a0-0366-40c4-b916-f0b61346cf77' class='xr-var-data-in' type='checkbox'><label for='data-3bf400a0-0366-40c4-b916-f0b61346cf77' 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(103458531)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>raster:bands</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>{&#x27;nodata&#x27;: 0, &#x27;spatial_resolutio...</div><input id='attrs-47bb323e-203d-4798-9f31-6a174cae333a' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-47bb323e-203d-4798-9f31-6a174cae333a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-050e7942-9cdd-4d0a-969f-ff51e69f55e3' class='xr-var-data-in' type='checkbox'><label for='data-050e7942-9cdd-4d0a-969f-ff51e69f55e3' 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;nodata&#x27;: 0, &#x27;spatial_resolution&#x27;: 10}, 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'>32651</div><input id='attrs-d101b358-2ffc-485b-a762-c6571528ee75' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-d101b358-2ffc-485b-a762-c6571528ee75' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-6d7b045f-7202-43ed-8994-756c5019aef5' class='xr-var-data-in' type='checkbox'><label for='data-6d7b045f-7202-43ed-8994-756c5019aef5' 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(32651)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-3ac46ce6-d9f6-49e9-8601-b4a2c4fe58a5' class='xr-section-summary-in' type='checkbox'  checked><label for='section-3ac46ce6-d9f6-49e9-8601-b4a2c4fe58a5' 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=32651, bounds=(271760.0, 1599970.0, 322390.0, 1639210.0), resolutions_xy=(10.0, 10.0))</dd><dt><span>crs :</span></dt><dd>epsg:32651</dd><dt><span>transform :</span></dt><dd>| 10.00, 0.00, 271760.00|\n",
       "| 0.00,-10.00, 1639210.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-5c28b96685cb1458a25519df2b6bfa98' (time: 1,\n",
       "                                                                band: 1,\n",
       "                                                                y: 3924, x: 5063)>\n",
       "dask.array<fetch_raster_window, shape=(1, 1, 3924, 5063), dtype=uint8, chunksize=(1, 1, 1024, 1024), chunktype=numpy.ndarray>\n",
       "Coordinates: (12/20)\n",
       "  * time               (time) datetime64[ns] 2020-06-01\n",
       "    id                 (time) <U8 '51P-2020'\n",
       "  * band               (band) <U4 'data'\n",
       "  * x                  (x) float64 2.718e+05 2.718e+05 ... 3.224e+05 3.224e+05\n",
       "  * y                  (y) float64 1.639e+06 1.639e+06 ... 1.6e+06 1.6e+06\n",
       "    proj:epsg          int64 32651\n",
       "    ...                 ...\n",
       "    proj:transform     object {0.0, 137987.1216733564, 10.0, 1799862.99727356...\n",
       "    proj:bbox          object {908347.1216733564, 137987.1216733564, 790452.9...\n",
       "    label:type         <U6 'raster'\n",
       "    file:size          int64 103458531\n",
       "    raster:bands       object {'nodata': 0, 'spatial_resolution': 10}\n",
       "    epsg               int64 32651\n",
       "Attributes:\n",
       "    spec:        RasterSpec(epsg=32651, bounds=(271760.0, 1599970.0, 322390.0...\n",
       "    crs:         epsg:32651\n",
       "    transform:   | 10.00, 0.00, 271760.00|\\n| 0.00,-10.00, 1639210.00|\\n| 0.0...\n",
       "    resolution:  10.0"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# The STAC metadata contains some information we'll want to use when creating\n",
    "# our merged dataset. Get the EPSG code of the first item and the nodata value.\n",
    "item = items[0]\n",
    "epsg = proj.ext(item).epsg\n",
    "nodata = item.assets[\"data\"].extra_fields[\"raster:bands\"][0][\"nodata\"]\n",
    "bounds_latlon = rasterio.features.bounds(area_of_interest)\n",
    "\n",
    "# Create a single DataArray from out multiple resutls with the corresponding\n",
    "# rasters projected to a single CRS. Note that we set the dtype to ubyte, which\n",
    "# matches our data, since stackstac will use float64 by default.\n",
    "stack = stackstac.stack(\n",
    "    items, epsg=epsg, dtype=np.ubyte, fill_value=nodata, bounds_latlon=bounds_latlon\n",
    ")\n",
    "\n",
    "stack"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "2c8b22c9-9260-48af-8dc9-bb46fb531673",
   "metadata": {},
   "source": [
    "Start up a local Dask cluster to allow us to do parallel reads. Use the following URL to open a dashboard in the Hub's Dask Extension. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "2161565e-163c-4996-a5cd-0c6b14390eee",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/srv/conda/envs/notebook/lib/python3.10/site-packages/distributed/node.py:183: UserWarning: Port 8787 is already in use.\n",
      "Perhaps you already have a cluster running?\n",
      "Hosting the HTTP server on port 46119 instead\n",
      "  warnings.warn(\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "/proxy/46119/status\n"
     ]
    }
   ],
   "source": [
    "client = dask.distributed.Client(processes=False)\n",
    "print(f\"/proxy/{client.scheduler_info()['services']['dashboard']}/status\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "a0e8e4e2-c477-4a43-885e-59fc1f5a86fb",
   "metadata": {},
   "source": [
    "### Mosaic and clip the raster\n",
    "\n",
    "So far, we haven't read in any data. Stackstac has used the STAC metadata to construct a DataArray that will contain our Item data. Let's mosaic the rasters across the `time` dimension (remember, they're all from a single synthesized \"time\" from 2020) and drop the single `band` dimension. Finally, we ask Dask to read the actual data by calling `.compute()`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "5108587f-10d1-4b52-9861-3aeb161bce46",
   "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-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.DataArray &#x27;stackstac-5c28b96685cb1458a25519df2b6bfa98&#x27; (y: 3924, x: 5063)&gt;\n",
       "array([[5, 5, 5, ..., 2, 2, 2],\n",
       "       [5, 5, 5, ..., 2, 2, 2],\n",
       "       [5, 5, 5, ..., 2, 2, 2],\n",
       "       ...,\n",
       "       [7, 7, 7, ..., 2, 2, 2],\n",
       "       [7, 7, 7, ..., 2, 2, 2],\n",
       "       [7, 7, 7, ..., 2, 2, 2]], dtype=uint8)\n",
       "Coordinates: (12/18)\n",
       "    band               &lt;U4 &#x27;data&#x27;\n",
       "  * x                  (x) float64 2.718e+05 2.718e+05 ... 3.224e+05 3.224e+05\n",
       "  * y                  (y) float64 1.639e+06 1.639e+06 ... 1.6e+06 1.6e+06\n",
       "    proj:epsg          int64 32651\n",
       "    label:classes      object {&#x27;name&#x27;: &#x27;&#x27;, &#x27;classes&#x27;: [&#x27;nodata&#x27;, &#x27;water&#x27;, &#x27;tr...\n",
       "    start_datetime     &lt;U20 &#x27;2020-01-01T00:00:00Z&#x27;\n",
       "    ...                 ...\n",
       "    proj:transform     object {0.0, 137987.1216733564, 10.0, 1799862.99727356...\n",
       "    proj:bbox          object {908347.1216733564, 137987.1216733564, 790452.9...\n",
       "    label:type         &lt;U6 &#x27;raster&#x27;\n",
       "    file:size          int64 103458531\n",
       "    raster:bands       object {&#x27;nodata&#x27;: 0, &#x27;spatial_resolution&#x27;: 10}\n",
       "    epsg               int64 32651\n",
       "Attributes:\n",
       "    spec:        RasterSpec(epsg=32651, bounds=(271760.0, 1599970.0, 322390.0...\n",
       "    crs:         epsg:32651\n",
       "    transform:   | 10.00, 0.00, 271760.00|\\n| 0.00,-10.00, 1639210.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-5c28b96685cb1458a25519df2b6bfa98'</div><ul class='xr-dim-list'><li><span class='xr-has-index'>y</span>: 3924</li><li><span class='xr-has-index'>x</span>: 5063</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-dcda6672-8259-4bd4-b51a-67992b15ea4d' class='xr-array-in' type='checkbox' checked><label for='section-dcda6672-8259-4bd4-b51a-67992b15ea4d' 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>5 5 5 5 5 5 5 1 1 1 1 1 1 1 1 1 1 ... 2 2 2 2 2 2 2 2 6 6 6 6 6 2 2 2</span></div><div class='xr-array-data'><pre>array([[5, 5, 5, ..., 2, 2, 2],\n",
       "       [5, 5, 5, ..., 2, 2, 2],\n",
       "       [5, 5, 5, ..., 2, 2, 2],\n",
       "       ...,\n",
       "       [7, 7, 7, ..., 2, 2, 2],\n",
       "       [7, 7, 7, ..., 2, 2, 2],\n",
       "       [7, 7, 7, ..., 2, 2, 2]], dtype=uint8)</pre></div></div></li><li class='xr-section-item'><input id='section-7d93dbfc-1068-46fc-806c-3233a09e078e' class='xr-section-summary-in' type='checkbox'  checked><label for='section-7d93dbfc-1068-46fc-806c-3233a09e078e' class='xr-section-summary' >Coordinates: <span>(18)</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>band</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U4</div><div class='xr-var-preview xr-preview'>&#x27;data&#x27;</div><input id='attrs-27b589be-69b3-45a9-9a55-5b73b6aa01c6' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-27b589be-69b3-45a9-9a55-5b73b6aa01c6' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-2d295420-6143-4bfd-b868-ffa102d4ec27' class='xr-var-data-in' type='checkbox'><label for='data-2d295420-6143-4bfd-b868-ffa102d4ec27' 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;data&#x27;, dtype=&#x27;&lt;U4&#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'>2.718e+05 2.718e+05 ... 3.224e+05</div><input id='attrs-08f4bab8-0875-4281-9b26-fc17f91da493' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-08f4bab8-0875-4281-9b26-fc17f91da493' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-d8d31502-ae11-4531-8e0f-87143b359ddd' class='xr-var-data-in' type='checkbox'><label for='data-d8d31502-ae11-4531-8e0f-87143b359ddd' 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([271760., 271770., 271780., ..., 322360., 322370., 322380.])</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'>1.639e+06 1.639e+06 ... 1.6e+06</div><input id='attrs-d849c3a2-8f6b-4350-a53b-bbbb0a347793' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-d849c3a2-8f6b-4350-a53b-bbbb0a347793' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e34de17c-67cb-4361-86fc-0548e51f4d7e' class='xr-var-data-in' type='checkbox'><label for='data-e34de17c-67cb-4361-86fc-0548e51f4d7e' 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([1639210., 1639200., 1639190., ..., 1600000., 1599990., 1599980.])</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'>32651</div><input id='attrs-6fcb0827-c0cf-4081-9412-932c5aaf18c3' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-6fcb0827-c0cf-4081-9412-932c5aaf18c3' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-7d406cf3-50c2-447a-85c6-a154ef653729' class='xr-var-data-in' type='checkbox'><label for='data-7d406cf3-50c2-447a-85c6-a154ef653729' 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(32651)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>label:classes</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>{&#x27;name&#x27;: &#x27;&#x27;, &#x27;classes&#x27;: [&#x27;nodata...</div><input id='attrs-48c2109a-67d5-4dca-9819-17ce81cb8fdc' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-48c2109a-67d5-4dca-9819-17ce81cb8fdc' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4f876885-c188-484a-870f-23f0c4d18873' class='xr-var-data-in' type='checkbox'><label for='data-4f876885-c188-484a-870f-23f0c4d18873' 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;name&#x27;: &#x27;&#x27;, &#x27;classes&#x27;: [&#x27;nodata&#x27;, &#x27;water&#x27;, &#x27;trees&#x27;, &#x27;grass&#x27;, &#x27;flooded veg&#x27;, &#x27;crops&#x27;, &#x27;scrub&#x27;, &#x27;built area&#x27;, &#x27;bare&#x27;, &#x27;snow/ice&#x27;, &#x27;clouds&#x27;]},\n",
       "      dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>start_datetime</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U20</div><div class='xr-var-preview xr-preview'>&#x27;2020-01-01T00:00:00Z&#x27;</div><input id='attrs-7e678eeb-27f6-4ca5-a39d-179f1caad12f' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-7e678eeb-27f6-4ca5-a39d-179f1caad12f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-85c8d791-ba31-41cd-a39a-fa7096fe9999' class='xr-var-data-in' type='checkbox'><label for='data-85c8d791-ba31-41cd-a39a-fa7096fe9999' 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;2020-01-01T00:00:00Z&#x27;, dtype=&#x27;&lt;U20&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>end_datetime</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U20</div><div class='xr-var-preview xr-preview'>&#x27;2021-01-01T00:00:00Z&#x27;</div><input id='attrs-02cf929d-a2ac-49fb-9a39-30083cc5d4aa' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-02cf929d-a2ac-49fb-9a39-30083cc5d4aa' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-d720a793-1b55-45e1-97a2-bc1c100647b7' class='xr-var-data-in' type='checkbox'><label for='data-d720a793-1b55-45e1-97a2-bc1c100647b7' 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;2021-01-01T00:00:00Z&#x27;, dtype=&#x27;&lt;U20&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>label:description</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U4</div><div class='xr-var-preview xr-preview'>&#x27;lulc&#x27;</div><input id='attrs-892b22f4-30a6-4cc5-970d-44223ec8acb0' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-892b22f4-30a6-4cc5-970d-44223ec8acb0' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-fb7eb678-b869-45d5-a3c1-663cb3ce679b' class='xr-var-data-in' type='checkbox'><label for='data-fb7eb678-b869-45d5-a3c1-663cb3ce679b' 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;lulc&#x27;, dtype=&#x27;&lt;U4&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>io:tile_id</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U3</div><div class='xr-var-preview xr-preview'>&#x27;51P&#x27;</div><input id='attrs-5bef5c67-5f12-4e45-a025-3c0ec88c427c' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-5bef5c67-5f12-4e45-a025-3c0ec88c427c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-2de699ef-c369-4df5-a325-c7323f9fddd0' class='xr-var-data-in' type='checkbox'><label for='data-2de699ef-c369-4df5-a325-c7323f9fddd0' 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;51P&#x27;, dtype=&#x27;&lt;U3&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>label:properties</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>None</div><input id='attrs-47b653fa-0eec-49b7-ba9f-530e3516e863' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-47b653fa-0eec-49b7-ba9f-530e3516e863' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-622c00df-9e07-4c6c-b4ad-0bc613610780' class='xr-var-data-in' type='checkbox'><label for='data-622c00df-9e07-4c6c-b4ad-0bc613610780' 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(None, dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>proj:shape</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>{77036, 100941}</div><input id='attrs-cea12768-e405-4d94-be7e-7f7c0e50fc9d' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-cea12768-e405-4d94-be7e-7f7c0e50fc9d' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-17d03858-973e-47ee-b5da-1d6d963b4e08' class='xr-var-data-in' type='checkbox'><label for='data-17d03858-973e-47ee-b5da-1d6d963b4e08' 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({77036, 100941}, dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>io:supercell_id</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U3</div><div class='xr-var-preview xr-preview'>&#x27;51P&#x27;</div><input id='attrs-35cfc1cf-6c60-498d-9a7c-e679737b159a' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-35cfc1cf-6c60-498d-9a7c-e679737b159a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-9d5b32de-1733-4a9b-b8de-13cbef3c8aeb' class='xr-var-data-in' type='checkbox'><label for='data-9d5b32de-1733-4a9b-b8de-13cbef3c8aeb' 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;51P&#x27;, dtype=&#x27;&lt;U3&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>proj:transform</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>{0.0, 137987.1216733564, 10.0, 1...</div><input id='attrs-4cce8a4d-041f-4084-b009-da1dd660ddc9' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-4cce8a4d-041f-4084-b009-da1dd660ddc9' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-eb62ddfc-fb4d-4ad6-91f4-3c514f2ea87f' class='xr-var-data-in' type='checkbox'><label for='data-eb62ddfc-fb4d-4ad6-91f4-3c514f2ea87f' 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.0, 137987.1216733564, 10.0, 1799862.9972735613, -10.0},\n",
       "      dtype=object)</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'>{908347.1216733564, 137987.12167...</div><input id='attrs-842bcbfa-6fac-4d68-bd71-7cc562ceb5df' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-842bcbfa-6fac-4d68-bd71-7cc562ceb5df' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-cdcf3ef2-08ed-45f0-96ea-b8fcea348d71' class='xr-var-data-in' type='checkbox'><label for='data-cdcf3ef2-08ed-45f0-96ea-b8fcea348d71' 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({908347.1216733564, 137987.1216733564, 790452.9972735613, 1799862.9972735613},\n",
       "      dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>label:type</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U6</div><div class='xr-var-preview xr-preview'>&#x27;raster&#x27;</div><input id='attrs-88089592-3e57-4c2d-8106-0810a039ed46' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-88089592-3e57-4c2d-8106-0810a039ed46' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-bad417ad-891b-4c66-ba25-53894b8efcf3' class='xr-var-data-in' type='checkbox'><label for='data-bad417ad-891b-4c66-ba25-53894b8efcf3' 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;raster&#x27;, dtype=&#x27;&lt;U6&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>file:size</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>103458531</div><input id='attrs-12305f99-4ed9-4fa7-9c48-f4be814e630b' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-12305f99-4ed9-4fa7-9c48-f4be814e630b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-a5e816d3-97ec-4a44-a1fe-8fac7f9fe552' class='xr-var-data-in' type='checkbox'><label for='data-a5e816d3-97ec-4a44-a1fe-8fac7f9fe552' 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(103458531)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>raster:bands</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>{&#x27;nodata&#x27;: 0, &#x27;spatial_resolutio...</div><input id='attrs-1fcf148f-2169-4476-871a-5d7f9419aba1' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-1fcf148f-2169-4476-871a-5d7f9419aba1' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-2026abb1-9fcb-403a-a36e-8a3e1d992b54' class='xr-var-data-in' type='checkbox'><label for='data-2026abb1-9fcb-403a-a36e-8a3e1d992b54' 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;nodata&#x27;: 0, &#x27;spatial_resolution&#x27;: 10}, 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'>32651</div><input id='attrs-98403b82-4135-41ca-9060-152ddc1b8a91' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-98403b82-4135-41ca-9060-152ddc1b8a91' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b62e2f20-eddc-4132-85d7-cf5b2593f258' class='xr-var-data-in' type='checkbox'><label for='data-b62e2f20-eddc-4132-85d7-cf5b2593f258' 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(32651)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-b3d2a601-c571-4bf0-bd32-1a7a8abdd652' class='xr-section-summary-in' type='checkbox'  checked><label for='section-b3d2a601-c571-4bf0-bd32-1a7a8abdd652' 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=32651, bounds=(271760.0, 1599970.0, 322390.0, 1639210.0), resolutions_xy=(10.0, 10.0))</dd><dt><span>crs :</span></dt><dd>epsg:32651</dd><dt><span>transform :</span></dt><dd>| 10.00, 0.00, 271760.00|\n",
       "| 0.00,-10.00, 1639210.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-5c28b96685cb1458a25519df2b6bfa98' (y: 3924, x: 5063)>\n",
       "array([[5, 5, 5, ..., 2, 2, 2],\n",
       "       [5, 5, 5, ..., 2, 2, 2],\n",
       "       [5, 5, 5, ..., 2, 2, 2],\n",
       "       ...,\n",
       "       [7, 7, 7, ..., 2, 2, 2],\n",
       "       [7, 7, 7, ..., 2, 2, 2],\n",
       "       [7, 7, 7, ..., 2, 2, 2]], dtype=uint8)\n",
       "Coordinates: (12/18)\n",
       "    band               <U4 'data'\n",
       "  * x                  (x) float64 2.718e+05 2.718e+05 ... 3.224e+05 3.224e+05\n",
       "  * y                  (y) float64 1.639e+06 1.639e+06 ... 1.6e+06 1.6e+06\n",
       "    proj:epsg          int64 32651\n",
       "    label:classes      object {'name': '', 'classes': ['nodata', 'water', 'tr...\n",
       "    start_datetime     <U20 '2020-01-01T00:00:00Z'\n",
       "    ...                 ...\n",
       "    proj:transform     object {0.0, 137987.1216733564, 10.0, 1799862.99727356...\n",
       "    proj:bbox          object {908347.1216733564, 137987.1216733564, 790452.9...\n",
       "    label:type         <U6 'raster'\n",
       "    file:size          int64 103458531\n",
       "    raster:bands       object {'nodata': 0, 'spatial_resolution': 10}\n",
       "    epsg               int64 32651\n",
       "Attributes:\n",
       "    spec:        RasterSpec(epsg=32651, bounds=(271760.0, 1599970.0, 322390.0...\n",
       "    crs:         epsg:32651\n",
       "    transform:   | 10.00, 0.00, 271760.00|\\n| 0.00,-10.00, 1639210.00|\\n| 0.0...\n",
       "    resolution:  10.0"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "merged = stackstac.mosaic(stack, dim=\"time\", axis=None, nodata=0).squeeze().compute()\n",
    "merged"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b5f31c78-fe4d-4b0d-bdcb-0289abc7e860",
   "metadata": {},
   "source": [
    "Now a quick plot to check that we've got the data we want."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "534e9b80-3c86-4029-a236-14653df64333",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<img src=\"https://ai4edatasetspublicassets.blob.core.windows.net/assets/notebook-output/datasets-io-lulc-io-lulc-example.ipynb/8.png\"/>"
      ],
      "text/plain": [
       "<Figure size 640x480 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "merged.plot()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "cc3c5100-12f2-4488-a19a-cc174efeee6e",
   "metadata": {},
   "source": [
    "It looks good, but it doesn't look like a land cover map. The source GeoTIFFs contain a colormap and the STAC metadata contains the class names. We'll open one of the source files just to read this metadata and construct the right colors and names for our plot."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "17570d98-c60a-4aa9-a6fa-9164cb9ba733",
   "metadata": {},
   "outputs": [],
   "source": [
    "class_names = merged.coords[\"label:classes\"].item()[\"classes\"]\n",
    "class_count = len(class_names)\n",
    "\n",
    "with rasterio.open(item.assets[\"data\"].href) as src:\n",
    "    colormap_def = src.colormap(1)  # get metadata colormap for band 1\n",
    "    colormap = [\n",
    "        np.array(colormap_def[i]) / 255 for i in range(class_count)\n",
    "    ]  # transform to matplotlib color format\n",
    "\n",
    "cmap = ListedColormap(colormap)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "b9d365ca-a063-4640-a00f-731f340f40d3",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<img src=\"https://ai4edatasetspublicassets.blob.core.windows.net/assets/notebook-output/datasets-io-lulc-io-lulc-example.ipynb/10.png\"/>"
      ],
      "text/plain": [
       "<Figure size 1500x750 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots(\n",
    "    figsize=(12, 6), dpi=125, subplot_kw=dict(projection=ccrs.epsg(epsg)), frameon=False\n",
    ")\n",
    "p = merged.plot(\n",
    "    ax=ax,\n",
    "    transform=ccrs.epsg(epsg),\n",
    "    cmap=cmap,\n",
    "    add_colorbar=False,\n",
    "    vmin=0,\n",
    "    vmax=class_count,\n",
    ")\n",
    "ax.set_title(\"Manila 2020 LULC\")\n",
    "\n",
    "cbar = plt.colorbar(p)\n",
    "cbar.set_ticks(range(class_count))\n",
    "cbar.set_ticklabels(class_names)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e742d47e-dda6-4bdf-a441-854f23b3f752",
   "metadata": {},
   "source": [
    "That looks better. Let's also plot a histogram of the pixel values to see the distribution of land cover types within our area of interest. We can reuse the colormap we generated to help tie the two visualizations together."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "746155c1-1a24-4696-b2f4-ae2deef5990e",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<img src=\"https://ai4edatasetspublicassets.blob.core.windows.net/assets/notebook-output/datasets-io-lulc-io-lulc-example.ipynb/11.png\"/>"
      ],
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "colors = list(cmap.colors)\n",
    "\n",
    "ax = (\n",
    "    pd.value_counts(merged.data.ravel(), sort=False)\n",
    "    .sort_index()\n",
    "    .reindex(range(len(colors)), fill_value=0)\n",
    "    .rename(dict(enumerate(class_names)))\n",
    "    .plot.barh(color=colors, rot=0, width=0.9)\n",
    ")\n",
    "ax.set(\n",
    "    title=\"Distribution of Land Cover classes\",\n",
    "    ylabel=\"Landcover class\",\n",
    "    xlabel=\"Pixel count\",\n",
    ");"
   ]
  }
 ],
 "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.6"
  },
  "widgets": {
   "application/vnd.jupyter.widget-state+json": {
    "state": {
     "32d4bb24f892449da58da91bb003e193": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.5.0",
      "model_name": "VBoxModel",
      "state": {
       "layout": "IPY_MODEL_5d442c2da9f24247b080b5a6ea045460"
      }
     },
     "5d442c2da9f24247b080b5a6ea045460": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "1.2.0",
      "model_name": "LayoutModel",
      "state": {}
     }
    },
    "version_major": 2,
    "version_minor": 0
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}