Python3.x 爬取历史探空数据并绘图

kmdqlb

楼主，你好，请问 一、如何得到 修改SkewT项目下面的SkewT.py后优化的图像：

二、下载好数据用skewt.py 无法绘图

数据如下：

-----------------------------------------------------------------------------
   PRES   HGHT   TEMP   DWPT   RELH   MIXR   DRCT   SKNT   THTA   THTE   THTV
    hPa     m      C      C      %    g/kg    deg   knot     K      K      K
-----------------------------------------------------------------------------
  814.0   1892   12.6    4.6     58   6.57    225      4  303.1  323.1  304.3
  807.0   1965   14.6    3.6     48   6.18    227      5  305.9  325.0  307.1
  703.0   3110    4.4    0.2     74   5.55    254     15  306.9  324.2  308.0
  700.0   3145    4.0   -0.5     72   5.30    255     16  306.9  323.4  307.9
  673.0   3465    2.6   -9.4     41   2.80    261     25  308.8  317.9  309.3
  655.0   3685    3.4  -24.9     10   0.78    265     31  312.1  314.9  312.2
  651.0   3734    3.6  -28.4      8   0.57    264     32  312.9  314.9  313.0
  586.0   4580   -0.6  -37.5      4   0.26    255     45  317.5  318.5  317.6
  573.0   4761   -1.5  -39.5      4   0.22    254     44  318.5  319.4  318.5
  500.0   5830   -9.9  -41.9      5   0.19    250     43  320.9  321.7  320.9
  452.0   6593  -15.2  -46.1      5   0.14    240     49  323.6  324.2  323.7
  418.0   7184  -19.3  -49.3      5   0.10    228     65  325.7  326.1  325.7
  411.0   7309  -19.9  -41.5     13   0.25    225     68  326.5  327.5  326.5
  400.0   7510  -20.9  -28.9     49   0.88    225     70  327.7  331.0  327.9
  393.0   7640  -22.1  -28.4     57   0.94    225     72  327.8  331.3  328.0
  392.0   7659  -22.3  -28.3     58   0.95    225     72  327.8  331.3  328.0
  362.0   8238  -27.7  -31.0     73   0.80    228     77  328.1  331.1  328.3
  361.0   8258  -26.5  -38.5     31   0.38    228     77  330.0  331.5  330.1
  356.0   8359  -26.1  -44.1     17   0.21    229     78  331.9  332.7  331.9
  300.0   9580  -33.5  -61.5      4   0.03    235     89  338.0  338.2  338.1
  268.0  10363  -35.9  -62.9      5   0.03    241    100  345.6  345.8  345.6
  250.0  10840  -39.7  -65.7      5   0.02    245    107  346.9  347.0  346.9
  221.0  11663  -46.1  -70.5      5   0.01    245    109  349.5  349.5  349.5
  200.0  12330  -51.3  -74.3      5   0.01    240    101  351.4  351.4  351.4
  165.0  13565  -57.9  -79.9      4   0.00    237     85  360.2  360.2  360.2
  150.0  14160  -63.1                         235     78  361.2         361.2
  147.0  14281  -63.7                         235     78  362.2         362.2
  105.0  16288  -73.9                         235     41  379.4         379.4
  100.0  16570  -75.5                         235     41  381.6         381.6
   93.3  16969  -78.5                         232     39  383.3         383.3
   91.2  17100  -77.1                         231     38  388.6         388.6
   88.0  17304  -78.3                         230     37  390.2         390.2
   83.4  17610  -80.1                         235     35  392.6         392.6
   73.0  18362  -81.3                         250     31  405.2         405.2
   70.0  18600  -78.9                         255     29  415.3         415.3
   65.0  19024  -77.3                         259     24  427.6         427.6
   64.0  19115  -75.7                         260     23  433.2         433.2
   59.4  19554  -67.7                         251     18  460.3         460.3
   59.0  19595  -67.7                         250     17  461.1         461.1
   51.0  20471  -68.4                         225     19  479.1         479.1
   50.0  20590  -68.5                         225     21  481.6         481.6
   47.8  20858  -68.9                         229     23  486.9         486.9
   44.0  21364  -64.7                         235     27  508.8         508.8
   41.0  21796  -61.1                         235     23  528.1         528.1
   38.6  22164  -58.1                         233     20  545.0         545.0
   36.0  22603  -57.7                         230     16  557.1         557.1
   30.0  23750  -56.5                         240      8  590.0         590.0
   29.0  23966  -56.4                         250      8  596.0         596.0
   27.0  24421  -56.2                         265      4  609.0         609.0
   26.4  24564  -56.1                         286      2  613.1         613.1
   26.0  24662  -55.9                         300      0  616.5         616.5
   24.0  25177  -54.6                          80      2  634.4         634.4
   21.0  26036  -52.5                          60      6  665.5         665.5
   20.0  26350  -51.7                          60      6  677.2         677.2
   18.0  27035  -50.7                                     701.0         701.0

                         Station identifier: ZPPP
                             Station number: 56778
                           Observation time: 170101/1200
                           Station latitude: 25.01
                          Station longitude: 102.68
                          Station elevation: 1892.0
                               Lifted index: 4.05
    LIFT computed using virtual temperature: 3.67
      Convective Available Potential Energy: 0.00
             CAPE using virtual temperature: 0.00
                      Convective Inhibition: 0.00
             CINS using virtual temperature: 0.00
                     Bulk Richardson Number: 0.00
          Bulk Richardson Number using CAPV: 0.00
  Temp [K] of the Lifted Condensation Level: 274.21
