Merge pull request #145 from jlarsen-usgs/main

fix rasterio index() call

Author: mnfienen

notebooks/part0_python_intro/10_Rasterio.ipynb

@@ -156,17 +156,16 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "list(zip(*points))"
+    "xpts, ypts = list(zip(*points))\n",
+    "# points"
    ]
   },
   {
    "cell_type": "markdown",
    "id": "404a972a",
    "metadata": {},
    "source": [
-    "#### Now let's sample elevations at these points\n",
-    "\n",
-    "If we do the unpacking in a function call (e.g. `src.index(*zip(*points))`), we need to add a leading `*` to unpack the tuple into two separate arguments"
+    "#### Now let's sample elevations at these points\n"
    ]
   },
   {
@@ -176,10 +175,14 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "with rasterio.open(input_rasters[2015]) as src:       \n",
-    "        meta = src.meta\n",
-    "        data = src.read(1)\n",
-    "        i, j = src.index(*zip(*points))\n",
+    "with rasterio.open(input_rasters[2015]) as src:  \n",
+    "    meta = src.meta\n",
+    "    data = src.read(1)\n",
+    "    i, j = [], []\n",
+    "    for ix in range(len(xpts)):\n",
+    "        ti, tj = src.index(xpts[ix], ypts[ix])\n",
+    "        i.append(ti)\n",
+    "        j.append(tj)\n",
     "print(data[i, j])"
    ]
   },
@@ -993,7 +996,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.11.7"
+   "version": "3.11.10"
   }
  },
  "nbformat": 4,