标签:
MERRA是NOAA的一种再分析资料,HDF数据遵循COARDS协议,读取比较简单。
脚本程序:
#Add data file folder = ‘D:/Temp/hdf/‘ fns = ‘MERRA300.prod.assim.tavg1_2d_slv_Nx.20150101.SUB.hdf‘ fn = folder + fns f = addfile(fn) vname = ‘t2m‘ t = f[vname] t2m = t[0,:,:] #Plot axesm() mlayer = shaperead(‘D:/Temp/map/country1.shp‘) geoshow(mlayer) layer = imshowm(t2m, 20) colorbar(layer, orientation=‘horizontal‘) title(vname + ‘ at time=0‘) axism()
标签:
原文地址:http://www.cnblogs.com/yaqiang/p/4629389.html