added the example mapping with cartopy

Author: grallewellyn

docs/python-usage.ipynb

@@ -36,10 +36,24 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": null,
    "id": "5a65137f",
    "metadata": {},
-   "outputs": [],
+   "outputs": [
+    {
+     "ename": "",
+     "evalue": "",
+     "output_type": "error",
+     "traceback": [
+      "\u001b[1;31mJupyter cannot be started. Error attempting to locate Jupyter: Running cells with 'Python 3.13.0 64-bit' requires jupyter and notebook package.\n",
+      "\u001b[1;31mRun the following command to install 'jupyter and notebook' into the Python environment. \n",
+      "\u001b[1;31mCommand: 'python -m pip install jupyter notebook -U\n",
+      "\u001b[1;31mor\n",
+      "\u001b[1;31mconda install jupyter notebook -U'\n",
+      "\u001b[1;31mClick <a href='https://aka.ms/installJupyterForVSCode'>here</a> for more info."
+     ]
+    }
+   ],
    "source": [
     "import os\n",
     "from io import BytesIO\n",
@@ -57,7 +71,9 @@
     "import xml.etree.ElementTree as xmlet\n",
     "import lxml.etree as xmltree\n",
     "from PIL import Image as plimg\n",
+    "from PIL import ImageDraw\n",
     "import numpy as np\n",
+    "import pandas as pd\n",
     "from owslib.wms import WebMapService\n",
     "from IPython.display import Image, display\n",
     "import rasterio\n",
@@ -67,6 +83,7 @@
     "from rasterio.mask import mask\n",
     "from rasterio.warp import calculate_default_transform, reproject, Resampling\n",
     "import fiona\n",
+    "from datetime import datetime, timedelta\n",
     "%matplotlib inline"
    ]
   },
@@ -989,6 +1006,68 @@
     "m"
    ]
   },
+  {
+   "cell_type": "markdown",
+   "id": "8b567f79",
+   "metadata": {},
+   "source": [
+    "### Animated Web Map with WMS\n",
+    "\n",
+    "<DIV align=\"left\" style=\"line-height:1.5em;\">\n",
+    "<p>\n",
+    "\n",
+    "The next example shows how to display IMERG_Precipitation_Rate layer in an animated web map (may require additional Python libraries). \n",
+    "\n",
+    "</p>\n",
+    "</DIV>"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "2b7d7657",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "wms = WebMapService('https://gibs.earthdata.nasa.gov/wms/epsg4326/best/wms.cgi?', version='1.3.0')\n",
+    "\n",
+    "layers = ['MODIS_Aqua_CorrectedReflectance_TrueColor',\n",
+    "          'IMERG_Precipitation_Rate',\n",
+    "          'Reference_Features',\n",
+    "          'Reference_Labels']\n",
+    "color = 'rgb(255,255,255)'\n",
+    "frames = []\n",
+    "\n",
+    "start_date = datetime(2022, 9, 25)\n",
+    "end_date = datetime(2022, 10, 1)\n",
+    "dates = pd.date_range(start_date,end_date-timedelta(days=1),freq='d')\n",
+    "\n",
+    "for day in dates:\n",
+    "    datatime = day.strftime(\"%Y-%m-%d\")\n",
+    "    img = wms.getmap(layers=layers,  # Layers\n",
+    "                     srs='epsg:4326',  # Map projection\n",
+    "                     bbox=(-87, 18, -72, 35),  # Bounds\n",
+    "                     size=(600,600),  # Image size\n",
+    "                     time=datatime,  # Time of data\n",
+    "                     format='image/png',  # Image format\n",
+    "                     transparent=True)  # Nodata transparency\n",
+    "\n",
+    "    image = plimg.open(img)\n",
+    "    draw = ImageDraw.Draw(image)\n",
+    "    (x, y) = (50, 20)\n",
+    "    draw.text((x, y), f'IMERG Precipitation Rate   -   {datatime}', fill=color)\n",
+    "    frames.append(image)\n",
+    "    \n",
+    "frames[0].save('IMERG_Precipitation_Rate_Ian.gif',\n",
+    "               format='GIF',\n",
+    "               append_images=frames,\n",
+    "               save_all=True,\n",
+    "               duration=1000,\n",
+    "               loop=0)\n",
+    "\n",
+    "Image('IMERG_Precipitation_Rate_Ian.gif')"
+   ]
+  },
   {
    "attachments": {},
    "cell_type": "markdown",
@@ -1912,6 +1991,47 @@
     "print('')"
    ]
   },
+  {
+   "cell_type": "markdown",
+   "id": "d5e46225",
+   "metadata": {},
+   "source": [
+    "<a id=\"visualize_wmts_raster_data\"></a>\n",
+    "\n",
+    "### Visualize WMTS Raster Data By Cartopy\n",
+    "\n",
+    "<DIV align=\"left\" style=\"line-height:1.5em;\">\n",
+    "<p>\n",
+    "    \n",
+    "This example shows how to display WMTS raster data on a map using Cartopy.\n",
+    "    \n",
+    "</p>\n",
+    "</DIV>"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "50131451",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Define the WMTS URL\n",
+    "wmts_url = \"https://gibs.earthdata.nasa.gov/wmts/epsg4326/best/wmts.cgi\"\n",
+    "\n",
+    "# Create a map with PlateCarree projection\n",
+    "ax = plt.axes(projection=ccrs.PlateCarree())\n",
+    "ax.set_extent([-180, 180, -90, 90])  # World extent\n",
+    "\n",
+    "# Add WMTS layer with a specific date\n",
+    "layer = \"VIIRS_SNPP_SurfaceReflectance_BandsM11-M7-M5\"\n",
+    "time = \"2023-01-01\"\n",
+    "ax.add_wmts(wmts_url, layer_name=layer, wmts_kwargs={\"time\": time})\n",
+    "\n",
+    "plt.title('Land Surface Reflectance (Bands M11-M7-M5, Best Available, VIIRS, Suomi NPP)')\n",
+    "plt.show()"
+   ]
+  },
   {
    "attachments": {},
    "cell_type": "markdown",
@@ -3694,11 +3814,9 @@
  ],
  "metadata": {
   "celltoolbar": "Raw Cell Format",
-  "interpreter": {
-   "hash": "aee8b7b246df8f9039afb4144a1f6fd8d2ca17a180786b69acc140d282b71a49"
-  },
   "kernelspec": {
-   "display_name": "Python 3.10.4 64-bit",
+   "display_name": "Python 3.13.0 ('gibsPythonNotebook2')",
+   "language": "python",
    "name": "python3"
   },
   "language_info": {
@@ -3711,7 +3829,12 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.10.4 (v3.10.4:9d38120e33, Mar 23 2022, 17:29:05) [Clang 13.0.0 (clang-1300.0.29.30)]"
+   "version": "3.13.0"
+  },
+  "vscode": {
+   "interpreter": {
+    "hash": "cd6b7991ad4bc00fdd3e8f0f0395227ccb9271f102bce6443e69587a933de936"
+   }
   }
  },
  "nbformat": 4,