examples/swot_cmr_sim.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0b7e6c08-fba9-405d-a278-767ec692efb4",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "# Imports\n",
    "from sliderule import sliderule, earthdata, swot\n",
    "import matplotlib.pyplot as plt\n",
    "import geopandas as gpd\n",
    "import numpy"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "f4155c77-d1bb-42c2-b47a-1089823a8c89",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "# Query EarthData for SWOT data at Area of Interest\n",
    "region = sliderule.toregion(\"../data/grandmesa.geojson\")\n",
    "granules = earthdata.cmr(short_name=\"SWOT_SIMULATED_L2_KARIN_SSH_ECCO_LLC4320_CALVAL_V1\", polygon=region[\"poly\"], time_start=None)\n",
    "granules"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "6f27772e-630a-4e50-bcd3-aa7e1d6acff4",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "# Boundary information from EarthData for SWOT_L2_LR_SSH_Expert_238_002_20120705T114746_20120705T123852_DG10_01.nc \n",
    "track = {\n",
    "    \"points\": [\n",
    "        {\n",
    "          \"longitude\": -180,\n",
    "          \"latitude\": 77.32443\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -180,\n",
    "          \"latitude\": 75.96161\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -162.09458,\n",
    "          \"latitude\": 72.95769\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -148.71771,\n",
    "          \"latitude\": 68.02475\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -139.57473,\n",
    "          \"latitude\": 61.46376\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -132.76541,\n",
    "          \"latitude\": 52.52883\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -128.23535,\n",
    "          \"latitude\": 42.33934\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -124.53764,\n",
    "          \"latitude\": 28.97745\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -116.14332,\n",
    "          \"latitude\": -22.21489\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -112.54563,\n",
    "          \"latitude\": -38.63235\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -107.59535,\n",
    "          \"latitude\": -52.02395\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -101.02934,\n",
    "          \"latitude\": -61.71245\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -96.51192,\n",
    "          \"latitude\": -65.80947\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -91.33566,\n",
    "          \"latitude\": -69.12003\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -77.82439,\n",
    "          \"latitude\": -74.1454\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -59.19556,\n",
    "          \"latitude\": -77.20514\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -35.45207,\n",
    "          \"latitude\": -78.29215\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -35.45501,\n",
    "          \"latitude\": -77.0331\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -59.11783,\n",
    "          \"latitude\": -75.83566\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -77.9601,\n",
    "          \"latitude\": -72.38965\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -90.38687,\n",
    "          \"latitude\": -67.38657\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -99.59434,\n",
    "          \"latitude\": -60.19465\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -106.37693,\n",
    "          \"latitude\": -50.25634\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -111.51882,\n",
    "          \"latitude\": -36.48096\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -115.3008,\n",
    "          \"latitude\": -19.27607\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -123.63896,\n",
    "          \"latitude\": 31.6175\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -127.38371,\n",
    "          \"latitude\": 45.02024\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -131.90094,\n",
    "          \"latitude\": 54.92655\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -139.12559,\n",
    "          \"latitude\": 63.9081\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -148.94375,\n",
    "          \"latitude\": 70.29256\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -162.17736,\n",
    "          \"latitude\": 74.64083\n",
    "        },\n",
    "        {\n",
    "          \"longitude\": -180,\n",
    "          \"latitude\": 77.32443\n",
    "        }\n",
    "    ]\n",
    "}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "6ff89d99-3045-48db-b084-8757521e2462",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Build GeoDataFrame of Boundary\n",
    "lat = [p['latitude'] for p in track['points']]\n",
    "lon = [p['longitude'] for p in track['points']]\n",
    "data = [x for x in range(len(lon))]\n",
    "a = {\"latitude\": lat, \"longitude\": lon, \"time\": numpy.array(data)}\n",
    "gdf = sliderule.todataframe(a)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "8abcc3a0-2a86-4c3f-9438-7e9befa320c4",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "# Plot Boundary and Area of Interest on World Map\n",
    "worldmap = gpd.read_file(gpd.datasets.get_path(\"naturalearth_lowres\"))\n",
    "fig, ax = plt.subplots(figsize=(12, 6))\n",
    "worldmap.plot(color=\"lightgrey\", ax=ax)\n",
    "gdf.plot(ax=ax)\n",
    "region['gdf'].buffer(1.0).plot(ax=ax, color='green')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3ba51afd-98fe-4ac3-a885-d9a7975bc009",
   "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.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}