added reprojection use case

Author: grallewellyn

docs/python-usage.ipynb

@@ -65,6 +65,7 @@
     "from shapely.geometry import box\n",
     "import urllib.request\n",
     "from rasterio.mask import mask\n",
+    "from rasterio.warp import calculate_default_transform, reproject, Resampling\n",
     "import fiona\n",
     "%matplotlib inline"
    ]
@@ -700,6 +701,73 @@
     "print('')"
    ]
   },
+  {
+   "cell_type": "markdown",
+   "id": "4b6ca9ca",
+   "metadata": {},
+   "source": [
+    "### Visualize WMS Raster Data In Any Projection\n",
+    "\n",
+    "<DIV align=\"left\" style=\"line-height:1.5em;\">\n",
+    "<p>\n",
+    "    \n",
+    "This example shows how to reproject WMS Raster Data to the projection of your choice using rasterio\n",
+    "    \n",
+    "</p>\n",
+    "</DIV>"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "3641a183",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Here you can set the new projection to whatever you like \n",
+    "dst_crs = 'EPSG:2261'\n",
+    "\n",
+    "# Save a global extents tiff file\n",
+    "# Connect to GIBS WMS Service\n",
+    "wms = WebMapService('https://gibs.earthdata.nasa.gov/wms/epsg4326/best/wms.cgi?', version='1.1.1')\n",
+    "\n",
+    "# Configure request for MODIS_Terra_SurfaceReflectance_Bands143\n",
+    "img = wms.getmap(layers=['MODIS_Terra_SurfaceReflectance_Bands143'],  # Layers\n",
+    "                 srs='epsg:4326',  # Map projection\n",
+    "                 bbox=(-180,-90,180,90),  # Bounds\n",
+    "                 size=(1200, 600),  # Image size\n",
+    "                 time='2021-11-25',  # Time of data\n",
+    "                 format='image/tiff',  # Image format\n",
+    "                 transparent=True)  # Nodata transparency\n",
+    "\n",
+    "# Save output TIFF to a file\n",
+    "out = open('global_extents.tiff', 'wb')\n",
+    "out.write(img.read())\n",
+    "out.close()\n",
+    "\n",
+    "with rasterio.open('global_extents.tiff') as src:\n",
+    "    transform, width, height = calculate_default_transform(\n",
+    "        src.crs, dst_crs, src.width, src.height, *src.bounds)\n",
+    "    kwargs = src.meta.copy()\n",
+    "    kwargs.update({\n",
+    "        'crs': dst_crs,\n",
+    "        'transform': transform,\n",
+    "        'width': width,\n",
+    "        'height': height\n",
+    "    })\n",
+    "\n",
+    "    with rasterio.open(os.getcwd()+'/reprojectedImage.byte.tif', 'w', **kwargs) as dst:\n",
+    "        for i in range(1, src.count + 1):\n",
+    "            reproject(\n",
+    "                source=rasterio.band(src, i),\n",
+    "                destination=rasterio.band(dst, i),\n",
+    "                src_transform=src.transform,\n",
+    "                src_crs=src.crs,\n",
+    "                dst_transform=transform,\n",
+    "                dst_crs=dst_crs,\n",
+    "                resampling=Resampling.nearest)"
+   ]
+  },
   {
    "attachments": {},
    "cell_type": "markdown",