1

我正在分析天气模型的输出，使用 python 模块 netCDF4 并使用 mpl_toolkits Basemap 制作图像。

我做了一些气压的等高线图，这很容易。现在我想为风添加倒钩，最后得到一个只有海风的情节。经过一番摆弄，我发现底图函数 fillcontinents() 实际上覆盖了这些风倒钩，即使稍后调用了 barbs 函数。因此，删除 fillcontinents() 函数后，所有的风倒钩都得到了想要的结果。

有没有办法填满大陆，但情节中仍然有我所有的倒钩？

我的代码看起来有点像这样（我遗漏了一些图形参数，它们对我的问题并不重要）：

from netCDF4 import Dataset
from mpl_toolkits.basemap import Basemap
import numpy as np

file = 'W:/projects/wrfout_d01_2013-10-27_120000'
data = Dataset(file)

# Load grid
longrid = data.variables['XLONG'][:].squeeze() # Removes time dimension
latgrid = data.variables['XLAT'][:].squeeze()

# Load surface pressure
p_surf = data.variables['PSFC'][:].squeeze()

# Load wind components
u = data.variables['U10'][:].squeeze()
v = data.variables['V10'][:].squeeze()

m = Basemap(projection='cyl', resolution='l',
                    llcrnrlat= 42., llcrnrlon= -12.,
                    urcrnrlat= 57., urcrnrlon= 12.)
m.drawcountries()
m.drawcoastlines()
m.fillcontinents('#FAF0E6') ### Blanking out this line gives barbs over land
m.drawmeridians(np.arange(-180,180,10),labels=[0,0,1,1])
m.drawparallels(np.arange(-90,90,10),labels=[1,0,0,0])

x,y = m(longrid,latgrid)
ct = m.contour(x,y,p_surf, colors='#1C8CFC') ### So far, so good
ct.clabel(fmt='%1i')

m.barbs(x,y,u,v)

plt.savefig('p_surf-barbs-map.png')
4

1 回答 1

0

在 Paul 的帮助下，我找到了解决方案：

绘制倒钩时，添加一个具有高值（我使用 100）的变量 zorder，如下所示：

m.barbs(x, y, u, v, zorder=100)
于 2016-11-01T08:43:41.967 回答