如何使用Cartopy将等值线图中心定位至太平洋区域

核心问题分析 #

使用ccrs.PlateCarree(central_longitude=180)时，若你的经度数据是-180~180范围，会导致太平洋区域被切割，出现异常。需要先调整数据的经度范围或添加循环点，匹配投影中心的设置。

具体解决方案 #

1. 处理经度数据，避免区域断裂 #

如果你的经度数组是-180到180，可以直接转换为0到360范围：

# 假设你的经度变量为lon
lon = np.where(lon < 0, lon + 360, lon)

或者使用Cartopy自带的工具添加循环点，修复180°经线附近的断裂问题：

from cartopy.util import add_cyclic_point

# 假设sst_anomaly是海表温度距平数据，lon是经度数组
sst_anomaly_cyclic, lon_cyclic = add_cyclic_point(sst_anomaly, coord=lon)

2. 正确配置投影与绘图 #

使用ccrs.PlateCarree(central_longitude=180)作为地图投影，绘图时明确指定数据的原始坐标系为ccrs.PlateCarree()。完整示例代码：

import matplotlib.pyplot as plt
import cartopy.crs as ccrs
from cartopy.util import add_cyclic_point
import numpy as np

# 替换为你的真实数据
lon = np.arange(-180, 180, 2)
lat = np.arange(-90, 90, 2)
lon_grid, lat_grid = np.meshgrid(lon, lat)
sst_anomaly = np.sin(lat_grid/30) * np.cos(lon_grid/60)

# 处理循环点，消除太平洋区域的断裂
sst_anomaly_cyclic, lon_cyclic = add_cyclic_point(sst_anomaly, coord=lon)

# 创建太平洋中心的绘图对象
fig, ax = plt.subplots(figsize=(10, 6), subplot_kw={'projection': ccrs.PlateCarree(central_longitude=180)})

# 添加基础地图元素
ax.coastlines()
ax.add_feature(ccrs.cartopy.feature.LAND, facecolor='lightgray')

# 绘制等值线填色图
contour = ax.contourf(lon_cyclic, lat_grid, sst_anomaly_cyclic, 
                      transform=ccrs.PlateCarree(),
                      cmap='RdBu_r', levels=20)

# 设置易读的经纬度标签
ax.set_xticks([120, 160, 180, -160, -120], crs=ccrs.PlateCarree())
ax.set_yticks([-60, -30, 0, 30, 60], crs=ccrs.PlateCarree())
ax.set_xticklabels(['120°E', '160°E', '180°', '160°W', '120°W'])
ax.set_yticklabels(['60°S', '30°S', '0°', '30°N', '60°N'])

# 添加色标
plt.colorbar(contour, ax=ax, orientation='horizontal', pad=0.05, label='Sea Surface Temperature Anomaly (°C)')

plt.title('Pacific-Centered Sea Surface Temperature Anomaly')
plt.show()

3. 常见异常排查 #

• 若等值线在180°经线附近断裂：必须用add_cyclic_point处理数据，补上循环点。
• 若投影后数据位置偏移：确保绘图时指定了transform=ccrs.PlateCarree()，明确数据的原始坐标系与地图投影的区别。

内容的提问来源于stack exchange，提问作者jose paulo Oliveira