Investigating EMIT L2A GLT in Python.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "from os.path import join, expanduser\n",
    "import numpy as np\n",
    "from affine import Affine\n",
    "from rasters import MultiRaster, RasterGrid\n",
    "from emit_tools import load_EMIT_swath, extract_GLT, apply_GLT, ortho_xr\n",
    "from emit_tools import FILL_VALUE as fill_value\n",
    "from emit_tools import GLT_NODATA_VALUE as GLT_nodata_value"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "filename = expanduser(join(\"~\", \"data\", \"EMIT_L2A_RFL_001_20240521T232728_2414216_003.nc\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body[data-theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block !important;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-index-preview {\n",
       "  grid-column: 2 / 5;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data,\n",
       ".xr-index-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data,\n",
       ".xr-index-data-in:checked ~ .xr-index-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-index-name div,\n",
       ".xr-index-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data,\n",
       ".xr-index-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2,\n",
       ".xr-no-icon {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.Dataset&gt;\n",
       "Dimensions:           (downtrack: 1280, crosstrack: 1242, wavelengths: 285,\n",
       "                       ortho_y: 1893, ortho_x: 2327)\n",
       "Coordinates:\n",
       "  * downtrack         (downtrack) int64 0 1 2 3 4 5 ... 1275 1276 1277 1278 1279\n",
       "  * crosstrack        (crosstrack) int64 0 1 2 3 4 ... 1237 1238 1239 1240 1241\n",
       "    lon               (downtrack, crosstrack) float64 ...\n",
       "    lat               (downtrack, crosstrack) float64 ...\n",
       "    elev              (downtrack, crosstrack) float64 ...\n",
       "    glt_x             (ortho_y, ortho_x) float64 ...\n",
       "    glt_y             (ortho_y, ortho_x) float64 ...\n",
       "  * wavelengths       (wavelengths) float32 381.0 388.4 ... 2.486e+03 2.493e+03\n",
       "    fwhm              (wavelengths) float32 ...\n",
       "    good_wavelengths  (wavelengths) float32 ...\n",
       "Dimensions without coordinates: ortho_y, ortho_x\n",
       "Data variables:\n",
       "    reflectance       (downtrack, crosstrack, wavelengths) float32 ...\n",
       "Attributes: (12/39)\n",
       "    ncei_template_version:             NCEI_NetCDF_Swath_Template_v2.0\n",
       "    summary:                           The Earth Surface Mineral Dust Source ...\n",
       "    keywords:                          Imaging Spectroscopy, minerals, EMIT, ...\n",
       "    Conventions:                       CF-1.63\n",
       "    sensor:                            EMIT (Earth Surface Mineral Dust Sourc...\n",
       "    instrument:                        EMIT\n",
       "    ...                                ...\n",
       "    spatialResolution:                 0.000542232520256367\n",
       "    spatial_ref:                       GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHER...\n",
       "    geotransform:                      [-1.16103646e+02  5.42232520e-04 -0.00...\n",
       "    day_night_flag:                    Day\n",
       "    title:                             EMIT L2A Estimated Surface Reflectance...\n",
       "    granule_id:                        EMIT_L2A_RFL_001_20240521T232728_24142...</pre><div class='xr-wrap' style='display:none'><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-00bac9e1-ccd7-4eee-a039-6454eb227e09' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-00bac9e1-ccd7-4eee-a039-6454eb227e09' class='xr-section-summary'  title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'><ul class='xr-dim-list'><li><span class='xr-has-index'>downtrack</span>: 1280</li><li><span class='xr-has-index'>crosstrack</span>: 1242</li><li><span class='xr-has-index'>wavelengths</span>: 285</li><li><span>ortho_y</span>: 1893</li><li><span>ortho_x</span>: 2327</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-2a3f6b74-8615-4fce-952e-090c5e8dfe51' class='xr-section-summary-in' type='checkbox'  checked><label for='section-2a3f6b74-8615-4fce-952e-090c5e8dfe51' 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'>downtrack</span></div><div class='xr-var-dims'>(downtrack)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>0 1 2 3 4 ... 1276 1277 1278 1279</div><input id='attrs-4d0704df-053d-4612-8b12-e55d7a9bbf37' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-4d0704df-053d-4612-8b12-e55d7a9bbf37' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5e4f62b2-7df9-4ceb-ae01-91d8a4c78694' class='xr-var-data-in' type='checkbox'><label for='data-5e4f62b2-7df9-4ceb-ae01-91d8a4c78694' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([   0,    1,    2, ..., 1277, 1278, 1279])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>crosstrack</span></div><div class='xr-var-dims'>(crosstrack)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>0 1 2 3 4 ... 1238 1239 1240 1241</div><input id='attrs-8ca1a3f6-cb1d-4be4-814f-049e88bd6e25' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-8ca1a3f6-cb1d-4be4-814f-049e88bd6e25' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-57179206-f452-487f-948b-64a4de45720d' class='xr-var-data-in' type='checkbox'><label for='data-57179206-f452-487f-948b-64a4de45720d' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([   0,    1,    2, ..., 1239, 1240, 1241])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>lon</span></div><div class='xr-var-dims'>(downtrack, crosstrack)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-7c4917ea-09e4-4193-b8ca-f7393663922a' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-7c4917ea-09e4-4193-b8ca-f7393663922a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-bcd4594d-dacb-4478-8b7c-12b57594f6a5' class='xr-var-data-in' type='checkbox'><label for='data-bcd4594d-dacb-4478-8b7c-12b57594f6a5' 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>long_name :</span></dt><dd>Longitude (WGS-84)</dd><dt><span>units :</span></dt><dd>degrees east</dd></dl></div><div class='xr-var-data'><pre>[1589760 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>lat</span></div><div class='xr-var-dims'>(downtrack, crosstrack)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-bd54c876-2466-4f26-a5af-70d5d93e6e06' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-bd54c876-2466-4f26-a5af-70d5d93e6e06' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5d70416d-b90e-489d-a41a-717a0bb3da92' class='xr-var-data-in' type='checkbox'><label for='data-5d70416d-b90e-489d-a41a-717a0bb3da92' 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>long_name :</span></dt><dd>Latitude (WGS-84)</dd><dt><span>units :</span></dt><dd>degrees north</dd></dl></div><div class='xr-var-data'><pre>[1589760 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>elev</span></div><div class='xr-var-dims'>(downtrack, crosstrack)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-d372038e-dfa5-49f6-af1a-8428e14d4688' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-d372038e-dfa5-49f6-af1a-8428e14d4688' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-75dd26a1-13c0-4234-a280-c26dffe4b271' class='xr-var-data-in' type='checkbox'><label for='data-75dd26a1-13c0-4234-a280-c26dffe4b271' 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>long_name :</span></dt><dd>Surface Elevation</dd><dt><span>units :</span></dt><dd>m</dd></dl></div><div class='xr-var-data'><pre>[1589760 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>glt_x</span></div><div class='xr-var-dims'>(ortho_y, ortho_x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-a9daf9c6-37d8-46f1-9d36-a63cdef63899' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-a9daf9c6-37d8-46f1-9d36-a63cdef63899' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-57852105-54a4-46c5-85c6-d9adbbe59182' class='xr-var-data-in' type='checkbox'><label for='data-57852105-54a4-46c5-85c6-d9adbbe59182' 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>long_name :</span></dt><dd>GLT Sample Lookup</dd><dt><span>units :</span></dt><dd>pixel location</dd></dl></div><div class='xr-var-data'><pre>[4405011 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>glt_y</span></div><div class='xr-var-dims'>(ortho_y, ortho_x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-06f8f36a-0709-419f-b89e-d2e9a82155e2' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-06f8f36a-0709-419f-b89e-d2e9a82155e2' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-12ba0039-e63e-410e-b56a-546d3d8e519c' class='xr-var-data-in' type='checkbox'><label for='data-12ba0039-e63e-410e-b56a-546d3d8e519c' 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>long_name :</span></dt><dd>GLT Line Lookup</dd><dt><span>units :</span></dt><dd>pixel location</dd></dl></div><div class='xr-var-data'><pre>[4405011 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>wavelengths</span></div><div class='xr-var-dims'>(wavelengths)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>381.0 388.4 ... 2.486e+03 2.493e+03</div><input id='attrs-977587e6-8eb0-49e3-a555-7ebc45392e17' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-977587e6-8eb0-49e3-a555-7ebc45392e17' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5302441f-9c83-4b17-a648-b2798c88e6f2' class='xr-var-data-in' type='checkbox'><label for='data-5302441f-9c83-4b17-a648-b2798c88e6f2' 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>long_name :</span></dt><dd>Wavelength Centers</dd><dt><span>units :</span></dt><dd>nm</dd></dl></div><div class='xr-var-data'><pre>array([ 381.00558,  388.4092 ,  395.81583, ..., 2478.153  , 2485.5386 ,\n",
       "       2492.9238 ], dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>fwhm</span></div><div class='xr-var-dims'>(wavelengths)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-4393e961-de82-4baf-be15-b108cb12f676' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-4393e961-de82-4baf-be15-b108cb12f676' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4599cf19-30b4-4f8b-a9eb-9c7a4b7353ee' class='xr-var-data-in' type='checkbox'><label for='data-4599cf19-30b4-4f8b-a9eb-9c7a4b7353ee' 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>long_name :</span></dt><dd>Full Width at Half Max</dd><dt><span>units :</span></dt><dd>nm</dd></dl></div><div class='xr-var-data'><pre>[285 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>good_wavelengths</span></div><div class='xr-var-dims'>(wavelengths)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-08ef697f-223e-4005-8ea1-555d433d887e' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-08ef697f-223e-4005-8ea1-555d433d887e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4a762731-2ea6-42e7-9226-142c50f70340' class='xr-var-data-in' type='checkbox'><label for='data-4a762731-2ea6-42e7-9226-142c50f70340' 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>long_name :</span></dt><dd>Wavelengths where reflectance is useable: 1 = good data, 0 = bad data</dd><dt><span>units :</span></dt><dd>unitless</dd></dl></div><div class='xr-var-data'><pre>[285 values with dtype=float32]</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-7b8dccc9-6ed5-44cc-b835-4026abef62da' class='xr-section-summary-in' type='checkbox'  checked><label for='section-7b8dccc9-6ed5-44cc-b835-4026abef62da' class='xr-section-summary' >Data variables: <span>(1)</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>reflectance</span></div><div class='xr-var-dims'>(downtrack, crosstrack, wavelengths)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-f9bf580d-6c0f-4c7e-a942-58b69244567b' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-f9bf580d-6c0f-4c7e-a942-58b69244567b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-2f76134b-a0d7-4b35-8d49-603b3303c608' class='xr-var-data-in' type='checkbox'><label for='data-2f76134b-a0d7-4b35-8d49-603b3303c608' 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>long_name :</span></dt><dd>Surface Reflectance</dd><dt><span>units :</span></dt><dd>unitless</dd></dl></div><div class='xr-var-data'><pre>[453081600 values with dtype=float32]</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-790c64da-ba41-4937-b275-ff03669af119' class='xr-section-summary-in' type='checkbox'  ><label for='section-790c64da-ba41-4937-b275-ff03669af119' class='xr-section-summary' >Indexes: <span>(3)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-index-name'><div>downtrack</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-360e5515-c8a0-4fb1-ac0b-b790177da19b' class='xr-index-data-in' type='checkbox'/><label for='index-360e5515-c8a0-4fb1-ac0b-b790177da19b' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Index([   0,    1,    2,    3,    4,    5,    6,    7,    8,    9,\n",
       "       ...\n",
       "       1270, 1271, 1272, 1273, 1274, 1275, 1276, 1277, 1278, 1279],\n",
       "      dtype=&#x27;int64&#x27;, name=&#x27;downtrack&#x27;, length=1280))</pre></div></li><li class='xr-var-item'><div class='xr-index-name'><div>crosstrack</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-013060a9-54da-4cee-bc36-142b0e6089c4' class='xr-index-data-in' type='checkbox'/><label for='index-013060a9-54da-4cee-bc36-142b0e6089c4' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Index([   0,    1,    2,    3,    4,    5,    6,    7,    8,    9,\n",
       "       ...\n",
       "       1232, 1233, 1234, 1235, 1236, 1237, 1238, 1239, 1240, 1241],\n",
       "      dtype=&#x27;int64&#x27;, name=&#x27;crosstrack&#x27;, length=1242))</pre></div></li><li class='xr-var-item'><div class='xr-index-name'><div>wavelengths</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-5e469e47-8377-48f3-a538-9c31b74901ca' class='xr-index-data-in' type='checkbox'/><label for='index-5e469e47-8377-48f3-a538-9c31b74901ca' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Index([ 381.0055847167969,  388.4092102050781,  395.8158264160156,\n",
       "       403.22540283203125, 410.63800048828125,  418.0535888671875,\n",
       "        425.4721374511719,  432.8927001953125,  440.3172607421875,\n",
       "        447.7427978515625,\n",
       "       ...\n",
       "        2426.440185546875,  2433.830322265625,   2441.21826171875,\n",
       "          2448.6064453125,  2455.994384765625,  2463.381591796875,\n",
       "        2470.767822265625,  2478.153076171875,   2485.53857421875,\n",
       "           2492.923828125],\n",
       "      dtype=&#x27;float32&#x27;, name=&#x27;wavelengths&#x27;, length=285))</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-f96c90dd-9d47-434d-8cd8-f2b678cbe31a' class='xr-section-summary-in' type='checkbox'  ><label for='section-f96c90dd-9d47-434d-8cd8-f2b678cbe31a' class='xr-section-summary' >Attributes: <span>(39)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>ncei_template_version :</span></dt><dd>NCEI_NetCDF_Swath_Template_v2.0</dd><dt><span>summary :</span></dt><dd>The Earth Surface Mineral Dust Source Investigation (EMIT) is an Earth Ventures-Instrument (EVI-4) Mission that maps the surface mineralogy of arid dust source regions via imaging spectroscopy in the visible and short-wave infrared (VSWIR). Installed on the International Space Station (ISS), the EMIT instrument is a Dyson imaging spectrometer that uses contiguous spectroscopic measurements from 410 to 2450 nm to resolve absoprtion features of iron oxides, clays, sulfates, carbonates, and other dust-forming minerals. During its one-year mission, EMIT will observe the sunlit Earth&#x27;s dust source regions that occur within +/-52° latitude and produce maps of the source regions that can be used to improve forecasts of the role of mineral dust in the radiative forcing (warming or cooling) of the atmosphere.\\n\\nThis file contains L2A estimated surface reflectances and geolocation data. Reflectance estimates are created using an Optimal Estimation technique - see ATBD for details. Reflectance values are reported as fractions (relative to 1). Geolocation data (latitude, longitude, height) and a lookup table to project the data are also included.</dd><dt><span>keywords :</span></dt><dd>Imaging Spectroscopy, minerals, EMIT, dust, radiative forcing</dd><dt><span>Conventions :</span></dt><dd>CF-1.63</dd><dt><span>sensor :</span></dt><dd>EMIT (Earth Surface Mineral Dust Source Investigation)</dd><dt><span>instrument :</span></dt><dd>EMIT</dd><dt><span>platform :</span></dt><dd>ISS</dd><dt><span>institution :</span></dt><dd>NASA Jet Propulsion Laboratory/California Institute of Technology</dd><dt><span>license :</span></dt><dd>https://science.nasa.gov/earth-science/earth-science-data/data-information-policy/</dd><dt><span>naming_authority :</span></dt><dd>LPDAAC</dd><dt><span>date_created :</span></dt><dd>2024-05-23T15:58:15Z</dd><dt><span>keywords_vocabulary :</span></dt><dd>NASA Global Change Master Directory (GCMD) Science Keywords</dd><dt><span>stdname_vocabulary :</span></dt><dd>NetCDF Climate and Forecast (CF) Metadata Convention</dd><dt><span>creator_name :</span></dt><dd>Jet Propulsion Laboratory/California Institute of Technology</dd><dt><span>creator_url :</span></dt><dd>https://earth.jpl.nasa.gov/emit/</dd><dt><span>project :</span></dt><dd>Earth Surface Mineral Dust Source Investigation</dd><dt><span>project_url :</span></dt><dd>https://earth.jpl.nasa.gov/emit/</dd><dt><span>publisher_name :</span></dt><dd>NASA LPDAAC</dd><dt><span>publisher_url :</span></dt><dd>https://lpdaac.usgs.gov</dd><dt><span>publisher_email :</span></dt><dd>lpdaac@usgs.gov</dd><dt><span>identifier_product_doi_authority :</span></dt><dd>https://doi.org</dd><dt><span>flight_line :</span></dt><dd>emit20240521t232728_o14216_s000</dd><dt><span>time_coverage_start :</span></dt><dd>2024-05-21T23:27:28+0000</dd><dt><span>time_coverage_end :</span></dt><dd>2024-05-21T23:27:40+0000</dd><dt><span>software_build_version :</span></dt><dd>010618</dd><dt><span>software_delivery_version :</span></dt><dd>010618</dd><dt><span>product_version :</span></dt><dd>V001</dd><dt><span>history :</span></dt><dd>PGE Run Command: {python /beegfs/store/emit/ops/repos/emit-sds-l2a/spectrum_quality.py /tmp/emit/ops/emit20240521t232728_emit.L2AReflectance_20240523t023036/output/emit20240521t232728_rfl /tmp/emit/ops/emit20240521t232728_emit.L2AReflectance_20240523t023036/output/emit20240521t232728_rfl_quality.txt}, PGE Input Files: {radiance_file=/beegfs/store/emit/ops/data/acquisitions/20240521/emit20240521t232728/l1b/emit20240521t232728_o14216_s000_l1b_rdn_b0106_v01.img, pixel_locations_file=/beegfs/store/emit/ops/data/acquisitions/20240521/emit20240521t232728/l1b/emit20240521t232728_o14216_s000_l1b_loc_b0106_v01.img, observation_parameters_file=/beegfs/store/emit/ops/data/acquisitions/20240521/emit20240521t232728/l1b/emit20240521t232728_o14216_s000_l1b_obs_b0106_v01.img, surface_model_config=/beegfs/store/emit/ops/repos/emit-sds-l2a/surface/surface_20221020.json}</dd><dt><span>crosstrack_orientation :</span></dt><dd>as seen on ground</dd><dt><span>easternmost_longitude :</span></dt><dd>-114.84187059769744</dd><dt><span>northernmost_latitude :</span></dt><dd>33.9515655017347</dd><dt><span>westernmost_longitude :</span></dt><dd>-116.103645672334</dd><dt><span>southernmost_latitude :</span></dt><dd>32.925119340889395</dd><dt><span>spatialResolution :</span></dt><dd>0.000542232520256367</dd><dt><span>spatial_ref :</span></dt><dd>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;]],AXIS[&quot;Latitude&quot;,NORTH],AXIS[&quot;Longitude&quot;,EAST],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]]</dd><dt><span>geotransform :</span></dt><dd>[-1.16103646e+02  5.42232520e-04 -0.00000000e+00  3.39515655e+01\n",
       " -0.00000000e+00 -5.42232520e-04]</dd><dt><span>day_night_flag :</span></dt><dd>Day</dd><dt><span>title :</span></dt><dd>EMIT L2A Estimated Surface Reflectance 60 m V001</dd><dt><span>granule_id :</span></dt><dd>EMIT_L2A_RFL_001_20240521T232728_2414216_003</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset>\n",
       "Dimensions:           (downtrack: 1280, crosstrack: 1242, wavelengths: 285,\n",
       "                       ortho_y: 1893, ortho_x: 2327)\n",
       "Coordinates:\n",
       "  * downtrack         (downtrack) int64 0 1 2 3 4 5 ... 1275 1276 1277 1278 1279\n",
       "  * crosstrack        (crosstrack) int64 0 1 2 3 4 ... 1237 1238 1239 1240 1241\n",
       "    lon               (downtrack, crosstrack) float64 ...\n",
       "    lat               (downtrack, crosstrack) float64 ...\n",
       "    elev              (downtrack, crosstrack) float64 ...\n",
       "    glt_x             (ortho_y, ortho_x) float64 ...\n",
       "    glt_y             (ortho_y, ortho_x) float64 ...\n",
       "  * wavelengths       (wavelengths) float32 381.0 388.4 ... 2.486e+03 2.493e+03\n",
       "    fwhm              (wavelengths) float32 ...\n",
       "    good_wavelengths  (wavelengths) float32 ...\n",
       "Dimensions without coordinates: ortho_y, ortho_x\n",
       "Data variables:\n",
       "    reflectance       (downtrack, crosstrack, wavelengths) float32 ...\n",
       "Attributes: (12/39)\n",
       "    ncei_template_version:             NCEI_NetCDF_Swath_Template_v2.0\n",
       "    summary:                           The Earth Surface Mineral Dust Source ...\n",
       "    keywords:                          Imaging Spectroscopy, minerals, EMIT, ...\n",
       "    Conventions:                       CF-1.63\n",
       "    sensor:                            EMIT (Earth Surface Mineral Dust Sourc...\n",
       "    instrument:                        EMIT\n",
       "    ...                                ...\n",
       "    spatialResolution:                 0.000542232520256367\n",
       "    spatial_ref:                       GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHER...\n",
       "    geotransform:                      [-1.16103646e+02  5.42232520e-04 -0.00...\n",
       "    day_night_flag:                    Day\n",
       "    title:                             EMIT L2A Estimated Surface Reflectance...\n",
       "    granule_id:                        EMIT_L2A_RFL_001_20240521T232728_24142..."
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "swath_ds = load_EMIT_swath(filename)\n",
    "swath_ds"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[[0, 0],\n",
       "        [0, 0],\n",
       "        [0, 0],\n",
       "        ...,\n",
       "        [0, 0],\n",
       "        [0, 0],\n",
       "        [0, 0]],\n",
       "\n",
       "       [[0, 0],\n",
       "        [0, 0],\n",
       "        [0, 0],\n",
       "        ...,\n",
       "        [0, 0],\n",
       "        [0, 0],\n",
       "        [0, 0]],\n",
       "\n",
       "       [[0, 0],\n",
       "        [0, 0],\n",
       "        [0, 0],\n",
       "        ...,\n",
       "        [0, 0],\n",
       "        [0, 0],\n",
       "        [0, 0]],\n",
       "\n",
       "       ...,\n",
       "\n",
       "       [[0, 0],\n",
       "        [0, 0],\n",
       "        [0, 0],\n",
       "        ...,\n",
       "        [0, 0],\n",
       "        [0, 0],\n",
       "        [0, 0]],\n",
       "\n",
       "       [[0, 0],\n",
       "        [0, 0],\n",
       "        [0, 0],\n",
       "        ...,\n",
       "        [0, 0],\n",
       "        [0, 0],\n",
       "        [0, 0]],\n",
       "\n",
       "       [[0, 0],\n",
       "        [0, 0],\n",
       "        [0, 0],\n",
       "        ...,\n",
       "        [0, 0],\n",
       "        [0, 0],\n",
       "        [0, 0]]])"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "GLT_array = extract_GLT(swath_ds)\n",
    "GLT_array"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "numpy.ndarray"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "type(GLT_array)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(1893, 2327, 2)"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "GLT_array.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body[data-theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block !important;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-index-preview {\n",
       "  grid-column: 2 / 5;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data,\n",
       ".xr-index-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data,\n",
       ".xr-index-data-in:checked ~ .xr-index-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-index-name div,\n",
       ".xr-index-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data,\n",
       ".xr-index-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2,\n",
       ".xr-no-icon {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.Dataset&gt;\n",
       "Dimensions:           (latitude: 1893, longitude: 2327, wavelengths: 285)\n",
       "Coordinates:\n",
       "  * wavelengths       (wavelengths) float32 381.0 388.4 ... 2.486e+03 2.493e+03\n",
       "    fwhm              (wavelengths) float32 ...\n",
       "    good_wavelengths  (wavelengths) float32 ...\n",
       "  * latitude          (latitude) float64 33.95 33.95 33.95 ... 32.93 32.93 32.93\n",
       "  * longitude         (longitude) float64 -116.1 -116.1 -116.1 ... -114.8 -114.8\n",
       "    elev              (latitude, longitude) float32 -9.999e+03 ... -9.999e+03\n",
       "    spatial_ref       int64 0\n",
       "Data variables:\n",
       "    reflectance       (latitude, longitude, wavelengths) float32 -9.999e+03 ....\n",
       "Attributes: (12/39)\n",
       "    ncei_template_version:             NCEI_NetCDF_Swath_Template_v2.0\n",
       "    summary:                           The Earth Surface Mineral Dust Source ...\n",
       "    keywords:                          Imaging Spectroscopy, minerals, EMIT, ...\n",
       "    Conventions:                       CF-1.63\n",
       "    sensor:                            EMIT (Earth Surface Mineral Dust Sourc...\n",
       "    instrument:                        EMIT\n",
       "    ...                                ...\n",
       "    spatialResolution:                 0.000542232520256367\n",
       "    spatial_ref:                       GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHER...\n",
       "    geotransform:                      [-1.16103646e+02  5.42232520e-04 -0.00...\n",
       "    day_night_flag:                    Day\n",
       "    title:                             EMIT L2A Estimated Surface Reflectance...\n",
       "    granule_id:                        EMIT_L2A_RFL_001_20240521T232728_24142...</pre><div class='xr-wrap' style='display:none'><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-d0cb90d3-28ea-4cd4-aa8d-5f2c7a275f05' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-d0cb90d3-28ea-4cd4-aa8d-5f2c7a275f05' class='xr-section-summary'  title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'><ul class='xr-dim-list'><li><span class='xr-has-index'>latitude</span>: 1893</li><li><span class='xr-has-index'>longitude</span>: 2327</li><li><span class='xr-has-index'>wavelengths</span>: 285</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-aa7b7024-5786-4cf3-96cc-b5771f736d0e' class='xr-section-summary-in' type='checkbox'  checked><label for='section-aa7b7024-5786-4cf3-96cc-b5771f736d0e' class='xr-section-summary' >Coordinates: <span>(7)</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'>wavelengths</span></div><div class='xr-var-dims'>(wavelengths)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>381.0 388.4 ... 2.486e+03 2.493e+03</div><input id='attrs-3a8b308a-2670-4f60-bd8d-baea4d532d1c' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-3a8b308a-2670-4f60-bd8d-baea4d532d1c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4d6dcc70-e4a2-4db5-b142-feca7c72eea0' class='xr-var-data-in' type='checkbox'><label for='data-4d6dcc70-e4a2-4db5-b142-feca7c72eea0' 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>long_name :</span></dt><dd>Wavelength Centers</dd><dt><span>units :</span></dt><dd>nm</dd></dl></div><div class='xr-var-data'><pre>array([ 381.00558,  388.4092 ,  395.81583, ..., 2478.153  , 2485.5386 ,\n",
       "       2492.9238 ], dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>fwhm</span></div><div class='xr-var-dims'>(wavelengths)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-f659a5d0-a6c4-466d-931a-8854e21caa41' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-f659a5d0-a6c4-466d-931a-8854e21caa41' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-25b6480f-8349-4ee2-93cb-515f8b00f15b' class='xr-var-data-in' type='checkbox'><label for='data-25b6480f-8349-4ee2-93cb-515f8b00f15b' 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>long_name :</span></dt><dd>Full Width at Half Max</dd><dt><span>units :</span></dt><dd>nm</dd></dl></div><div class='xr-var-data'><pre>[285 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>good_wavelengths</span></div><div class='xr-var-dims'>(wavelengths)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-3ebf1fc8-bc65-44ac-969b-af7142e13c35' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-3ebf1fc8-bc65-44ac-969b-af7142e13c35' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-2ab5f0ed-5071-4a57-8d98-b73bd5a48ccb' class='xr-var-data-in' type='checkbox'><label for='data-2ab5f0ed-5071-4a57-8d98-b73bd5a48ccb' 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>long_name :</span></dt><dd>Wavelengths where reflectance is useable: 1 = good data, 0 = bad data</dd><dt><span>units :</span></dt><dd>unitless</dd></dl></div><div class='xr-var-data'><pre>[285 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>latitude</span></div><div class='xr-var-dims'>(latitude)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>33.95 33.95 33.95 ... 32.93 32.93</div><input id='attrs-2fbbb072-080b-4315-8b99-a91a324e7144' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-2fbbb072-080b-4315-8b99-a91a324e7144' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-92ba021b-6756-4cfb-b442-a5fcd407fe5e' class='xr-var-data-in' type='checkbox'><label for='data-92ba021b-6756-4cfb-b442-a5fcd407fe5e' 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>long_name :</span></dt><dd>Latitude (WGS-84)</dd><dt><span>units :</span></dt><dd>degrees north</dd></dl></div><div class='xr-var-data'><pre>array([33.951294, 33.950752, 33.95021 , ..., 32.926475, 32.925933, 32.92539 ])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>longitude</span></div><div class='xr-var-dims'>(longitude)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-116.1 -116.1 ... -114.8 -114.8</div><input id='attrs-bc2330d9-0603-400d-ac00-3e64a934099f' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-bc2330d9-0603-400d-ac00-3e64a934099f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-6ce646e1-d115-4292-bd3c-b20fd7c8699f' class='xr-var-data-in' type='checkbox'><label for='data-6ce646e1-d115-4292-bd3c-b20fd7c8699f' 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>long_name :</span></dt><dd>Longitude (WGS-84)</dd><dt><span>units :</span></dt><dd>degrees east</dd></dl></div><div class='xr-var-data'><pre>array([-116.103375, -116.102832, -116.10229 , ..., -114.843226, -114.842684,\n",
       "       -114.842142])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>elev</span></div><div class='xr-var-dims'>(latitude, longitude)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>-9.999e+03 ... -9.999e+03</div><input id='attrs-657c3610-6d1d-4b04-a500-79c5a34891d2' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-657c3610-6d1d-4b04-a500-79c5a34891d2' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-57019a5a-15f1-4e01-b1ed-3ad36bf7d309' class='xr-var-data-in' type='checkbox'><label for='data-57019a5a-15f1-4e01-b1ed-3ad36bf7d309' 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>long_name :</span></dt><dd>Surface Elevation</dd><dt><span>units :</span></dt><dd>m</dd></dl></div><div class='xr-var-data'><pre>array([[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "       [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "       [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "       ...,\n",
       "       [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "       [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "       [-9999., -9999., -9999., ..., -9999., -9999., -9999.]],\n",
       "      dtype=float32)</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-06711c32-de77-4bf2-a649-6c899465b150' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-06711c32-de77-4bf2-a649-6c899465b150' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-bcc145ab-4cd8-434d-88c2-8a7ad9826878' class='xr-var-data-in' type='checkbox'><label for='data-bcc145ab-4cd8-434d-88c2-8a7ad9826878' 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>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;]],AXIS[&quot;Latitude&quot;,NORTH],AXIS[&quot;Longitude&quot;,EAST],AUTHORITY[&quot;EPSG&quot;,&quot;4326&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>grid_mapping_name :</span></dt><dd>latitude_longitude</dd><dt><span>spatial_ref :</span></dt><dd>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;]],AXIS[&quot;Latitude&quot;,NORTH],AXIS[&quot;Longitude&quot;,EAST],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]]</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-38e931b2-4fe3-48c5-9d20-66cd51b8e9ce' class='xr-section-summary-in' type='checkbox'  checked><label for='section-38e931b2-4fe3-48c5-9d20-66cd51b8e9ce' class='xr-section-summary' >Data variables: <span>(1)</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>reflectance</span></div><div class='xr-var-dims'>(latitude, longitude, wavelengths)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>-9.999e+03 ... -9.999e+03</div><input id='attrs-3c2efaa0-8be8-4bd3-ba3d-2c821bad3ee8' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-3c2efaa0-8be8-4bd3-ba3d-2c821bad3ee8' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-79a08087-a645-40f1-9865-8a8fc2dc8883' class='xr-var-data-in' type='checkbox'><label for='data-79a08087-a645-40f1-9865-8a8fc2dc8883' 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>long_name :</span></dt><dd>Surface Reflectance</dd><dt><span>units :</span></dt><dd>unitless</dd></dl></div><div class='xr-var-data'><pre>array([[[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        ...,\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.]],\n",
       "\n",
       "       [[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        ...,\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.]],\n",
       "\n",
       "       [[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        ...,\n",
       "...\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.]],\n",
       "\n",
       "       [[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        ...,\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.]],\n",
       "\n",
       "       [[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        ...,\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.]]],\n",
       "      dtype=float32)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-6f19ed9b-2e19-40e5-aaf5-dbad0907036d' class='xr-section-summary-in' type='checkbox'  ><label for='section-6f19ed9b-2e19-40e5-aaf5-dbad0907036d' class='xr-section-summary' >Indexes: <span>(3)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-index-name'><div>wavelengths</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-9cecf3c4-8fe3-40d7-92e6-2fb180bae432' class='xr-index-data-in' type='checkbox'/><label for='index-9cecf3c4-8fe3-40d7-92e6-2fb180bae432' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Index([ 381.0055847167969,  388.4092102050781,  395.8158264160156,\n",
       "       403.22540283203125, 410.63800048828125,  418.0535888671875,\n",
       "        425.4721374511719,  432.8927001953125,  440.3172607421875,\n",
       "        447.7427978515625,\n",
       "       ...\n",
       "        2426.440185546875,  2433.830322265625,   2441.21826171875,\n",
       "          2448.6064453125,  2455.994384765625,  2463.381591796875,\n",
       "        2470.767822265625,  2478.153076171875,   2485.53857421875,\n",
       "           2492.923828125],\n",
       "      dtype=&#x27;float32&#x27;, name=&#x27;wavelengths&#x27;, length=285))</pre></div></li><li class='xr-var-item'><div class='xr-index-name'><div>latitude</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-d0af9f3a-0e8b-48a5-b4e1-37722d4ea5fa' class='xr-index-data-in' type='checkbox'/><label for='index-d0af9f3a-0e8b-48a5-b4e1-37722d4ea5fa' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Index([ 33.95129438547457, 33.950752152954315, 33.950209920434055,\n",
       "         33.9496676879138,  33.94912545539354,  33.94858322287329,\n",
       "        33.94804099035303, 33.947498757832776, 33.946956525312515,\n",
       "        33.94641429279226,\n",
       "       ...\n",
       "        32.93027054983183,  32.92972831731157,  32.92918608479132,\n",
       "        32.92864385227106, 32.928101619750805, 32.927559387230545,\n",
       "        32.92701715471029,  32.92647492219003,  32.92593268966978,\n",
       "       32.925390457149526],\n",
       "      dtype=&#x27;float64&#x27;, name=&#x27;latitude&#x27;, length=1893))</pre></div></li><li class='xr-var-item'><div class='xr-index-name'><div>longitude</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-3d874034-880b-41a6-92b4-70f167edb5be' class='xr-index-data-in' type='checkbox'/><label for='index-3d874034-880b-41a6-92b4-70f167edb5be' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Index([-116.10337455607387, -116.10283232355361, -116.10229009103337,\n",
       "       -116.10174785851311, -116.10120562599285, -116.10066339347259,\n",
       "       -116.10012116095234, -116.09957892843208, -116.09903669591182,\n",
       "       -116.09849446339156,\n",
       "       ...\n",
       "       -114.84702180663987, -114.84647957411961, -114.84593734159935,\n",
       "        -114.8453951090791, -114.84485287655885, -114.84431064403859,\n",
       "       -114.84376841151834, -114.84322617899808, -114.84268394647782,\n",
       "       -114.84214171395756],\n",
       "      dtype=&#x27;float64&#x27;, name=&#x27;longitude&#x27;, length=2327))</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-31e31a53-b7e8-4e48-9832-85909c0d0937' class='xr-section-summary-in' type='checkbox'  ><label for='section-31e31a53-b7e8-4e48-9832-85909c0d0937' class='xr-section-summary' >Attributes: <span>(39)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>ncei_template_version :</span></dt><dd>NCEI_NetCDF_Swath_Template_v2.0</dd><dt><span>summary :</span></dt><dd>The Earth Surface Mineral Dust Source Investigation (EMIT) is an Earth Ventures-Instrument (EVI-4) Mission that maps the surface mineralogy of arid dust source regions via imaging spectroscopy in the visible and short-wave infrared (VSWIR). Installed on the International Space Station (ISS), the EMIT instrument is a Dyson imaging spectrometer that uses contiguous spectroscopic measurements from 410 to 2450 nm to resolve absoprtion features of iron oxides, clays, sulfates, carbonates, and other dust-forming minerals. During its one-year mission, EMIT will observe the sunlit Earth&#x27;s dust source regions that occur within +/-52° latitude and produce maps of the source regions that can be used to improve forecasts of the role of mineral dust in the radiative forcing (warming or cooling) of the atmosphere.\\n\\nThis file contains L2A estimated surface reflectances and geolocation data. Reflectance estimates are created using an Optimal Estimation technique - see ATBD for details. Reflectance values are reported as fractions (relative to 1). Geolocation data (latitude, longitude, height) and a lookup table to project the data are also included.</dd><dt><span>keywords :</span></dt><dd>Imaging Spectroscopy, minerals, EMIT, dust, radiative forcing</dd><dt><span>Conventions :</span></dt><dd>CF-1.63</dd><dt><span>sensor :</span></dt><dd>EMIT (Earth Surface Mineral Dust Source Investigation)</dd><dt><span>instrument :</span></dt><dd>EMIT</dd><dt><span>platform :</span></dt><dd>ISS</dd><dt><span>institution :</span></dt><dd>NASA Jet Propulsion Laboratory/California Institute of Technology</dd><dt><span>license :</span></dt><dd>https://science.nasa.gov/earth-science/earth-science-data/data-information-policy/</dd><dt><span>naming_authority :</span></dt><dd>LPDAAC</dd><dt><span>date_created :</span></dt><dd>2024-05-23T15:58:15Z</dd><dt><span>keywords_vocabulary :</span></dt><dd>NASA Global Change Master Directory (GCMD) Science Keywords</dd><dt><span>stdname_vocabulary :</span></dt><dd>NetCDF Climate and Forecast (CF) Metadata Convention</dd><dt><span>creator_name :</span></dt><dd>Jet Propulsion Laboratory/California Institute of Technology</dd><dt><span>creator_url :</span></dt><dd>https://earth.jpl.nasa.gov/emit/</dd><dt><span>project :</span></dt><dd>Earth Surface Mineral Dust Source Investigation</dd><dt><span>project_url :</span></dt><dd>https://earth.jpl.nasa.gov/emit/</dd><dt><span>publisher_name :</span></dt><dd>NASA LPDAAC</dd><dt><span>publisher_url :</span></dt><dd>https://lpdaac.usgs.gov</dd><dt><span>publisher_email :</span></dt><dd>lpdaac@usgs.gov</dd><dt><span>identifier_product_doi_authority :</span></dt><dd>https://doi.org</dd><dt><span>flight_line :</span></dt><dd>emit20240521t232728_o14216_s000</dd><dt><span>time_coverage_start :</span></dt><dd>2024-05-21T23:27:28+0000</dd><dt><span>time_coverage_end :</span></dt><dd>2024-05-21T23:27:40+0000</dd><dt><span>software_build_version :</span></dt><dd>010618</dd><dt><span>software_delivery_version :</span></dt><dd>010618</dd><dt><span>product_version :</span></dt><dd>V001</dd><dt><span>history :</span></dt><dd>PGE Run Command: {python /beegfs/store/emit/ops/repos/emit-sds-l2a/spectrum_quality.py /tmp/emit/ops/emit20240521t232728_emit.L2AReflectance_20240523t023036/output/emit20240521t232728_rfl /tmp/emit/ops/emit20240521t232728_emit.L2AReflectance_20240523t023036/output/emit20240521t232728_rfl_quality.txt}, PGE Input Files: {radiance_file=/beegfs/store/emit/ops/data/acquisitions/20240521/emit20240521t232728/l1b/emit20240521t232728_o14216_s000_l1b_rdn_b0106_v01.img, pixel_locations_file=/beegfs/store/emit/ops/data/acquisitions/20240521/emit20240521t232728/l1b/emit20240521t232728_o14216_s000_l1b_loc_b0106_v01.img, observation_parameters_file=/beegfs/store/emit/ops/data/acquisitions/20240521/emit20240521t232728/l1b/emit20240521t232728_o14216_s000_l1b_obs_b0106_v01.img, surface_model_config=/beegfs/store/emit/ops/repos/emit-sds-l2a/surface/surface_20221020.json}</dd><dt><span>crosstrack_orientation :</span></dt><dd>as seen on ground</dd><dt><span>easternmost_longitude :</span></dt><dd>-114.84187059769744</dd><dt><span>northernmost_latitude :</span></dt><dd>33.9515655017347</dd><dt><span>westernmost_longitude :</span></dt><dd>-116.103645672334</dd><dt><span>southernmost_latitude :</span></dt><dd>32.925119340889395</dd><dt><span>spatialResolution :</span></dt><dd>0.000542232520256367</dd><dt><span>spatial_ref :</span></dt><dd>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;]],AXIS[&quot;Latitude&quot;,NORTH],AXIS[&quot;Longitude&quot;,EAST],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]]</dd><dt><span>geotransform :</span></dt><dd>[-1.16103646e+02  5.42232520e-04 -0.00000000e+00  3.39515655e+01\n",
       " -0.00000000e+00 -5.42232520e-04]</dd><dt><span>day_night_flag :</span></dt><dd>Day</dd><dt><span>title :</span></dt><dd>EMIT L2A Estimated Surface Reflectance 60 m V001</dd><dt><span>granule_id :</span></dt><dd>EMIT_L2A_RFL_001_20240521T232728_2414216_003</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset>\n",
       "Dimensions:           (latitude: 1893, longitude: 2327, wavelengths: 285)\n",
       "Coordinates:\n",
       "  * wavelengths       (wavelengths) float32 381.0 388.4 ... 2.486e+03 2.493e+03\n",
       "    fwhm              (wavelengths) float32 ...\n",
       "    good_wavelengths  (wavelengths) float32 ...\n",
       "  * latitude          (latitude) float64 33.95 33.95 33.95 ... 32.93 32.93 32.93\n",
       "  * longitude         (longitude) float64 -116.1 -116.1 -116.1 ... -114.8 -114.8\n",
       "    elev              (latitude, longitude) float32 -9.999e+03 ... -9.999e+03\n",
       "    spatial_ref       int64 0\n",
       "Data variables:\n",
       "    reflectance       (latitude, longitude, wavelengths) float32 -9.999e+03 ....\n",
       "Attributes: (12/39)\n",
       "    ncei_template_version:             NCEI_NetCDF_Swath_Template_v2.0\n",
       "    summary:                           The Earth Surface Mineral Dust Source ...\n",
       "    keywords:                          Imaging Spectroscopy, minerals, EMIT, ...\n",
       "    Conventions:                       CF-1.63\n",
       "    sensor:                            EMIT (Earth Surface Mineral Dust Sourc...\n",
       "    instrument:                        EMIT\n",
       "    ...                                ...\n",
       "    spatialResolution:                 0.000542232520256367\n",
       "    spatial_ref:                       GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHER...\n",
       "    geotransform:                      [-1.16103646e+02  5.42232520e-04 -0.00...\n",
       "    day_night_flag:                    Day\n",
       "    title:                             EMIT L2A Estimated Surface Reflectance...\n",
       "    granule_id:                        EMIT_L2A_RFL_001_20240521T232728_24142..."
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ortho_ds = ortho_xr(swath_ds)\n",
    "ortho_ds"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[[0.04310642, 0.04336452, 0.04150592, ..., 0.10212765,\n",
       "         0.0974834 , 0.09012822],\n",
       "        [0.06145171, 0.06501866, 0.06859049, ..., 0.10978349,\n",
       "         0.10137388, 0.09884916],\n",
       "        [0.06124046, 0.0649093 , 0.0685825 , ..., 0.10782409,\n",
       "         0.09501854, 0.0907679 ],\n",
       "        ...,\n",
       "        [0.00507217, 0.00802279, 0.01125222, ..., 0.1837036 ,\n",
       "         0.18001686, 0.17860676],\n",
       "        [0.0060593 , 0.00890038, 0.00998863, ..., 0.19765685,\n",
       "         0.19258572, 0.18802734],\n",
       "        [0.01689515, 0.01768825, 0.01971469, ..., 0.18779285,\n",
       "         0.17850284, 0.17630193]],\n",
       "\n",
       "       [[0.05653568, 0.06005755, 0.06357399, ..., 0.10565821,\n",
       "         0.10205421, 0.09661389],\n",
       "        [0.06207323, 0.06573673, 0.0694032 , ..., 0.10960843,\n",
       "         0.09955885, 0.09995052],\n",
       "        [0.06266887, 0.06635007, 0.07002756, ..., 0.12251341,\n",
       "         0.10631712, 0.10210592],\n",
       "        ...,\n",
       "        [0.0090206 , 0.01041509, 0.01469167, ..., 0.19406983,\n",
       "         0.18753886, 0.18678044],\n",
       "        [0.00804483, 0.00854951, 0.01287901, ..., 0.2105122 ,\n",
       "         0.19998693, 0.19420823],\n",
       "        [0.01344281, 0.02022264, 0.02332068, ..., 0.21165165,\n",
       "         0.20166427, 0.1968817 ]],\n",
       "\n",
       "       [[0.05645982, 0.06005775, 0.06365536, ..., 0.10396454,\n",
       "         0.09737607, 0.09334052],\n",
       "        [0.0636641 , 0.06723313, 0.07079227, ..., 0.11946846,\n",
       "         0.10748741, 0.10533743],\n",
       "        [0.0503701 , 0.05209369, 0.05276025, ..., 0.13338569,\n",
       "         0.11541577, 0.10610302],\n",
       "        ...,\n",
       "        [0.01045959, 0.01191389, 0.01809295, ..., 0.18949942,\n",
       "         0.18770707, 0.1849792 ],\n",
       "        [0.01119158, 0.01632101, 0.01772896, ..., 0.20461799,\n",
       "         0.19886766, 0.19580711],\n",
       "        [0.01879507, 0.02066785, 0.02619297, ..., 0.20391232,\n",
       "         0.19330026, 0.18949336]],\n",
       "\n",
       "       ...,\n",
       "\n",
       "       [[0.04757445, 0.04710032, 0.04930563, ..., 0.15936494,\n",
       "         0.15221749, 0.14613159],\n",
       "        [0.05635054, 0.05588401, 0.05641483, ..., 0.17072305,\n",
       "         0.14819695, 0.14232107],\n",
       "        [0.05913376, 0.06211373, 0.0591522 , ..., 0.18061386,\n",
       "         0.15313587, 0.14717843],\n",
       "        ...,\n",
       "        [0.00169001, 0.01058776, 0.01867289, ..., 0.29611236,\n",
       "         0.2422675 , 0.23527871],\n",
       "        [0.0066209 , 0.01140852, 0.02137864, ..., 0.30386373,\n",
       "         0.25074586, 0.21783502],\n",
       "        [0.01381256, 0.01634106, 0.02840667, ..., 0.32458004,\n",
       "         0.25021455, 0.21957657]],\n",
       "\n",
       "       [[0.04608515, 0.04636445, 0.04815491, ..., 0.15685537,\n",
       "         0.14791432, 0.14348774],\n",
       "        [0.05720487, 0.05567257, 0.0558444 , ..., 0.17613028,\n",
       "         0.14585048, 0.14365467],\n",
       "        [0.06053159, 0.06280502, 0.06347847, ..., 0.19072337,\n",
       "         0.16295457, 0.15611555],\n",
       "        ...,\n",
       "        [0.00047705, 0.00738141, 0.01732354, ..., 0.30685014,\n",
       "         0.25620067, 0.24252747],\n",
       "        [0.00173604, 0.00858373, 0.0177619 , ..., 0.31598717,\n",
       "         0.24594146, 0.21550393],\n",
       "        [0.01271567, 0.01607181, 0.02918276, ..., 0.3471389 ,\n",
       "         0.27922907, 0.24179955]],\n",
       "\n",
       "       [[0.04805192, 0.04812167, 0.04807356, ..., 0.15694118,\n",
       "         0.14717287, 0.14105979],\n",
       "        [0.05594174, 0.05660632, 0.05843432, ..., 0.17907235,\n",
       "         0.15017025, 0.1443422 ],\n",
       "        [0.05993449, 0.0628731 , 0.06273907, ..., 0.1795007 ,\n",
       "         0.15194707, 0.15029328],\n",
       "        ...,\n",
       "        [0.00143897, 0.00733537, 0.01272336, ..., 0.29998982,\n",
       "         0.24480523, 0.22862697],\n",
       "        [0.00134862, 0.00810808, 0.01724624, ..., 0.31884253,\n",
       "         0.2547975 , 0.23689198],\n",
       "        [0.01563575, 0.02371439, 0.03452612, ..., 0.34024465,\n",
       "         0.28424233, 0.24582131]]], dtype=float32)"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "swath_array = swath_ds[\"reflectance\"].data\n",
    "swath_array"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(1280, 1242, 285)"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "swath_array.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "numpy.ndarray"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "type(swath_array)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        ...,\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.]],\n",
       "\n",
       "       [[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        ...,\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.]],\n",
       "\n",
       "       [[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        ...,\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.]],\n",
       "\n",
       "       ...,\n",
       "\n",
       "       [[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        ...,\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.]],\n",
       "\n",
       "       [[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        ...,\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.]],\n",
       "\n",
       "       [[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        ...,\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.]]],\n",
       "      dtype=float32)"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ortho_array = apply_GLT(swath_array, GLT_array)\n",
    "ortho_array"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "numpy.ndarray"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "type(ortho_array)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(1893, 2327, 285)"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ortho_array.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        ...,\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.]],\n",
       "\n",
       "       [[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        ...,\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.]],\n",
       "\n",
       "       [[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        ...,\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.]],\n",
       "\n",
       "       ...,\n",
       "\n",
       "       [[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        ...,\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.]],\n",
       "\n",
       "       [[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        ...,\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.]],\n",
       "\n",
       "       [[-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        ...,\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.],\n",
       "        [-9999., -9999., -9999., ..., -9999., -9999., -9999.]]],\n",
       "      dtype=float32)"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Build Output Dataset\n",
    "if swath_array.ndim == 2:\n",
    "    swath_array = swath_array[:, :, np.newaxis]\n",
    "\n",
    "# ortho_array_shape = (GLT_array.shape[0], GLT_array.shape[1], swath_array.shape[-1])\n",
    "\n",
    "latitude_length = ortho_array.shape[0]\n",
    "longitude_length = ortho_array.shape[1]\n",
    "band_length = swath_array.shape[-1]\n",
    "ortho_array_shape = (latitude_length, longitude_length, band_length)\n",
    "\n",
    "ortho_array = np.full(\n",
    "    ortho_array_shape,\n",
    "    fill_value,\n",
    "    dtype=np.float32,\n",
    ")\n",
    "\n",
    "valid_GLT = np.all(GLT_array != GLT_nodata_value, axis=-1)\n",
    "\n",
    "# Adjust for One based Index - make a copy to prevent decrementing multiple times inside ortho_xr when applying the glt to elev\n",
    "glt_array_copy = GLT_array.copy()\n",
    "glt_array_copy[valid_GLT] -= 1\n",
    "ortho_array[valid_GLT, :] = swath_array[\n",
    "    glt_array_copy[valid_GLT, 1], glt_array_copy[valid_GLT, 0], :\n",
    "]\n",
    "\n",
    "ortho_array"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(1893, 2327, 285)"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ortho_array.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(1893, 2327, 285)"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "latitude_length = GLT_array.shape[0]\n",
    "longitude_length = GLT_array.shape[1]\n",
    "band_length = swath_array.shape[-1]\n",
    "ortho_array_shape = (latitude_length, longitude_length, band_length)\n",
    "ortho_array_shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([-1.16103646e+02,  5.42232520e-04, -0.00000000e+00,  3.39515655e+01,\n",
       "       -0.00000000e+00, -5.42232520e-04])"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "swath_ds.geotransform"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "xarray.core.dataset.Dataset"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "type(swath_ds)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Affine(0.000542232520256367, -0.0, -116.103645672334,\n",
       "       -0.0, -0.000542232520256367, 33.9515655017347)"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "affine = Affine.from_gdal(*swath_ds.geotransform)\n",
    "affine"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{\n",
       "  \"dimensions\": {\n",
       "    \"rows\": 2327,\n",
       "    \"cols\": 1893\n",
       "  },\n",
       "  \"bbox\": {\n",
       "    \"xmin\": -116.103645672334,\n",
       "    \"ymin\": 32.68979042709813,\n",
       "    \"xmax\": -115.0771995114887,\n",
       "    \"ymax\": 33.9515655017347\n",
       "  },\n",
       "  \"crs\": \"EPSG:4326\",\n",
       "  \"resolution\": {\n",
       "    \"cell_width\": 0.000542232520256367,\n",
       "    \"cell_height\": -0.000542232520256367\n",
       "  }\n",
       "}"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "grid = RasterGrid.from_affine(affine, longitude_length, latitude_length)\n",
    "grid"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "emit_tools",
   "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.11.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}