GDF.simplify 弄乱了几何图形

Question

我正在尝试在地图上绘制河流流域。为了减小生成的矢量图的大小，我申请了GeoSeries.simplify()。

import cartopy
import cartopy.crs as ccrs

import matplotlib.pyplot as plt
import geopandas as gpd

# %%% Earth
fig = plt.figure()

latlon_proj = ccrs.PlateCarree()
axis_proj = ccrs.Orthographic()

ax = plt.axes(
    projection=axis_proj
    )

# %%% Major River Basins

mrb_basins = gpd.read_file('mrb_basins.json') # 520 entries

mrb_basins['geometry'] = mrb_basins['geometry'].simplify(0.1)

for shape in mrb_basins['geometry']:

    feat = cartopy.feature.ShapelyFeature(
        [shape],
        latlon_proj,
        facecolor='red',
        )

    ax.add_feature(feat)

mrb_basins.plot()

问题是，生成的地球地图完全被红色形状覆盖。

如果我删除mrb_basins['geometry'] = mrb_basins['geometry'].simplify(0.1)这一行，情况并非如此。

如何在保持完整性的同时简化几何？

主要流域数据集可在here获取。

Answer 1

GeoSeries.simplify() 并不总是返回有效的几何图形，因为 GEOS 使用了底层简化算法。而且 cartopy 无法绘制无效的几何图形。

您需要先修复几何图形，然后再将它们传递给 cartopy。简单的技巧是调用buffer(0)。

mrb_basins['geometry'] = mrb_basins['geometry'].simplify(0.1).buffer(0)

那么你的代码就可以正常工作了。