极地立体图

Question

是否可以绘制最新版本的极地立体地图上的等高线？ 我想看看这样做的一个例子，因为我正在努力解决它自己！

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极立体轮廓图。如果您的数据有问题，请查看open issues tagged "Geometry transforms"，看看您是否收到类似的东西，如果没有，请拨打go ahead and open an issue，我们可以查看它。

注意：此答案与cartopy v0.5.x（即刚刚发布到v0.5之前）有关，此处提及的许多错误应该有望在未来的版本中被压缩。

希望帮助，