Python 3.x 如何使用cartopy在地磁坐标中绘制数据？

我是python初学者，我使用cartopy将数据指向地理坐标。然而，我想用地磁坐标来绘制我的数据。我看到cartopy中有一个旋转杆选项。我不确定的是：如果我使用磁极的位置（纬度和经度）旋转磁极，这会给我正确的结果吗？因为南北磁极彼此不对称（因为太阳风的影响），但地理磁极是对称的。据我所知，cartopy中的旋转极点选项根据给定的新极点位置对称地移动地理极点

此外，地球上磁北极的纬度和经度每年都在变化。所以，我想，对于我们感兴趣的每个日期，磁北极的位置应该分别更新。这些值可以从多个地方获得，因此这不是问题

接下来是我的实际问题：我应该在下面的代码中做些什么来正确地绘制地磁坐标中的数据？在这种情况下，cartopy中的旋转杆选项是否有效？任何关于处理此类问题的其他方法的建议也将不胜感激

这是我用于在地理坐标中绘制的内容：

# set the map coverage:
extent = [-90, -60, 30, 60]
central_lon = np.mean(extent[:2])
central_lat = np.mean(extent[2:])

fig2 = plt.figure(figsize=(7, 7))
ax2 = plt.axes(projection=ccrs.Orthographic(central_longitude=central_lon, central_latitude=central_lat))
ax2.set_extent(extent, ccrs.PlateCarree())
ax2.add_feature(cartopy.feature.OCEAN, color='white', alpha=1, zorder=0)
ax2.add_feature(cartopy.feature.LAND, edgecolor='white', 
                color='silver', alpha=0.3, zorder=10)
ax2.add_feature(cartopy.feature.LAKES, color='white', alpha=1, zorder=0)

# plot the pointing direction.
for i in range (0,ind1):
# downsampling to every 10th arrow ([::10])
    L=i*10
# Orientation (X body vector)
    if 360+lon[L] > 360+OLon: # if spacecraft is in the east 
        ax2.quiver(np.array([lon[L]]), np.array([lat[L]]), 
                np.array([-XV[L][2]]), np.array([XV[L][1]]), 
                color='black', scale=7, width=.006, axes=ax2,
                transform=ccrs.PlateCarree(), angles = "xy", zorder=20)
ax2.text(-0.09, 0.5, 'Geographic Latitude', fontsize=14, va='bottom', ha='center',
        rotation='vertical', rotation_mode='anchor',
        transform=ax2.transAxes)
ax2.text(0.5, -0.12, 'Geographic Longitude', fontsize=14, va='bottom', ha='center',
        rotation='horizontal', rotation_mode='anchor',
        transform=ax2.transAxes)

---

如果在v小尺度上绘制数据进行可视化，南北磁极不对称的事实并不重要。但是，当我想要绘制全局图时，这种差异会产生影响，对吗@斯沃特柴