0

我正在尝试使用 python 3.6 在底图上绘制来自 netCDF4 文件的降水数据。我正在使用的代码在非常相似的数据集上运行良好。唯一的区别是之前的沉淀数据以“cm”为单位，而我目前尝试绘制的数据以“kg m-2 s-1”为单位。在此文件中找到的变量是 time、time_bnds、lat、lat_bnds、lon、lon_bnds 和 pr。pr 是降水变量，也是我对绘图感兴趣的变量。

这是我的代码

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*2592000,range(0,1000,10),cmap=cm.s3pcpn,latlon=True)

#data is multiplied by the amount of seconds in a month

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

plt.show()

如果您想尝试使用我正在使用的相同文件运行此代码，请从 CMIP5 网站下载 netCDF4 文件： http ://pcmdi9.llnl.gov/ 我的文件名：(/data/CCSM4/pr_Amon_CCSM4_historical_r1i1p1_185001-200512 .nc) 但是原始数据是 netCDF3_CLASSIC 格式，因此您必须将其更改为 netCDF4。

4

1 回答 1

1

如果唯一的变化是降水量的单位，那么如何将数据转换为厘米，这样您就知道它会起作用并且您已经熟悉如何解释数据。你甚至可以编写一个函数来转换它。

如果您转换它，并且它的格式与您期望的相同，但它仍然不起作用，您会发现您的数据存在其他问题，因此尝试这种方法还有另一个好处。

于 2017-06-05T20:34:46.553 回答