溶解在geopandas中后删除多边形的内线

Question

我正在使用一个城市的地理数据框，其中每个单位都是一个地区（行政部门）。 它的情节是这样的：

import geopandas as gpd
df = gpd.read_file('districts_lima.geojson')
df.plot()

然后，我使用名为zone的变量将一些地理单位合并为更大的组。 结果是：

df2 = df.dissolve(by='zone', aggfunc='sum')
df2.plot(column='population', legend=True, cmap='Blues')

我唯一的问题是，当我用较暗的边界重现同一个图时，很明显一些合并的多边形（区域）有内线，这些线是原始地理数据框的内区边界。 这在此图中清楚地显示：

df2.plot(column='population', legend=True, cmap='Blues', edgecolor='black')

有没有办法使用 geopandas 删除多边形的内线，这样它们就不会出现在最后一个图中？

我的数据可以在这里找到。

Answer 1

我实际上找到了一个很好的解决方案，它特别适用于我的问题是在应用geopandas dissolve()属性后geopandas 。 显然，问题是由相邻内部单元的边界线中不明显的差异产生的，这阻止了折叠以删除生成的多边形的内部线。

这是通过添加一个小缓冲区来加宽多边形线来解决的，这样那些不明显的差异就被消除了，并且初始多边形的每条内线实际上都重叠了。 具体来说，我需要添加：

df2['geometry'] = df2['geometry'].buffer(0.0001)

前

df2 = df.dissolve(by='zone', aggfunc='sum')

所以现在 plot 命令

df2.plot(column='population', legend=True, cmap='Blues', edgecolor='Black')

产量：

Answer 2

使用轴 (ax1) 启用在公共轴上绘制多个图层。 还可以使用 zorder 来排列图层，较高的值将图层置于较低值的图层之上。

fig, ax1 = plt.subplots(1, 1)
...
# bottom layer (has thick black lines)
df2.plot(column='population', legend=True, cmap='Blues', edgecolor='black', ax=ax1, zorder=5, lw=6)
# top layers (thin boundary lines)
df2.plot(column='population', legend=True, cmap='Blues', ax=ax1, zorder=6)

希望顶层将隐藏几乎所有底层，但部分外部边界线。