Updated rasterio syntax and other small syntax and comment changes.

Author: kevinlacaille

jupyter-notebooks/coastal-erosion-example/1_rasterio_firstlook.ipynb

@@ -119,6 +119,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
+    "# Open our image\n",
     "satdat = rasterio.open(image_file)\n",
     "\n",
     "# satdat is our open dataset object\n",
@@ -170,8 +171,9 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "# PlanetScope 4-band band order: BGRN\n",
+    "# PlanetScope 3-band band order: BGR\n",
     "\n",
+    "# Extract the rasterbands\n",
     "with rasterio.open(image_file) as src:\n",
     "    blue = src.read(1)\n",
     "    green = src.read(2)\n",

jupyter-notebooks/coastal-erosion-example/2_rasterbands.ipynb

@@ -110,6 +110,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
+    "# Open our image\n",
     "satdat = rasterio.open(image_file)\n",
     "\n",
     "# satdat is our open dataset object\n",
@@ -130,8 +131,7 @@
    "outputs": [],
    "source": [
     "# Get the image's coordinate reference system (WGS84 / UTM 15N)\n",
-    "\n",
-    "satdat.crs"
+    "print(satdat.crs)"
    ]
   },
   {
@@ -148,7 +148,6 @@
    "outputs": [],
    "source": [
     "# Minimum bounding box in projected units (meters), in our zone\n",
-    "\n",
     "print(satdat.bounds)"
    ]
   },
@@ -159,7 +158,6 @@
    "outputs": [],
    "source": [
     "# Get dimensions, in projected units (using the example GeoTIFF, that's meters)\n",
-    "\n",
     "width_in_projected_units = satdat.bounds.right - satdat.bounds.left\n",
     "height_in_projected_units = satdat.bounds.top - satdat.bounds.bottom\n",
     "\n",
@@ -173,7 +171,6 @@
    "outputs": [],
    "source": [
     "# Number of rows and columns (pixels)\n",
-    "\n",
     "print(\"Rows: {}, Columns: {}\".format(satdat.height, satdat.width))"
    ]
   },
@@ -183,7 +180,8 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "# This dataset's projection uses meters as projected units.  What are the dimensions of a single pixel in meters?\n",
+    "# This dataset's projection uses meters as projected units.  \n",
+    "# What are the dimensions of a single pixel in meters?\n",
     "\n",
     "xres = width_in_projected_units / satdat.width\n",
     "yres = height_in_projected_units / satdat.height\n",
@@ -223,8 +221,8 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "# All of the metadata required to create an image of the same dimensions, datatype, format, etc. is stored in\n",
-    "# the dataset's profile:\n",
+    "# All of the metadata required to create an image of the same dimensions, \n",
+    "# datatype, format, etc. is stored in the dataset's profile:\n",
     "\n",
     "satdat.profile"
    ]
@@ -263,7 +261,6 @@
    "outputs": [],
    "source": [
     "# PlanetScope 3-band band order: BGR\n",
-    "\n",
     "with rasterio.open(image_file) as src:\n",
     "    blue = src.read(1)\n",
     "    green = src.read(2)\n",
@@ -277,7 +274,6 @@
    "outputs": [],
    "source": [
     "# To construct a visual image, we will need the red, green, and blue bands\n",
-    "\n",
     "visual_image = np.dstack((blue, green, red))"
    ]
   },
