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()
时间: 2024-10-18 06:02:03