Geopandas 0.10.1中使用overlay how='difference'且keep_geom_type=True时返回GeometryCollection是否为Bug？

Hey there, this definitely looks like an edge-case bug that wasn't fully resolved by PR #1582, even though that PR was meant to fix issues with keep_geom_type not enforcing the expected geometry type. Let's break this down:

What's Going Wrong #

The keep_geom_type=True parameter is supposed to guarantee that the output GeoDataFrame only contains geometry types matching the left input (in your case, Polygons/MultiPolygons). But when your left MultiPolygon shares boundaries with multiple right polygons, the difference operation is leaving behind line string "cuts" that get bundled into a GeometryCollection instead of being filtered out or merged into a valid MultiPolygon.

This is likely an oversight in the PR #1582 fix—while it addressed common cases, it didn't account for shared boundary scenarios where leftover linear elements end up in the result.

Quick Workarounds #

Until an official fix is released, you can clean up the invalid geometries manually:

1. Extract Polygons from GeometryCollections
   Write a helper function to strip out non-polygon elements from any stray GeometryCollections:

   from shapely.geometry import GeometryCollection, Polygon, MultiPolygon
   
   def clean_geometry(geom):
       if isinstance(geom, GeometryCollection):
           # Grab only polygon/multipolygon components
           valid_polys = [g for g in geom.geoms if isinstance(g, (Polygon, MultiPolygon))]
           if len(valid_polys) == 1:
               return valid_polys[0]
           elif len(valid_polys) > 1:
               return MultiPolygon(valid_polys)
           else:
               return None  # No valid polygons to keep
       return geom
   
   # Apply to your result
   res_difference['geometry'] = res_difference['geometry'].apply(clean_geometry)
   res_difference = res_difference.dropna(subset=['geometry'])
   

2. Re-Validate Geometries
   Even with make_valid=True in the overlay call, running an extra validation pass can help ensure clean geometries:

   res_difference['geometry'] = res_difference['geometry'].make_valid()
   

Next Steps for Confirming the Bug #

• Check newer Geopandas versions: The fix from PR #1582 might have been refined in releases after 0.10.1 (try 0.11.x or later to see if the issue is resolved).
• Report the bug: Head to the Geopandas GitHub repo and open a new issue. Include a minimal reproducible example (simplified versions of your left/right dataframes) and the exact steps you took—this will help maintainers replicate and fix the edge case.

内容的提问来源于stack exchange，提问作者banchan