0

我正在使用 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()
4

0 回答 0