python - 如何将接触的多边形与 geopandas 合并

Question

我需要将 shapefile 中的所有接触多边形与 geopandas 合并。那么有没有解决办法。例如，下图必须是一个多边形，而不是多个多边形。

谢谢

Answer 1

这取决于您希望如何处理附加到这些多边形的数据。如果您不关心它们并且只对几何图形感兴趣，一种选择是解散整个 GeoDataFrame，然后将其分解。

combined_polygons = gdf.dissolve().explode()

但是，这可能不是最有效的解决方案。最好的方法是确定邻接组件并根据这些组件进行分解。您可以使用 libpysal 轻松做到这一点。

import libpysal

# create spatial weights matrix
W = libpysal.weights.Queen.from_dataframe(gdf)

# get component labels
components = W.component_labels

combined_polygons = gdf.dissolve(by=components)

后者将允许您在 dissolve 中指定 aggfunc 以管理其他属性。