Meteoinfo Script Example: Reading fy3a AOD HDF file

FY3A Satellite has AOD product data, HDF format, here is an example to read and display this type of data using the Meteoinfo script program.

The script is as follows:

#-----------------------------------------------------#Author:yaqiang Wang#date:2015-3-18#purpose:read fy3a AOD hdf5 data#Note:sample#----------------------------------------------------- fromOrg.meteoinfo.data.meteodataImportMeteodatainfo fromOrg.meteoinfo.data.meteodataImportDimension fromOrg.meteoinfo.data.meteodataImportDimensiontype fromOrg.meteoinfo.data.meteodataImportDrawmeteodata fromOrg.meteoinfo.legendImportLegendscheme fromOrg.meteoinfo.shapeImportShapetypes fromUcar.nc2ImportNetcdffileImportOs.pathImportJarray#Set Data DirectoryDataDir ='D:/TEMP/HDF'#Create Meteodatainfo ObjectMDI =Meteodatainfo ()#Read HDF Data Filefn = Os.path.join (DataDir,'FY3A_MERSI_GBAL_L3_ASL_MLT_GLL_20140331_AOAM_5000M_MS.HDF')ifOs.path.isfile (FN):Printfn Mdi.opennetcdfdata (FN) Datainfo=Mdi.getdatainfo ()#print Mdi.getinfotext ()Xmin =-30.0ymin=-90.0Xnum= 7200Ynum= 3600Xdelt= 0.05Ydelt= 0.05xlist=[] ylist= []     forIinchRange (0,xnum): Xlist.append (xmin+ Xdelt *i) forIinchRange (0,ynum): Ylist.append (ymin+ Ydelt *i) X= Jarray.array (Xlist,'D') Y= Jarray.array (Ylist,'D') XDim=Dimension (dimensiontype.x) xdim.setvalues (X) datainfo.setxdimension (xDim) Ydim=Dimension (DIMENSIONTYPE.Y) ydim.setvalues (Y) datainfo.setydimension (Ydim) var= Datainfo.getvariable ('aerosol_optical_thickness_of_mersi_550nm')    Printvar.getname () dimlist=[Ydim, XDim] var.setdimensions (dimlist) GData=Mdi.getgriddata (Var.getname ()) Gdata.yreverse () GData= Gdata.div (10000.0) Gdata.missingvalue=-3.2767#ALS = Legendmanage.createlegendschemefromgriddata (GData, Legendtype.graduatedcolor, Shapetypes.polygon)ALS =legendscheme (Shapetypes.polygon) als.importfromxmlfile (Os.path.join (DataDir,'Aod.lgs')) Alayer= Drawmeteodata.createrasterlayer (GData,"TEST_HDF", ALS) MF=miapp.getmapdocument (). Getactivemapframe () Mf.addlayer (Alayer) mf.movelayer (alayer, 0)Print 'finished!'

Meteoinfo Script Example: Reading fy3a AOD HDF file