Flesh out demo w/ JGIS

Author: mfisher87

modules/06-geojupyter/demo.ipynb

@@ -49,7 +49,36 @@
     "\n",
     "DATA_DIR = Path().cwd() / \"data\"\n",
     "INEQUALITY_GEOJSON_FILE = DATA_DIR / \"redlining_newhaven_ct.geojson\"\n",
-    "NDVI_FILE = DATA_DIR / \"ndvi.tif\""
+    "NDVI_GEOTIFF_FILE = DATA_DIR / \"ndvi.tif\"\n",
+    "INEQUALITY_NDVI_GEOJSON_FILE = DATA_DIR / \"redlining_ndvi_newhaven_ct.geojson\""
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "caaed0ec-07be-4675-bf2f-c626d500d443",
+   "metadata": {},
+   "source": [
+    "### Set up a JupyterGIS project\n",
+    "\n",
+    "We can build a JupyterGIS project from scratch in Python. Let's start by adding an OpenStreetMap basemap and displaying the widget in a side panel."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "4eb60a7f-dd77-4043-8812-32cac1bdd800",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "from jupytergis import GISDocument\n",
+    "\n",
+    "jgis_project = GISDocument()\n",
+    "jgis_project.add_raster_layer(\n",
+    "    url=\"https://tile.openstreetmap.org/{z}/{x}/{y}.png\",\n",
+    "    name=\"Basemap\",\n",
+    ")\n",
+    "\n",
+    "jgis_project.sidecar()"
    ]
   },
   {
@@ -82,6 +111,8 @@
     ")\n",
     "\n",
     "new_haven_redlining =  redlines.execute().set_crs(\"EPSG:4326\")\n",
+    "# TODO: Set a numeric grade for graduated symbology\n",
+    "\n",
     "new_haven_redlining.to_file(INEQUALITY_GEOJSON_FILE, engine=\"fiona\")\n",
     "\n",
     "new_haven_bbox = new_haven_redlining.total_bounds"
@@ -92,7 +123,9 @@
    "id": "63f92fa9-49cd-4e9e-a19d-7d6a60747b65",
    "metadata": {},
    "source": [
-    "Let's explore the data a little bit. Hover over the polygons after running the cell below! Does anything jump out at you?"
+    "#### Explore the data\n",
+    "\n",
+    "Let's explore the data a little bit. After running the cell below, **right-click the \"New Haven neighborhood redlining\" layer** in the JupyterGIS interface, and **select \"Zoom to layer\"**. "
    ]
   },
   {
@@ -102,15 +135,30 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "new_haven_redlining.explore(column=\"grade\", cmap=\"inferno\")"
+    "jgis_project.add_geojson_layer(\n",
+    "    path=INEQUALITY_GEOJSON_FILE,\n",
+    "    name=\"New Haven neighborhood redlining\",\n",
+    ");"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "7663ab54-2d69-479b-b8ce-cbdc0bc13143",
+   "metadata": {},
+   "source": [
+    "Now, with the \"New Haven neighborhood redlining\" layer selected, **click the `i` (identify) icon in the toolbar** at the top of the JupyterGIS interface.\n",
+    "\n",
+    "Select some neighborhoods and view their \"Grade\" and \"Category\" attributes.\n",
+    "\n",
+    "TODO: Symbologize on numeric grade"
    ]
   },
   {
    "cell_type": "markdown",
    "id": "61803898-9832-44a9-a04f-23bdb1a4ac3c",
    "metadata": {},
    "source": [
-    "### Calculating NDVI\n",
+    "### Calculate NDVI\n",
     "\n",
     "We're going to calculate NDVI from Sentinel-2 data."
    ]
@@ -202,11 +250,43 @@
     "ndvi.rio.reproject(\n",
     "    \"EPSG:4326\",\n",
     ").rio.to_raster(\n",
-    "    raster_path=NDVI_FILE, \n",
+    "    raster_path=NDVI_GEOTIFF_FILE, \n",
     "    driver=\"COG\",\n",
     ")"
    ]
   },
+  {
+   "cell_type": "markdown",
+   "id": "14a184a5-00ab-4347-b576-bc384b715c92",
+   "metadata": {},
+   "source": [
+    "#### Explore the data\n",
+    "\n",
+    "Let's explore the data in JupyterGIS again.\n",
+    "This time, we'll add the layer with the GUI.\n",
+    "\n",
+    "If the \"identify\" tool is still active, click the `i` icon in the toolbar again to disable it.\n",
+    "\n",
+    "Now, **click the `+` icon in the toolbar** to open the new layer interface.\n",
+    "**Select \"Add Raster Layer\", then \"New GeoTiff Layer\"**.\n",
+    "\n",
+    "**Select \"Browse Server Files\"** and then **navigate to `ndvi.tif` in the module 6 data directory** (`workshop-open-source-geospatial/modules/06-geojupyter/data/ndvi.tif`).\n",
+    "\n",
+    "Click **Select**.\n",
+    "\n",
+    "**Set the \"Min\" field to `0` and \"Max\" to `1`**.\n",
+    "\n",
+    "**Uncheck \"Normalize\"**.\n",
+    "\n",
+    "Scroll down to **input the layer name as \"NDVI\"**.\n",
+    "\n",
+    "Finally, **right-click the \"NDVI\" layer** and **select \"Edit Symbology\"**. The symbology menu may take a moment to load. Be patient! **Select \"Classify\" then click \"OK\".**\n",
+    "\n",
+    "The brighter areas have a higher NDVI value, and the darker areas have a lower one.\n",
+    "\n",
+    "We can use the identify tool (`i` icon in the toolbar) to explore the raw values."
+   ]
+  },
   {
    "cell_type": "markdown",
    "id": "41a3f550-04da-4170-9e69-b4131d428fdd",
@@ -227,15 +307,30 @@
     "from exactextract import exact_extract\n",
     "\n",
     "new_haven_redlining_and_ndvi = exact_extract(\n",
-    "    NDVI_FILE,\n",
+    "    NDVI_GEOTIFF_FILE,\n",
     "    new_haven_redlining,\n",
     "    \"mean_ndvi=mean\",\n",
     "    include_geom = True,\n",
     "    include_cols=[\"label\", \"grade\", \"city\", \"fill\"],\n",
     "    output=\"pandas\",\n",
     ")\n",
     "\n",
-    "new_haven_redlining_and_ndvi.explore(column=\"mean_ndvi\")"
+    "new_haven_redlining_and_ndvi.set_crs(\n",
+    "    \"EPSG:4326\"\n",
+    ").to_file(INEQUALITY_NDVI_GEOJSON_FILE, engine=\"fiona\")"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "8244a541-59d2-4372-a8b6-9f7f8b6483de",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# jgis_project.add_geojson_layer(\n",
+    "#     path=INEQUALITY_NDVI_GEOJSON_FILE,\n",
+    "#     name=\"New Haven neighborhood redlining w/ NDVI\",\n",
+    "# )  # TODO"
    ]
   },
   {