我想这就是你的想法

import geopandas as gpd

# load your geodataframe ..

# self join on geodataframe to get all polygon intersections
intersects = gdf.sjoin(gdf, how="left", predicate="intersects")

# dissolve intersections on right index indices using the minimum value
intersects_diss = intersects.dissolve("id_right",aggfunc="min")

# dissolve again on left index using minimum
intersects_diss = intersects_diss.reset_index().dissolve("id_left",aggfunc="min")

任何遇到这个问题的人都想知道如何解决这个问题，如果我们有“多边形链”，例如：

• 多边形A与B相交
• 多边形B与A和C相交
• 多边形C与B相交

这就是你如何解决它（灵感来自@NielsFlohr的答案）：

import geopandas as gpd

# load your geodataframe ..

# self join on geodataframe to get all polygon intersections
intersects = gdf.sjoin(gdf[["geometry"]], how="left", predicate="intersects").reset_index()

# dissolve intersections on right index indices aggregating to list
intersects_diss = intersects.dissolve(
    "index_right",
    aggfunc=lambda x: x.tolist(),
)
intersects_diss["id"] = gdf["id"]
intersects_diss.geometry = intersects_diss.geometry.buffer(0)

# determine which polygons are connected
intersects_diss["polygon_group"] = None
for i in range(len(polygons)):
    if intersects_diss.loc[i, "polygon_group"] is None:
        intersects_diss.loc[i, "polygon_group"] = i
    poly_group = intersects_diss.loc[i, "polygon_group"]
    intersects_diss.loc[intersects_diss["index"].apply(lambda x: i in x), "polygon_group"] = poly_group

# dissolve on the polygon_group
intersects_diss = intersects_diss.dissolve(
    "polygon_group",
    aggfunc="min",
)
intersects_diss.geometry = intersects_diss.geometry.buffer(0)

可能不是最好的方法，但它有效。