merge a few cells in 01

Author: simonusher

01_introduction_to_geospatial.ipynb

@@ -433,7 +433,7 @@
    "metadata": {
     "scrolled": false,
     "slideshow": {
-     "slide_type": "fragment"
+     "slide_type": "subslide"
     }
    },
    "outputs": [],
@@ -449,64 +449,23 @@
     }
    },
    "source": [
-    "## Using GeoDataFrames in SRAI"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {
-    "slideshow": {
-     "slide_type": "-"
-    }
-   },
-   "outputs": [],
-   "source": [
-    "poland_gdf = countries[countries[\"NAME\"] == \"Poland\"]\n",
-    "poland_gdf"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {
-    "slideshow": {
-     "slide_type": "subslide"
-    }
-   },
-   "outputs": [],
-   "source": [
-    "from srai.regionalizers import AdministrativeBoundaryRegionalizer\n",
-    "\n",
-    "regionalizer = AdministrativeBoundaryRegionalizer(admin_level=4)\n",
-    "regions_gdf = regionalizer.transform(poland_gdf)\n",
-    "regions_gdf.head(5)"
+    "## Let's go deeper - [Shapely](https://shapely.readthedocs.io/en/stable/manual.html) objects\n",
+    "- GeoPandas uses Shapely - geometry column\n",
+    "- geometric operations"
    ]
   },
   {
    "cell_type": "code",
    "execution_count": null,
    "metadata": {
     "slideshow": {
-     "slide_type": "subslide"
+     "slide_type": "fragment"
     }
    },
    "outputs": [],
    "source": [
-    "regions_gdf.explore()"
-   ]
-  },
-  {
-   "cell_type": "markdown",
-   "metadata": {
-    "slideshow": {
-     "slide_type": "slide"
-    }
-   },
-   "source": [
-    "## Let's go deeper - [Shapely](https://shapely.readthedocs.io/en/stable/manual.html) objects\n",
-    "- GeoPandas uses Shapely - geometry column\n",
-    "- geometric operations"
+    "pl_de_gdf = countries[countries[\"NAME\"].isin([\"Poland\", \"Germany\"])]\n",
+    "pl_de_gdf"
    ]
   },
   {
@@ -519,7 +478,8 @@
    },
    "outputs": [],
    "source": [
-    "type(regions_gdf.iloc[2].geometry)"
+    "merged_geom = pl_de_gdf.unary_union\n",
+    "merged_geom"
    ]
   },
   {
@@ -532,10 +492,7 @@
    },
    "outputs": [],
    "source": [
-    "print(regions_gdf.iloc[2].name)\n",
-    "voivodeship_region = regions_gdf.iloc[2:3]\n",
-    "voivodeship_geom = voivodeship_region.geometry[0]\n",
-    "voivodeship_geom"
+    "type(merged_geom)"
    ]
   },
   {
@@ -559,7 +516,7 @@
    },
    "outputs": [],
    "source": [
-    "voivodeship_geom.area"
+    "merged_geom.area"
    ]
   },
   {
@@ -573,7 +530,7 @@
    "outputs": [],
    "source": [
     "## minimum bounding region\n",
-    "voivodeship_geom.bounds"
+    "merged_geom.bounds"
    ]
   },
   {
@@ -584,7 +541,7 @@
     }
    },
    "source": [
-    "## Let's create a geometry object"
+    "## Creating a geometry manually"
    ]
   },
   {
@@ -598,7 +555,7 @@
    "outputs": [],
    "source": [
     "from shapely.geometry import LineString\n",
-    "bounds = voivodeship_geom.bounds\n",
+    "bounds = merged_geom.bounds\n",
     "line = LineString(\n",
     "    [(bounds[0], bounds[1]),\n",
     "    (bounds[2], bounds[3]),]\n",
@@ -627,7 +584,7 @@
    },
    "outputs": [],
    "source": [
-    "gpd.GeoSeries([line, voivodeship_geom]).plot(cmap='tab10')"
+    "gpd.GeoSeries([line, merged_geom]).plot(cmap='tab10')"
    ]
   },
   {
@@ -651,7 +608,7 @@
    },
    "outputs": [],
    "source": [
-    "line.within(voivodeship_geom)"
+    "line.within(merged_geom)"
    ]
   },
   {
@@ -664,7 +621,7 @@
    },
    "outputs": [],
    "source": [
-    "line.intersects(voivodeship_geom)"
+    "line.intersects(merged_geom)"
    ]
   },
   {
@@ -679,6 +636,15 @@
     "You can use the same spatial operations as in Shapely, on entire GeoDataFrames."
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "pl_de_gdf"
+   ]
+  },
   {
    "cell_type": "markdown",
    "metadata": {
@@ -752,7 +718,7 @@
    },
    "outputs": [],
    "source": [
-    "gpd.GeoSeries([line, voivodeship_geom]).plot(cmap='tab10')"
+    "gpd.GeoSeries([line, merged_geom]).plot(cmap='tab10')"
    ]
   },
   {
@@ -768,7 +734,7 @@
     "from srai.regionalizers import geocode_to_region_gdf, H3Regionalizer\n",
     "from utils import CB_SAFE_PALLETE\n",
     "\n",
-    "regionized = H3Regionalizer(resolution=6).transform(voivodeship_region)\n",
+    "regionized = H3Regionalizer(resolution=3).transform(pl_de_gdf)\n",
     "regionized[\"intersects\"] = regionized.intersects(line)\n",
     "regionized.explore(\"intersects\")"
    ]
@@ -816,7 +782,11 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {},
+   "metadata": {
+    "slideshow": {
+     "slide_type": "fragment"
+    }
+   },
    "outputs": [],
    "source": [
     "from srai.loaders import OSMOnlineLoader\n",
@@ -926,7 +896,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.10.11"
+   "version": "3.10.12"
   },
   "rise": {
    "controls": false,