Basemap中如何用海岸线级高分辨率填充海洋以掩盖低分辨率等高线？

• 核心思路：直接复用Basemap加载的高分辨率海岸线多边形，通过多边形布尔运算生成海洋区域，完全避免高分辨率栅格掩码的资源浪费，同时匹配海岸线的细节精度。

实现步骤及代码 #

1. 初始化Basemap并加载高分辨率海岸线

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
from shapely.geometry import Polygon, MultiPolygon
from shapely.ops import unary_union

# 初始化Basemap（根据业务需求调整投影和范围）
m = Basemap(projection='merc', llcrnrlat=-60, urcrnrlat=60,
            llcrnrlon=-180, urcrnrlon=180, resolution='h')

# 加载高分辨率海岸线，先隐藏线条后续再单独绘制
coast_paths = m.drawcoastlines(resolution='h', linewidth=0).get_paths()

2. 提取并合并陆地多边形

# 将海岸线路径转换为Shapely多边形并去重合并
land_polys = []
for path in coast_paths:
    coords = path.vertices
    land_polys.append(Polygon(coords))
land_union = unary_union(land_polys)

3. 生成海洋区域并填充

# 创建覆盖整个地图范围的多边形
map_bounds = Polygon([(m.xmin, m.ymin), (m.xmax, m.ymin),
                      (m.xmax, m.ymax), (m.xmin, m.ymax)])

# 计算海洋区域：地图范围减去所有陆地
ocean_region = map_bounds.difference(land_union)

# 填充海洋区域
if isinstance(ocean_region, MultiPolygon):
    for poly in ocean_region.geoms:
        x, y = poly.exterior.xy
        m.fill(x, y, color='#2081C3', zorder=2)
else:
    x, y = ocean_region.exterior.xy
    m.fill(x, y, color='#2081C3', zorder=2)

4. 叠加业务图层

# 绘制低分辨率等高线（zorder高于海洋，确保显示在上方）
# m.contourf(your_lon, your_lat, your_data, zorder=3)

# 重新绘制高分辨率海岸线（zorder最高，保证边界清晰）
m.drawcoastlines(resolution='h', linewidth=0.5, color='k', zorder=4)

扩展：包含湖泊填充 #

如果需要处理湖泊，提取湖泊路径后从海洋区域中剔除，再单独填充湖泊：

# 加载湖泊路径
lake_paths = m.drawlakes(linewidth=0).get_paths()
lake_polys = [Polygon(path.vertices) for path in lake_paths]
lake_union = unary_union(lake_polys)

# 海洋区域减去湖泊，得到纯海洋区域
ocean_only = ocean_region.difference(lake_union)
# 填充湖泊
for poly in lake_union.geoms:
    x, y = poly.exterior.xy
    m.fill(x, y, color='#87CEEB', zorder=2)

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