MeteoInfoLab脚本示例：AOD

MeteoInfo

MODIS的气溶胶光学厚度（AOD）产品应用很广，数据可以在Giovanni上下载：http://disc.sci.gsfc.nasa.gov/giovanni/overview/index.html。有HDF, netCDF和文本格式，这里示例打开netCDF格式数据，提取AOD数据并绘图。


脚本程序如下：

1. f = addfile('D:/Temp/nc/MOD08_D3.A2015121.051.2015122103938.pscs_000500931513.Optical_Depth_Land_And_Ocean_Mean.G3.nc')
   
2. v_aod = f['Optical_Depth_Land_And_Ocean_Mean']
   
3. aod = v_aod['33:41','113:122']
   
4. worldmap()
   
5. china = shaperead('D:/Temp/map/bou2_4p.shp')
   
6. geoshow(china, linecolor='gray')
   
7. world = shaperead('D:/Temp/map/country1.shp')
   
8. geoshow(world)
   
9. levels = arange(0, 1.0, 0.02)
   
10. layer = contourfm(aod, levels, interpolate=True)
    
11. title('Aerosol Optical Depth at 550 nm')
    
12. colorbar(layer)

f = addfile('D:/Temp/nc/MOD08_D3.A2015121.051.2015122103938.pscs_000500931513.Optical_Depth_Land_And_Ocean_Mean.G3.nc')<br /> v_aod = f['Optical_Depth_Land_And_Ocean_Mean']<br /> aod = v_aod['33:41','113:122']<br /> worldmap()<br /> china = shaperead('D:/Temp/map/bou2_4p.shp')<br /> geoshow(china, linecolor='gray')<br /> world = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(world)<br /> levels = arange(0, 1.0, 0.02)<br /> layer = contourfm(aod, levels, interpolate=True)<br /> title('Aerosol Optical Depth at 550 nm')<br /> colorbar(layer)

要和网站上的图形一致的话需要在等值线追踪前将格点数据插值成2倍的分辨率（layer = contourfm(aod, levels, interpolate=True)）。

HDF格式数据文件中的经纬度维有问题，需要在脚本中指定，示例如下：

1. f = addfile('D:/Temp/hdf/MOD08_D3.A2015121.051.2015122103938.pscs_000500931513.Optical_Depth_Land_And_Ocean_Mean.G3.hdf')
   
2. v_aod = f['Optical_Depth_Land_And_Ocean_Mean']
   
3. #Set x/y
   
4. x = linspace(-180.0,180.0, 360)
   
5. y = linspace(-90.0,90.0, 180)
   
6. #Set x/y dimensions
   
7. v_aod.setdim('Y', y)
   
8. v_aod.setdim('X', x)
   
9. #Get data
   
10. aod = v_aod['33:41','113:122']
    
11. #Plot
    
12. axesm()
    
13. china = shaperead('D:/Temp/map/bou2_4p.shp')
    
14. geoshow(china, linecolor='gray')
    
15. world = shaperead('D:/Temp/map/country1.shp')
    
16. geoshow(world)
    
17. levels = arange(0, 1.0, 0.02)
    
18. layer = contourfm(aod, levels, interpolate=True)
    
19. title('Aerosol Optical Depth at 550 nm')
    
20. colorbar(layer)

f = addfile('D:/Temp/hdf/MOD08_D3.A2015121.051.2015122103938.pscs_000500931513.Optical_Depth_Land_And_Ocean_Mean.G3.hdf')<br /> v_aod = f['Optical_Depth_Land_And_Ocean_Mean']<br /> #Set x/y<br /> x = linspace(-180.0,180.0, 360)<br /> y = linspace(-90.0,90.0, 180)<br /> #Set x/y dimensions<br /> v_aod.setdim('Y', y)<br /> v_aod.setdim('X', x)<br /> #Get data<br /> aod = v_aod['33:41','113:122']<br /> #Plot<br /> axesm()<br /> china = shaperead('D:/Temp/map/bou2_4p.shp')<br /> geoshow(china, linecolor='gray')<br /> world = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(world)<br /> levels = arange(0, 1.0, 0.02)<br /> layer = contourfm(aod, levels, interpolate=True)<br /> title('Aerosol Optical Depth at 550 nm')<br /> colorbar(layer)

931404656

不能进入数据网站更新了新版本

stu_data

与Matlab好像啊

houxiaogang

王老师，我运行程序怎么是错误的。

MeteoInfo

houxiaogang 发表于 2015-6-20 00:18
王老师，我运行程序怎么是错误的。

之前的例子是读netCDF格式数据文件的，对于HDF格式文件需要不同的处理，已经在1楼中加入了处理HDF文件的脚本（需要最新版MeteoInfo 1.2.6）。

houxiaogang

MeteoInfo 发表于 2015-6-20 16:46
之前的例子是读netCDF格式数据文件的，对于HDF格式文件需要不同的处理，已经在1楼中加入了处理HDF文件的 ...

王老师，端午节安康！非常感谢您的热心帮助！百度我们新疆不能用，如果您有时间给我邮箱发一发可以吗？再次感谢您！我邮箱：hxg510487306@qq.com

houxiaogang

MeteoInfo 发表于 2015-6-20 16:46
之前的例子是读netCDF格式数据文件的，对于HDF格式文件需要不同的处理，已经在1楼中加入了处理HDF文件的 ...

或者除了百度网盘哪儿可以下载到最新版MeteoInfo 1.2.6？

MeteoInfo

houxiaogang 发表于 2015-6-22 12:40
或者除了百度网盘哪儿可以下载到最新版MeteoInfo 1.2.6？

MeteoInfo网站：http://www.meteothinker.com/

houxiaogang

MeteoInfo 发表于 2015-6-22 17:43
MeteoInfo网站：http://www.meteothinker.com/

非常感谢王老师及时回复，我把MeteoInfo 1.2.6下载下来了，但是我的还是出错了。我跟踪了一下，好像运行到这个语句：
#Get data
aod = v_aod[[33,41],[113,122]]
就出错了，请王老师指点！

MeteoInfo

houxiaogang 发表于 2015-6-22 22:58
非常感谢王老师及时回复，我把MeteoInfo 1.2.6下载下来了，但是我的还是出错了。我跟踪了一下，好像运行 ...

错误信息是什么？把你用的数据文件和脚本都贴出来看看。