Meteoinfolab Script Example: MODIS AOD

Modis Aerosol Optical thickness (AOD) products are widely used, and data can be downloaded on Giovanni: http://disc.sci.gsfc.nasa.gov/giovanni/overview/index.html. There are HDF, netcdf and text formats, where examples open NETCDF format data, extract AOD data and draw.

The script is as follows:

f = AddFile ('d:/temp/nc/mod08_d3. A2015121.051.2015122103938.PSCS_000500931513.OPTICAL_DEPTH_LAND_AND_OCEAN_MEAN.G3.NC') V_aod= f['Optical_depth_land_and_ocean_mean']aod= v_aod[[33,41],[113,122]]axesm () China= Shaperead ('d:/temp/map/bou2_4p.shp') Geoshow (China, LineColor='Gray') World= Shaperead ('d:/temp/map/country1.shp') Geoshow (world) levels= Arange (0, 1.0, 0.02) Layer= CONTOURFM (AOD, Levels, interpolate=True) title ('aerosol Optical Depth at 550 nm') Colorbar (layer) show ()

Meteoinfolab Script Example: MODIS AOD