Meteoinfolab Script Example: AIRS Swath HDF Data

Examples of airs Swath HDF data in polar stereographic (South Pole) projection close to the rectangle, you need to read from the data latitude and longitude and related data array, using the Surfacem function to draw Swath data (scatter), In the Surfacem function, the latitude data is projected into the destination projected coordinate system (defined in the AXESM function), and surface data (that is, lattice data) is generated from the 2-D scatter data using the so-called surface interpolation method. 2-D Scatter data can be made into grids (non-rectangular grids), traversing each lattice point of the target grid point data, and the grid's scatter value is assigned to that grid point, which falls in which grid in the grid. Some of the awkward, in fact, the algorithm is not very complex, and MATLAB in the function of the SURFACEM function similar.

Script Program:

#Add Data Filefolder ='d:/temp/hdf/'FNS='AIRS.2002.12.31.001.L2.CC_H.V4.0.21.0.G06100185050.HDF'fn= Folder +Fnsf=AddFile (FN) Lon_v= f['Longitude']lat_v= f['Latitude']lon=Lon_v[:,:]lat=Lat_v[:,:]vname='radiances'Rad_v=F[vname]rad= rad_v[:,:,567]#PlotPlot,proj = Axesm (proj='Stere', lat_0=-90, Gridline=true, griddx=30, griddy=30) Mlayer= Shaperead ('d:/temp/map/country1.shp') geoshow (mlayer) Levs= Arange (40, 90, 1)#slayer = Scatterm (lon, LAT, Rad, Levs, edge=true)Layer =Surfacem (Lon, LAT, Rad, levs) colorbar (layer, orientation='Horizontal') title ('{0}\n {1}'. Format (FNS, VName)) axism ()

Meteoinfolab Script Example: AIRS Swath HDF Data