基线图未在0度纵向左侧绘制

2024-10-04 01:25:04 发布

您现在位置:Python中文网/ 问答频道 /正文

我对Python相当陌生,我正试图绘制非洲大陆的长期平均降雨量数据。我有各种NetCDF文件,它们已经被剪切成只包含长期平均值-我只需要绘制它。 我的问题是数据只在0度经线的右侧绘制。我想这是因为基础地图需要-180到180个坐标,我的数据是0到360。然而,我所做的一切似乎都不管用。在

下面是代码(它给出了正确的绘图,只在0度的左边截断):

nc = dataset(GISS-E2-H_MAM_plots.nc)

prcp = nc.variables['pr'][0,:,:]
pr = 86400*prcp[:]
lon=nc.variables['lon']
lat=nc.variables['lat']
[lonall, latall] = np.meshgrid(lon, lat)
fig = plt.figure()
m = Basemap(projection='cyl', llcrnrlat=-25, urcrnrlat=15, llcrnrlon=-20, urcrnrlon=60)
m.drawcoastlines()
m.drawcountries()
m.drawparallels(np.arange(-90.,90.,10.), labels = [1,0,0,0], fontsize = 10)
m.drawmeridians(np.arange(-180., 180., 10.), labels = [0,0,0,1], fontsize = 10)
levels=np.arange(2, 11.6, 0.8)
mymapf = plt.contourf(lonall, latall, pr, levels, cmap=plt.cm.gist_rainbow_r)

我试着用下面的方法将数据移动180,然后np.滚动一直移动它。在

^{pr2}$

在类似的情况下,这对一位同事有效,但对我没有效果。 任何帮助将不胜感激!在


Tags: 数据labelsnp绘制pltprvariableslon