merge main

momepy/functional/_elements.py

@@ -230,8 +230,13 @@
     """
 
     if simplify and not SHPLY_GE_250:
         raise ImportError("Coverage simplification requires shapely 2.5 or higher.")
 
+    if isinstance(geometry.index, MultiIndex):
+        raise ValueError(
+            "MultiIndex is not supported in `momepy.enclosed_tessellation`."
+        )
+
     # convert to GeoDataFrame and add position (we will need it later)
     enclosures = enclosures.geometry.to_frame()
     enclosures["position"] = range(len(enclosures))
@@ -320,10 +325,10 @@
             as_gdf=True,
         )
         if to_simplify:
             simpl_collection = shapely.coverage_simplify(
                 tess.geometry, tolerance=1e-1, simplify_boundary=False
             )
             tess.geometry = gpd.GeoSeries(simpl_collection.geoms).values
         tess[enclosure_id] = ix
         return tess
 

momepy/functional/tests/test_elements.py

@@ -334,47 +334,78 @@
 
     @pytest.mark.skipif(not SHPLY_GE_250, reason="coverage_simplify required")
     def test_simplified_tesselations(self):
         n_workers = -1
         tessellations = mm.enclosed_tessellation(
             self.df_buildings, self.enclosures.geometry, n_jobs=n_workers
         )
         simplified_tessellations = mm.enclosed_tessellation(
             self.df_buildings, self.enclosures.geometry, simplify=True, n_jobs=n_workers
         )
         ## empty enclosures should be unmodified
         assert_geodataframe_equal(
             tessellations[tessellations.index < 0],
             simplified_tessellations[simplified_tessellations.index < 0],
         )
         ## simplification should result in less total points
         orig_points = shapely.get_coordinates(
             tessellations[tessellations.index >= 0].geometry
         ).shape
         simpl_points = shapely.get_coordinates(
             simplified_tessellations[simplified_tessellations.index >= 0].geometry
         ).shape
         assert orig_points > simpl_points
 
         ## simplification should not modify the external borders of tesselation cells\
         orig_grouper = tessellations.groupby("enclosure_index")
         simpl_grouper = simplified_tessellations.groupby("enclosure_index")
         for idx in np.union1d(
             tessellations["enclosure_index"].unique(),
             simplified_tessellations["enclosure_index"].unique(),
         ):
             orig_group = orig_grouper.get_group(idx).dissolve().boundary
             enclosure = self.enclosures.loc[[idx]].dissolve().boundary
 
             simpl_group = simpl_grouper.get_group(idx).dissolve().boundary
 
             ## simplified is not different to enclosure
             ## this needs to be redone
             assert np.isclose(simpl_group.difference(enclosure).area, 0)
 
             # simplified is not different to original tess
             ## this needs to be redone
             assert np.isclose(simpl_group.difference(orig_group).area, 0)
 
+    def test_multi_index(self):
+        buildings = self.df_buildings.set_index(["uID", "uID"])
+        with pytest.raises(
+            ValueError,
+            match="MultiIndex is not supported in `momepy.morphological_tessellation`.",
+        ):
+            mm.morphological_tessellation(buildings)
+        with pytest.raises(
+            ValueError,
+            match="MultiIndex is not supported in `momepy.enclosed_tessellation`.",
+        ):
+            mm.enclosed_tessellation(buildings, self.enclosures)
+        with pytest.raises(
+            ValueError,
+            match="MultiIndex is not supported in `momepy.verify_tessellation`.",
+        ):
+            mm.verify_tessellation(buildings, self.enclosures)
+
+        with pytest.raises(
+            ValueError,
+            match="MultiIndex is not supported in `momepy.get_nearest_node`.",
+        ):
+            mm.get_nearest_node(
+                buildings, self.enclosures, self.enclosures, self.enclosures
+            )
+
+        with pytest.raises(
+            ValueError, match="MultiIndex is not supported in `momepy.generate_blocks`"
+        ):
+            mm.generate_blocks(buildings, self.enclosures, self.enclosures)
+
     def test_tess_single_building_edge_case(self):
         tessellations = mm.enclosed_tessellation(
             self.df_buildings, self.enclosures.geometry, n_jobs=-1