python处理探空

lovechang1314

今天需要处理探空资料，看了半天发现不能直接读取，写了代码分享一下。

1. data = pd.read_csv(filein, delim_whitespace=True, skiprows=2, names=[
   
2.                    'hpa', 'high', 't', 'td', 'wd', 'ws'])
   
3. print(filein, data is None)
   
4. wsna = data.index[data.ws.isnull()]
   
5. sta = data.hpa[wsna].values
   
6. alldata = []
   
7. for inum, num in enumerate(wsna):
   
8.     if sta[inum] in stas:
   
9.         if inum != len(wsna) - 1:
   
10.             tmpdata = data[num + 1:wsna[inum + 1]]
    
11.             tmpdata['sta'] = sta[inum]
    
12.             alldata.append(tmpdata)
    
13.         else:
    
14.             tmpdata = data[num + 1:]
    
15.             tmpdata['sta'] = sta[inum]
    
16.             alldata.append(tmpdata)
    
17. alldata = pd.concat(alldata)
    
18. # print(filein.split('.')[0][-8:])
    
19. alldata['date'] = pd.to_datetime(
    
20.     filein.split('.')[0][-8:], format='%y%m%d%H')

复制代码

得到结果如下：

weinihou

楼主是用pandas来搞啊，有个包叫  SkewT ，可以在github或者pip上找到，专门处理探空sounding的
比如可以在这里找找https://github.com/tjlang/SkewT

自含绘图和一些简单计算

然后metpy没用过，借楼主推荐，学习下~！
谢谢分享~~~

fjqxyyj

很给力，学习了

lgz

很给力，学习了

schliezer

可以出skew-plot图吗？

lovechang1314

schliezer 发表于 2017-4-17 11:53
可以出skew-plot图吗？

直接调用metpy可以出

schliezer

lovechang1314 发表于 2017-4-17 13:28
直接调用metpy可以出

好强大.......

zq200866

楼主，能把数据和完整代码分享出来吗？

chongzika

能不能给个简单的例子文件，让大家实际操作下

lovechang1314

zq200866 发表于 2017-4-20 20:58
楼主，能把数据和完整代码分享出来吗？

随便找一个micaps的探空文件测试就行啦，我的是做别的用，对你可能没啥大用，所以只贴了读取