MeteoInfoLab脚本示例：TOMS HDF数据

MeteoInfo

TOMS (Total Ozone Mapping Spectrometer)数据是全球臭氧观测。

脚本程序：

1. #Add data file
   
2. folder = 'D:/Temp/hdf/'
   
3. fns = 'TOMS-EP_L3-TOMSEPL3_2000m0101_v8.HDF'
   
4. fn = folder + fns
   
5. f = addfile(fn)
   
6. vname = 'Ozone'
   
7. v = f[vname]
   
8. ozone = v[::-1,:]
   
9. #Plot
   
10. axesm()
    
11. mlayer = shaperead('D:/Temp/map/country1.shp')
    
12. geoshow(mlayer)
    
13. layer = imshowm(ozone, 20)
    
14. colorbar(layer, orientation='horizontal')
    
15. title('TOMS - Column Amount Ozone')
    
16. axism([-180,180,-90,90])

#Add data file <br /> folder = 'D:/Temp/hdf/' <br /> fns = 'TOMS-EP_L3-TOMSEPL3_2000m0101_v8.HDF' <br /> fn = folder + fns <br /> f = addfile(fn) <br /> vname = 'Ozone' <br /> v = f[vname] <br /> ozone = v[::-1,:] <br /> #Plot <br /> axesm() <br /> mlayer = shaperead('D:/Temp/map/country1.shp') <br /> geoshow(mlayer) <br /> layer = imshowm(ozone, 20) <br /> colorbar(layer, orientation='horizontal') <br /> title('TOMS - Column Amount Ozone') <br /> axism([-180,180,-90,90])