MeteoInfoLab脚本示例：AOD

风格而才

王老师您好！我想问一下这个怎么通过MYD04_3k的气溶胶数据画这个图，不知道有没有脚本参考。麻烦王老师了。

MeteoInfo

风格而才 发表于 2018-10-17 09:54
王老师您好！我想问一下这个怎么通过MYD04_3k的气溶胶数据画这个图，不知道有没有脚本参考。麻烦王老师了。

你可以加入MeteoInfo QQ群，提供一个示例数据文件以便写参考脚本。

风格而才

王老师您好！我提供了一个2018年1月1日的MYD04_3K的数据。我想东北地区某一天的AOD空间分布，用1°的数据我会画，但是1°的数据太粗。所以我想用10公里或3公里的数据画一下，这里面涉及到数据的预处理，所以想麻烦王老师问一下meteoinfo能否利用这个数据直接用脚本画这种空间分布图图。

MeteoInfo

风格而才 发表于 2018-10-17 15:33
王老师您好！我提供了一个2018年1月1日的MYD04_3K的数据。我想东北地区某一天的AOD空间分布，用1°的数据我 ...

参考这个：

1. fn = 'C:/Temp/aod/MYD04_3K.A2018001.0450.061.2018003204236.hdf'
   
2. f = addfile(fn)
   
3. lon = f['Longitude'][:,:]
   
4. lat = f['Latitude'][:,:]
   
5. aod = f['Optical_Depth_Land_And_Ocean'][:,:]
   
6. aod[aod<0] = nan
   
7. 
   
8. #Plot
   
9. axesm()
   
10. geoshow('cn_province')
    
11. geoshow('country', edgecolor='k')
    
12. levs = arange(0, 2, 0.1)
    
13. layer = pcolorm(lon, lat, aod, levs)
    
14. colorbar(layer, orientation='horizontal', aspect=40)
    
15. title('AOD')

fn = 'C:/Temp/aod/MYD04_3K.A2018001.0450.061.2018003204236.hdf'<br /> f = addfile(fn)<br /> lon = f['Longitude'][:,:]<br /> lat = f['Latitude'][:,:]<br /> aod = f['Optical_Depth_Land_And_Ocean'][:,:]<br /> aod[aod<0] = nan<br /> <br /> #Plot<br /> axesm()<br /> geoshow('cn_province')<br /> geoshow('country', edgecolor='k')<br /> levs = arange(0, 2, 0.1)<br /> layer = pcolorm(lon, lat, aod, levs)<br /> colorbar(layer, orientation='horizontal', aspect=40)<br /> title('AOD')

风格而才

MeteoInfo 发表于 2018-10-17 18:52
参考这个：

多谢老师！ 我再问一下有没有办法不对数据做几何校正和拼接，然后把这一天的几个数据通过脚本画一个日平均图。

tianxiao

牛，值得学习。。。

tobewealthy

感谢王老师！

tobewealthy

请问王老师可不可以将处理后的AOD数据输出成excel或其他格式的文件，输出该用什么语句

MeteoInfo

tobewealthy 发表于 2020-3-8 10:42
请问王老师可不可以将处理后的AOD数据输出成excel或其他格式的文件，输出该用什么语句

参考这里：http://www.meteothink.org/docs/meteoinfolab/dataset/midata.html

wyytwyyt

王老师，您好，请问一下，为什么画出来那么多小圆点~麻烦王老师啦~