我想绘制风向量。我使用Basemap和这个例子http://basemaptutorial.readthedocs.org/en/latest/plotting_data.html#quiver。 在我的文件中,经度从0到360.我使用latlon = True将其转换为-180:180。但是,当我使用“点”来定义点密度时,矢量仅绘制在我图的一半上。 这是代码: 当然,我使用这段代码:
m =底图(projection ='cyl',llcrnrlat = 30,urcrnrlat = 80,llcrnrlon = -40,urcrnrlon = 40,resolution ='c')
latvar = nc.variables ['lat']
lat = latvar [:]
lon = nc.variables ['lon'] [:]
X,Y = m(lon,lat)
lons,lats = meshgrid(lon,lat)
X4,Y4 = m(离子,拉特)
varU = ncU.variables ['var1'] [0,0,:,:]
varV = ncV.variables ['var2'] [0,0,:,:]
speed = np.sqrt(varUvarU + varVvarV)
yy = np.arange(0,len(Y [:]),3)
xx = np.arange(0,len(X [:]),3)
points = np.meshgrid(yy,xx)
m.quiver(X4 [分],Y4 [分],varU [分],天花病毒[分],速度[分],CMAP = CMAP,latlon =真)
有人可以帮我解决这个问题吗?谢谢。
答案 0 :(得分:1)
你有很多事情发生,我不确定你是否需要。我已经尝试了一些微小的修改代码,它可以工作:
md-is-locked-open