python - 极地立体图

Question

是否可以使用最新版本的 cartopy 在极地立体地图上绘制等高线？我想看一个例子来说明这是如何完成的，因为我自己正在努力解决这个问题！

Answer 1

立体投影引起了一些头痛，并且可能是对 cartopy 的多边形转换代码提出最多问题的投影。

以下示例显示了如何使用 cartopy 生成极坐标立体图。请注意：即使使用此代码，也可以调整示例数据分辨率并发现实际渲染该图需要大约 30 分钟（这是一个错误，我们需要尽快进行排序）。

import cartopy.crs as ccrs
import matplotlib.pyplot as plt

from cartopy.examples.waves import sample_data

ax = plt.axes(projection=ccrs.NorthPolarStereo())

x, y, z = sample_data((100, 200))
cs = ax.contourf(x, y, z, 50,
                 transform=ccrs.PlateCarree(),
                 cmap='gist_ncar')
ax.coastlines()

# without the set_global, currently, the plot is tiny because the limits
# are being erroneously being set (opened issue for that)
ax.set_global()

plt.show()

希望这将向您展示如何在 cartopy 中制作极地立体等高线图。如果您的数据有问题，请查看标记为“几何变换”的未解决问题，看看您是否得到类似的东西，如果没有，请继续打开一个问题，我们可以进行调查。

注意：这个答案与 cartopy v0.5.x（即在 v0.5 版本之前）有关，这里提到的许多错误应该有望在未来的版本中得到解决。

希望有帮助，