MeteoInfoLab脚本示例：读取绘制自动气象站数据

MeteoInfo

国家级和区域级自动气象站数据文件有统一的格式，都是文本文件，文件名类似AWS_WEA_NAT_MUL_FTM-YYYYMMDDHH.TXT，这个提供一个示例脚本来读取文件中所有站点的气象要素数据并绘制图形。

1. fn = r'U:\data\micaps\auto\nat\2016\AWS_WEA_NAT_MUL_FTM-2016120100.TXT'
   
2. colnames = ['year','month','day','hour','stid','lat','lon','alt', \
   
3.     'p_alt','wd_2min','ws_2min','wd_10min','ws_10min','wd_max_10min', \
   
4.     'ws_max_10min', 'time_max_10min','wd','ws','wd_max','ws_max', \
   
5.     'time_max','rain','t','t_max','time_t_max','t_min','time_t_min', \
   
6.     't_max_24h','t_min_24h','rh','rh_min','time_rh_min','press', \
   
7.     'press_max','time_press_max','press_min','time_press_min', \
   
8.     'press_sl']
   
9. table = readtable(fn, format='%4i%s%10f%i%4f%i%3f%i%f%i%4f%i%2f%i%f%i%f', \
   
10.     headerlines=-1, colnames=colnames)
    
11. ws = table['ws_10min']
    
12. ws[ws>100000] = nan
    
13. wd = table['wd_10min']
    
14. wd[wd>100000] = nan
    
15. lon = table['lon']
    
16. lat = table['lat']
    
17. #Plot
    
18. axesm()
    
19. lchina = shaperead('U:/map/bou2_4p.shp')
    
20. geoshow(lchina, edgecolor='k')
    
21. #layer = scatterm(lon, lat, ws, edge=False)
    
22. layer = quiverm(lon, lat, wd, ws, ws, isuv=False, size=10)
    
23. quiverkey(layer, 0.78, 0.18, 10, bbox={'edge':False, 'fill':True})
    
24. colorbar(layer)
    
25. xlim(70,138)
    
26. ylim(15,55)
    
27. title('Wind speed')

fn = r'U:\data\micaps\auto\nat\2016\AWS_WEA_NAT_MUL_FTM-2016120100.TXT'<br /> colnames = ['year','month','day','hour','stid','lat','lon','alt', \<br />     'p_alt','wd_2min','ws_2min','wd_10min','ws_10min','wd_max_10min', \<br />     'ws_max_10min', 'time_max_10min','wd','ws','wd_max','ws_max', \<br />     'time_max','rain','t','t_max','time_t_max','t_min','time_t_min', \<br />     't_max_24h','t_min_24h','rh','rh_min','time_rh_min','press', \<br />     'press_max','time_press_max','press_min','time_press_min', \<br />     'press_sl']<br /> table = readtable(fn, format='%4i%s%10f%i%4f%i%3f%i%f%i%4f%i%2f%i%f%i%f', \<br />     headerlines=-1, colnames=colnames)<br /> ws = table['ws_10min']<br /> ws[ws>100000] = nan<br /> wd = table['wd_10min']<br /> wd[wd>100000] = nan<br /> lon = table['lon']<br /> lat = table['lat']<br /> #Plot<br /> axesm()<br /> lchina = shaperead('U:/map/bou2_4p.shp')<br /> geoshow(lchina, edgecolor='k')<br /> #layer = scatterm(lon, lat, ws, edge=False)<br /> layer = quiverm(lon, lat, wd, ws, ws, isuv=False, size=10)<br /> quiverkey(layer, 0.78, 0.18, 10, bbox={'edge':False, 'fill':True})<br /> colorbar(layer)<br /> xlim(70,138)<br /> ylim(15,55)<br /> title('Wind speed')

topmad

1. colnames = ['year','month','day','hour','stid','lat','lon','alt', \
   
2. 'p_alt','wd_2min','ws_2min','wd_10min','ws_10min','wd_max_10min', \
   
3. 'ws_max_10min', 'time_max_10min','wd','ws','wd_max','ws_max', \
   
4. 'time_max','rain','t','t_max','time_t_max','t_min','time_t_min', \
   
5. 't_max_24h','t_min_24h','rh','rh_min','time_rh_min','press', \
   
6. 'press_max','time_press_max','press_min','time_press_min', \
   
7. 'press_sl']
   

复制代码

王老师
  这里面关于要素的描述，您用的是英文的简称，这个是您自己随便定义的 还是有什么参考的标准呢

如果有 能否指点一二   CIMISS元数据和数据元  是一个坑啊

MeteoInfo

topmad 发表于 2016-12-27 09:44
王老师
  这里面关于要素的描述，您用的是英文的简称，这个是您自己随便定义的 还是有什么参考的标准呢
...

根据数据格式说明里的中文字段名称自己弄成英文简写

topmad

MeteoInfo 发表于 2016-12-27 09:53
根据数据格式说明里的中文字段名称自己弄成英文简写

好吧
看来 还是要找一套标准 实在不行 就只能自己写一个标准了

MeteoInfo

topmad 发表于 2016-12-27 11:07
好吧
看来 还是要找一套标准 实在不行 就只能自己写一个标准了

标准控呀

topmad

MeteoInfo 发表于 2016-12-27 11:15
标准控呀

我也不想控...
关键这个业务要求啊...找了半天  现在弄了个中美结合的 不伦不类  

到时候要是我弄个啥标准，请王老师来评审

无尽的兰博

好东西  谢谢分享

ddd212

学习学习

独孤酒见

王老师，您这个例子中的风向wd数据是0-360的值吗（除了那个大于100000），我的是wd在0-360之间，最后运行不出图啊。

独孤酒见

王老师及各位高人，我想请教两个问题。具体描述见附件。
第一，能不能在读入不同文件时，出来的图的风速色标都共用一个色标。
第二，为什么我换一个输入的文本文件后，我后面的geoshow那个地方运行出错。