我正在尝试在matplotlib的填充等高线图上叠加轮廓以获取某些大气数据。但是,我的轮廓标签并不总是出现在屏幕上。以下是一个例子:
如您所见,轮廓标签仅出现在最里面的几个轮廓上。
知道我的轮廓范围早先定义为
list(range(950,1052,4))
我有以下代码来实际绘图:
parallels = np.arange(0.,90,5.)
basem.drawparallels(parallels,labels=[1,0,0,0],fontsize=10)
# draw meridians
meridians = np.arange(180.,360.,5.)
basem.drawmeridians(meridians,labels=[0,0,0,1],fontsize=10)
basem.drawstates()
basem.drawcountries()
if clevs != 0:
cs = basem.contourf(x,y, plotted_var, clevs)
cl = basem.contour(x,y, plotted_var, clevsl, colors='k')
plt.clabel(cl, fmt="%1.0f", fontsize=8)
else:
cs = basem.contourf(x,y, plotted_var, cmap=plt.get_cmap(colorbar),
vmin = vmin, vmax = vmax)
cbar = basem.colorbar(cs, location='bottom', pad = "5%")
cbar.set_label(units)
此外,我的底图定义是:
basem = Basemap(width=5800000,height=3000000,
rsphere=(6378137.00,6356752.3142),\
resolution='h',area_thresh=1000.,projection='lcc',\
lat_1=45.,lat_2=55,lat_0=40,lon_0=-102.)
这是一个错误还是只有我失踪的东西?如果可以的话,我试图避免使用手册。
输入数据是全局数据集(GFS天气模型)。 x,y由下式获得:
lons2, lats2 = np.meshgrid(lons, lats)
x,y = basem(lons2, lats2)
lons,lats是:
lons = [0.0, 0.25, 0.5, 0.75, 1.0, ..., 359.75, 360.0]
lats = [-90, -89.75, ..., 89.75, 90]
答案 0 :(得分:1)
我似乎已经解决了这个问题:
我必须做的是将数据绑定到底图上可以看到的内容。我的输入数据集是一个全局数据集,当我绑定它时(见下面的白色区域),轮廓标签主要出现在地图边界内。这似乎仍然是一个错误 - 我选择的显示区域几乎改变了所有其他气象图形程序(如GrADS)中轮廓标签的位置,但我现在要对此进行标记。