summary.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "legislative-triumph",
   "metadata": {},
   "outputs": [],
   "source": [
    "#!pip install -q sat-search stackstac stac-vrt rioxarray==0.3.1"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "committed-fence",
   "metadata": {},
   "source": [
    "## STAC + Dask\n",
    "\n",
    "STAC and Dask pair nicely. Dask works best when it can build up a large chain of tasks and then execute that computation later on. When you're operating on the xarray DataArray the actual computation is lazy. But we *do* need some kind of \"schema\" to build up the original DataArray in the first place. Things like the sizes, coordinates, dimension names, etc.\n",
    "\n",
    "You *can* do this with a bunch of COGs in blob storage, it's just slow. Whether your're stitching together a bunch of COGs are stacking them over time, you'll need to open each COG to read its metadata (CRS, bands, bounding box, transform, etc.). Opening each COG from Blob Storage involves *at least* one HTTP request, often more.\n",
    "\n",
    "The neat opportunity afforded by STAC is to get all of that metadata up front with just one or a handful of HTTP requests. We can construct the DataArray using just information gleaned from STAC, without opening a single COG. Then when it's time to do the actual computation, we start opening COGs, reading data, and operating on it.\n",
    "\n",
    "So that's the opportunity: **efficiently create DataArrays from metadata in a STAC ItemCollection, delaying opening the COGs until we read the actual data.**."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "moderate-softball",
   "metadata": {},
   "source": [
    "## Comparison between stackstac and stac-vrt\n",
    "\n",
    "Both [`stackstac`](https://stackstac.readthedocs.io/) and [`stac-vrt`](https://stac-vrt.readthedocs.io/en/latest/) are exploiting this opportunity. This notebook explores some of their similarities and differences.\n",
    "\n",
    "### stackstac example\n",
    "\n",
    "`stackstac` exposes a single function, `stackstac.stack`, which -- you guessed it -- stacks a STAC ItemCollection. It's currently able to stack a *timeseries* of STAC items covering the roughly the same spatial extent.\n",
    "\n",
    "![](https://www.esri.com/arcgis-blog/wp-content/uploads/2021/02/timeseries_imagestack.jpg)\n",
    "\n",
    "(image from https://www.esri.com/arcgis-blog/products/arcgis-pro/imagery/change-detection-of-time-series-imagery/)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "authentic-bedroom",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "13"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# https://stackstac.readthedocs.io/en/latest/basic.html\n",
    "import stackstac\n",
    "import satsearch\n",
    "\n",
    "lon, lat = -105.78, 35.79\n",
    "\n",
    "items = satsearch.Search(\n",
    "    url=\"https://earth-search.aws.element84.com/v0\",\n",
    "    intersects=dict(type=\"Point\", coordinates=[lon, lat]),\n",
    "    collections=[\"sentinel-s2-l2a-cogs\"],\n",
    "    datetime=\"2020-04-01/2020-05-01\"\n",
    ").items()\n",
    "len(items)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "interested-nightmare",
   "metadata": {},
   "source": [
    "There are 13 observations (time periods), each with 17 bands. They can be stacked into a 4-D cube with `stackstac.stack`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "improved-relationship",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block;\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-9918e037162926aaa5ec5d240d1a9174&#x27; (time: 13, band: 17, y: 10980, x: 10980)&gt;\n",
       "dask.array&lt;fetch_raster_window, shape=(13, 17, 10980, 10980), dtype=float64, chunksize=(1, 1, 1024, 1024), chunktype=numpy.ndarray&gt;\n",
       "Coordinates: (12/23)\n",
       "  * time                        (time) datetime64[ns] 2020-04-01T18:04:04 ......\n",
       "    id                          (time) &lt;U24 &#x27;S2B_13SDV_20200401_0_L2A&#x27; ... &#x27;S...\n",
       "  * band                        (band) &lt;U8 &#x27;overview&#x27; &#x27;visual&#x27; ... &#x27;WVP&#x27; &#x27;SCL&#x27;\n",
       "  * x                           (x) float64 4e+05 4e+05 ... 5.097e+05 5.098e+05\n",
       "  * y                           (y) float64 4e+06 4e+06 ... 3.89e+06 3.89e+06\n",
       "    created                     (time) &lt;U24 &#x27;2020-09-05T06:23:47.836Z&#x27; ... &#x27;2...\n",
       "    ...                          ...\n",
       "    sentinel:valid_cloud_cover  bool True\n",
       "    sentinel:data_coverage      (time) float64 33.85 100.0 33.9 ... 100.0 34.29\n",
       "    constellation               &lt;U10 &#x27;sentinel-2&#x27;\n",
       "    platform                    (time) &lt;U11 &#x27;sentinel-2b&#x27; ... &#x27;sentinel-2b&#x27;\n",
       "    proj:epsg                   int64 32613\n",
       "    title                       (band) object None ... &#x27;Scene Classification ...\n",
       "Attributes:\n",
       "    spec:        RasterSpec(epsg=32613, bounds=(399960.0, 3890220.0, 509760.0...\n",
       "    crs:         epsg:32613\n",
       "    transform:   | 10.00, 0.00, 399960.00|\\n| 0.00,-10.00, 4000020.00|\\n| 0.0...\n",
       "    resolution:  10.0</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'>'stackstac-9918e037162926aaa5ec5d240d1a9174'</div><ul class='xr-dim-list'><li><span class='xr-has-index'>time</span>: 13</li><li><span class='xr-has-index'>band</span>: 17</li><li><span class='xr-has-index'>y</span>: 10980</li><li><span class='xr-has-index'>x</span>: 10980</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-1e6744b4-3522-42cf-b245-54a2b3ca839f' class='xr-array-in' type='checkbox' checked><label for='section-1e6744b4-3522-42cf-b245-54a2b3ca839f' 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><td> </td><th> Array </th><th> Chunk </th></tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr><th> Bytes </th><td> 213.15 GB </td> <td> 8.39 MB </td></tr>\n",
       "    <tr><th> Shape </th><td> (13, 17, 10980, 10980) </td> <td> (1, 1, 1024, 1024) </td></tr>\n",
       "    <tr><th> Count </th><td> 27304 Tasks </td><td> 26741 Chunks </td></tr>\n",
       "    <tr><th> Type </th><td> float64 </td><td> numpy.ndarray </td></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=\"1\" y1=\"0\" x2=\"1\" y2=\"25\" />\n",
       "  <line x1=\"3\" y1=\"0\" x2=\"3\" y2=\"25\" />\n",
       "  <line x1=\"5\" y1=\"0\" x2=\"5\" y2=\"25\" />\n",
       "  <line x1=\"7\" y1=\"0\" x2=\"7\" y2=\"25\" />\n",
       "  <line x1=\"9\" y1=\"0\" x2=\"9\" y2=\"25\" />\n",
       "  <line x1=\"11\" y1=\"0\" x2=\"11\" y2=\"25\" />\n",
       "  <line x1=\"13\" y1=\"0\" x2=\"13\" y2=\"25\" />\n",
       "  <line x1=\"15\" y1=\"0\" x2=\"15\" y2=\"25\" />\n",
       "  <line x1=\"17\" y1=\"0\" x2=\"17\" y2=\"25\" />\n",
       "  <line x1=\"19\" y1=\"0\" x2=\"19\" y2=\"25\" />\n",
       "  <line x1=\"21\" y1=\"0\" x2=\"21\" y2=\"25\" />\n",
       "  <line x1=\"23\" y1=\"0\" x2=\"23\" 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\" >13</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=\"11\" x2=\"109\" y2=\"26\" />\n",
       "  <line x1=\"95\" y1=\"22\" x2=\"109\" y2=\"37\" />\n",
       "  <line x1=\"95\" y1=\"33\" x2=\"109\" y2=\"48\" />\n",
       "  <line x1=\"95\" y1=\"44\" x2=\"109\" y2=\"59\" />\n",
       "  <line x1=\"95\" y1=\"55\" x2=\"109\" y2=\"70\" />\n",
       "  <line x1=\"95\" y1=\"67\" x2=\"109\" y2=\"82\" />\n",
       "  <line x1=\"95\" y1=\"78\" x2=\"109\" y2=\"93\" />\n",
       "  <line x1=\"95\" y1=\"89\" x2=\"109\" y2=\"104\" />\n",
       "  <line x1=\"95\" y1=\"100\" x2=\"109\" y2=\"115\" />\n",
       "  <line x1=\"95\" y1=\"111\" x2=\"109\" y2=\"126\" />\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=\"95\" y1=\"0\" x2=\"95\" y2=\"120\" />\n",
       "  <line x1=\"96\" y1=\"1\" x2=\"96\" y2=\"121\" />\n",
       "  <line x1=\"97\" y1=\"2\" x2=\"97\" y2=\"122\" />\n",
       "  <line x1=\"98\" y1=\"3\" x2=\"98\" y2=\"123\" />\n",
       "  <line x1=\"99\" y1=\"4\" x2=\"99\" y2=\"124\" />\n",
       "  <line x1=\"100\" y1=\"5\" x2=\"100\" y2=\"125\" />\n",
       "  <line x1=\"101\" y1=\"6\" x2=\"101\" y2=\"126\" />\n",
       "  <line x1=\"102\" y1=\"7\" x2=\"102\" y2=\"127\" />\n",
       "  <line x1=\"102\" y1=\"7\" x2=\"102\" y2=\"127\" />\n",
       "  <line x1=\"103\" y1=\"8\" x2=\"103\" y2=\"128\" />\n",
       "  <line x1=\"104\" y1=\"9\" x2=\"104\" y2=\"129\" />\n",
       "  <line x1=\"105\" y1=\"10\" x2=\"105\" y2=\"130\" />\n",
       "  <line x1=\"106\" y1=\"11\" x2=\"106\" y2=\"131\" />\n",
       "  <line x1=\"107\" y1=\"12\" x2=\"107\" y2=\"132\" />\n",
       "  <line x1=\"108\" y1=\"13\" x2=\"108\" y2=\"133\" />\n",
       "  <line x1=\"109\" y1=\"14\" x2=\"109\" y2=\"134\" />\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=\"95\" y1=\"0\" x2=\"215\" y2=\"0\" />\n",
       "  <line x1=\"96\" y1=\"1\" x2=\"216\" y2=\"1\" />\n",
       "  <line x1=\"97\" y1=\"2\" x2=\"217\" y2=\"2\" />\n",
       "  <line x1=\"98\" y1=\"3\" x2=\"218\" y2=\"3\" />\n",
       "  <line x1=\"99\" y1=\"4\" x2=\"219\" y2=\"4\" />\n",
       "  <line x1=\"100\" y1=\"5\" x2=\"220\" y2=\"5\" />\n",
       "  <line x1=\"101\" y1=\"6\" x2=\"221\" y2=\"6\" />\n",
       "  <line x1=\"102\" y1=\"7\" x2=\"222\" y2=\"7\" />\n",
       "  <line x1=\"102\" y1=\"7\" x2=\"222\" y2=\"7\" />\n",
       "  <line x1=\"103\" y1=\"8\" x2=\"223\" y2=\"8\" />\n",
       "  <line x1=\"104\" y1=\"9\" x2=\"224\" y2=\"9\" />\n",
       "  <line x1=\"105\" y1=\"10\" x2=\"225\" y2=\"10\" />\n",
       "  <line x1=\"106\" y1=\"11\" x2=\"226\" y2=\"11\" />\n",
       "  <line x1=\"107\" y1=\"12\" x2=\"227\" y2=\"12\" />\n",
       "  <line x1=\"108\" y1=\"13\" x2=\"228\" y2=\"13\" />\n",
       "  <line x1=\"109\" y1=\"14\" x2=\"229\" y2=\"14\" />\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=\"106\" y1=\"0\" x2=\"121\" y2=\"14\" />\n",
       "  <line x1=\"117\" y1=\"0\" x2=\"132\" y2=\"14\" />\n",
       "  <line x1=\"128\" y1=\"0\" x2=\"143\" y2=\"14\" />\n",
       "  <line x1=\"139\" y1=\"0\" x2=\"154\" y2=\"14\" />\n",
       "  <line x1=\"150\" y1=\"0\" x2=\"165\" y2=\"14\" />\n",
       "  <line x1=\"162\" y1=\"0\" x2=\"177\" y2=\"14\" />\n",
       "  <line x1=\"173\" y1=\"0\" x2=\"188\" y2=\"14\" />\n",
       "  <line x1=\"184\" y1=\"0\" x2=\"199\" y2=\"14\" />\n",
       "  <line x1=\"195\" y1=\"0\" x2=\"210\" y2=\"14\" />\n",
       "  <line x1=\"206\" y1=\"0\" x2=\"221\" 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=\"26\" x2=\"229\" y2=\"26\" />\n",
       "  <line x1=\"109\" y1=\"37\" x2=\"229\" y2=\"37\" />\n",
       "  <line x1=\"109\" y1=\"48\" x2=\"229\" y2=\"48\" />\n",
       "  <line x1=\"109\" y1=\"59\" x2=\"229\" y2=\"59\" />\n",
       "  <line x1=\"109\" y1=\"70\" x2=\"229\" y2=\"70\" />\n",
       "  <line x1=\"109\" y1=\"82\" x2=\"229\" y2=\"82\" />\n",
       "  <line x1=\"109\" y1=\"93\" x2=\"229\" y2=\"93\" />\n",
       "  <line x1=\"109\" y1=\"104\" x2=\"229\" y2=\"104\" />\n",
       "  <line x1=\"109\" y1=\"115\" x2=\"229\" y2=\"115\" />\n",
       "  <line x1=\"109\" y1=\"126\" x2=\"229\" y2=\"126\" />\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=\"121\" y1=\"14\" x2=\"121\" y2=\"134\" />\n",
       "  <line x1=\"132\" y1=\"14\" x2=\"132\" y2=\"134\" />\n",
       "  <line x1=\"143\" y1=\"14\" x2=\"143\" y2=\"134\" />\n",
       "  <line x1=\"154\" y1=\"14\" x2=\"154\" y2=\"134\" />\n",
       "  <line x1=\"165\" y1=\"14\" x2=\"165\" y2=\"134\" />\n",
       "  <line x1=\"177\" y1=\"14\" x2=\"177\" y2=\"134\" />\n",
       "  <line x1=\"188\" y1=\"14\" x2=\"188\" y2=\"134\" />\n",
       "  <line x1=\"199\" y1=\"14\" x2=\"199\" y2=\"134\" />\n",
       "  <line x1=\"210\" y1=\"14\" x2=\"210\" y2=\"134\" />\n",
       "  <line x1=\"221\" y1=\"14\" x2=\"221\" y2=\"134\" />\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\" >10980</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)\">10980</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)\">17</text>\n",
       "</svg>\n",
       "</td>\n",
       "</tr>\n",
       "</table></div></div></li><li class='xr-section-item'><input id='section-e168544a-9f45-47ac-86ce-fe1f60fcdfae' class='xr-section-summary-in' type='checkbox'  checked><label for='section-e168544a-9f45-47ac-86ce-fe1f60fcdfae' class='xr-section-summary' >Coordinates: <span>(23)</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-04-01T18:04:04 ... 2020-05-...</div><input id='attrs-51c6dfd3-6825-43de-86e8-51b6b810af3b' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-51c6dfd3-6825-43de-86e8-51b6b810af3b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f5c06d4c-6030-4029-86f7-518405064c4d' class='xr-var-data-in' type='checkbox'><label for='data-f5c06d4c-6030-4029-86f7-518405064c4d' 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-04-01T18:04:04.000000000&#x27;, &#x27;2020-04-03T17:54:07.000000000&#x27;,\n",
       "       &#x27;2020-04-06T18:04:04.000000000&#x27;, &#x27;2020-04-08T17:54:08.000000000&#x27;,\n",
       "       &#x27;2020-04-11T18:04:03.000000000&#x27;, &#x27;2020-04-13T17:54:10.000000000&#x27;,\n",
       "       &#x27;2020-04-16T18:04:06.000000000&#x27;, &#x27;2020-04-18T17:54:06.000000000&#x27;,\n",
       "       &#x27;2020-04-21T18:04:01.000000000&#x27;, &#x27;2020-04-23T17:54:12.000000000&#x27;,\n",
       "       &#x27;2020-04-26T18:04:09.000000000&#x27;, &#x27;2020-04-28T17:54:06.000000000&#x27;,\n",
       "       &#x27;2020-05-01T18:04:03.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;U24</div><div class='xr-var-preview xr-preview'>&#x27;S2B_13SDV_20200401_0_L2A&#x27; ... &#x27;...</div><input id='attrs-72f494a0-a288-41dd-9f6a-4cb3c27a8678' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-72f494a0-a288-41dd-9f6a-4cb3c27a8678' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-a2e443af-18e9-43ef-bc91-808855a58b9c' class='xr-var-data-in' type='checkbox'><label for='data-a2e443af-18e9-43ef-bc91-808855a58b9c' 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_13SDV_20200401_0_L2A&#x27;, &#x27;S2A_13SDV_20200403_0_L2A&#x27;,\n",
       "       &#x27;S2A_13SDV_20200406_0_L2A&#x27;, &#x27;S2B_13SDV_20200408_0_L2A&#x27;,\n",
       "       &#x27;S2B_13SDV_20200411_0_L2A&#x27;, &#x27;S2A_13SDV_20200413_0_L2A&#x27;,\n",
       "       &#x27;S2A_13SDV_20200416_0_L2A&#x27;, &#x27;S2B_13SDV_20200418_0_L2A&#x27;,\n",
       "       &#x27;S2B_13SDV_20200421_0_L2A&#x27;, &#x27;S2A_13SDV_20200423_0_L2A&#x27;,\n",
       "       &#x27;S2A_13SDV_20200426_0_L2A&#x27;, &#x27;S2B_13SDV_20200428_0_L2A&#x27;,\n",
       "       &#x27;S2B_13SDV_20200501_0_L2A&#x27;], dtype=&#x27;&lt;U24&#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;U8</div><div class='xr-var-preview xr-preview'>&#x27;overview&#x27; &#x27;visual&#x27; ... &#x27;WVP&#x27; &#x27;SCL&#x27;</div><input id='attrs-3dd218b0-32e7-478c-99e4-aebc95cb7afe' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-3dd218b0-32e7-478c-99e4-aebc95cb7afe' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e75a28a1-9eab-4e65-8a94-7ccd19f67cc4' class='xr-var-data-in' type='checkbox'><label for='data-e75a28a1-9eab-4e65-8a94-7ccd19f67cc4' 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;overview&#x27;, &#x27;visual&#x27;, &#x27;B01&#x27;, &#x27;B02&#x27;, &#x27;B03&#x27;, &#x27;B04&#x27;, &#x27;B05&#x27;, &#x27;B06&#x27;, &#x27;B07&#x27;,\n",
       "       &#x27;B08&#x27;, &#x27;B8A&#x27;, &#x27;B09&#x27;, &#x27;B11&#x27;, &#x27;B12&#x27;, &#x27;AOT&#x27;, &#x27;WVP&#x27;, &#x27;SCL&#x27;], dtype=&#x27;&lt;U8&#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'>4e+05 4e+05 ... 5.097e+05 5.098e+05</div><input id='attrs-690ffdbb-bada-497c-91a9-e04ff21edff3' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-690ffdbb-bada-497c-91a9-e04ff21edff3' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f35c0231-0fb9-4bfd-911c-328276d38f8e' class='xr-var-data-in' type='checkbox'><label for='data-f35c0231-0fb9-4bfd-911c-328276d38f8e' 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([399965., 399975., 399985., ..., 509735., 509745., 509755.])</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'>4e+06 4e+06 ... 3.89e+06 3.89e+06</div><input id='attrs-8dd163f1-aafc-4738-b5ff-bc7da0cf7647' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-8dd163f1-aafc-4738-b5ff-bc7da0cf7647' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b7197cbb-8f07-4fd8-a4b7-057549a4e3b1' class='xr-var-data-in' type='checkbox'><label for='data-b7197cbb-8f07-4fd8-a4b7-057549a4e3b1' 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([4000025., 4000015., 4000005., ..., 3890255., 3890245., 3890235.])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>created</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>&lt;U24</div><div class='xr-var-preview xr-preview'>&#x27;2020-09-05T06:23:47.836Z&#x27; ... &#x27;...</div><input id='attrs-a8361830-41e5-4f56-8969-88935214f463' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-a8361830-41e5-4f56-8969-88935214f463' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-908193ae-bb10-4812-a5ff-170a470a5448' class='xr-var-data-in' type='checkbox'><label for='data-908193ae-bb10-4812-a5ff-170a470a5448' 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-09-05T06:23:47.836Z&#x27;, &#x27;2020-09-23T15:43:10.502Z&#x27;,\n",
       "       &#x27;2020-09-05T11:18:48.722Z&#x27;, &#x27;2020-08-31T15:23:56.855Z&#x27;,\n",
       "       &#x27;2020-09-23T16:16:43.693Z&#x27;, &#x27;2020-09-24T06:31:45.567Z&#x27;,\n",
       "       &#x27;2020-09-24T05:31:29.634Z&#x27;, &#x27;2020-09-19T07:13:25.290Z&#x27;,\n",
       "       &#x27;2020-09-05T11:20:20.745Z&#x27;, &#x27;2020-08-23T09:50:31.532Z&#x27;,\n",
       "       &#x27;2020-09-05T13:56:08.881Z&#x27;, &#x27;2020-08-31T13:13:14.956Z&#x27;,\n",
       "       &#x27;2020-09-24T07:38:40.841Z&#x27;], dtype=&#x27;&lt;U24&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>data_coverage</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>33.85 100 33.9 ... 32.84 100 34.29</div><input id='attrs-3848931f-c5aa-4804-a470-333bd5d8c81d' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-3848931f-c5aa-4804-a470-333bd5d8c81d' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-d789fa44-8838-4fe2-ae0c-0e1a9eeeb6f7' class='xr-var-data-in' type='checkbox'><label for='data-d789fa44-8838-4fe2-ae0c-0e1a9eeeb6f7' 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([33.85, 100, 33.9, 100, 33.87, 100, 33.37, 100, 34.09, None, 32.84,\n",
       "       100, 34.29], dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>sentinel:sequence</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U1</div><div class='xr-var-preview xr-preview'>&#x27;0&#x27;</div><input id='attrs-04277e97-0c15-4888-bd28-b5a3d90318dd' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-04277e97-0c15-4888-bd28-b5a3d90318dd' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-af4d7c99-ad8f-45b0-b780-9a6a541648aa' class='xr-var-data-in' type='checkbox'><label for='data-af4d7c99-ad8f-45b0-b780-9a6a541648aa' 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;0&#x27;, dtype=&#x27;&lt;U1&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>updated</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>&lt;U24</div><div class='xr-var-preview xr-preview'>&#x27;2020-09-05T06:23:47.836Z&#x27; ... &#x27;...</div><input id='attrs-532ade57-efee-46fa-a609-ffd2bb6b3308' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-532ade57-efee-46fa-a609-ffd2bb6b3308' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-dfb77f2a-ffb9-4554-b907-b16410b3b39d' class='xr-var-data-in' type='checkbox'><label for='data-dfb77f2a-ffb9-4554-b907-b16410b3b39d' 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-09-05T06:23:47.836Z&#x27;, &#x27;2020-09-23T15:43:10.502Z&#x27;,\n",
       "       &#x27;2020-09-05T11:18:48.722Z&#x27;, &#x27;2020-08-31T15:23:56.855Z&#x27;,\n",
       "       &#x27;2020-09-23T16:16:43.693Z&#x27;, &#x27;2020-09-24T06:31:45.567Z&#x27;,\n",
       "       &#x27;2020-09-24T05:31:29.634Z&#x27;, &#x27;2020-09-19T07:13:25.290Z&#x27;,\n",
       "       &#x27;2020-09-05T11:20:20.745Z&#x27;, &#x27;2020-08-23T09:50:31.532Z&#x27;,\n",
       "       &#x27;2020-09-05T13:56:08.881Z&#x27;, &#x27;2020-08-31T13:13:14.956Z&#x27;,\n",
       "       &#x27;2020-09-24T07:38:40.841Z&#x27;], dtype=&#x27;&lt;U24&#x27;)</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'>29.24 1.16 27.26 ... 87.33 5.41</div><input id='attrs-4eabfe80-26bc-4b12-b40e-a50c32c302ed' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-4eabfe80-26bc-4b12-b40e-a50c32c302ed' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-07373d5e-2024-4c67-b08e-9d348466330f' class='xr-var-data-in' type='checkbox'><label for='data-07373d5e-2024-4c67-b08e-9d348466330f' 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([ 29.24,   1.16,  27.26,   1.15,   9.37, 100.  ,   4.53,  35.  ,\n",
       "        14.97,  31.78,  99.42,  87.33,   5.41])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>sentinel:latitude_band</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U1</div><div class='xr-var-preview xr-preview'>&#x27;S&#x27;</div><input id='attrs-72487742-2db8-444d-aa0b-67febc2b280f' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-72487742-2db8-444d-aa0b-67febc2b280f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-02c77422-a4f9-48ab-864d-222a238c41f1' class='xr-var-data-in' type='checkbox'><label for='data-02c77422-a4f9-48ab-864d-222a238c41f1' 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;S&#x27;, dtype=&#x27;&lt;U1&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>sentinel:utm_zone</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>13</div><input id='attrs-29286eba-106d-4dc6-8722-4fb51773928d' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-29286eba-106d-4dc6-8722-4fb51773928d' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-959efeef-b072-49f3-bc7a-31b20d411472' class='xr-var-data-in' type='checkbox'><label for='data-959efeef-b072-49f3-bc7a-31b20d411472' 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(13)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>sentinel:product_id</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>&lt;U60</div><div class='xr-var-preview xr-preview'>&#x27;S2B_MSIL2A_20200401T174909_N021...</div><input id='attrs-51dd4707-c8f5-45bb-aac4-1a67018b5a04' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-51dd4707-c8f5-45bb-aac4-1a67018b5a04' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4ccc8464-aa70-4536-bcc7-36f856645178' class='xr-var-data-in' type='checkbox'><label for='data-4ccc8464-aa70-4536-bcc7-36f856645178' 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_20200401T174909_N0214_R141_T13SDV_20200401T220155&#x27;,\n",
       "       &#x27;S2A_MSIL2A_20200403T173901_N0214_R098_T13SDV_20200403T220105&#x27;,\n",
       "       &#x27;S2A_MSIL2A_20200406T174901_N0214_R141_T13SDV_20200406T221027&#x27;,\n",
       "       &#x27;S2B_MSIL2A_20200408T173859_N0214_R098_T13SDV_20200408T215856&#x27;,\n",
       "       &#x27;S2B_MSIL2A_20200411T174909_N0214_R141_T13SDV_20200411T220443&#x27;,\n",
       "       &#x27;S2A_MSIL2A_20200413T173901_N0214_R098_T13SDV_20200413T235616&#x27;,\n",
       "       &#x27;S2A_MSIL2A_20200416T174911_N0214_R141_T13SDV_20200416T221101&#x27;,\n",
       "       &#x27;S2B_MSIL2A_20200418T173859_N0214_R098_T13SDV_20200418T214953&#x27;,\n",
       "       &#x27;S2B_MSIL2A_20200421T174859_N0214_R141_T13SDV_20200421T214716&#x27;,\n",
       "       &#x27;S2A_MSIL2A_20200423T173911_N0214_R098_T13SDV_20200423T234141&#x27;,\n",
       "       &#x27;S2A_MSIL2A_20200426T174911_N0214_R141_T13SDV_20200426T221124&#x27;,\n",
       "       &#x27;S2B_MSIL2A_20200428T173859_N0214_R098_T13SDV_20200428T214309&#x27;,\n",
       "       &#x27;S2B_MSIL2A_20200501T174909_N0214_R141_T13SDV_20200501T220131&#x27;],\n",
       "      dtype=&#x27;&lt;U60&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>gsd</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>10</div><input id='attrs-5d002de8-bc14-4764-80a5-ee1db2655bca' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-5d002de8-bc14-4764-80a5-ee1db2655bca' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-88ad39fc-fd52-4b23-bc0d-a2e2558c0795' class='xr-var-data-in' type='checkbox'><label for='data-88ad39fc-fd52-4b23-bc0d-a2e2558c0795' 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)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>view:off_nadir</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>0</div><input id='attrs-67bbce9a-d751-407d-97f3-d2ff71cbf549' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-67bbce9a-d751-407d-97f3-d2ff71cbf549' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-ea374bc7-6fac-4147-9e63-8858b46bedb1' class='xr-var-data-in' type='checkbox'><label for='data-ea374bc7-6fac-4147-9e63-8858b46bedb1' 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>sentinel:grid_square</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U2</div><div class='xr-var-preview xr-preview'>&#x27;DV&#x27;</div><input id='attrs-e66f4a3e-14c4-401b-a3d9-37a7e9c6950e' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-e66f4a3e-14c4-401b-a3d9-37a7e9c6950e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c05d3a3f-9e67-478d-b2c1-46dea0458da8' class='xr-var-data-in' type='checkbox'><label for='data-c05d3a3f-9e67-478d-b2c1-46dea0458da8' 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;DV&#x27;, dtype=&#x27;&lt;U2&#x27;)</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-583bcbff-4f28-4ee9-85c2-fdb20a1b0007' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-583bcbff-4f28-4ee9-85c2-fdb20a1b0007' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-bf23ce22-c12d-4c8a-bfe8-882348dfe0f7' class='xr-var-data-in' type='checkbox'><label for='data-bf23ce22-c12d-4c8a-bfe8-882348dfe0f7' 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>sentinel:valid_cloud_cover</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>bool</div><div class='xr-var-preview xr-preview'>True</div><input id='attrs-47e99674-6e50-4f00-97fa-378112772cd1' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-47e99674-6e50-4f00-97fa-378112772cd1' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-9a622a27-6189-42fb-a2a3-3bbd62ec23ba' class='xr-var-data-in' type='checkbox'><label for='data-9a622a27-6189-42fb-a2a3-3bbd62ec23ba' 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(True)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>sentinel:data_coverage</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>33.85 100.0 33.9 ... 100.0 34.29</div><input id='attrs-1d8c450f-ab79-4103-8f8b-1755f68cf691' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-1d8c450f-ab79-4103-8f8b-1755f68cf691' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-583404f1-2e3e-4caa-8ddd-b89bb3931f60' class='xr-var-data-in' type='checkbox'><label for='data-583404f1-2e3e-4caa-8ddd-b89bb3931f60' 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([ 33.85, 100.  ,  33.9 , 100.  ,  33.87, 100.  ,  33.37, 100.  ,\n",
       "        34.09, 100.  ,  32.84, 100.  ,  34.29])</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-ed0117b7-bf79-4b23-a022-2de858f636b6' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-ed0117b7-bf79-4b23-a022-2de858f636b6' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c66172f8-6be8-4bab-b28e-19313bf05a51' class='xr-var-data-in' type='checkbox'><label for='data-c66172f8-6be8-4bab-b28e-19313bf05a51' 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>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-2b&#x27;</div><input id='attrs-f662dc17-3ff6-4b30-bd96-e56d0c915117' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-f662dc17-3ff6-4b30-bd96-e56d0c915117' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-14335e17-b1a5-42d5-95e5-4d1c871f069f' class='xr-var-data-in' type='checkbox'><label for='data-14335e17-b1a5-42d5-95e5-4d1c871f069f' 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;, &#x27;sentinel-2a&#x27;, &#x27;sentinel-2b&#x27;,\n",
       "       &#x27;sentinel-2b&#x27;, &#x27;sentinel-2a&#x27;, &#x27;sentinel-2a&#x27;, &#x27;sentinel-2b&#x27;,\n",
       "       &#x27;sentinel-2b&#x27;, &#x27;sentinel-2a&#x27;, &#x27;sentinel-2a&#x27;, &#x27;sentinel-2b&#x27;,\n",
       "       &#x27;sentinel-2b&#x27;], dtype=&#x27;&lt;U11&#x27;)</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'>32613</div><input id='attrs-3e84f3b1-0d6c-49c5-bb24-17a8942ca8df' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-3e84f3b1-0d6c-49c5-bb24-17a8942ca8df' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-394932e5-f912-4700-b42f-4b59b37c1937' class='xr-var-data-in' type='checkbox'><label for='data-394932e5-f912-4700-b42f-4b59b37c1937' 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(32613)</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'>object</div><div class='xr-var-preview xr-preview'>None ... &#x27;Scene Classification M...</div><input id='attrs-8a1dc4c7-5554-42a1-9fc1-a6c3eebe2700' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-8a1dc4c7-5554-42a1-9fc1-a6c3eebe2700' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-1ef5e37e-14df-4d6b-91ad-6e989d6cbf54' class='xr-var-data-in' type='checkbox'><label for='data-1ef5e37e-14df-4d6b-91ad-6e989d6cbf54' 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, &#x27;True color image&#x27;, &#x27;Band 1 (coastal)&#x27;, &#x27;Band 2 (blue)&#x27;,\n",
       "       &#x27;Band 3 (green)&#x27;, &#x27;Band 4 (red)&#x27;, &#x27;Band 5&#x27;, &#x27;Band 6&#x27;, &#x27;Band 7&#x27;,\n",
       "       &#x27;Band 8 (nir)&#x27;, &#x27;Band 8A&#x27;, &#x27;Band 9&#x27;, &#x27;Band 11 (swir16)&#x27;,\n",
       "       &#x27;Band 12 (swir22)&#x27;, &#x27;Aerosol Optical Thickness (AOT)&#x27;,\n",
       "       &#x27;Water Vapour (WVP)&#x27;, &#x27;Scene Classification Map (SCL)&#x27;],\n",
       "      dtype=object)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-2bb016ca-32dc-4de9-b0ab-cbcddfddfb50' class='xr-section-summary-in' type='checkbox'  checked><label for='section-2bb016ca-32dc-4de9-b0ab-cbcddfddfb50' 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=32613, bounds=(399960.0, 3890220.0, 509760.0, 4000020.0), resolutions_xy=(10.0, 10.0))</dd><dt><span>crs :</span></dt><dd>epsg:32613</dd><dt><span>transform :</span></dt><dd>| 10.00, 0.00, 399960.00|\n",
       "| 0.00,-10.00, 4000020.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-9918e037162926aaa5ec5d240d1a9174' (time: 13, band: 17, y: 10980, x: 10980)>\n",
       "dask.array<fetch_raster_window, shape=(13, 17, 10980, 10980), dtype=float64, chunksize=(1, 1, 1024, 1024), chunktype=numpy.ndarray>\n",
       "Coordinates: (12/23)\n",
       "  * time                        (time) datetime64[ns] 2020-04-01T18:04:04 ......\n",
       "    id                          (time) <U24 'S2B_13SDV_20200401_0_L2A' ... 'S...\n",
       "  * band                        (band) <U8 'overview' 'visual' ... 'WVP' 'SCL'\n",
       "  * x                           (x) float64 4e+05 4e+05 ... 5.097e+05 5.098e+05\n",
       "  * y                           (y) float64 4e+06 4e+06 ... 3.89e+06 3.89e+06\n",
       "    created                     (time) <U24 '2020-09-05T06:23:47.836Z' ... '2...\n",
       "    ...                          ...\n",
       "    sentinel:valid_cloud_cover  bool True\n",
       "    sentinel:data_coverage      (time) float64 33.85 100.0 33.9 ... 100.0 34.29\n",
       "    constellation               <U10 'sentinel-2'\n",
       "    platform                    (time) <U11 'sentinel-2b' ... 'sentinel-2b'\n",
       "    proj:epsg                   int64 32613\n",
       "    title                       (band) object None ... 'Scene Classification ...\n",
       "Attributes:\n",
       "    spec:        RasterSpec(epsg=32613, bounds=(399960.0, 3890220.0, 509760.0...\n",
       "    crs:         epsg:32613\n",
       "    transform:   | 10.00, 0.00, 399960.00|\\n| 0.00,-10.00, 4000020.00|\\n| 0.0...\n",
       "    resolution:  10.0"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ds = stackstac.stack(items)\n",
    "ds"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "insured-mayor",
   "metadata": {},
   "source": [
    "### stac-vrt example\n",
    "\n",
    "`stac-vrt` also exposes a single function, `build_vrt` in this case. It takes a STAC ItemCollection and builds a GDAL VRT. It's currently focused on stitching together many images taken at the same time.\n",
    "\n",
    "![](https://desktop.arcgis.com/en/arcmap/10.3/manage-data/raster-and-images/GUID-BCA5B031-B811-424B-9F54-BAB2224FBAD0-web.gif)\n",
    "\n",
    "(image from https://desktop.arcgis.com/en/arcmap/10.3/manage-data/raster-and-images/what-is-a-mosaic.htm)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "labeled-dressing",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "10"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import requests\n",
    "\n",
    "stac_items = requests.get(\n",
    "    url=\"https://pct-pqe-staging.westeurope.cloudapp.azure.com/stac/v1/collections/naip/items\"\n",
    ").json()[\"features\"]\n",
    "len(stac_items)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "native-fireplace",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "CPU times: user 26.5 ms, sys: 0 ns, total: 26.5 ms\n",
      "Wall time: 25.9 ms\n"
     ]
    }
   ],
   "source": [
    "import stac_vrt\n",
    "\n",
    "%time vrt = stac_vrt.build_vrt(stac_items, data_type=\"Byte\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "graduate-syndrome",
   "metadata": {},
   "source": [
    "`vrt` is a string that's a valid GDAL VRT."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "together-contractor",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<VRTDataset rasterXSize=\"83790\" rasterYSize=\"196870\">\n",
      "  <SRS dataAxisToSRSAxisMapping=\"1,2\">PROJCS[\"NAD83 / UTM zone 17N\",GEOGCS[\"NAD83\",DATUM[\"North_American_Datum_1983\",SPHEROID[\"GRS 1980\",6378137,298.257222101,AUTHORITY[\"EPSG\",\"7019\"]],AUTHORITY[\"EPSG\",\"6269\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4269\"]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"latitude_of_origin\",0],PARAMETER[\"central_meridian\",-81],PARAMETER[\"scale_factor\",0.9996],PARAMETER[\"false_easting\",500000],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"Easting\",EAST],AXIS[\"Northing\",NORTH],AUTHORITY[\"EPSG\",\"26917\"]]</SRS>\n",
      "  <GeoTransform>524382.0, 0.6, 0.0, 3097458.0, 0.0, -0.6</GeoTransform>\n",
      "  <VRTRasterBand dataType=\"Byte\" band=\"1\">\n",
      "    <ColorInterp>Red</ColorInterp>\n",
      "    <SimpleSource>\n",
      "      <SourceFilename relativeToVRT=\"0\">/vsicurl/https://naipeuwest.blob.core.windows.net/naip/v002/fl/2019/fl_60cm_2\n"
     ]
    }
   ],
   "source": [
    "print(vrt[:1000])"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "worse-murder",
   "metadata": {},
   "source": [
    "Libraries like rasterio and rioxarray know how to work with VRTs."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "diverse-olive",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block;\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 (band: 4, y: 196870, x: 83790)&gt;\n",
       "dask.array&lt;open_rasterio-2a3a1a5b470aa0bcf12bc69837fce833&lt;this-array&gt;, shape=(4, 196870, 83790), dtype=uint8, chunksize=(1, 4096, 4096), chunktype=numpy.ndarray&gt;\n",
       "Coordinates:\n",
       "  * band         (band) int64 1 2 3 4\n",
       "  * y            (y) float64 3.097e+06 3.097e+06 ... 2.979e+06 2.979e+06\n",
       "  * x            (x) float64 5.244e+05 5.244e+05 ... 5.747e+05 5.747e+05\n",
       "    spatial_ref  int64 0\n",
       "Attributes:\n",
       "    scale_factor:  1.0\n",
       "    add_offset:    0.0\n",
       "    grid_mapping:  spatial_ref</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'></div><ul class='xr-dim-list'><li><span class='xr-has-index'>band</span>: 4</li><li><span class='xr-has-index'>y</span>: 196870</li><li><span class='xr-has-index'>x</span>: 83790</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-0767b795-6ebb-4816-b972-74a343b60e15' class='xr-array-in' type='checkbox' checked><label for='section-0767b795-6ebb-4816-b972-74a343b60e15' 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, 4096, 4096), meta=np.ndarray&gt;</span></div><div class='xr-array-data'><table>\n",
       "<tr>\n",
       "<td>\n",
       "<table>\n",
       "  <thead>\n",
       "    <tr><td> </td><th> Array </th><th> Chunk </th></tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr><th> Bytes </th><td> 65.98 GB </td> <td> 16.78 MB </td></tr>\n",
       "    <tr><th> Shape </th><td> (4, 196870, 83790) </td> <td> (1, 4096, 4096) </td></tr>\n",
       "    <tr><th> Count </th><td> 4117 Tasks </td><td> 4116 Chunks </td></tr>\n",
       "    <tr><th> Type </th><td> uint8 </td><td> numpy.ndarray </td></tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</td>\n",
       "<td>\n",
       "<svg width=\"126\" height=\"184\" style=\"stroke:rgb(0,0,0);stroke-width:1\" >\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"10\" y1=\"0\" x2=\"24\" y2=\"14\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"10\" y1=\"4\" x2=\"24\" y2=\"19\" />\n",
       "  <line x1=\"10\" y1=\"12\" x2=\"24\" y2=\"27\" />\n",
       "  <line x1=\"10\" y1=\"17\" x2=\"24\" y2=\"32\" />\n",
       "  <line x1=\"10\" y1=\"24\" x2=\"24\" y2=\"39\" />\n",
       "  <line x1=\"10\" y1=\"29\" x2=\"24\" y2=\"44\" />\n",
       "  <line x1=\"10\" y1=\"37\" x2=\"24\" y2=\"52\" />\n",
       "  <line x1=\"10\" y1=\"44\" x2=\"24\" y2=\"59\" />\n",
       "  <line x1=\"10\" y1=\"49\" x2=\"24\" y2=\"64\" />\n",
       "  <line x1=\"10\" y1=\"57\" x2=\"24\" y2=\"72\" />\n",
       "  <line x1=\"10\" y1=\"62\" x2=\"24\" y2=\"77\" />\n",
       "  <line x1=\"10\" y1=\"69\" x2=\"24\" y2=\"84\" />\n",
       "  <line x1=\"10\" y1=\"74\" x2=\"24\" y2=\"89\" />\n",
       "  <line x1=\"10\" y1=\"82\" x2=\"24\" y2=\"97\" />\n",
       "  <line x1=\"10\" y1=\"89\" x2=\"24\" y2=\"104\" />\n",
       "  <line x1=\"10\" y1=\"94\" x2=\"24\" y2=\"109\" />\n",
       "  <line x1=\"10\" y1=\"102\" x2=\"24\" y2=\"117\" />\n",
       "  <line x1=\"10\" y1=\"107\" x2=\"24\" y2=\"122\" />\n",
       "  <line x1=\"10\" y1=\"114\" x2=\"24\" y2=\"129\" />\n",
       "  <line x1=\"10\" y1=\"120\" x2=\"24\" y2=\"134\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"10\" y1=\"0\" x2=\"10\" y2=\"120\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"13\" y1=\"3\" x2=\"13\" y2=\"123\" />\n",
       "  <line x1=\"17\" y1=\"7\" x2=\"17\" y2=\"127\" />\n",
       "  <line x1=\"21\" y1=\"11\" x2=\"21\" y2=\"131\" />\n",
       "  <line x1=\"24\" y1=\"14\" x2=\"24\" y2=\"134\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"10.0,0.0 24.9485979497544,14.948597949754403 24.9485979497544,134.9485979497544 10.0,120.0\" style=\"fill:#8B4903A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"10\" y1=\"0\" x2=\"61\" y2=\"0\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"13\" y1=\"3\" x2=\"64\" y2=\"3\" />\n",
       "  <line x1=\"17\" y1=\"7\" x2=\"68\" y2=\"7\" />\n",
       "  <line x1=\"21\" y1=\"11\" x2=\"72\" y2=\"11\" />\n",
       "  <line x1=\"24\" y1=\"14\" x2=\"76\" y2=\"14\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"10\" y1=\"0\" x2=\"24\" y2=\"14\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"12\" y1=\"0\" x2=\"27\" y2=\"14\" />\n",
       "  <line x1=\"14\" y1=\"0\" x2=\"29\" y2=\"14\" />\n",
       "  <line x1=\"17\" y1=\"0\" x2=\"32\" y2=\"14\" />\n",
       "  <line x1=\"19\" y1=\"0\" x2=\"34\" y2=\"14\" />\n",
       "  <line x1=\"22\" y1=\"0\" x2=\"37\" y2=\"14\" />\n",
       "  <line x1=\"24\" y1=\"0\" x2=\"39\" y2=\"14\" />\n",
       "  <line x1=\"27\" y1=\"0\" x2=\"42\" y2=\"14\" />\n",
       "  <line x1=\"29\" y1=\"0\" x2=\"44\" y2=\"14\" />\n",
       "  <line x1=\"32\" y1=\"0\" x2=\"47\" y2=\"14\" />\n",
       "  <line x1=\"37\" y1=\"0\" x2=\"52\" y2=\"14\" />\n",
       "  <line x1=\"39\" y1=\"0\" x2=\"54\" y2=\"14\" />\n",
       "  <line x1=\"42\" y1=\"0\" x2=\"57\" y2=\"14\" />\n",
       "  <line x1=\"44\" y1=\"0\" x2=\"59\" y2=\"14\" />\n",
       "  <line x1=\"47\" y1=\"0\" x2=\"62\" y2=\"14\" />\n",
       "  <line x1=\"49\" y1=\"0\" x2=\"64\" y2=\"14\" />\n",
       "  <line x1=\"52\" y1=\"0\" x2=\"67\" y2=\"14\" />\n",
       "  <line x1=\"54\" y1=\"0\" x2=\"69\" y2=\"14\" />\n",
       "  <line x1=\"57\" y1=\"0\" x2=\"72\" y2=\"14\" />\n",
       "  <line x1=\"61\" y1=\"0\" x2=\"76\" y2=\"14\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"10.0,0.0 61.07329709960888,0.0 76.02189504936328,14.948597949754403 24.9485979497544,14.948597949754403\" style=\"fill:#8B4903A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"24\" y1=\"14\" x2=\"76\" y2=\"14\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"24\" y1=\"19\" x2=\"76\" y2=\"19\" />\n",
       "  <line x1=\"24\" y1=\"27\" x2=\"76\" y2=\"27\" />\n",
       "  <line x1=\"24\" y1=\"32\" x2=\"76\" y2=\"32\" />\n",
       "  <line x1=\"24\" y1=\"39\" x2=\"76\" y2=\"39\" />\n",
       "  <line x1=\"24\" y1=\"44\" x2=\"76\" y2=\"44\" />\n",
       "  <line x1=\"24\" y1=\"52\" x2=\"76\" y2=\"52\" />\n",
       "  <line x1=\"24\" y1=\"59\" x2=\"76\" y2=\"59\" />\n",
       "  <line x1=\"24\" y1=\"64\" x2=\"76\" y2=\"64\" />\n",
       "  <line x1=\"24\" y1=\"72\" x2=\"76\" y2=\"72\" />\n",
       "  <line x1=\"24\" y1=\"77\" x2=\"76\" y2=\"77\" />\n",
       "  <line x1=\"24\" y1=\"84\" x2=\"76\" y2=\"84\" />\n",
       "  <line x1=\"24\" y1=\"89\" x2=\"76\" y2=\"89\" />\n",
       "  <line x1=\"24\" y1=\"97\" x2=\"76\" y2=\"97\" />\n",
       "  <line x1=\"24\" y1=\"104\" x2=\"76\" y2=\"104\" />\n",
       "  <line x1=\"24\" y1=\"109\" x2=\"76\" y2=\"109\" />\n",
       "  <line x1=\"24\" y1=\"117\" x2=\"76\" y2=\"117\" />\n",
       "  <line x1=\"24\" y1=\"122\" x2=\"76\" y2=\"122\" />\n",
       "  <line x1=\"24\" y1=\"129\" x2=\"76\" y2=\"129\" />\n",
       "  <line x1=\"24\" y1=\"134\" x2=\"76\" y2=\"134\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"24\" y1=\"14\" x2=\"24\" y2=\"134\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"27\" y1=\"14\" x2=\"27\" y2=\"134\" />\n",
       "  <line x1=\"29\" y1=\"14\" x2=\"29\" y2=\"134\" />\n",
       "  <line x1=\"32\" y1=\"14\" x2=\"32\" y2=\"134\" />\n",
       "  <line x1=\"34\" y1=\"14\" x2=\"34\" y2=\"134\" />\n",
       "  <line x1=\"37\" y1=\"14\" x2=\"37\" y2=\"134\" />\n",
       "  <line x1=\"39\" y1=\"14\" x2=\"39\" y2=\"134\" />\n",
       "  <line x1=\"42\" y1=\"14\" x2=\"42\" y2=\"134\" />\n",
       "  <line x1=\"44\" y1=\"14\" x2=\"44\" y2=\"134\" />\n",
       "  <line x1=\"47\" y1=\"14\" x2=\"47\" y2=\"134\" />\n",
       "  <line x1=\"52\" y1=\"14\" x2=\"52\" y2=\"134\" />\n",
       "  <line x1=\"54\" y1=\"14\" x2=\"54\" y2=\"134\" />\n",
       "  <line x1=\"57\" y1=\"14\" x2=\"57\" y2=\"134\" />\n",
       "  <line x1=\"59\" y1=\"14\" x2=\"59\" y2=\"134\" />\n",
       "  <line x1=\"62\" y1=\"14\" x2=\"62\" y2=\"134\" />\n",
       "  <line x1=\"64\" y1=\"14\" x2=\"64\" y2=\"134\" />\n",
       "  <line x1=\"67\" y1=\"14\" x2=\"67\" y2=\"134\" />\n",
       "  <line x1=\"69\" y1=\"14\" x2=\"69\" y2=\"134\" />\n",
       "  <line x1=\"72\" y1=\"14\" x2=\"72\" y2=\"134\" />\n",
       "  <line x1=\"76\" y1=\"14\" x2=\"76\" y2=\"134\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"24.9485979497544,14.948597949754403 76.02189504936328,14.948597949754403 76.02189504936328,134.9485979497544 24.9485979497544,134.9485979497544\" style=\"fill:#8B4903A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Text -->\n",
       "  <text x=\"50.485246\" y=\"154.948598\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" >83790</text>\n",
       "  <text x=\"96.021895\" y=\"74.948598\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(-90,96.021895,74.948598)\">196870</text>\n",
       "  <text x=\"7.474299\" y=\"147.474299\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(45,7.474299,147.474299)\">4</text>\n",
       "</svg>\n",
       "</td>\n",
       "</tr>\n",
       "</table></div></div></li><li class='xr-section-item'><input id='section-37896420-b39f-456b-ba23-15ce893e2e49' class='xr-section-summary-in' type='checkbox'  checked><label for='section-37896420-b39f-456b-ba23-15ce893e2e49' class='xr-section-summary' >Coordinates: <span>(4)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>band</span></div><div class='xr-var-dims'>(band)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>1 2 3 4</div><input id='attrs-ca0e6375-f9bb-4db7-acd2-a4a3c2693628' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-ca0e6375-f9bb-4db7-acd2-a4a3c2693628' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-99f8d692-7b39-4914-9c74-b626f76d36ad' class='xr-var-data-in' type='checkbox'><label for='data-99f8d692-7b39-4914-9c74-b626f76d36ad' 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, 2, 3, 4])</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'>3.097e+06 3.097e+06 ... 2.979e+06</div><input id='attrs-461ae844-122c-42d7-92c9-4ae16b840d98' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-461ae844-122c-42d7-92c9-4ae16b840d98' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-7c6d5079-9db0-49de-a432-67db00093d86' class='xr-var-data-in' type='checkbox'><label for='data-7c6d5079-9db0-49de-a432-67db00093d86' 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([3097457.7, 3097457.1, 3097456.5, ..., 2979337.5, 2979336.9, 2979336.3])</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'>5.244e+05 5.244e+05 ... 5.747e+05</div><input id='attrs-12cebff0-4d6d-46b4-9ae5-7308425b3ec4' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-12cebff0-4d6d-46b4-9ae5-7308425b3ec4' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-50901a0e-7b56-4322-80d1-1395a3dec06f' class='xr-var-data-in' type='checkbox'><label for='data-50901a0e-7b56-4322-80d1-1395a3dec06f' 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([524382.3, 524382.9, 524383.5, ..., 574654.5, 574655.1, 574655.7])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>spatial_ref</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>0</div><input id='attrs-aad38d3b-5252-47e5-b185-d2883e690132' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-aad38d3b-5252-47e5-b185-d2883e690132' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-ac935919-c60f-4c28-9364-d8ee3c3d1450' class='xr-var-data-in' type='checkbox'><label for='data-ac935919-c60f-4c28-9364-d8ee3c3d1450' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>crs_wkt :</span></dt><dd>PROJCS[&quot;NAD83 / UTM zone 17N&quot;,GEOGCS[&quot;NAD83&quot;,DATUM[&quot;North_American_Datum_1983&quot;,SPHEROID[&quot;GRS 1980&quot;,6378137,298.257222101,AUTHORITY[&quot;EPSG&quot;,&quot;7019&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6269&quot;]],PRIMEM[&quot;Greenwich&quot;,0,AUTHORITY[&quot;EPSG&quot;,&quot;8901&quot;]],UNIT[&quot;degree&quot;,0.0174532925199433,AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;4269&quot;]],PROJECTION[&quot;Transverse_Mercator&quot;],PARAMETER[&quot;latitude_of_origin&quot;,0],PARAMETER[&quot;central_meridian&quot;,-81],PARAMETER[&quot;scale_factor&quot;,0.9996],PARAMETER[&quot;false_easting&quot;,500000],PARAMETER[&quot;false_northing&quot;,0],UNIT[&quot;metre&quot;,1,AUTHORITY[&quot;EPSG&quot;,&quot;9001&quot;]],AXIS[&quot;Easting&quot;,EAST],AXIS[&quot;Northing&quot;,NORTH],AUTHORITY[&quot;EPSG&quot;,&quot;26917&quot;]]</dd><dt><span>semi_major_axis :</span></dt><dd>6378137.0</dd><dt><span>semi_minor_axis :</span></dt><dd>6356752.314140356</dd><dt><span>inverse_flattening :</span></dt><dd>298.257222101</dd><dt><span>reference_ellipsoid_name :</span></dt><dd>GRS 1980</dd><dt><span>longitude_of_prime_meridian :</span></dt><dd>0.0</dd><dt><span>prime_meridian_name :</span></dt><dd>Greenwich</dd><dt><span>geographic_crs_name :</span></dt><dd>NAD83</dd><dt><span>horizontal_datum_name :</span></dt><dd>North American Datum 1983</dd><dt><span>projected_crs_name :</span></dt><dd>NAD83 / UTM zone 17N</dd><dt><span>grid_mapping_name :</span></dt><dd>transverse_mercator</dd><dt><span>latitude_of_projection_origin :</span></dt><dd>0.0</dd><dt><span>longitude_of_central_meridian :</span></dt><dd>-81.0</dd><dt><span>false_easting :</span></dt><dd>500000.0</dd><dt><span>false_northing :</span></dt><dd>0.0</dd><dt><span>scale_factor_at_central_meridian :</span></dt><dd>0.9996</dd><dt><span>spatial_ref :</span></dt><dd>PROJCS[&quot;NAD83 / UTM zone 17N&quot;,GEOGCS[&quot;NAD83&quot;,DATUM[&quot;North_American_Datum_1983&quot;,SPHEROID[&quot;GRS 1980&quot;,6378137,298.257222101,AUTHORITY[&quot;EPSG&quot;,&quot;7019&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6269&quot;]],PRIMEM[&quot;Greenwich&quot;,0,AUTHORITY[&quot;EPSG&quot;,&quot;8901&quot;]],UNIT[&quot;degree&quot;,0.0174532925199433,AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;4269&quot;]],PROJECTION[&quot;Transverse_Mercator&quot;],PARAMETER[&quot;latitude_of_origin&quot;,0],PARAMETER[&quot;central_meridian&quot;,-81],PARAMETER[&quot;scale_factor&quot;,0.9996],PARAMETER[&quot;false_easting&quot;,500000],PARAMETER[&quot;false_northing&quot;,0],UNIT[&quot;metre&quot;,1,AUTHORITY[&quot;EPSG&quot;,&quot;9001&quot;]],AXIS[&quot;Easting&quot;,EAST],AXIS[&quot;Northing&quot;,NORTH],AUTHORITY[&quot;EPSG&quot;,&quot;26917&quot;]]</dd><dt><span>GeoTransform :</span></dt><dd>524382.0 0.6 0.0 3097458.0 0.0 -0.6</dd></dl></div><div class='xr-var-data'><pre>array(0)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-28244ff5-9bc4-4fe0-ba62-c259e754d023' class='xr-section-summary-in' type='checkbox'  checked><label for='section-28244ff5-9bc4-4fe0-ba62-c259e754d023' class='xr-section-summary' >Attributes: <span>(3)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>scale_factor :</span></dt><dd>1.0</dd><dt><span>add_offset :</span></dt><dd>0.0</dd><dt><span>grid_mapping :</span></dt><dd>spatial_ref</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.DataArray (band: 4, y: 196870, x: 83790)>\n",
       "dask.array<open_rasterio-2a3a1a5b470aa0bcf12bc69837fce833<this-array>, shape=(4, 196870, 83790), dtype=uint8, chunksize=(1, 4096, 4096), chunktype=numpy.ndarray>\n",
       "Coordinates:\n",
       "  * band         (band) int64 1 2 3 4\n",
       "  * y            (y) float64 3.097e+06 3.097e+06 ... 2.979e+06 2.979e+06\n",
       "  * x            (x) float64 5.244e+05 5.244e+05 ... 5.747e+05 5.747e+05\n",
       "    spatial_ref  int64 0\n",
       "Attributes:\n",
       "    scale_factor:  1.0\n",
       "    add_offset:    0.0\n",
       "    grid_mapping:  spatial_ref"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import rioxarray\n",
    "\n",
    "ds = rioxarray.open_rasterio(vrt, chunks=(1, 4096, 4096))\n",
    "ds"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "optional-findings",
   "metadata": {},
   "source": [
    "## Comparision\n",
    "\n",
    "1. **focus**: timeseries at a location vs. single timestep at many locations\n",
    "\n",
    "`stackstac` currently focuses on stacking a stack of STAC items along the `time` dimension. `stac-vrt` focuses on stitching together many slightly overlapping images from approximately the same time.\n",
    "\n",
    "I suspect that `stackstac` can grow to handle the \"stiching many images at a single timestep\" case. I'm not sure if `stac-vrt` can grow to handle the \"stack of images along time\" case, mainly because I don't know if GDAL's VRT can (nested VRTs, especially without local files?).\n",
    "\n",
    "2. **implementation**: build a `VRT` vs. build a `DataArray`\n",
    "\n",
    "`stac-vrt` is a bit simpler. It's solely focused on building a VRT from an `ItemCollection`. It knows nothing about I/O, Dask, etc. It inherits all the limitations of VRTs.\n",
    "\n",
    "`stackstac` is more complicated, but I think more flexible. It implements the code to do \"this chunk of this Dask array comes from reading that chunk of that source file\". I suspect that that's more flexible than `stack-vrt`."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}