Cartopy遇含NaN数据时投影范围自动变更的问题求助

嘿，我当初从Basemap转Cartopy的时候也踩过这个坑！默认情况下Cartopy确实会根据你绘制的数据自动缩放地图范围，但要固定住绘图域其实很简单，核心就是提前手动锁定地图的显示范围，不管数据里有没有NaN都不会改变。

下面给你两种最常用的解决方法，附代码示例：

1. 用ax.set_extent()固定范围（推荐） #

这是最直接的方式，就像你在Basemap初始化时指定llcrnrlon、llcrnrlat等参数一样。你只需要在绘制任何数据之前，调用set_extent()方法，传入你想要的固定区域的经纬度（或投影坐标）范围，同时指定对应的坐标参考系（CRS）。

举个完整的例子：

import cartopy.crs as ccrs
import matplotlib.pyplot as plt
import numpy as np

# 模拟带NaN的数据集
lon = np.linspace(-120, -60, 100)
lat = np.linspace(20, 50, 100)
data = np.random.rand(100, 100)
data[20:40, 20:40] = np.nan  # 制造一片无效值区域

# 创建地图子图（用Mercator投影为例）
fig, ax = plt.subplots(subplot_kw={'projection': ccrs.Mercator()})

# 关键步骤：固定绘图范围（这里用经纬度指定，所以CRS用PlateCarree）
ax.set_extent([-120, -60, 20, 50], crs=ccrs.PlateCarree())

# 先绘制地图背景（海岸线、地形等）
ax.coastlines(linewidth=0.8)
ax.stock_img()  # 用内置的地形底图，也可以替换成你自己的地形数据

# 再叠加数据，即使有NaN，地图范围也不会变
im = ax.imshow(
    data,
    extent=[lon.min(), lon.max(), lat.min(), lat.max()],
    origin='lower',
    transform=ccrs.PlateCarree(),  # 数据的坐标系统
    cmap='viridis'
)

plt.colorbar(im, shrink=0.8)
plt.show()

注意要点： #

• set_extent()的第一个参数是[xmin, xmax, ymin, ymax]，顺序别搞反；
• 第二个参数crs必须和你传入的范围的坐标系统一致：如果是经纬度范围，就用ccrs.PlateCarree()；如果是投影后的坐标（比如米单位），就用对应的投影CRS（比如ccrs.Mercator()）；
• 一定要先调用set_extent()再绘制数据，这样后续的绘图操作就不会触发自动缩放了。

2. 绘制数据时指定extent参数（适用于imshow/pcolormesh） #

如果你习惯在绘制数据时直接控制范围，也可以在调用imshow()或pcolormesh()时传入extent参数，同时关闭Cartopy的自动缩放功能：

# 在绘制数据时添加这两个参数
im = ax.imshow(
    data,
    extent=[-120, -60, 20, 50],  # 直接指定数据显示的范围
    origin='lower',
    transform=ccrs.PlateCarree(),
    cmap='viridis',
    aspect='auto'  # 避免地图被拉伸变形
)
# 额外锁定坐标轴范围
ax.set_xlim(-120, -60)
ax.set_ylim(20, 50)

不过这种方法不如第一种灵活，尤其是当你同时绘制多个图层时，还是推荐用set_extent()统一控制。

总的来说，只要提前用set_extent()固定好你想要的区域，Cartopy就会像Basemap一样保持绘图域不变，NaN区域自然会显示为白色（或你设置的NaN颜色），完全符合你的需求～

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