MeteoInfoLab脚本示例：AOD

houxiaogang

MeteoInfo 发表于 2015-6-22 23:30
错误信息是什么？把你用的数据文件和脚本都贴出来看看。

王老师回复这么快，我贴出来了，谢谢您！

MeteoInfo

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

你的数据不是全球的，需要改x, y的范围：

1. f = addfile('D:/Temp/hdf/MOD08_D3.A20150501.hdf')
   
2. v_aod = f['Optical_Depth_Land_And_Ocean_Mean']
   
3. #Set x/y
   
4. x = linspace(112.0,123.0, 11)
   
5. y = linspace(32.0,42.0, 10)
   
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. aod = v_aod[:,:]
    
12. #Plot
    
13. axesm()
    
14. china = shaperead('D:/Temp/map/bou2_4p.shp')
    
15. geoshow(china, linecolor='gray')
    
16. world = shaperead('D:/Temp/map/country1.shp')
    
17. geoshow(world)
    
18. levels = arange(0, 1.0, 0.02)
    
19. layer = contourfm(aod, levels, interpolate=True)
    
20. title('Aerosol Optical Depth at 550 nm')
    
21. colorbar(layer)
    
22. show()

f = addfile('D:/Temp/hdf/MOD08_D3.A20150501.hdf')<br /> v_aod = f['Optical_Depth_Land_And_Ocean_Mean']<br /> #Set x/y<br /> x = linspace(112.0,123.0, 11)<br /> y = linspace(32.0,42.0, 10)<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 /> aod = v_aod[:,:]<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)<br /> show()

houxiaogang

MeteoInfo 发表于 2015-6-22 23:47
你的数据不是全球的，需要改x, y的范围：

王老师太牛了！结果出来了，非常非常感谢！谢谢王老师这么晚回复我！我还有如下几个问题：
1、这个图怎么存成bmp（或其他格式）图；
2、怎么把经纬网格加上去（如图例）；
3、Colorbar怎么可以设置成如图例一样。
王老师早点休息，这个几个问题可以以后再帮助我。

houxiaogang

MeteoInfo真的很牛，一个小小的建议就是错误提示如果能比较清楚点就很好了！

MeteoInfo

houxiaogang 发表于 2015-6-22 23:55
王老师太牛了！结果出来了，非常非常感谢！谢谢王老师这么晚回复我！我还有如下几个问题：
1、这个图怎 ...

1、输出图片用savefig函数，参加此帖：http://bbs.06climate.com/forum.php?mod=viewthread&tid=34995

2、加入grid(True)语句，目前网格线是画在最底层，以后会加上画在最顶层的选项。

3、重新定义levels和cmap，cmap的取值可以看MeteoInfo路径中colormaps目录中的文件名，不用写后缀。

1. f = addfile('D:/Temp/hdf/MOD08_D3.A20150501.hdf')
   
2. v_aod = f['Optical_Depth_Land_And_Ocean_Mean']
   
3. #Set x/y
   
4. x = linspace(112.0,123.0, 11)
   
5. y = linspace(32.0,42.0, 10)
   
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. aod = v_aod[:,:]
    
12. #Plot
    
13. axesm()
    
14. china = shaperead('D:/Temp/map/bou2_4p.shp')
    
15. geoshow(china, linecolor='gray')
    
16. world = shaperead('D:/Temp/map/country1.shp')
    
17. geoshow(world)
    
18. #levels = arange(0, 1.0, 0.02)
    
19. levels = arange(0.1, 0.98, 0.08)
    
20. layer = contourfm(aod, levels, cmap='grads_rainbow', interpolate=True)
    
21. title('Aerosol Optical Depth at 550 nm')
    
22. grid(True)
    
23. colorbar(layer)
    
24. show()
    
25. savefig('D:/Temp/test/aod.png', 500, 600)

f = addfile('D:/Temp/hdf/MOD08_D3.A20150501.hdf')<br /> v_aod = f['Optical_Depth_Land_And_Ocean_Mean']<br /> #Set x/y<br /> x = linspace(112.0,123.0, 11)<br /> y = linspace(32.0,42.0, 10)<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 /> aod = v_aod[:,:]<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 /> levels = arange(0.1, 0.98, 0.08) <br /> layer = contourfm(aod, levels, cmap='grads_rainbow', interpolate=True)<br /> title('Aerosol Optical Depth at 550 nm')<br /> grid(True)<br /> colorbar(layer)<br /> show()<br /> savefig('D:/Temp/test/aod.png', 500, 600)

houxiaogang

MeteoInfo 发表于 2015-6-23 09:32
1、输出图片用savefig函数，参加此帖：http://bbs.06climate.com/forum.php?mod=viewthread&tid=34995

...

谢谢王老师的指导与帮助！结果出来了，以后好好学学这个软件，我也多推广！

MeteoInfo

houxiaogang 发表于 2015-6-23 10:45
谢谢王老师的指导与帮助！结果出来了，以后好好学学这个软件，我也多推广！

嗯，软件有用户才能有生命力。MeteoInfoLab还在很初级的阶段，也希望用户能提出功能需求并汇报bug。

z黑色青春

王老师，如果用MATLAB处理这个数据呢，哪个更简单些

MeteoInfo

z黑色青春 发表于 2015-7-13 15:45
王老师，如果用MATLAB处理这个数据呢，哪个更简单些

你可以自己尝试一下，给出结果。

z黑色青春

z黑色青春 发表于 2015-7-13 15:45
王老师，如果用MATLAB处理这个数据呢，哪个更简单些

好的 谢谢您