MODIS的气溶胶光学厚度（AOD）产品应用很广，数据可以在Giovanni上下载：http://disc.sci.gsfc.nasa.gov/giovanni/overview/index.html。有HDF, netCDF和文本格式，这里示例打开netCDF格式数据，提取AOD数据并绘图。
脚本程序如下：
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()
原文：http://www.cnblogs.com/yaqiang/p/4585754.html