基于Siphon和MetPy获取探空数据并绘图

hzwjy

Siphon（链接：https://unidata.github.io/siphon/latest/index.html）
Siphon 可用于获取THREDDS数据服务提供的多种观测数据，也可以获取怀俄明大学的探空数据（http://weather.uwyo.edu/upperair/seasia.html）和IGRA2探空数据集。返回的结果为pandas的DataFrame类型，变量选取和处理都比较方便。
可用 conda 安装：

1. conda install -c conda-forge siphon

复制代码

MetPy （链接：https://unidata.github.io/MetPy/latest/index.html）
MetPy 可以用于气象数据的读取、变量计算和数据可视化，在站点图和2维图像绘制上都有独特的优势。
conda 安装命令：

1. conda install -c conda-forge metpy

复制代码

获取探空数据和绘制探空图的py脚本：

1. import numpy as np
   
2. import matplotlib.pyplot as plt
   
3. import metpy.calc as mpcalc
   
4. from metpy.plots import SkewT
   
5. from metpy.units import units
   
6. from datetime import datetime
   
7. from siphon.simplewebservice.wyoming import WyomingUpperAir
   
8. 
   
9. date = datetime(2017, 8, 15, 0)
   
10. site = '56187'
    
11. df = WyomingUpperAir.request_data(date, site)
    
12. 
    
13. p = df['pressure'].values * units(df.units['pressure'])
    
14. T = df['temperature'].values * units(df.units['temperature'])
    
15. Td = df['dewpoint'].values * units(df.units['dewpoint'])
    
16. u = df['u_wind'].values * units(df.units['u_wind'])
    
17. v = df['v_wind'].values * units(df.units['v_wind'])
    
18. 
    
19. # 避免风向标过于密集
    
20. interval = np.logspace(2, 3, num=25) * units('hPa')
    
21. idx = mpcalc.resample_nn_1d(p, interval)
    
22. 
    
23. fig = plt.figure(figsize=(9, 9))
    
24. 
    
25. skew = SkewT(fig, rotation=30)
    
26. 
    
27. skew.ax.set_ylim(1000, 100)
    
28. skew.ax.set_xlim(-40, 40)
    
29. 
    
30. skew.plot(p, T, 'r')
    
31. skew.plot(p, Td, 'g')
    
32. skew.plot_barbs(p[idx], u[idx], v[idx])
    
33. 
    
34. skew.plot_dry_adiabats()
    
35. skew.plot_moist_adiabats()
    
36. skew.plot_mixing_lines()

复制代码

得到的结果

探空绘图

qxjzkx@163.com

提示如下错误：
Warning (from warnings module):
  File "D:\Python37\lib\site-packages\pint\quantity.py", line 1377
    warnings.warn("The unit of the quantity is stripped.", UnitStrippedWarning)
UnitStrippedWarning: The unit of the quantity is stripped.
pint安装好了啊。

风往北吹

qxjzkx@163.com 发表于 2019-4-16 17:47
提示如下错误：
Warning (from warnings module):
  File "D:\Python37\lib\site-packages\pint\quantity ...

不是报错，只是警告！ 解决办法：https://github.com/dopplershift/ ... aed093126b975d84830

931404656

速度单位跟国标不一致哦  要更正下哈

谢斌

速度单位跟国标不一致哦  要更正下哈

大罗宾

请问为什么读出来没有相对湿度呢

cibaouc

No module named 'siphon'

林头儿

plt.show()

arctichare

好好好！现在网站直接点不出来了，用这个方法还可以下！

beiruoyiqiao

安装siphon依赖protobuf库，好像很复杂还没搞明白怎么安装protobuf库