如何去除Cartopy绘图中湖泊内部的河流线段？

你遇到的这个问题其实挺常见的——Cartopy默认的LAKES和RIVERS是两个独立的矢量图层，压根没做空间上的自动裁剪，所以河流就会“穿”过湖泊显示出来。你说试了zorder没效果，大概率是图层顺序或者zorder的设置搞反了，先试试下面这个最简单的办法：

方法1：调整图层绘制顺序/设置zorder #

只要让湖泊图层位于河流图层的上方，湖泊的填充色就会直接覆盖掉内部的河流线段。你可以直接调整添加图层的顺序，或者用zorder参数明确指定层级：

# 方式A：调整添加顺序（先加河流，再加湖泊，湖泊自然显示在上方）
ax.add_feature(cfeature.RIVERS, edgecolor='blue')
ax.add_feature(cfeature.LAKES, edgecolor='black', facecolor=cfeature.COLORS['water'])

# 方式B：用zorder强制设置层级（给湖泊更高的zorder值，确保它在河流上面）
ax.add_feature(cfeature.LAKES, edgecolor='black', facecolor=cfeature.COLORS['water'], zorder=2)
ax.add_feature(cfeature.RIVERS, edgecolor='blue', zorder=1)

如果这个方法能解决你的问题，那绝对是最省事的。但如果你的需求是湖泊需要透明显示，或者想真正移除河流在湖泊内的部分（而不是单纯被覆盖），那可以试试下面的空间裁剪方法：

方法2：用Geopandas+Shapely裁剪河流数据 #

这个方法需要先安装额外的依赖包：

pip install geopandas shapely

核心思路是：加载更高精度的河流和湖泊数据，裁剪到你的绘图范围，然后把河流中与湖泊重叠的部分剪掉，最后再添加到地图上。

完整代码示例：

import cartopy.crs as ccrs
import cartopy.feature as cfeature
import matplotlib.pyplot as plt
import geopandas as gpd
from shapely.ops import unary_union, difference
from cartopy.io import shapereader

# 定义你的绘图范围
llcrnrlon = 102.37615967
llcrnrlat = 47.20949936
urcrnrlon = 110.55716705
urcrnrlat = 55.39049911
proj = ccrs.PlateCarree()

# 加载Natural Earth 10m分辨率的河流和湖泊数据（Cartopy会自动下载）
river_shp = shapereader.natural_earth(resolution='10m', category='physical', name='rivers_lake_centerlines')
lake_shp = shapereader.natural_earth(resolution='10m', category='physical', name='lakes')

# 读取为GeoDataFrame格式，方便空间操作
rivers = gpd.read_file(river_shp)
lakes = gpd.read_file(lake_shp)

# 裁剪数据到我们需要的绘图范围
bbox = gpd.GeoDataFrame(
    {'geometry': [gpd.points_from_xy([llcrnrlon, urcrnrlon], [llcrnrlat, urcrnrlat]).envelope]},
    crs=proj
)
lakes_cropped = gpd.clip(lakes, bbox)
rivers_cropped = gpd.clip(rivers, bbox)

# 合并所有湖泊的几何形状，方便后续批量裁剪河流
lakes_union = unary_union(lakes_cropped.geometry)

# 处理每条河流，移除与湖泊重叠的部分
processed_rivers = []
for _, row in rivers_cropped.iterrows():
    if row.geometry.intersects(lakes_union):
        # 计算河流与湖泊的差集，只保留湖泊外的河流部分
        new_geom = difference(row.geometry, lakes_union)
        if not new_geom.is_empty:
            processed_rivers.append({'geometry': new_geom})
    else:
        processed_rivers.append({'geometry': row.geometry})

processed_rivers_gdf = gpd.GeoDataFrame(processed_rivers, crs=proj)

# 创建绘图并添加元素
fig, ax = plt.subplots(figsize=(10, 10), subplot_kw=dict(projection=proj))
ax.set_extent([llcrnrlon, urcrnrlon, llcrnrlat, urcrnrlat])

# 添加湖泊
ax.add_feature(cfeature.LAKES, edgecolor='black', facecolor=cfeature.COLORS['water'])

# 添加处理后的河流（已经去掉湖泊内的线段）
processed_rivers_gdf.plot(ax=ax, edgecolor='blue', linewidth=1, transform=proj)

# 添加经纬度网格
gl = ax.gridlines(crs=proj, draw_labels=True,
                  linewidth=1, color='gray', alpha=0.5, linestyle='--')
gl.top_labels = False
gl.right_labels = False
gl.xlabel_style = {'size': 12, 'color': 'black'}
gl.ylabel_style = {'size': 12, 'color': 'black'}

plt.show()

这个方法会真正移除河流在湖泊内的线段，而不是单纯覆盖，适合有特殊显示需求的场景。

备注：内容来源于stack exchange，提问作者Vyacheslav Gradov