Pres [hPa] of the Lifted Condensation Level: 682.16
     Mean mixed layer potential temperature: 305.90
              Mean mixed layer mixing ratio: 6.11
Precipitable water [mm] for entire sounding: 9.43

zdx-1982

必须赞一个

曦尘

谢谢分享，学习学习！

zdx_1982

No module named 'UserDict'不知道问题出在哪

zdx_1982

求教下，下载的日期格式需要改成什么样的呢，ValueError: time data '-----------------------------------------------------------------------------' does not match format '%HZ %d %b %Y'

简单的世界

跺碎那片云 发表于 2019-3-21 08:28
昨天检查了下，是模块加载出错，可能没装成功

前辈我遇到的问题提示和你一模一样，请问怎么才能正确安装模块

sreveryone

非常感谢这个教程的分享

简单的世界

Traceback (most recent call last):
  File "E:\tankong\SkewT-plot.py", line 7, in <module>
    S=SkewT.Sounding("./56146.txt")
  File "E:\python\lib\site-packages\skewt\SkewT.py", line 232, in __init__
    self.uwyofile(filename)
  File "E:\python\lib\site-packages\skewt\SkewT.py", line 499, in uwyofile
    dstr, "%HZ %d %b %Y").strftime("%Y-%m-%d_%H:%M:%S")
  File "E:\python\lib\_strptime.py", line 577, in _strptime_datetime
    tt, fraction, gmtoff_fraction = _strptime(data_string, format)
  File "E:\python\lib\_strptime.py", line 359, in _strptime
    (data_string, format))
ValueError: time data '-----------------------------------------------------------------------------' does not match format '%HZ %d %b %Y'
解决了安装问题。又出现这个什么鬼{:5_275:}

hzwjy

1. #!/usr/bin/env python
   
2. # coding: utf-8
   
3. 
   
4. from datetime import datetime
   
5. from siphon.simplewebservice.wyoming import WyomingUpperAir
   
6. 
   
7. date = datetime(2019,3,28,12)
   
8. df = WyomingUpperAir.request_data(date, "56187")
   

复制代码

下载数据这一步有更简洁的方法，siphon 获取探空数据
MetPy 出的教程推荐大家看看
MetPy 教程

cz51

不错