02_EMIT_Animations.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "import holoviews as hv\n",
    "import hvplot.xarray  # noqa: F401\n",
    "from fsspec.implementations.http import HTTPFileSystem\n",
    "\n",
    "from emit_tools import emit_xarray\n",
    "from utils import get_earthdata_token, load_emit_granule"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Band animations\n",
    "\n",
    "This notebook will create an interactive visualisation explorer or animation that\n",
    "steps through the bands.\n",
    "\n",
    "Since the hyperspectral sensor splits the spectrum into 285 bands, there are a lot\n",
    "of combinations that could be used if we stick to RGB or single-band visualisation\n",
    "as demonstrated in the first notebook. This notebook will allow us to quickly\n",
    "explore the scene across the spectra. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "%%time\n",
    "# Loading data can take around 3-4 minutes on a 100 Mbps connection\n",
    "\n",
    "# See README.md for instructions on how to get an Earthdata token\n",
    "token = get_earthdata_token()\n",
    "\n",
    "# Refer to the README.md for instructions on how to find different granule IDs\n",
    "granule = \"EMIT_L2A_RFL_001_20230316T045211_2307503_006\" # Canberra\n",
    "\n",
    "ds = load_emit_granule(granule, token)\n",
    "ds"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "# Render an image for each band and build an interactive explorer\n",
    "# This takes a few minutes to render if you ask for many images\n",
    "\n",
    "# All wavelengths... generally doesn't work because there's too many!\n",
    "# wavelengths = ds.wavelengths.values\n",
    "        \n",
    "# A subset of wavelengths. Change the third value to change the step size.\n",
    "wavelengths = ds.wavelengths.values[range(0, len(ds.wavelengths.values), 20)]\n",
    "\n",
    "animation = hv.HoloMap(\n",
    "    {\n",
    "        band: hv.Image(ds.sel(bands=band).reflectance).opts(\n",
    "            cmap=\"bokeh\", colorbar=True, frame_width=600, aspect=\"equal\", xaxis='bare', yaxis='bare', invert_axes=True\n",
    "        )\n",
    "        for band in wavelengths\n",
    "    },\n",
    "    kdims=\"band\",\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "# View the interactive explorer, it gets slow with lots of images\n",
    "animation\n",
    "\n",
    "# These need selenium installed...\n",
    "# # Inline animated output\n",
    "# hv.output(animation, holomap=\"gif\", fps=4)\n",
    "\n",
    "# # Save animation as a file\n",
    "# hv.save(animation, 'spectral_animation.gif', fps=4, toolbar=False)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}