Meteoinfolab Script Example: fy-3a AOD HDF Data

The FY3A satellite has AOD product data in HDF format, with a global range of 0.05 degrees resolution. To reset the X/y dimension after reading the data file variable, the data is reversed on the Y axis and needs to be divided by 10000 to get the AOD value.

Script Program:

#Add Data Filefn ='D:/TEMP/HDF/FY3A_MERSI_GBAL_L3_ASL_MLT_GLL_20140331_AOAM_5000M_MS.HDF'F=AddFile (FN)#Get Data Variablev = f['aerosol_optical_thickness_of_mersi_550nm']#Set x/yx = Arange1 ( -30.0,7200, 0.05) y= Arange1 ( -90.0,3600, 0.05)#Set x/y dimensionsV.setdim ('Y', y, 0) V.setdim ('X', X, 1)#Get Data ArrayAOD =v[:,:]#AOD = v[[15, 55],[70, []]AOD = aod[::-1,:]/10000AOD= aod[[15, 55],[70, 140]]#Plotaxesm () World= Shaperead ('d:/temp/map/country1.shp') China= Shaperead ('d:/temp/map/bou2_4p.shp') Geoshow (China, Edgecolor='Gray') geoshow (world) Levs= Arange1 (0, 20, 0.05) cols= Makecolors (cmap='Matlab_jet') Cols[0]='W'Layer= IMSHOWM (AOD, Levs, colors=cols) Colorbar (layer) title ('fy-3a AOD (2014-03-31)')

Meteoinfolab Script Example: fy-3a AOD HDF Data