参考此脚本:
- import glob
- year = 2015
- sday = 121
- eday = 122
- aods = []
- while sday <= eday:
- fn = glob.glob('D:/Temp/hdf/MOD08_D3.A%i%i*.hdf' % (year, sday))[0]
- print fn
- f = addfile(fn)
- v_aod = f['Optical_Depth_Land_And_Ocean_Mean']
- #Get data
- aod = v_aod[:,:]
- aods.append(aod)
- sday += 1
- #Mean
- aod = mean(aods)
- #Set x/y
- x = linspace(-180.0,180.0, 360)
- y = linspace(-90.0,90.0, 180)
- #Plot
- axesm()
- geoshow('country', linecolor='gray')
- levels = arange(0, 1.0, 0.02)
- layer = imshowm(x, y, aod, levels)
- title('Aerosol Optical Depth at 550 nm')
- colorbar(layer)
|