cartopy - 使用 cartopy Orthographic 和 RotatedPole 将经度偏移 180 度

Question

我试图将来自北极的蓝点放置在本初子午线下方（经度=0），但看到这些点沿着日期线（经度=180）向下移动。编码：

#!/usr/bin/env python
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
Lon = 0
ax = plt.axes( projection=ccrs.Orthographic( central_latitude=70,
                                             central_longitude=Lon) )
ax.set_global()
vector_crs = ccrs.RotatedPole( pole_latitude=90, pole_longitude=Lon )
ax.plot([Lon, Lon, Lon, Lon, Lon, Lon],   # longitude
        [ 90,  80,  70,  60,  50,  40],   # latitude
        'bo', markersize=5, transform=vector_crs)
ax.stock_img()
plt.show()

可能与变换有关，但我还没弄清楚是什么。Cartopy 版本 0.14.2，Python 3.6。

Answer 1

我认为问题来自您定义的变换投影：

ax = plt.axes(projection=vector_crs)
ax.coastlines()
plt.show()

请注意，变换投影中的这个简单的海岸线图看起来像中心经度为 180° 的 Plate Carree 图。考虑到这一点，让我们看一下使用 Plate Carree 投影在绘图上绘制样本数据点，以便尝试简化要绘制的地图：

ax = plt.axes(projection=ccrs.PlateCarree()) 
ax.plot([Lon, Lon, Lon, Lon, Lon, Lon],
        [ 90,  80,  70,  60,  50,  40],
        'bo', markersize=5, transform=vector_crs)
ax.stock_img()
plt.show()

与您的示例一样，这些点并未出现在我们预期的位置。最后，让我们在将点绘制到正交地图上时尝试使用 Plate Carree 投影作为点的变换：

ax = plt.axes(projection=ccrs.Orthographic(central_latitude=70,
                                           central_longitude=Lon))
ax.set_global()
ax.plot([Lon, Lon, Lon, Lon, Lon, Lon],
        [ 90,  80,  70,  60,  50,  40],
        'bo', markersize=5, transform=ccrs.PlateCarree())
ax.stock_img()
plt.show()

这似乎提供了更多您正在寻找的情节。