cartopy - NorthPolarStereo 的中心经度

Question

我想用 180 在图的底部绘制北半球的极地立体图，以便我可以强调太平洋地区。我正在使用来自 git 的最新 cartopy，并且可以制作极坐标立体图没有问题，但我无法弄清楚如何更改图底部的经度。我尝试将经度范围设置为 [-180, 180] 但这没有帮助，而且 NorthPolarStereo() 不接受任何关键字参数，如 central_longitude。目前这可能吗？

Answer 1

此功能现已在 Cartopy (v0.6.x) 中实现。以下示例在北半球极地立体投影中生成两个子图，一个使用默认设置，另一个更改中心经度:

"""Stereographic plot with adjusted central longitude."""
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
from cartopy.examples.waves import sample_data


# read sample data
x, y, z = sample_data(shape=(73, 145))

fig = plt.figure(figsize=(8, 4))

# first plot with default settings
ax1 = fig.add_subplot(121, projection=ccrs.NorthPolarStereo())
cs1 = ax1.contourf(x, y, z, 50, transform=ccrs.PlateCarree(),
                   cmap='gist_ncar')
ax1.set_extent([0, 360, 0, 90], crs=ccrs.PlateCarree())
ax1.coastlines()
ax1.set_title('Centred on 0$^\circ$ (default)')

# second plot with 90W at the bottom of the plot
ax2 = fig.add_subplot(
    122, projection=ccrs.NorthPolarStereo(central_longitude=-90))
cs2 = ax2.contourf(x, y, z, 50, transform=ccrs.PlateCarree(),
                   cmap='gist_ncar')
ax2.set_extent([0, 360, 0, 90], crs=ccrs.PlateCarree())
ax2.coastlines()
ax2.set_title('Centred on 90$^\circ$W')

plt.show()

这个脚本的输出是: