0

我正在使用 cartopy 制作北极立体图,描绘极地有风的情况,使用 regrid_shape 提供更均匀的矢量密度。该场看起来大部分是正确的,但在看极点时有点紧张,具有一些似乎不在基础数据中的收敛和发散行为。在绘图之前,我已经从风中删除了 lat=90 值。

看我的阴谋

这是一个用不同软件制作的流线图。我认为这可能是实际风的合理表示。

可能是正确的风场

这是我的代码:

h=np.shape(lats)[0]//2+1  #to start from lat > 0

ax = plt.axes(projection=ccrs.NorthPolarStereo())
ax.set_extent([-180,180,30,90],crs=ccrs.PlateCarree())

ax.quiver(lons, lats[h:-2], u10[h:-2,:],v10[h:-2,:],
       transform=ccrs.PlateCarree(),regrid_shape=20)   

当我问的时候,我似乎也无法让 NorthPolarStereo(central_longitude=27) 与 set_extent 一起工作。

谢谢!

4

1 回答 1

1

我不认为这是一个错误,只是你的向量没有很好地缩放,使它们看起来像它们在极点附近会聚(或交叉)。你有一些强风朝向极点,然后迅速改变方向,矢量都被绘制成直线,所以如果在转弯之前绘制一个,它将超出实际轨迹。这是绘制向量的一般问题。

如果你玩缩放,你可以让它看起来更好,或者你可以尝试使用cartopy 的 streamplot方法来制作一个更像你的目标的绘图(使用 streamplot 的示例在这里:http:// scitools.org.uk/cartopy/docs/latest/examples/streamplot.html

于 2017-02-16T19:38:47.310 回答