examples/gedi_l4a.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0e7a3594-2353-476d-a42a-b6bf1b279c61",
   "metadata": {},
   "outputs": [],
   "source": [
    "import time\n",
    "import matplotlib.pyplot as plt\n",
    "from sliderule import gedi\n",
    "import sliderule"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "54acbd2a-23d8-4cb9-ac06-b3893a3b45db",
   "metadata": {},
   "outputs": [],
   "source": [
    "# initialize client (notebook only processes one granule, so one node is sufficient)\n",
    "gedi.init(\"slideruleearth.io\", verbose=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e1421300-b0d1-4139-928e-f8bc612a8ac9",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Specify region of interest from geojson\n",
    "poly_fn = 'grandmesa.geojson'\n",
    "region = sliderule.toregion(poly_fn)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "b2e2f316-35e2-4345-bd05-8327e031ba0e",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Build GEDI L4A Request Parameters\n",
    "parms = {\n",
    "    \"poly\": region[\"poly\"],\n",
    "    \"degrade_flag\": 0,\n",
    "    \"l2_quality_flag\": 1,\n",
    "    \"beam\": 0\n",
    "}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "02682045-e626-4864-b4eb-1ed04d77d778",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Latch Start Time\n",
    "perf_start = time.perf_counter()\n",
    "\n",
    "# Request GEDI Data\n",
    "gedi04a = gedi.gedi04ap(parms)\n",
    "                        \n",
    "# Latch Stop Time\n",
    "perf_stop = time.perf_counter()\n",
    "\n",
    "# Display Statistics\n",
    "perf_duration = perf_stop - perf_start\n",
    "print(\"Completed in {:.3f} seconds of wall-clock time\".format(perf_duration))\n",
    "print(\"Received {} footprints\".format(gedi04a.shape[0]))\n",
    "if len(gedi04a) > 0:\n",
    "    print(\"Beams: {}\".format(gedi04a[\"beam\"].unique()))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "92d482e1-08e2-49a4-9ff7-155b7952e178",
   "metadata": {},
   "outputs": [],
   "source": [
    "gedi04a"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "eeade03e-0dc7-42fe-a879-8e59cce8e6af",
   "metadata": {},
   "outputs": [],
   "source": [
    "# plot elevations and vegetation density\n",
    "f, ax = plt.subplots(1, 2, figsize=[12,8])\n",
    "ax[0].set_title(\"Elevations\")\n",
    "ax[0].set_aspect('equal')\n",
    "gedi04a.plot(ax=ax[0], column='elevation', cmap='inferno', s=0.1)\n",
    "ax[1].set_title(\"Vegetation Density\")\n",
    "ax[1].set_aspect('equal')\n",
    "gedi04a.plot(ax=ax[1], column='agbd', cmap='inferno', s=0.1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "b7243ae5-f3bb-4219-8486-bde773f4effa",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}