tutorials/rtc-qualitative-assessment.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "db8d71aa",
   "metadata": {},
   "source": [
    "## Radiometric terrain correction qualitative assessment \n",
    "\n",
    "[sarsen](https://github.com/bopen/sarsen) radiometric correction implements the gamma flatting algorithm proposed by Small, David (2011). [Flattening Gamma: Radiometric Terrain Correction for SAR Imagery. Geoscience and Remote Sensing, IEEE Transactions on. 49. 3081 - 3093. 10.1109/TGRS.2011.2120616](https://www.researchgate.net/publication/224230688_Flattening_Gamma_Radiometric_Terrain_Correction_for_SAR_Imagery)\n",
    "\n",
    "In this notebook, we perform the same qualitative test done by the authors to evaluate the goodness of the gamma flattening correction. For an introduction to terrain correction, see [Customizable radiometric terrain corrections of Sentinel-1 products\n",
    "](https://planetarycomputer.microsoft.com/docs/tutorials/customizable-rtc-sentinel1/)."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "6d604a80-9222-46f2-97cb-afe5bfb9bdbd",
   "metadata": {},
   "source": [
    "### Install Dependencies and Imports\n",
    "Additional dependencies: `sarsen==0.9.2`"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "064a137e-1437-40d8-a20b-1123b0b70c94",
   "metadata": {},
   "outputs": [],
   "source": [
    "# !pip install sarsen==0.9.2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "9f82185a-60a4-4651-8bdf-6fb1ee82c884",
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "%config InlineBackend.figure_format = 'retina'\n",
    "\n",
    "\n",
    "import matplotlib.pyplot as plt\n",
    "from matplotlib.colors import Normalize\n",
    "from matplotlib import cm\n",
    "\n",
    "plt.rcParams[\"figure.figsize\"] = (10, 7)\n",
    "plt.rcParams[\"font.size\"] = 14"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "1391c47f-d76e-4840-8eed-d547b4784331",
   "metadata": {},
   "outputs": [],
   "source": [
    "import os\n",
    "import tempfile\n",
    "\n",
    "import numpy as np\n",
    "import xarray as xr\n",
    "from scipy.stats import gaussian_kde as kde\n",
    "from sarsen import apps, geocoding, orbit, scene\n",
    "\n",
    "import adlfs\n",
    "import planetary_computer\n",
    "import pystac_client\n",
    "import stackstac\n",
    "import rioxarray  # noqa: F401\n",
    "\n",
    "plt.rcParams[\"figure.figsize\"] = (10, 7)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "86970ee8",
   "metadata": {},
   "source": [
    "### Processing definitions"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "107608ba-51f5-4a3d-ab70-bc9fb34eafed",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'/tmp'"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# create a temporary directory where to store downloaded data\n",
    "tmp_dir = tempfile.gettempdir()\n",
    "\n",
    "# DEM path\n",
    "dem_path = os.path.join(tmp_dir, \"South-of-Redmond-10m.tif\")\n",
    "\n",
    "# path to Sentinel-1 input product in the Planetary Computer\n",
    "product_folder = \"GRD/2021/12/17/IW/DV/S1B_IW_GRDH_1SDV_20211217T141304_20211217T141329_030066_039705_9048\"  # noqa: E501\n",
    "\n",
    "# band to be processed\n",
    "measurement_group = \"IW/VV\"\n",
    "\n",
    "tmp_dir"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "4dd0a0f2-4476-48fa-bceb-6f5a872c59de",
   "metadata": {},
   "source": [
    "#### Area of interest definition: South-of-Redmond (Seattle, US)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "8ec8652a-7736-40cc-b8d9-1134493e5888",
   "metadata": {},
   "outputs": [],
   "source": [
    "lon, lat = [-121.95, 47.04]\n",
    "buffer = 0.2\n",
    "bbox = [lon - buffer, lat - buffer, lon + buffer, lat + buffer]"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b144d7a4-839b-4e41-95ca-6d20418432d9",
   "metadata": {},
   "source": [
    "#### DEMs discovery\n",
    "\n",
    "Here we use the DEM with a 10-meter ground sample distance (GDS) available on the Planetary Computer. Note that **any DEM supported by GDAL/Proj can be used**.\n",
    "\n",
    "Using `pystac_client` we can search the Planetary Computer's STAC endpoint for items matching our query parameters.  \n",
    "As multiple DEMs acquired at different times are available in this area, we select the DEMs with 10-meter GDS and perform the average of the remaining DEMs along the time dimension."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "e08b28b1-a706-4005-b1e8-633a748e7f69",
   "metadata": {},
   "outputs": [],
   "source": [
    "catalog = pystac_client.Client.open(\n",
    "    \"https://planetarycomputer.microsoft.com/api/stac/v1\",\n",
    ")\n",
    "search = catalog.search(\n",
    "    collections=\"3dep-seamless\", bbox=bbox, query={\"gsd\": {\"eq\": 10}}\n",
    ")\n",
    "items = planetary_computer.sign(search.get_all_items())"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "7503fa77-ba5a-492c-9662-ae2f9b2ba5f4",
   "metadata": {},
   "source": [
    "Here we load the data into an xarray `DataArray` using stackstac."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "625eb3d3-90b0-429b-addc-e0ac59e60942",
   "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-69a5cbaa571c4ef68cbeaba52a270af4&#x27; (time: 4,\n",
       "                                                                y: 4321, x: 4321)&gt;\n",
       "dask.array&lt;getitem, shape=(4, 4321, 4321), dtype=float64, chunksize=(1, 1024, 1024), chunktype=numpy.ndarray&gt;\n",
       "Coordinates: (12/13)\n",
       "  * time             (time) datetime64[ns] 2018-02-02 2018-02-08 ... 2020-01-07\n",
       "    id               (time) &lt;U10 &#x27;n48w122-13&#x27; &#x27;n47w122-13&#x27; ... &#x27;n48w123-13&#x27;\n",
       "    band             &lt;U4 &#x27;data&#x27;\n",
       "  * x                (x) float64 -122.2 -122.1 -122.1 ... -121.8 -121.8 -121.8\n",
       "  * y                (y) float64 47.24 47.24 47.24 47.24 ... 46.84 46.84 46.84\n",
       "    end_datetime     (time) &lt;U20 &#x27;2016-12-31T00:00:00Z&#x27; ... &#x27;2019-04-25T00:00...\n",
       "    ...               ...\n",
       "    proj:shape       object {10812}\n",
       "    start_datetime   (time) &lt;U20 &#x27;1952-01-01T00:00:00Z&#x27; ... &#x27;2017-09-14T00:00...\n",
       "    proj:epsg        int64 5498\n",
       "    gsd              int64 10\n",
       "    description      &lt;U1849 &#x27;This tile of the 3D Elevation Program (3DEP) sea...\n",
       "    epsg             int64 5498\n",
       "Attributes:\n",
       "    spec:        RasterSpec(epsg=5498, bounds=(-122.15000053746, 46.839907613...\n",
       "    crs:         epsg:5498\n",
       "    transform:   | 0.00, 0.00,-122.15|\\n| 0.00,-0.00, 47.24|\\n| 0.00, 0.00, 1...\n",
       "    resolution:  9.2592593e-05</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-69a5cbaa571c4ef68cbeaba52a270af4'</div><ul class='xr-dim-list'><li><span class='xr-has-index'>time</span>: 4</li><li><span class='xr-has-index'>y</span>: 4321</li><li><span class='xr-has-index'>x</span>: 4321</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-46726180-20c8-475e-836f-d01852a242bc' class='xr-array-in' type='checkbox' checked><label for='section-46726180-20c8-475e-836f-d01852a242bc' 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, 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> 569.79 MiB </td>\n",
       "                        <td> 8.00 MiB </td>\n",
       "                    </tr>\n",
       "                    \n",
       "                    <tr>\n",
       "                        <th> Shape </th>\n",
       "                        <td> (4, 4321, 4321) </td>\n",
       "                        <td> (1, 1024, 1024) </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <th> Count </th>\n",
       "                        <td> 208 Tasks </td>\n",
       "                        <td> 100 Chunks </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                    <th> Type </th>\n",
       "                    <td> float64 </td>\n",
       "                    <td> numpy.ndarray </td>\n",
       "                    </tr>\n",
       "                </tbody>\n",
       "            </table>\n",
       "        </td>\n",
       "        <td>\n",
       "        <svg width=\"194\" 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=\"28\" x2=\"24\" y2=\"43\" />\n",
       "  <line x1=\"10\" y1=\"56\" x2=\"24\" y2=\"71\" />\n",
       "  <line x1=\"10\" y1=\"85\" x2=\"24\" y2=\"100\" />\n",
       "  <line x1=\"10\" y1=\"113\" x2=\"24\" y2=\"128\" />\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:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"10\" y1=\"0\" x2=\"130\" y2=\"0\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"13\" y1=\"3\" x2=\"133\" y2=\"3\" />\n",
       "  <line x1=\"17\" y1=\"7\" x2=\"137\" y2=\"7\" />\n",
       "  <line x1=\"21\" y1=\"11\" x2=\"141\" y2=\"11\" />\n",
       "  <line x1=\"24\" y1=\"14\" x2=\"144\" 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=\"38\" y1=\"0\" x2=\"53\" y2=\"14\" />\n",
       "  <line x1=\"66\" y1=\"0\" x2=\"81\" y2=\"14\" />\n",
       "  <line x1=\"95\" y1=\"0\" x2=\"110\" y2=\"14\" />\n",
       "  <line x1=\"123\" y1=\"0\" x2=\"138\" y2=\"14\" />\n",
       "  <line x1=\"130\" y1=\"0\" x2=\"144\" y2=\"14\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"10.0,0.0 130.0,0.0 144.9485979497544,14.948597949754403 24.9485979497544,14.948597949754403\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"24\" y1=\"14\" x2=\"144\" y2=\"14\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"24\" y1=\"43\" x2=\"144\" y2=\"43\" />\n",
       "  <line x1=\"24\" y1=\"71\" x2=\"144\" y2=\"71\" />\n",
       "  <line x1=\"24\" y1=\"100\" x2=\"144\" y2=\"100\" />\n",
       "  <line x1=\"24\" y1=\"128\" x2=\"144\" y2=\"128\" />\n",
       "  <line x1=\"24\" y1=\"134\" x2=\"144\" 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=\"53\" y1=\"14\" x2=\"53\" y2=\"134\" />\n",
       "  <line x1=\"81\" y1=\"14\" x2=\"81\" y2=\"134\" />\n",
       "  <line x1=\"110\" y1=\"14\" x2=\"110\" y2=\"134\" />\n",
       "  <line x1=\"138\" y1=\"14\" x2=\"138\" y2=\"134\" />\n",
       "  <line x1=\"144\" y1=\"14\" x2=\"144\" y2=\"134\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"24.9485979497544,14.948597949754403 144.9485979497544,14.948597949754403 144.9485979497544,134.9485979497544 24.9485979497544,134.9485979497544\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Text -->\n",
       "  <text x=\"84.948598\" y=\"154.948598\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" >4321</text>\n",
       "  <text x=\"164.948598\" y=\"74.948598\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(-90,164.948598,74.948598)\">4321</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-3522941c-06cc-4cad-bc88-0df9b60617bf' class='xr-section-summary-in' type='checkbox'  checked><label for='section-3522941c-06cc-4cad-bc88-0df9b60617bf' class='xr-section-summary' >Coordinates: <span>(13)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>time</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>datetime64[ns]</div><div class='xr-var-preview xr-preview'>2018-02-02 ... 2020-01-07</div><input id='attrs-5dea2b79-10e8-4135-9cd7-0348c7c6dedf' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-5dea2b79-10e8-4135-9cd7-0348c7c6dedf' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-a7545b5a-62a2-440a-8b52-4373e819849f' class='xr-var-data-in' type='checkbox'><label for='data-a7545b5a-62a2-440a-8b52-4373e819849f' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([&#x27;2018-02-02T00:00:00.000000000&#x27;, &#x27;2018-02-08T00:00:00.000000000&#x27;,\n",
       "       &#x27;2018-02-08T00:00:00.000000000&#x27;, &#x27;2020-01-07T00:00:00.000000000&#x27;],\n",
       "      dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>id</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>&lt;U10</div><div class='xr-var-preview xr-preview'>&#x27;n48w122-13&#x27; ... &#x27;n48w123-13&#x27;</div><input id='attrs-c5f3cfb2-23c6-45ef-b378-72f53d6ebe01' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-c5f3cfb2-23c6-45ef-b378-72f53d6ebe01' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-aaf07900-ea7b-4d56-bde7-9b9c9e481a13' class='xr-var-data-in' type='checkbox'><label for='data-aaf07900-ea7b-4d56-bde7-9b9c9e481a13' 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;n48w122-13&#x27;, &#x27;n47w122-13&#x27;, &#x27;n47w123-13&#x27;, &#x27;n48w123-13&#x27;],\n",
       "      dtype=&#x27;&lt;U10&#x27;)</pre></div></li><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-c78d0e06-f823-4882-a805-885971a5e333' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-c78d0e06-f823-4882-a805-885971a5e333' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-d001e7bf-f41e-460f-bfc3-4971f6eee70b' class='xr-var-data-in' type='checkbox'><label for='data-d001e7bf-f41e-460f-bfc3-4971f6eee70b' 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'>-122.2 -122.1 ... -121.8 -121.8</div><input id='attrs-512eca3c-6d39-4553-af4c-46087466ea4b' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-512eca3c-6d39-4553-af4c-46087466ea4b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-52ad432c-6154-49c6-9fda-e33ce43d8ce9' class='xr-var-data-in' type='checkbox'><label for='data-52ad432c-6154-49c6-9fda-e33ce43d8ce9' 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([-122.150001, -122.149908, -122.149815, ..., -121.750186, -121.750093,\n",
       "       -121.750001])</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'>47.24 47.24 47.24 ... 46.84 46.84</div><input id='attrs-43c1872b-3e68-482c-857a-eb63237b605e' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-43c1872b-3e68-482c-857a-eb63237b605e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3ba34fb2-7bab-4132-a8f1-9885975e33b3' class='xr-var-data-in' type='checkbox'><label for='data-3ba34fb2-7bab-4132-a8f1-9885975e33b3' 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([47.24    , 47.239908, 47.239815, ..., 46.840185, 46.840093, 46.84    ])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>end_datetime</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>&lt;U20</div><div class='xr-var-preview xr-preview'>&#x27;2016-12-31T00:00:00Z&#x27; ... &#x27;2019...</div><input id='attrs-e0909be2-6578-460f-8cee-90c32fb6781a' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-e0909be2-6578-460f-8cee-90c32fb6781a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5a246dc5-af43-4540-a4b1-11497acce4f6' class='xr-var-data-in' type='checkbox'><label for='data-5a246dc5-af43-4540-a4b1-11497acce4f6' 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;2016-12-31T00:00:00Z&#x27;, &#x27;2016-12-31T00:00:00Z&#x27;,\n",
       "       &#x27;2016-12-31T00:00:00Z&#x27;, &#x27;2019-04-25T00:00:00Z&#x27;], dtype=&#x27;&lt;U20&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>threedep:region</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U7</div><div class='xr-var-preview xr-preview'>&#x27;n40w130&#x27;</div><input id='attrs-ac6daee4-ee1c-4a6a-9a54-6a1126b8b3a8' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-ac6daee4-ee1c-4a6a-9a54-6a1126b8b3a8' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b1fd0a94-dcfc-45a9-b386-f600714ffba4' class='xr-var-data-in' type='checkbox'><label for='data-b1fd0a94-dcfc-45a9-b386-f600714ffba4' 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;n40w130&#x27;, dtype=&#x27;&lt;U7&#x27;)</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'>{10812}</div><input id='attrs-c882c400-0684-44de-96d9-e88b0a319bf5' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-c882c400-0684-44de-96d9-e88b0a319bf5' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3f493795-2744-4a0f-81ff-0edc3a81840c' class='xr-var-data-in' type='checkbox'><label for='data-3f493795-2744-4a0f-81ff-0edc3a81840c' 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({10812}, 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'>(time)</div><div class='xr-var-dtype'>&lt;U20</div><div class='xr-var-preview xr-preview'>&#x27;1952-01-01T00:00:00Z&#x27; ... &#x27;2017...</div><input id='attrs-5c965fb0-531b-4dc5-9026-58160f06ced3' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-5c965fb0-531b-4dc5-9026-58160f06ced3' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-ed6c2d32-6903-4a69-b46c-ea6327921665' class='xr-var-data-in' type='checkbox'><label for='data-ed6c2d32-6903-4a69-b46c-ea6327921665' 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;1952-01-01T00:00:00Z&#x27;, &#x27;1956-01-01T00:00:00Z&#x27;,\n",
       "       &#x27;1947-01-01T00:00:00Z&#x27;, &#x27;2017-09-14T00:00:00Z&#x27;], dtype=&#x27;&lt;U20&#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'>5498</div><input id='attrs-9864e4ee-1f95-4afe-ac03-f0c5aee0dc30' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-9864e4ee-1f95-4afe-ac03-f0c5aee0dc30' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-667cd256-1240-4e2d-bba4-3a995774fddc' class='xr-var-data-in' type='checkbox'><label for='data-667cd256-1240-4e2d-bba4-3a995774fddc' 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(5498)</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-e0ee4727-f65d-48b3-9fd0-9db5297098be' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-e0ee4727-f65d-48b3-9fd0-9db5297098be' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-69fa6f8e-1406-4805-bc78-6f12987e0074' class='xr-var-data-in' type='checkbox'><label for='data-69fa6f8e-1406-4805-bc78-6f12987e0074' 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>description</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U1849</div><div class='xr-var-preview xr-preview'>&#x27;This tile of the 3D Elevation P...</div><input id='attrs-89804ac6-4978-4d80-a6e7-f4c81d45e167' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-89804ac6-4978-4d80-a6e7-f4c81d45e167' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b3a99643-599a-4257-a779-9b15a813a953' class='xr-var-data-in' type='checkbox'><label for='data-b3a99643-599a-4257-a779-9b15a813a953' 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;This tile of the 3D Elevation Program (3DEP) seamless products is 1/3 arc-second resolution.  3DEP data serve as the elevation layer of The National Map, and provide basic elevation information for Earth science studies and mapping applications in the United States. Scientists and resource managers use 3DEP data for global change research, hydrologic modeling, resource monitoring, mapping and visualization, and many other applications. 3DEP data compose an elevation dataset that consists of seamless layers and a high resolution layer.  Each of these layers consists of the best available raster elevation data of the conterminous United States, Alaska, Hawaii, territorial islands, Mexico and Canada.  3DEP data are updated continually as new data become available. Seamless 3DEP data are derived from diverse source data that are processed to a common coordinate system and unit of vertical measure. These data are distributed in geographic coordinates in units of decimal degrees, and in conformance with the North American Datum of 1983 (NAD 83). All elevation values are in meters and, over the conterminous United States, are referenced to the North American Vertical Datum of 1988 (NAVD 88). The vertical reference will vary in other areas. Seamless 3DEP data are available nationally (except for Alaska) at resolutions of 1 arc-second (approximately 30 meters) and 1/3 arc-second (approximately 10 meters). In most of Alaska, only lower resolution source data are available. As a result, most seamless 3DEP data for Alaska are at 2 arc-second (approximately 60 meters) grid spacing. Part of Alaska is available at the 1 and 1/3 arc-second resolutions from interferometric synthetic aperture radar (ifsar) collections starting in 2010.  Plans are in place for collection of statewide ifsar in Alaska. All 3DEP products are public domain.&#x27;,\n",
       "      dtype=&#x27;&lt;U1849&#x27;)</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'>5498</div><input id='attrs-2ed79729-7e5c-4e88-bfcd-05b033879932' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-2ed79729-7e5c-4e88-bfcd-05b033879932' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4886e9e0-6c5c-4ed4-8be3-bccc1a0b3e49' class='xr-var-data-in' type='checkbox'><label for='data-4886e9e0-6c5c-4ed4-8be3-bccc1a0b3e49' 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(5498)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-ad1ccfa7-df7a-4270-a3ef-5adc92750cca' class='xr-section-summary-in' type='checkbox'  checked><label for='section-ad1ccfa7-df7a-4270-a3ef-5adc92750cca' 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=5498, bounds=(-122.15000053746, 46.839907613503, -121.749907943107, 47.240000207856), resolutions_xy=(9.2592593e-05, 9.2592593e-05))</dd><dt><span>crs :</span></dt><dd>epsg:5498</dd><dt><span>transform :</span></dt><dd>| 0.00, 0.00,-122.15|\n",
       "| 0.00,-0.00, 47.24|\n",
       "| 0.00, 0.00, 1.00|</dd><dt><span>resolution :</span></dt><dd>9.2592593e-05</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.DataArray 'stackstac-69a5cbaa571c4ef68cbeaba52a270af4' (time: 4,\n",
       "                                                                y: 4321, x: 4321)>\n",
       "dask.array<getitem, shape=(4, 4321, 4321), dtype=float64, chunksize=(1, 1024, 1024), chunktype=numpy.ndarray>\n",
       "Coordinates: (12/13)\n",
       "  * time             (time) datetime64[ns] 2018-02-02 2018-02-08 ... 2020-01-07\n",
       "    id               (time) <U10 'n48w122-13' 'n47w122-13' ... 'n48w123-13'\n",
       "    band             <U4 'data'\n",
       "  * x                (x) float64 -122.2 -122.1 -122.1 ... -121.8 -121.8 -121.8\n",
       "  * y                (y) float64 47.24 47.24 47.24 47.24 ... 46.84 46.84 46.84\n",
       "    end_datetime     (time) <U20 '2016-12-31T00:00:00Z' ... '2019-04-25T00:00...\n",
       "    ...               ...\n",
       "    proj:shape       object {10812}\n",
       "    start_datetime   (time) <U20 '1952-01-01T00:00:00Z' ... '2017-09-14T00:00...\n",
       "    proj:epsg        int64 5498\n",
       "    gsd              int64 10\n",
       "    description      <U1849 'This tile of the 3D Elevation Program (3DEP) sea...\n",
       "    epsg             int64 5498\n",
       "Attributes:\n",
       "    spec:        RasterSpec(epsg=5498, bounds=(-122.15000053746, 46.839907613...\n",
       "    crs:         epsg:5498\n",
       "    transform:   | 0.00, 0.00,-122.15|\\n| 0.00,-0.00, 47.24|\\n| 0.00, 0.00, 1...\n",
       "    resolution:  9.2592593e-05"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# select DEMs with resolution 10 meters\n",
    "dem_raster_all = stackstac.stack(items, bounds=bbox).squeeze()\n",
    "dem_raster_all"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "19adb6a3-26a4-4fdb-862a-1f6aa4f2abf6",
   "metadata": {},
   "source": [
    "#### DEMs average along the time dimension"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "65ae53fb",
   "metadata": {},
   "outputs": [],
   "source": [
    "dem_raster_geo = dem_raster_all.compute().mean(\"time\")\n",
    "dem_raster_geo.rio.set_crs(dem_raster_all.rio.crs);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c508e760",
   "metadata": {},
   "source": [
    "#### Convert the DEM in UTM coordinates\n",
    "\n",
    "In order to facilitate the comparison between the RTC computed by sarsen with the RTC available on the Planetery Computer, here we convert the DEM in UTM.\n",
    "\n",
    "#### RESAMPLING UTM DEM"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "f6ad2b6e",
   "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-69a5cbaa571c4ef68cbeaba52a270af4&#x27; (y: 4487, x: 3100)&gt;\n",
       "array([[1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308],\n",
       "       [1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308],\n",
       "       [1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308],\n",
       "       ...,\n",
       "       [1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308],\n",
       "       [1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308],\n",
       "       [1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308]])\n",
       "Coordinates:\n",
       "  * x                (x) float64 5.643e+05 5.643e+05 ... 5.953e+05 5.953e+05\n",
       "  * y                (y) float64 5.233e+06 5.233e+06 ... 5.188e+06 5.188e+06\n",
       "    threedep:region  &lt;U7 &#x27;n40w130&#x27;\n",
       "    description      &lt;U1849 &#x27;This tile of the 3D Elevation Program (3DEP) sea...\n",
       "    epsg             int64 5498\n",
       "    proj:shape       object {10812}\n",
       "    proj:epsg        int64 5498\n",
       "    band             &lt;U4 &#x27;data&#x27;\n",
       "    gsd              int64 10\n",
       "    spatial_ref      int64 0\n",
       "Attributes:\n",
       "    _FillValue:  1.7976931348623157e+308</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-69a5cbaa571c4ef68cbeaba52a270af4'</div><ul class='xr-dim-list'><li><span class='xr-has-index'>y</span>: 4487</li><li><span class='xr-has-index'>x</span>: 3100</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-141edfb6-704e-46f8-b0db-9b4d62a3b17a' class='xr-array-in' type='checkbox' checked><label for='section-141edfb6-704e-46f8-b0db-9b4d62a3b17a' 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>1.798e+308 1.798e+308 1.798e+308 ... 1.798e+308 1.798e+308 1.798e+308</span></div><div class='xr-array-data'><pre>array([[1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308],\n",
       "       [1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308],\n",
       "       [1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308],\n",
       "       ...,\n",
       "       [1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308],\n",
       "       [1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308],\n",
       "       [1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308]])</pre></div></div></li><li class='xr-section-item'><input id='section-3dbbc33e-4edc-46b4-834a-9f95746b0e5d' class='xr-section-summary-in' type='checkbox'  checked><label for='section-3dbbc33e-4edc-46b4-834a-9f95746b0e5d' class='xr-section-summary' >Coordinates: <span>(10)</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'>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.643e+05 5.643e+05 ... 5.953e+05</div><input id='attrs-64530a57-f5a1-402f-afc1-18718f862753' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-64530a57-f5a1-402f-afc1-18718f862753' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-728f3fb1-efcf-4952-b1d1-c549d0829873' class='xr-var-data-in' type='checkbox'><label for='data-728f3fb1-efcf-4952-b1d1-c549d0829873' 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>axis :</span></dt><dd>X</dd><dt><span>long_name :</span></dt><dd>x coordinate of projection</dd><dt><span>standard_name :</span></dt><dd>projection_x_coordinate</dd><dt><span>units :</span></dt><dd>metre</dd></dl></div><div class='xr-var-data'><pre>array([564332.162431, 564342.162431, 564352.162431, ..., 595302.162431,\n",
       "       595312.162431, 595322.162431])</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'>5.233e+06 5.233e+06 ... 5.188e+06</div><input id='attrs-407e479e-f390-4b85-8156-b2f0530ea81f' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-407e479e-f390-4b85-8156-b2f0530ea81f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-34e11f9e-68ff-4fcf-945c-d5622628b18f' class='xr-var-data-in' type='checkbox'><label for='data-34e11f9e-68ff-4fcf-945c-d5622628b18f' 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>axis :</span></dt><dd>Y</dd><dt><span>long_name :</span></dt><dd>y coordinate of projection</dd><dt><span>standard_name :</span></dt><dd>projection_y_coordinate</dd><dt><span>units :</span></dt><dd>metre</dd></dl></div><div class='xr-var-data'><pre>array([5232593.548183, 5232583.548183, 5232573.548183, ..., 5187753.548183,\n",
       "       5187743.548183, 5187733.548183])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>threedep:region</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U7</div><div class='xr-var-preview xr-preview'>&#x27;n40w130&#x27;</div><input id='attrs-24620492-a65a-463d-b78f-1a5dfb43c5c9' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-24620492-a65a-463d-b78f-1a5dfb43c5c9' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5421bc2b-54a9-4e36-9b7f-b40e6aefd3c0' class='xr-var-data-in' type='checkbox'><label for='data-5421bc2b-54a9-4e36-9b7f-b40e6aefd3c0' 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;n40w130&#x27;, dtype=&#x27;&lt;U7&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>description</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U1849</div><div class='xr-var-preview xr-preview'>&#x27;This tile of the 3D Elevation P...</div><input id='attrs-c1a171cd-0f6c-4968-9946-4c9ee891825d' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-c1a171cd-0f6c-4968-9946-4c9ee891825d' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3441b2c8-4e47-49e6-8baf-ace185aac302' class='xr-var-data-in' type='checkbox'><label for='data-3441b2c8-4e47-49e6-8baf-ace185aac302' 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;This tile of the 3D Elevation Program (3DEP) seamless products is 1/3 arc-second resolution.  3DEP data serve as the elevation layer of The National Map, and provide basic elevation information for Earth science studies and mapping applications in the United States. Scientists and resource managers use 3DEP data for global change research, hydrologic modeling, resource monitoring, mapping and visualization, and many other applications. 3DEP data compose an elevation dataset that consists of seamless layers and a high resolution layer.  Each of these layers consists of the best available raster elevation data of the conterminous United States, Alaska, Hawaii, territorial islands, Mexico and Canada.  3DEP data are updated continually as new data become available. Seamless 3DEP data are derived from diverse source data that are processed to a common coordinate system and unit of vertical measure. These data are distributed in geographic coordinates in units of decimal degrees, and in conformance with the North American Datum of 1983 (NAD 83). All elevation values are in meters and, over the conterminous United States, are referenced to the North American Vertical Datum of 1988 (NAVD 88). The vertical reference will vary in other areas. Seamless 3DEP data are available nationally (except for Alaska) at resolutions of 1 arc-second (approximately 30 meters) and 1/3 arc-second (approximately 10 meters). In most of Alaska, only lower resolution source data are available. As a result, most seamless 3DEP data for Alaska are at 2 arc-second (approximately 60 meters) grid spacing. Part of Alaska is available at the 1 and 1/3 arc-second resolutions from interferometric synthetic aperture radar (ifsar) collections starting in 2010.  Plans are in place for collection of statewide ifsar in Alaska. All 3DEP products are public domain.&#x27;,\n",
       "      dtype=&#x27;&lt;U1849&#x27;)</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'>5498</div><input id='attrs-455c01da-8a4c-4219-a4ef-61981e84204b' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-455c01da-8a4c-4219-a4ef-61981e84204b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-202cc3f1-f532-4f80-a63e-9a8efc894340' class='xr-var-data-in' type='checkbox'><label for='data-202cc3f1-f532-4f80-a63e-9a8efc894340' 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(5498)</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'>{10812}</div><input id='attrs-4c58fd64-a480-49f6-b8fd-ec848a5baace' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-4c58fd64-a480-49f6-b8fd-ec848a5baace' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-86c99d14-8a8a-46ca-919f-e836cbee93b6' class='xr-var-data-in' type='checkbox'><label for='data-86c99d14-8a8a-46ca-919f-e836cbee93b6' 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({10812}, dtype=object)</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'>5498</div><input id='attrs-dd145cc4-93b8-4d85-9853-7404475b84b0' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-dd145cc4-93b8-4d85-9853-7404475b84b0' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-1a802eb5-0cf6-424a-bb72-554b85e1ba13' class='xr-var-data-in' type='checkbox'><label for='data-1a802eb5-0cf6-424a-bb72-554b85e1ba13' 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(5498)</pre></div></li><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-cc623add-1bd2-4d24-9998-85793f4e1d12' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-cc623add-1bd2-4d24-9998-85793f4e1d12' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-d9649b8f-36a0-4253-bc48-c79512d9f1d0' class='xr-var-data-in' type='checkbox'><label for='data-d9649b8f-36a0-4253-bc48-c79512d9f1d0' 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>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-8844a24c-1568-4640-a2d7-1dfdf75dbd27' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-8844a24c-1568-4640-a2d7-1dfdf75dbd27' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-de92e1b8-ee67-479b-9b30-22fd16e45690' class='xr-var-data-in' type='checkbox'><label for='data-de92e1b8-ee67-479b-9b30-22fd16e45690' 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>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-a917a911-a94e-4552-a3de-423c9f0a05aa' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-a917a911-a94e-4552-a3de-423c9f0a05aa' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-aaac922b-b4b2-49be-92d3-60609a47973e' class='xr-var-data-in' type='checkbox'><label for='data-aaac922b-b4b2-49be-92d3-60609a47973e' 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;WGS 84 / UTM zone 10N&quot;,GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],PRIMEM[&quot;Greenwich&quot;,0,AUTHORITY[&quot;EPSG&quot;,&quot;8901&quot;]],UNIT[&quot;degree&quot;,0.0174532925199433,AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]],PROJECTION[&quot;Transverse_Mercator&quot;],PARAMETER[&quot;latitude_of_origin&quot;,0],PARAMETER[&quot;central_meridian&quot;,-123],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;32610&quot;]]</dd><dt><span>semi_major_axis :</span></dt><dd>6378137.0</dd><dt><span>semi_minor_axis :</span></dt><dd>6356752.314245179</dd><dt><span>inverse_flattening :</span></dt><dd>298.257223563</dd><dt><span>reference_ellipsoid_name :</span></dt><dd>WGS 84</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>WGS 84</dd><dt><span>horizontal_datum_name :</span></dt><dd>World Geodetic System 1984</dd><dt><span>projected_crs_name :</span></dt><dd>WGS 84 / UTM zone 10N</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>-123.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;WGS 84 / UTM zone 10N&quot;,GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],PRIMEM[&quot;Greenwich&quot;,0,AUTHORITY[&quot;EPSG&quot;,&quot;8901&quot;]],UNIT[&quot;degree&quot;,0.0174532925199433,AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]],PROJECTION[&quot;Transverse_Mercator&quot;],PARAMETER[&quot;latitude_of_origin&quot;,0],PARAMETER[&quot;central_meridian&quot;,-123],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;32610&quot;]]</dd><dt><span>GeoTransform :</span></dt><dd>564327.1624312083 10.0 0.0 5232598.548182782 0.0 -10.0</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-afec104e-4bc7-403f-8b53-af6b33403fe0' class='xr-section-summary-in' type='checkbox'  checked><label for='section-afec104e-4bc7-403f-8b53-af6b33403fe0' class='xr-section-summary' >Attributes: <span>(1)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>_FillValue :</span></dt><dd>1.7976931348623157e+308</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.DataArray 'stackstac-69a5cbaa571c4ef68cbeaba52a270af4' (y: 4487, x: 3100)>\n",
       "array([[1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308],\n",
       "       [1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308],\n",
       "       [1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308],\n",
       "       ...,\n",
       "       [1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308],\n",
       "       [1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308],\n",
       "       [1.79769313e+308, 1.79769313e+308, 1.79769313e+308, ...,\n",
       "        1.79769313e+308, 1.79769313e+308, 1.79769313e+308]])\n",
       "Coordinates:\n",
       "  * x                (x) float64 5.643e+05 5.643e+05 ... 5.953e+05 5.953e+05\n",
       "  * y                (y) float64 5.233e+06 5.233e+06 ... 5.188e+06 5.188e+06\n",
       "    threedep:region  <U7 'n40w130'\n",
       "    description      <U1849 'This tile of the 3D Elevation Program (3DEP) sea...\n",
       "    epsg             int64 5498\n",
       "    proj:shape       object {10812}\n",
       "    proj:epsg        int64 5498\n",
       "    band             <U4 'data'\n",
       "    gsd              int64 10\n",
       "    spatial_ref      int64 0\n",
       "Attributes:\n",
       "    _FillValue:  1.7976931348623157e+308"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# find the UTM zone and project in UTM\n",
    "t_srs = dem_raster_geo.rio.estimate_utm_crs()\n",
    "dem_raster = dem_raster_geo.rio.reproject(t_srs, resolution=(10, 10))\n",
    "dem_raster"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "5c4a8e3b",
   "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-69a5cbaa571c4ef68cbeaba52a270af4&#x27; (y: 1500, x: 1400)&gt;\n",
       "array([[ 840.86486816,  837.69940186,  836.6675415 , ..., 1222.60290527,\n",
       "        1224.68151855, 1229.60192871],\n",
       "       [ 847.59912109,  843.43511963,  841.19110107, ..., 1226.77722168,\n",
       "        1229.88330078, 1235.05480957],\n",
       "       [ 853.09313965,  849.24804688,  847.56011963, ..., 1231.57214355,\n",
       "        1235.27197266, 1240.26989746],\n",
       "       ...,\n",
       "       [1057.71154785, 1056.44726562, 1053.34716797, ..., 1967.21203613,\n",
       "        1966.72473145, 1967.17895508],\n",
       "       [1053.78833008, 1052.40124512, 1049.59545898, ..., 1970.66943359,\n",
       "        1969.1439209 , 1969.60314941],\n",
       "       [1049.06298828, 1048.03967285, 1045.55126953, ..., 1973.11730957,\n",
       "        1971.60437012, 1971.93334961]])\n",
       "Coordinates:\n",
       "  * x                (x) float64 5.79e+05 5.79e+05 ... 5.93e+05 5.93e+05\n",
       "  * y                (y) float64 5.21e+06 5.21e+06 ... 5.195e+06 5.195e+06\n",
       "    threedep:region  &lt;U7 &#x27;n40w130&#x27;\n",
       "    description      &lt;U1849 &#x27;This tile of the 3D Elevation Program (3DEP) sea...\n",
       "    epsg             int64 5498\n",
       "    proj:shape       object {10812}\n",
       "    proj:epsg        int64 5498\n",
       "    band             &lt;U4 &#x27;data&#x27;\n",
       "    gsd              int64 10\n",
       "    spatial_ref      int64 0\n",
       "Attributes:\n",
       "    _FillValue:  1.7976931348623157e+308</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-69a5cbaa571c4ef68cbeaba52a270af4'</div><ul class='xr-dim-list'><li><span class='xr-has-index'>y</span>: 1500</li><li><span class='xr-has-index'>x</span>: 1400</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-9ff28201-135d-49eb-9594-9efdc3d41398' class='xr-array-in' type='checkbox' checked><label for='section-9ff28201-135d-49eb-9594-9efdc3d41398' 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>840.9 837.7 836.7 835.7 ... 1.988e+03 1.973e+03 1.972e+03 1.972e+03</span></div><div class='xr-array-data'><pre>array([[ 840.86486816,  837.69940186,  836.6675415 , ..., 1222.60290527,\n",
       "        1224.68151855, 1229.60192871],\n",
       "       [ 847.59912109,  843.43511963,  841.19110107, ..., 1226.77722168,\n",
       "        1229.88330078, 1235.05480957],\n",
       "       [ 853.09313965,  849.24804688,  847.56011963, ..., 1231.57214355,\n",
       "        1235.27197266, 1240.26989746],\n",
       "       ...,\n",
       "       [1057.71154785, 1056.44726562, 1053.34716797, ..., 1967.21203613,\n",
       "        1966.72473145, 1967.17895508],\n",
       "       [1053.78833008, 1052.40124512, 1049.59545898, ..., 1970.66943359,\n",
       "        1969.1439209 , 1969.60314941],\n",
       "       [1049.06298828, 1048.03967285, 1045.55126953, ..., 1973.11730957,\n",
       "        1971.60437012, 1971.93334961]])</pre></div></div></li><li class='xr-section-item'><input id='section-af650600-58dc-49bc-b767-00c97976ad96' class='xr-section-summary-in' type='checkbox'  checked><label for='section-af650600-58dc-49bc-b767-00c97976ad96' class='xr-section-summary' >Coordinates: <span>(10)</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'>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.79e+05 5.79e+05 ... 5.93e+05</div><input id='attrs-981dd883-38d5-4c66-8b06-f0a7e097f62e' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-981dd883-38d5-4c66-8b06-f0a7e097f62e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e7a59a12-2583-4d9d-9be0-6c0272a7c2b9' class='xr-var-data-in' type='checkbox'><label for='data-e7a59a12-2583-4d9d-9be0-6c0272a7c2b9' 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>axis :</span></dt><dd>X</dd><dt><span>long_name :</span></dt><dd>x coordinate of projection</dd><dt><span>standard_name :</span></dt><dd>projection_x_coordinate</dd><dt><span>units :</span></dt><dd>metre</dd></dl></div><div class='xr-var-data'><pre>array([579002.162431, 579012.162431, 579022.162431, ..., 592972.162431,\n",
       "       592982.162431, 592992.162431])</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'>5.21e+06 5.21e+06 ... 5.195e+06</div><input id='attrs-7355acaf-77b2-48cf-ad27-eef3047dff56' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-7355acaf-77b2-48cf-ad27-eef3047dff56' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3f2f7477-6bde-407c-bf19-e5227afec9fe' class='xr-var-data-in' type='checkbox'><label for='data-3f2f7477-6bde-407c-bf19-e5227afec9fe' 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>axis :</span></dt><dd>Y</dd><dt><span>long_name :</span></dt><dd>y coordinate of projection</dd><dt><span>standard_name :</span></dt><dd>projection_y_coordinate</dd><dt><span>units :</span></dt><dd>metre</dd></dl></div><div class='xr-var-data'><pre>array([5209993.548183, 5209983.548183, 5209973.548183, ..., 5195023.548183,\n",
       "       5195013.548183, 5195003.548183])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>threedep:region</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U7</div><div class='xr-var-preview xr-preview'>&#x27;n40w130&#x27;</div><input id='attrs-c378d615-a942-445b-84f9-8445d4fde231' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-c378d615-a942-445b-84f9-8445d4fde231' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4145a844-3b13-4138-9dbf-d4241910a71e' class='xr-var-data-in' type='checkbox'><label for='data-4145a844-3b13-4138-9dbf-d4241910a71e' 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;n40w130&#x27;, dtype=&#x27;&lt;U7&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>description</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U1849</div><div class='xr-var-preview xr-preview'>&#x27;This tile of the 3D Elevation P...</div><input id='attrs-bc090140-f928-450e-a418-b7382e99c5ca' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-bc090140-f928-450e-a418-b7382e99c5ca' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-170ca1ba-676d-4785-815c-9b184bcb5ee8' class='xr-var-data-in' type='checkbox'><label for='data-170ca1ba-676d-4785-815c-9b184bcb5ee8' 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;This tile of the 3D Elevation Program (3DEP) seamless products is 1/3 arc-second resolution.  3DEP data serve as the elevation layer of The National Map, and provide basic elevation information for Earth science studies and mapping applications in the United States. Scientists and resource managers use 3DEP data for global change research, hydrologic modeling, resource monitoring, mapping and visualization, and many other applications. 3DEP data compose an elevation dataset that consists of seamless layers and a high resolution layer.  Each of these layers consists of the best available raster elevation data of the conterminous United States, Alaska, Hawaii, territorial islands, Mexico and Canada.  3DEP data are updated continually as new data become available. Seamless 3DEP data are derived from diverse source data that are processed to a common coordinate system and unit of vertical measure. These data are distributed in geographic coordinates in units of decimal degrees, and in conformance with the North American Datum of 1983 (NAD 83). All elevation values are in meters and, over the conterminous United States, are referenced to the North American Vertical Datum of 1988 (NAVD 88). The vertical reference will vary in other areas. Seamless 3DEP data are available nationally (except for Alaska) at resolutions of 1 arc-second (approximately 30 meters) and 1/3 arc-second (approximately 10 meters). In most of Alaska, only lower resolution source data are available. As a result, most seamless 3DEP data for Alaska are at 2 arc-second (approximately 60 meters) grid spacing. Part of Alaska is available at the 1 and 1/3 arc-second resolutions from interferometric synthetic aperture radar (ifsar) collections starting in 2010.  Plans are in place for collection of statewide ifsar in Alaska. All 3DEP products are public domain.&#x27;,\n",
       "      dtype=&#x27;&lt;U1849&#x27;)</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'>5498</div><input id='attrs-d1b9db21-5ee3-404e-8867-bae0e2f85f5a' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-d1b9db21-5ee3-404e-8867-bae0e2f85f5a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3223b92a-db5f-4456-87fc-be15791db8bf' class='xr-var-data-in' type='checkbox'><label for='data-3223b92a-db5f-4456-87fc-be15791db8bf' 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(5498)</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'>{10812}</div><input id='attrs-1cee795e-e3d0-441b-95ad-54172b85047a' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-1cee795e-e3d0-441b-95ad-54172b85047a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3942fb01-b766-46a1-ba42-28c4e96e8a33' class='xr-var-data-in' type='checkbox'><label for='data-3942fb01-b766-46a1-ba42-28c4e96e8a33' 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({10812}, dtype=object)</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'>5498</div><input id='attrs-a141b672-09a0-4e40-ada9-e59cc5799235' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-a141b672-09a0-4e40-ada9-e59cc5799235' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4e9fb56e-4aeb-4197-9c03-f796024bb9ba' class='xr-var-data-in' type='checkbox'><label for='data-4e9fb56e-4aeb-4197-9c03-f796024bb9ba' 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(5498)</pre></div></li><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-dad64b13-a4db-4ca5-8eb5-b4519e46ab2b' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-dad64b13-a4db-4ca5-8eb5-b4519e46ab2b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c3b2db17-6b76-4d25-8760-45f71691d7c4' class='xr-var-data-in' type='checkbox'><label for='data-c3b2db17-6b76-4d25-8760-45f71691d7c4' 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>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-c1aa331e-9f7c-4158-847f-80315a2bbb30' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-c1aa331e-9f7c-4158-847f-80315a2bbb30' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3b40d45e-2afe-4f8a-a4b9-f7287eef9bfa' class='xr-var-data-in' type='checkbox'><label for='data-3b40d45e-2afe-4f8a-a4b9-f7287eef9bfa' 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>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-b08c4134-8ced-446c-865d-a6513bb39600' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-b08c4134-8ced-446c-865d-a6513bb39600' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-30d207a0-8387-4655-9227-bfd783b50ca2' class='xr-var-data-in' type='checkbox'><label for='data-30d207a0-8387-4655-9227-bfd783b50ca2' 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;WGS 84 / UTM zone 10N&quot;,GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],PRIMEM[&quot;Greenwich&quot;,0,AUTHORITY[&quot;EPSG&quot;,&quot;8901&quot;]],UNIT[&quot;degree&quot;,0.0174532925199433,AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]],PROJECTION[&quot;Transverse_Mercator&quot;],PARAMETER[&quot;latitude_of_origin&quot;,0],PARAMETER[&quot;central_meridian&quot;,-123],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;32610&quot;]]</dd><dt><span>semi_major_axis :</span></dt><dd>6378137.0</dd><dt><span>semi_minor_axis :</span></dt><dd>6356752.314245179</dd><dt><span>inverse_flattening :</span></dt><dd>298.257223563</dd><dt><span>reference_ellipsoid_name :</span></dt><dd>WGS 84</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>WGS 84</dd><dt><span>horizontal_datum_name :</span></dt><dd>World Geodetic System 1984</dd><dt><span>projected_crs_name :</span></dt><dd>WGS 84 / UTM zone 10N</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>-123.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;WGS 84 / UTM zone 10N&quot;,GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],PRIMEM[&quot;Greenwich&quot;,0,AUTHORITY[&quot;EPSG&quot;,&quot;8901&quot;]],UNIT[&quot;degree&quot;,0.0174532925199433,AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]],PROJECTION[&quot;Transverse_Mercator&quot;],PARAMETER[&quot;latitude_of_origin&quot;,0],PARAMETER[&quot;central_meridian&quot;,-123],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;32610&quot;]]</dd><dt><span>GeoTransform :</span></dt><dd>564327.1624312083 10.0 0.0 5232598.548182782 0.0 -10.0</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-85044860-f745-45dd-946d-8b817708a16e' class='xr-section-summary-in' type='checkbox'  checked><label for='section-85044860-f745-45dd-946d-8b817708a16e' class='xr-section-summary' >Attributes: <span>(1)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>_FillValue :</span></dt><dd>1.7976931348623157e+308</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.DataArray 'stackstac-69a5cbaa571c4ef68cbeaba52a270af4' (y: 1500, x: 1400)>\n",
       "array([[ 840.86486816,  837.69940186,  836.6675415 , ..., 1222.60290527,\n",
       "        1224.68151855, 1229.60192871],\n",
       "       [ 847.59912109,  843.43511963,  841.19110107, ..., 1226.77722168,\n",
       "        1229.88330078, 1235.05480957],\n",
       "       [ 853.09313965,  849.24804688,  847.56011963, ..., 1231.57214355,\n",
       "        1235.27197266, 1240.26989746],\n",
       "       ...,\n",
       "       [1057.71154785, 1056.44726562, 1053.34716797, ..., 1967.21203613,\n",
       "        1966.72473145, 1967.17895508],\n",
       "       [1053.78833008, 1052.40124512, 1049.59545898, ..., 1970.66943359,\n",
       "        1969.1439209 , 1969.60314941],\n",
       "       [1049.06298828, 1048.03967285, 1045.55126953, ..., 1973.11730957,\n",
       "        1971.60437012, 1971.93334961]])\n",
       "Coordinates:\n",
       "  * x                (x) float64 5.79e+05 5.79e+05 ... 5.93e+05 5.93e+05\n",
       "  * y                (y) float64 5.21e+06 5.21e+06 ... 5.195e+06 5.195e+06\n",
       "    threedep:region  <U7 'n40w130'\n",
       "    description      <U1849 'This tile of the 3D Elevation Program (3DEP) sea...\n",
       "    epsg             int64 5498\n",
       "    proj:shape       object {10812}\n",
       "    proj:epsg        int64 5498\n",
       "    band             <U4 'data'\n",
       "    gsd              int64 10\n",
       "    spatial_ref      int64 0\n",
       "Attributes:\n",
       "    _FillValue:  1.7976931348623157e+308"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# crop DEM to our area of interest and save it\n",
    "dem_corners = dict(x=slice(579000, 593000), y=slice(5210000, 5195000))\n",
    "\n",
    "dem_raster = dem_raster.sel(**dem_corners)\n",
    "dem_raster.rio.to_raster(dem_path)\n",
    "dem_raster"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "7082c7d4",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<img src=\"https://ai4edatasetspublicassets.blob.core.windows.net/assets/notebook-output/tutorials-rtc-qualitative-assessment.ipynb/11.png\"/>"
      ],
      "text/plain": [
       "<Figure size 720x504 with 2 Axes>"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 469,
       "width": 666
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "dem_raster.plot()\n",
    "_ = plt.title(\"DEM in UTM coordinates\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "4445aa2a-abf4-4162-ab0b-b22ae069e4e4",
   "metadata": {},
   "source": [
    "### Define GRD parameters"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "0c35e4ca",
   "metadata": {},
   "outputs": [],
   "source": [
    "grd_account_name = \"sentinel1euwest\"\n",
    "grd_storage_container = \"s1-grd\"\n",
    "grd_product_folder = f\"{grd_storage_container}/{product_folder}\"\n",
    "grd_local_path = os.path.join(tmp_dir, product_folder)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ac695ed8-10bb-4efe-a43e-372c289605e3",
   "metadata": {},
   "source": [
    "### Retrieve Sentinel-1 GRD"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "81012b36",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['s1-grd/GRD/2021/12/17/IW/DV/S1B_IW_GRDH_1SDV_20211217T141304_20211217T141329_030066_039705_9048/manifest.safe']"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "grd_token = planetary_computer.sas.get_token(\n",
    "    grd_account_name, grd_storage_container\n",
    ").token\n",
    "grd_fs = adlfs.AzureBlobFileSystem(grd_account_name, credential=grd_token)\n",
    "grd_fs.ls(f\"{grd_product_folder}/manifest.safe\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "0b1daaa9",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "/tmp/GRD/2021/12/17/IW/DV/S1B_IW_GRDH_1SDV_20211217T141304_20211217T141329_030066_039705_9048\n"
     ]
    }
   ],
   "source": [
    "grd_fs.get(grd_product_folder, grd_local_path, recursive=True)\n",
    "!ls -d {grd_local_path}"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5dbe585a-5467-445e-a9cb-cd8c5b66fe8e",
   "metadata": {},
   "source": [
    "### Process"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "54a3d505-a01c-4a88-a939-0169264609d0",
   "metadata": {},
   "source": [
    "#### GTC\n",
    "\n",
    "Here we compute the geometric terrain correction.\n",
    "\n",
    "Input parameters:\n",
    "- `product_urlpath`: product path\n",
    "- `measurement_group`: band to be processed in the form {swath}/{polarization} (see [xarray-sentinel](https://pypi.org/project/xarray-sentinel/) for more details)\n",
    "- `dem_urlpath`: path to the input DEM. sarsen supports all DEMs supported by GDAL/Proj for ECEF-translation. \n",
    "- `interp_method`: interpolation method, sarsen supports all interpolation methods supported by [xarray.Dataset.interp](https://xarray.pydata.org/en/stable/generated/xarray.Dataset.interp.html)\n",
    "- `chunks`: dask chunks\n",
    "- `output_urlpath`: output path\n",
    "\n",
    "The output is the input SAR image resampled on DEM coordinates"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "41add520",
   "metadata": {},
   "outputs": [],
   "source": [
    "gtc_path = os.path.basename(product_folder) + \".GTC.tif\"\n",
    "\n",
    "gtc = apps.terrain_correction(\n",
    "    grd_local_path,\n",
    "    measurement_group,\n",
    "    dem_path,\n",
    "    interp_method=\"nearest\",\n",
    "    override_product_files=\"{dirname}/{prefix}{swath}-{polarization}{ext}\",\n",
    "    chunks=2048,\n",
    "    output_urlpath=gtc_path,\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "27a84435",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<img src=\"https://ai4edatasetspublicassets.blob.core.windows.net/assets/notebook-output/tutorials-rtc-qualitative-assessment.ipynb/16.png\"/>"
      ],
      "text/plain": [
       "<Figure size 720x504 with 2 Axes>"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 447,
       "width": 629
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "_ = gtc.plot(vmax=0.4)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c7399622",
   "metadata": {},
   "source": [
    "#### RTC\n",
    "`sarsen` implements the radiometric terrain-correction [Gamma Flattening](https://ieeexplore.ieee.org/document/5752845) algorithm.\n",
    "\n",
    "#### Input parameters\n",
    "- `correct_radiometry`: `correct_radiometry`: default `None`. If `correct_radiometry=None`the radiometric terrain correction is not applied. `correct_radiometry=gamma_bilinear` applies the gamma flattening classic algorithm using bilinear interpolation to compute the weights. `correct_radiometry=gamma_nearest` applies the gamma flattening using nearest neighbours instead of bilinear interpolation. 'gamma_nearest' significantly reduces the processing time.\n",
    "- `grouping_area_factor`: scaling factor for the size of the image pixel where the areas are summed. By default, the `grouping_area_factor` is `(1, 1)`, which corresponds to Sentinel-1 input product pixel size. The `grouping_area_factor` shall be increased if the DEM resolution is lower than the Sentinel-1 input product resolution to avoid gaps and distortions the normalization factor. It can be also used to to speed up the computation or  the DEM resolution is lower than the Sentinel-1 input product resolution.  \n",
    "\n",
    "\n",
    "**Note**: The `grouping_area_factor` can be increased (i) to speed up the processing or (ii) when the input DEM resolution is low. The Gamma Flattening usually works properly if the pixel size of the input DEM is much smaller than the pixel size of the input Sentinel-1 product. Otherwise, the output may have radiometric distortions. This problem can be avoided by increasing the `grouping_area_factor`. Be aware that `grouping_area_factor` too high may degrade the final result."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "8416dd32-afa0-434a-8c9c-5d44bdd06cf9",
   "metadata": {},
   "outputs": [],
   "source": [
    "import warnings\n",
    "\n",
    "warnings.filterwarnings(\"ignore\", \"divide by zero\", RuntimeWarning)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "id": "bfdd3e96",
   "metadata": {},
   "outputs": [],
   "source": [
    "rtc_path = os.path.basename(grd_local_path) + \".RTC.tif\"\n",
    "\n",
    "rtc = apps.terrain_correction(\n",
    "    grd_local_path,\n",
    "    measurement_group,\n",
    "    dem_path,\n",
    "    interp_method=\"nearest\",\n",
    "    override_product_files=\"{dirname}/{prefix}{swath}-{polarization}{ext}\",\n",
    "    correct_radiometry=\"gamma_nearest\",\n",
    "    output_urlpath=rtc_path,\n",
    "    grouping_area_factor=(1, 1),\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "849a9019",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<img src=\"https://ai4edatasetspublicassets.blob.core.windows.net/assets/notebook-output/tutorials-rtc-qualitative-assessment.ipynb/19.png\"/>"
      ],
      "text/plain": [
       "<Figure size 720x504 with 2 Axes>"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 447,
       "width": 629
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "_ = rtc.plot(vmax=0.4)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "646acd93",
   "metadata": {},
   "source": [
    "### Gamma vs Local incidence Angle\n",
    "\n",
    "The backscattering depends on the ground material, the roughness and the local incidence angle. Small values of the local incidence angle correspond usually to higher values of backscattering, while high values of radar incidence angle corresponds to zones with low backscattering. The aim of the radiometric terrain correction is to remove the contribution due to the local incidence angle. \n",
    "For evaluating the goodness of the radiometric correction we will compare the incidence angle versus backscattering before and after the correction."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3bc92d66",
   "metadata": {},
   "source": [
    "#### Compute local incidence angle and ellipsoid incidence angle"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f93af76d",
   "metadata": {},
   "source": [
    "![SNOWFALL](https://forum.step.esa.int/uploads/default/original/3X/e/3/e372bf0efa0133c63977111b15614bc0c7d96fd8.png)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "e8499527-c404-4c45-b266-c0ec077db4d9",
   "metadata": {},
   "outputs": [],
   "source": [
    "dem = scene.open_dem_raster(dem_path)\n",
    "\n",
    "\n",
    "orbit_ecef = xr.open_dataset(\n",
    "    grd_local_path,\n",
    "    engine=\"sentinel-1\",\n",
    "    group=f\"{measurement_group}/orbit\",\n",
    "    override_product_files=\"{dirname}/{prefix}{swath}-{polarization}{ext}\",\n",
    ")\n",
    "dem_ecef = scene.convert_to_dem_ecef(dem, source_crs=t_srs)\n",
    "dem_ecef_1 = scene.convert_to_dem_ecef(dem + 1, source_crs=t_srs)\n",
    "ellissoid_normal = dem_ecef - dem_ecef_1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "a5f3dd32-007b-4c1c-9cda-88592fe78dd4",
   "metadata": {},
   "outputs": [],
   "source": [
    "acquisition = apps.simulate_acquisition(\n",
    "    dem_ecef, orbit_ecef.position, coordinate_conversion=None, correct_radiometry=True\n",
    ")\n",
    "oriented_area = scene.compute_dem_oriented_area(dem_ecef)\n",
    "dem_normal = -oriented_area / np.sqrt(xr.dot(oriented_area, oriented_area, dims=\"axis\"))\n",
    "\n",
    "# compute direction target sensor\n",
    "orbit_interpolator = orbit.OrbitPolyfitIterpolator.from_position(orbit_ecef.position)\n",
    "position_ecef = orbit_interpolator.position()\n",
    "velocity_ecef = orbit_interpolator.velocity()\n",
    "acquisition = geocoding.backward_geocode(dem_ecef, orbit_ecef.position, velocity_ecef)\n",
    "slant_range = np.sqrt((acquisition.dem_distance**2).sum(dim=\"axis\"))\n",
    "dem_direction = acquisition.dem_distance / slant_range"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "id": "920ee292-f695-44d4-9178-950fb5703a4b",
   "metadata": {},
   "outputs": [],
   "source": [
    "angle = np.arccos(xr.dot(dem_normal, dem_direction, dims=\"axis\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "id": "a2797f82-5471-4a75-9878-317e9b23865d",
   "metadata": {},
   "outputs": [],
   "source": [
    "ellipsoid_incidence = np.arccos(xr.dot(ellissoid_normal, dem_direction, dims=\"axis\"))"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5c69085f",
   "metadata": {},
   "source": [
    "#### Load Data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "id": "a45cd942",
   "metadata": {},
   "outputs": [],
   "source": [
    "ellipsoid_incidence = ellipsoid_incidence.load()\n",
    "angle = angle.load()\n",
    "gtc = gtc.load()\n",
    "rtc = rtc.load()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5311f17d-aa35-4725-8e28-dd91719cf91c",
   "metadata": {},
   "source": [
    "#### GTC and RTC distribution of foreslope, backslope, and flat terrain \n",
    "\n",
    "We classify taking into account the local incidence angle and the ellipsoid incidence angle\n",
    " the image in three zone:\n",
    "- foreslope, small incidence angles that correspond to high backscattering zones.\n",
    "- flat terrain\n",
    "- backslope, big incidence angles that correspond to low backscattering zones\n",
    "\n",
    "\n",
    "We will show that the radiometric correction mitigate the effect of the foreslope and backslope."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "id": "fec9e2e8-0dc7-49e5-a9ea-a12bbce248e9",
   "metadata": {},
   "outputs": [],
   "source": [
    "# foreslope backslope and flat terrain\n",
    "gtc_foreslope = xr.where(angle < ellipsoid_incidence - 0.20, gtc, np.nan)\n",
    "gtc_backslope = xr.where(angle > ellipsoid_incidence + 0.20, gtc, np.nan)\n",
    "gtc_flat = xr.where(\n",
    "    (angle < ellipsoid_incidence + 0.20) & (angle > ellipsoid_incidence - 0.20),\n",
    "    gtc,\n",
    "    np.nan,\n",
    ")\n",
    "\n",
    "rtc_foreslope = xr.where(angle < ellipsoid_incidence - 0.20, rtc, np.nan)\n",
    "rtc_backslope = xr.where(angle > ellipsoid_incidence + 0.20, rtc, np.nan)\n",
    "rtc_flat = xr.where(\n",
    "    (angle < ellipsoid_incidence + 0.20) & (angle > ellipsoid_incidence - 0.20),\n",
    "    rtc,\n",
    "    np.nan,\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "id": "e3abd619-34e2-4fec-845d-89e85adf4049",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<img src=\"https://ai4edatasetspublicassets.blob.core.windows.net/assets/notebook-output/tutorials-rtc-qualitative-assessment.ipynb/26.png\"/>"
      ],
      "text/plain": [
       "<Figure size 2160x648 with 6 Axes>"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 631,
       "width": 2109
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "f, axes = plt.subplots(nrows=1, ncols=3, figsize=(30, 9))\n",
    "\n",
    "dem_raster.plot(ax=axes[0])\n",
    "axes[0].grid(c=\"red\")\n",
    "axes[0].set_title(\"DEM\")\n",
    "\n",
    "gtc.plot(ax=axes[1], vmax=0.4)\n",
    "axes[1].grid(c=\"red\")\n",
    "axes[1].set_title(\"GTC\")\n",
    "\n",
    "\n",
    "angle.plot(ax=axes[2])\n",
    "axes[2].grid(c=\"red\")\n",
    "axes[2].set_title(\"Incidence Angle\")\n",
    "\n",
    "plt.tight_layout()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "id": "8923dfd6-b740-4a72-9f88-a67e07d47df6",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<img src=\"https://ai4edatasetspublicassets.blob.core.windows.net/assets/notebook-output/tutorials-rtc-qualitative-assessment.ipynb/27.png\"/>"
      ],
      "text/plain": [
       "<Figure size 2160x648 with 6 Axes>"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 631,
       "width": 2114
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "f, axes = plt.subplots(nrows=1, ncols=3, figsize=(30, 9))\n",
    "\n",
    "gtc_backslope.plot(ax=axes[0], vmax=0.4)\n",
    "axes[0].grid(c=\"red\")\n",
    "axes[0].set_title(\"GTC basckslope\")\n",
    "\n",
    "gtc_foreslope.plot(ax=axes[1], vmax=0.4)\n",
    "axes[1].grid(c=\"red\")\n",
    "axes[1].set_title(\"GTC foreslope\")\n",
    "\n",
    "\n",
    "gtc_flat.plot(ax=axes[2], vmax=0.4)\n",
    "axes[2].grid(c=\"red\")\n",
    "axes[2].set_title(\"GTC flat terrain\")\n",
    "\n",
    "\n",
    "plt.tight_layout()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "10bcc74a-902c-4fb3-ac85-f0012e307f68",
   "metadata": {},
   "source": [
    "#### GTC and RTC distribution of foreslope backslope and flat terrain in db"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "id": "12ef90a7",
   "metadata": {},
   "outputs": [],
   "source": [
    "gtc_backslope_db = np.log(gtc_backslope)\n",
    "gtc_foreslope_db = np.log(gtc_foreslope)\n",
    "gtc_flat_db = np.log(gtc_flat)\n",
    "\n",
    "rtc_backslope_db = np.log(rtc_backslope)\n",
    "rtc_foreslope_db = np.log(rtc_foreslope)\n",
    "rtc_flat_db = np.log(rtc_flat)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "id": "68081745-5198-441f-87bf-b9b7e49d81ed",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0.5, 0, 'gamma_T DB')"
      ]
     },
     "execution_count": 29,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"https://ai4edatasetspublicassets.blob.core.windows.net/assets/notebook-output/tutorials-rtc-qualitative-assessment.ipynb/29.png\"/>"
      ],
      "text/plain": [
       "<Figure size 2160x864 with 2 Axes>"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 724,
       "width": 1717
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "f, axes = plt.subplots(nrows=1, ncols=2, figsize=(30, 12))\n",
    "\n",
    "gtc_backslope_db_h = xr.plot.hist(\n",
    "    gtc_backslope_db,\n",
    "    range=[-10, 5],\n",
    "    bins=50,\n",
    "    color=[0, 1, 0, 0.4],\n",
    "    density=True,\n",
    "    ax=axes[0],\n",
    "    label=\"backslope\",\n",
    ")\n",
    "gtc_foreslope_db_h = xr.plot.hist(\n",
    "    gtc_foreslope_db,\n",
    "    range=[-10, 5],\n",
    "    bins=50,\n",
    "    color=[1, 0, 0, 0.4],\n",
    "    density=True,\n",
    "    ax=axes[0],\n",
    "    label=\"foreslope\",\n",
    ")\n",
    "gtc_flat_db_h = xr.plot.hist(\n",
    "    gtc_flat_db,\n",
    "    range=[-10, 5],\n",
    "    bins=50,\n",
    "    color=[0, 0, 1, 0.4],\n",
    "    density=True,\n",
    "    ax=axes[0],\n",
    "    label=\"flat\",\n",
    ")\n",
    "\n",
    "rtc_backslope_db_h = xr.plot.hist(\n",
    "    rtc_backslope_db,\n",
    "    range=[-10, 5],\n",
    "    bins=50,\n",
    "    color=[0, 1, 0, 0.4],\n",
    "    density=True,\n",
    "    ax=axes[1],\n",
    "    label=\"backslope\",\n",
    ")\n",
    "rtc_foreslope_db_h = xr.plot.hist(\n",
    "    rtc_foreslope_db,\n",
    "    range=[-10, 5],\n",
    "    bins=50,\n",
    "    color=[1, 0, 0, 0.4],\n",
    "    density=True,\n",
    "    ax=axes[1],\n",
    "    label=\"foreslope\",\n",
    ")\n",
    "rtc_flat_db_h = xr.plot.hist(\n",
    "    rtc_flat_db,\n",
    "    range=[-10, 5],\n",
    "    bins=50,\n",
    "    color=[0, 0, 1, 0.4],\n",
    "    density=True,\n",
    "    ax=axes[1],\n",
    "    label=\"flat\",\n",
    ")\n",
    "axes[0].legend()\n",
    "axes[1].legend()\n",
    "\n",
    "axes[0].set_title(\"GTC distribution\")\n",
    "axes[0].set_xlabel(\"beta DB\")\n",
    "\n",
    "axes[1].set_title(\"RTC distribution\")\n",
    "axes[1].set_xlabel(\"gamma_T DB\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "cff38443-b393-4ace-bf3f-bc64318347a1",
   "metadata": {},
   "source": [
    "**Before the radiometic correction the distribution of backslope flat terrain and foreslope are separate, after they  overlap completely.**"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "456b4611",
   "metadata": {},
   "source": [
    "#### RTC and GTC scatterplots incidence angle vs backcattering\n",
    "\n",
    "We compare the incidence_angle-backscattering scatterplots of the GTC and the RTC to show that in the former there is a correlation between the incidence angle and backscattering, while in the RTC the correlation has been removed."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "id": "82786f69",
   "metadata": {},
   "outputs": [],
   "source": [
    "def makeColours(vals):\n",
    "    colours = np.zeros((len(vals), 3))\n",
    "    norm = Normalize(vmin=vals.min(), vmax=vals.max())\n",
    "\n",
    "    # Can put any colormap you like here.\n",
    "    colours = [cm.ScalarMappable(norm=norm, cmap=\"plasma\").to_rgba(val) for val in vals]\n",
    "\n",
    "    return colours\n",
    "\n",
    "\n",
    "def plot_density(x, y, ax=None):\n",
    "    if ax is None:\n",
    "        ax = plt\n",
    "    mask = np.logical_and(np.isfinite(y), np.isfinite(x))\n",
    "    samples = np.stack([x[mask], y[mask]])\n",
    "\n",
    "    densObj = kde(samples)\n",
    "    colours = makeColours(densObj.evaluate(samples))\n",
    "    fig = ax.scatter(samples[0], samples[1], color=colours)\n",
    "    return fig"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "id": "92dc79c9-099b-4c28-b41d-a46ce68cdb87",
   "metadata": {},
   "outputs": [],
   "source": [
    "gtc_db = np.log(gtc).load()\n",
    "rtc_db = np.log(rtc).load()\n",
    "angle = angle.load()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "id": "498a010d-6465-4cce-89c9-1feeb5cd5a0d",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0.5, 0, 'local incidence angle')"
      ]
     },
     "execution_count": 32,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"https://ai4edatasetspublicassets.blob.core.windows.net/assets/notebook-output/tutorials-rtc-qualitative-assessment.ipynb/32.png\"/>"
      ],
      "text/plain": [
       "<Figure size 2160x864 with 2 Axes>"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 724,
       "width": 1734
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "f, axes = plt.subplots(nrows=1, ncols=2, figsize=(30, 12))\n",
    "\n",
    "fig = plot_density(\n",
    "    angle[::10, ::10].data.ravel(), gtc_db[::10, ::10].data.ravel(), ax=axes[0]\n",
    ")\n",
    "fig = plot_density(\n",
    "    angle[::10, ::10].data.ravel(), rtc_db[::10, ::10].data.ravel(), ax=axes[1]\n",
    ")\n",
    "\n",
    "axes[0].set_title(\"GTC\")\n",
    "axes[0].set_ylabel(\"beta DB\")\n",
    "axes[0].set_xlabel(\"local incidence angle\")\n",
    "\n",
    "axes[1].set_title(\"RTC\")\n",
    "axes[1].set_ylabel(\"gamma_T DB\")\n",
    "axes[1].set_xlabel(\"local incidence angle\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e9a66c4d-ec35-49a9-8b71-1594615d6403",
   "metadata": {},
   "source": [
    "**Note that in the first plot there is a clear correlation between the angle of incidence and backscattering, while in the second graph there is no correlation.**"
   ]
  }
 ],
 "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.5"
  },
  "widgets": {
   "application/vnd.jupyter.widget-state+json": {
    "state": {
     "9aa2d4b3f58a4da480ec2bb37e4a4d6b": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "1.2.0",
      "model_name": "LayoutModel",
      "state": {}
     },
     "c64ed80821d9497abe19ad21f01b5e93": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.5.0",
      "model_name": "VBoxModel",
      "state": {
       "layout": "IPY_MODEL_9aa2d4b3f58a4da480ec2bb37e4a4d6b"
      }
     }
    },
    "version_major": 2,
    "version_minor": 0
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}