Py-Eddy-Tracker中地图投影与轨迹显示的同步调整技巧-CSDN博客

本文链接：https://blog.csdn.net/gitblog_07634/article/details/149012089

Py-Eddy-Tracker中地图投影与轨迹显示的同步调整技巧

py-eddy-tracker 项目地址: https://gitcode.com/gh_mirrors/py/py-eddy-tracker

在使用Py-Eddy-Tracker进行涡旋轨迹可视化时，经常会遇到地图投影与轨迹显示不匹配的问题。本文将以一个典型场景为例，详细介绍如何正确设置投影参数，确保地图背景与轨迹数据同步变换。

问题背景

当我们需要将地图中心从默认的大西洋区域调整到太平洋区域时（例如设置central_longitude=180），虽然地图背景能够正确变换，但涡旋轨迹数据却仍然保持原来的显示方式。这种现象在使用Cartopy进行地理数据可视化时较为常见。

核心概念

数据投影(transform)：原始数据的坐标参考系统
目标投影(projection)：最终地图显示的投影方式
参考经度(ref)：轨迹数据显示的经度参考系

解决方案

正确的做法是在绘制轨迹时明确指定数据的原始投影系统，而非目标投影系统。具体需要修改transform参数为ccrs.PlateCarree()，而不是使用与地图相同的投影对象。

# 正确设置方式
a_20y.plot(ax, ref=180, transform=ccrs.PlateCarree())

原理分析

transform参数：这个参数应该表示输入数据的坐标系统，而不是你想要转换到的坐标系统。大多数地理数据默认使用PlateCarree投影（等距圆柱投影）。
ref参数：用于调整轨迹显示的经度参考系，与投影变换配合使用可以实现地图中心的调整。
投影链：Cartopy内部会自动处理从数据投影到目标投影的转换，开发者只需要正确指定这两个端点即可。

完整示例代码

import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cfeature
from cartopy.mpl.gridliner import LongitudeFormatter, LatitudeFormatter

# 设置中心经度为180度（太平洋中心）
lonshift = 180
proj = ccrs.PlateCarree(central_longitude=lonshift)

# 创建图形和坐标轴
fig, ax = plt.subplots(1, 1, figsize=(40, 15), subplot_kw={'projection': proj})

# 配置网格线
gl = ax.gridlines(draw_labels=True, y_inline=True)
gl.top_labels = False
gl.xformatter = LongitudeFormatter()
gl.yformatter = LatitudeFormatter()

# 添加陆地特征
land = cfeature.NaturalEarthFeature(
    category='physical',
    name='land',
    scale='10m',
    facecolor='grey')
ax.add_feature(land)

# 绘制涡旋轨迹（注意transform的设置）
a_20y.plot(ax, ref=180, label="Anticyclonic", color="r", lw=1, transform=ccrs.PlateCarree())
c_20y.plot(ax, ref=180, label="Cyclonic", color="b", lw=1, transform=ccrs.PlateCarree())

ax.legend(fontsize=20, loc='lower right')