MeteoInfoLab脚本示例：吸收性气溶胶指数

MeteoInfo

OMI卫星吸收性气溶胶指数（Absorbing Aerosol Index）产品数据在此网页下载：http://www.temis.nl/airpollution/absaai/#OMI_AAI 。数据是文本文件，读取数据并绘图的示例脚本如下：

1. #Read data from file
   
2. tt = datetime.datetime(2016,4,6)
   
3. fn = 'W:/SDS_Asian/SDS_obs/aai_omi/' + tt.strftime('%y%m%d') + '.esr'
   
4. f = open(fn)
   
5. for i in range(3):
   
6.     f.readline()
   
7. data = []
   
8. for line in f:
   
9.     line = line.strip('\n').split(' lat')[0][1:]
   
10.     #print line
    
11.     for m in range(0, len(line), 3):
    
12.         d = line[m:m+3]
    
13.         data.append(float(d))
    
14. f.close()
    
15. #Reshape data
    
16. xn = 288
    
17. yn = 180
    
18. data = array(data)
    
19. data = data.reshape([yn,xn])
    
20. data[data>=996.0] = nan
    
21. data = (data - 450) * 0.1
    
22. #Set lon/lat
    
23. lon = arange1(-179.375, xn, 1.25)
    
24. lat = arange1(-89.5, yn, 1.0)
    
25. #Plot
    
26. axesm()
    
27. lworld = shaperead('U:/map/country1.shp')
    
28. geoshow(lworld, edgecolor=[50,50,50])
    
29. levs = arange(0, 3.1, 0.1)
    
30. layer = imshowm(lon, lat, data, levs, cmap='wh-bl-gr-ye-re', fill_color='lightgray')
    
31. colorbar(layer)
    
32. title('Absorbing Aerosol Index (' + tt.strftime('%Y-%m-%d') + ')')

#Read data from file<br /> tt = datetime.datetime(2016,4,6)<br /> fn = 'W:/SDS_Asian/SDS_obs/aai_omi/' + tt.strftime('%y%m%d') + '.esr'<br /> f = open(fn)<br /> for i in range(3):<br />     f.readline()<br /> data = []<br /> for line in f:<br />     line = line.strip('\n').split(' lat')[0][1:]<br />     #print line<br />     for m in range(0, len(line), 3):<br />         d = line[m:m+3]<br />         data.append(float(d))<br /> f.close()<br /> #Reshape data<br /> xn = 288<br /> yn = 180<br /> data = array(data)<br /> data = data.reshape([yn,xn])<br /> data[data>=996.0] = nan<br /> data = (data - 450) * 0.1<br /> #Set lon/lat<br /> lon = arange1(-179.375, xn, 1.25)<br /> lat = arange1(-89.5, yn, 1.0)<br /> #Plot<br /> axesm()<br /> lworld = shaperead('U:/map/country1.shp')<br /> geoshow(lworld, edgecolor=[50,50,50])<br /> levs = arange(0, 3.1, 0.1)<br /> layer = imshowm(lon, lat, data, levs, cmap='wh-bl-gr-ye-re', fill_color='lightgray')<br /> colorbar(layer)<br /> title('Absorbing Aerosol Index (' + tt.strftime('%Y-%m-%d') + ')')

突破的渴望

王老师您好，这个是meteoinfo的代码？ 具体怎么用啊？谢谢王老师，为什么我的 meteoinfo 打开会有这个错误呢？  UnicodeDecodeError: 'unicodeescape' codec can't decode bytes in position 8-10: truncated \xXX
打搅王老师了~~

MeteoInfo

突破的渴望 发表于 2017-5-11 11:21
王老师您好，这个是meteoinfo的代码？ 具体怎么用啊？谢谢王老师，为什么我的 meteoinfo 打开会有这个错误 ...

参考这个帖子：http://bbs.06climate.com/forum.p ... &extra=page%3D1

突破的渴望

王老师您好，为什么运行这个出来这个问题呢？ 这个代码 还需要修改什么吧？下载的文件 是TXT的 还是nc的啊？谢谢王老师~~

MeteoInfo

突破的渴望 发表于 2017-5-11 20:41
王老师您好，为什么运行这个出来这个问题呢？这个代码 还需要修改什么吧？下载的文件 是TXT的 还是nc的啊？ ...

看错误信息呀，文件不存在

突破的渴望

王老师您好 下载哪个文件呢？ 这个之前 没用过， 里面的文件下载下来了 但是好像 还是不对，是哪个呢？ 谢谢王老师~~~

MeteoInfo

突破的渴望 发表于 2017-5-11 20:58
王老师您好 下载哪个文件呢？这个之前 没用过， 里面的文件下载下来了 但是好像 还是不对，是哪个呢？ 谢谢 ...

“WARNING: We are moving from ASCII to NetCDF. For the level-3 data, this means that the gridded daily and monthly data are now only available in NetCDF format. The ASCII files are not produced any more.”

数据提供方修改数据格式了，现在是netCDF格式。读起来应该更简单，这里有很多例子可以参考：http://bbs.06climate.com/forum.p ... &extra=page%3D1