Python:底图降水数据未显示

时间:2017-06-02 21:24:10

标签: python matplotlib-basemap netcdf4

我使用python 3.6绘制CMIP5的降水数据,我下载的文件是netCDF4文件。我已经在另一个类似的文件上使用了这个代码,它运行得很好,所以我不确定是什么问题。我没有收到任何有关此代码的错误消息,它只显示一个世界地图,它应该是各种颜色的所有颜色。在此文件中找到的变量是time,time_bnds,lat,lat_bnds,lon,lon_bnds和prc。 prc是降水变量,也是我对绘图感兴趣的变量。任何想法都会有所帮助,谢谢!

这是我的代码

from mpl_toolkits.basemap import Basemap, cm
from netCDF4 import Dataset as NetCDFFile
import matplotlib.pyplot as plt

nc = NetCDFFile('filename.nc','r')

p = nc.variables['prc']
data = p[:,:,0]

fig = plt.figure(figsize=(8,8))
ax = fig.add_axes([0.1,0.1,0.8,0.8])

m = Basemap(projection='cyl',lon_0=180,lat_0=0,resolution='l')

m.drawcoastlines()
m.drawstates()
m.drawcountries()


ny = data.shape[0]; nx = data.shape[1]
lons, lats = m.makegrid(nx,ny) 
x,y = m(lons, lats) # compute map proj coordinates.

cs=plt.contourf(x,-y,data,range(0,1000,10),cmap=cm.s3pcpn,latlon=True)


cbar = m.colorbar(cs,location='bottom',pad="5%")
cbar.set_label('mm')

plt.show() 

0 个答案:

没有答案