@@ -295,7 +291,8 @@
    "outputs": [],
    "source": [
     "plt.imshow(visual_image)\n",
-    "plt.title(\"Visual Image\")"
+    "plt.title(\"Visual Image\")\n",
+    "plt.show()"
    ]
   },
   {

jupyter-notebooks/coastal-erosion-example/3_Compute_NDWI.ipynb

@@ -202,7 +202,8 @@
     "fig = plt.imshow(ndwi)\n",
     "fig.set_cmap('Blues')\n",
     "plt.title(\"NDWI\")\n",
-    "plt.colorbar()"
+    "plt.colorbar()\n",
+    "plt.show()"
    ]
   },
   {

jupyter-notebooks/coastal-erosion-example/4_masks_and_filters.ipynb

@@ -115,7 +115,12 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "blue, green, red, nir = satdat.read()"
+    "# Extract rasterbands\n",
+    "with rasterio.open(image_file) as src:\n",
+    "    blue = src.read(1)\n",
+    "    green = src.read(2)\n",
+    "    red = src.read(3)\n",
+    "    nir = src.read(4)"
    ]
   },
   {
@@ -134,6 +139,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
+    "# Generate the water and land masks\n",
     "water_mask = np.full(ndwi.shape, np.nan)\n",
     "land_mask = np.full(ndwi.shape, np.nan)\n",
     "\n",
@@ -167,7 +173,8 @@
     "\n",
     "plt.imshow(water_mask)\n",
     "plt.title(\"Houston water mask\")\n",
-    "plt.colorbar()"
+    "plt.colorbar()\n",
+    "plt.show()"
    ]
   },
   {
@@ -180,7 +187,8 @@
     "\n",
     "plt.imshow(land_mask)\n",
     "plt.title(\"Houston land mask\")\n",
-    "plt.colorbar()"
+    "plt.colorbar()\n",
+    "plt.show()"
    ]
   },
   {
@@ -232,8 +240,7 @@
     "mask_closed_opened = cv2.morphologyEx(mask_closed, cv2.MORPH_OPEN, opening_kernel)\n",
     "\n",
     "# Ensure the clipped areas remain clipped\n",
-    "mask_closed_opened[mask_closed_opened == 0] = np.nan\n",
-    "\n"
+    "mask_closed_opened[mask_closed_opened == 0] = np.nan"
    ]
   },
   {
@@ -243,7 +250,8 @@
    "outputs": [],
    "source": [
     "plt.imshow(mask_closed_opened)\n",
-    "plt.title(\"Houston filtered water mask\")"
+    "plt.title(\"Houston filtered water mask\")\n",
+    "plt.show()"
    ]
   },
   {

jupyter-notebooks/coastal-erosion-example/5_plotting_a_histogram.ipynb

@@ -96,8 +96,13 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "# Load the 4 bands into 2d arrays - recall that we previously learned PlanetScope band order is BGRN.\n",
-    "blue, green, red, nir = satdat.read()"
+    "# Load the 4 bands into 2d arrays\n",
+    "# recall that we previously learned 4-band PlanetScope band order is BGRN.\n",
+    "with rasterio.open(image_file) as src:\n",
+    "    blue = src.read(1)\n",
+    "    green = src.read(2)\n",
+    "    red = src.read(3)\n",
+    "    nir = src.read(4)"
    ]
   },
   {
@@ -108,10 +113,10 @@
    "outputs": [],
    "source": [
     "# Define a new figure\n",
-    "fig2 = plt.figure(figsize=(20,10))\n",
+    "fig = plt.figure(figsize=(20,10))\n",
     "\n",
     "# Give this new figure a subplot, which will contain the histogram itself\n",
-    "ax = fig2.add_subplot(111)\n",
+    "ax = fig.add_subplot(111)\n",
     "\n",
     "# Add a title & (x,y) labels to the plot\n",
     "plt.title(\"Histogram Example\", fontsize=18, fontweight='bold')\n",
@@ -133,7 +138,7 @@
     "ax.hist(x,bins,color=color)\n",
     "\n",
     "# Save the generated figure to an external image file\n",
-    "fig2.savefig(\"histogram.png\", dpi=200, bbox_inches='tight', pad_inches=0.7)\n",
+    "fig.savefig(\"histogram.png\", dpi=200, bbox_inches='tight', pad_inches=0.7)\n",
     "\n",
     "# Finally - let's take a look!\n",
     "plt.show()"