露点、温度和相对湿度

MeteoInfo

MICAPS第1类数据中有温度和露点，但没有相对湿度，需要计算获得。我不是学气象的，在网上查了一些资料，对这三者的关系有了简单的认识。在MeteoInfo类库的org.meteoinfo.data.analysis包中增加了一个MeteoMath类，希望能将气象相关算法陆续加入进去，目前实现了计算饱和水汽压（cal_Es(double tc)方法）、计算相对湿度（cal_RH(double t, double tdc）方法）、计算露点（cal_Tdc(double t, double rh）方法）。主要参考了“西风江南”提供的算法（http://bbs.06climate.com/forum.php?mod=viewthread&tid=15263&highlight=%C2%B6%B5%E3%2B%CF%E0%B6%D4%CA%AA%B6%C8），另外此网页有更简洁的算法（http://www.gorhamschaffler.com/humidity_formulas.htm），如果能有高手对这些算法进行一些讲解就更好了。和在线计算器的结果（http://go.vaisala.com/humiditycalculator/vaisala_humidity_calculator.html）进行了对比，结果很接近。需要用到最新的MeteoInfo Java 1.1.9R4版（http://www.meteothinker.com/Downloads.html）。

这是一个示例脚本：

1. from org.meteoinfo.data.analysis import MeteoMath
   
2. 
   
3. t = 20
   
4. #dt = 10
   
5. rh = 50
   
6. #rh = MeteoMath.cal_RH(t, dt)
   
7. dt = MeteoMath.cal_Tdc(t, rh)
   
8. print 'RH: ' + str(rh)
   
9. print 'Dewpoint: ' + str(dt)

from org.meteoinfo.data.analysis import MeteoMath<br /> <br /> t = 20<br /> #dt = 10<br /> rh = 50<br /> #rh = MeteoMath.cal_RH(t, dt)<br /> dt = MeteoMath.cal_Tdc(t, rh)<br /> print 'RH: ' + str(rh)<br /> print 'Dewpoint: ' + str(dt)

从MICAPS第1类数据中提取温度和露点，然后计算相对湿度的脚本：

1. #--------------------------------------------------------        
   
2. # Author: Yaqiang Wang                                          
   
3. # Date: 2015-1-14                                               
   
4. # Purpose: Extract station data from micaps1 data and calculate
   
5. #          relative humidity from temperature and dewpoint
   
6. # Note: Sample                                                   
   
7. #-----------------------------------------------------------
   
8. from org.meteoinfo.data.meteodata import MeteoDataInfo
   
9. from org.meteoinfo.table import DataTable, DataTypes
   
10. from org.meteoinfo.data.analysis import MeteoMath
    
11. import os
    
12. import datetime
    
13. 
    
14. #---- Set data directories
    
15. print 'Start script...'
    
16. mdir = 'U:/data/micaps/2014/plot'
    
17. outdir = 'D:/Temp/test'
    
18. 
    
19. #---- Set stations
    
20. stations = ['54511', '58362']
    
21. 
    
22. #---- Set variables
    
23. # Avialable variables in MICAPS 1 data:
    
24. #  "Altitude", "Grade", "CloudCover", "WindDirection", "WindSpeed", "Pressure",
    
25. #  "PressVar3h", "WeatherPast1", "WeatherPast2", "Precipitation6h", "LowCloudShape",
    
26. #  "LowCloudAmount", "LowCloudHeight", "DewPoint", "Visibility", "WeatherNow",
    
27. #  "Temperature", "MiddleCloudShape", "HighCloudShape"
    
28. variables = ['Temperature','DewPoint']
    
29. 
    
30. #---- Set start/end times
    
31. stime = datetime.datetime.now()
    
32. sdate = datetime.datetime(2014, 10, 1, 2)
    
33. edate = datetime.datetime(2014, 10, 2, 3)
    
34. #edate = datetime.datetime(2014, 12, 31, 24)
    
35. 
    
36. #---- Extract data
    
37. mdi = MeteoDataInfo()
    
38. rTable = DataTable()
    
39. n = 0
    
40. while sdate < edate:
    
41.     mfn = os.path.join(mdir, sdate.strftime('%y%m%d%H') + '.000')
    
42.     if os.path.exists(mfn):
    
43.         print mfn
    
44.         mdi.openMICAPSData(mfn)
    
45.         stdata = mdi.getStationInfoData()
    
46.         t = mdi.getDataInfo().getTime(0)
    
47.         result = stdata.extractStationData(stations, variables, t)
    
48.         result.addColumn('RH', DataTypes.String)
    
49.         for row in result.getRows():
    
50.             temp = float(row.getValue('Temperature'))
    
51.             dp = float(row.getValue('DewPoint'))
    
52.             rh = MeteoMath.cal_RH(temp, dp)
    
53.             row.setValue('RH', '%.2f'%rh)
    
54.         print result.toString()
    
55.         if n == 0:
    
56.             rTable = result
    
57.         else:
    
58.             for row in result.getRows():
    
59.                 rTable.appendRow(row)
    
60.         n += 1
    
61. 
    
62.     sdate = sdate + datetime.timedelta(hours=3)
    
63. 
    
64. #---- Output result
    
65. print '**********************************'
    
66. #print rTable.toString()
    
67. #outfn = os.path.join(outdir, 'stdata_test.csv')
    
68. #rTable.saveAsCSVFile(outfn)
    
69. #outfn = os.path.join(outdir, 'stdata_test.txt')
    
70. #rTable.saveAsASCIIFile(outfn)
    
71. 
    
72. print 'Run time: %s'%(datetime.datetime.now() - stime)
    
73. print 'Finished!'

#--------------------------------------------------------        <br /> # Author: Yaqiang Wang                                           <br /> # Date: 2015-1-14                                               <br /> # Purpose: Extract station data from micaps1 data and calculate<br /> #          relative humidity from temperature and dewpoint<br /> # Note: Sample                                                   <br /> #-----------------------------------------------------------<br /> from org.meteoinfo.data.meteodata import MeteoDataInfo<br /> from org.meteoinfo.table import DataTable, DataTypes<br /> from org.meteoinfo.data.analysis import MeteoMath<br /> import os<br /> import datetime<br /> <br /> #---- Set data directories<br /> print 'Start script...'<br /> mdir = 'U:/data/micaps/2014/plot'<br /> outdir = 'D:/Temp/test'<br /> <br /> #---- Set stations<br /> stations = ['54511', '58362']<br /> <br /> #---- Set variables<br /> # Avialable variables in MICAPS 1 data: <br /> #  "Altitude", "Grade", "CloudCover", "WindDirection", "WindSpeed", "Pressure",<br /> #  "PressVar3h", "WeatherPast1", "WeatherPast2", "Precipitation6h", "LowCloudShape",<br /> #  "LowCloudAmount", "LowCloudHeight", "DewPoint", "Visibility", "WeatherNow",<br /> #  "Temperature", "MiddleCloudShape", "HighCloudShape"<br /> variables = ['Temperature','DewPoint']<br /> <br /> #---- Set start/end times<br /> stime = datetime.datetime.now()<br /> sdate = datetime.datetime(2014, 10, 1, 2)<br /> edate = datetime.datetime(2014, 10, 2, 3)<br /> #edate = datetime.datetime(2014, 12, 31, 24)<br /> <br /> #---- Extract data<br /> mdi = MeteoDataInfo()<br /> rTable = DataTable()<br /> n = 0<br /> while sdate < edate:<br />     mfn = os.path.join(mdir, sdate.strftime('%y%m%d%H') + '.000')<br />     if os.path.exists(mfn):<br />         print mfn<br />         mdi.openMICAPSData(mfn)<br />         stdata = mdi.getStationInfoData()<br />         t = mdi.getDataInfo().getTime(0)<br />         result = stdata.extractStationData(stations, variables, t)<br />         result.addColumn('RH', DataTypes.String)<br />         for row in result.getRows():<br />             temp = float(row.getValue('Temperature'))<br />             dp = float(row.getValue('DewPoint'))<br />             rh = MeteoMath.cal_RH(temp, dp)<br />             row.setValue('RH', '%.2f'%rh)<br />         print result.toString()<br />         if n == 0:<br />             rTable = result<br />         else:<br />             for row in result.getRows():<br />                 rTable.appendRow(row)<br />         n += 1<br /> <br />     sdate = sdate + datetime.timedelta(hours=3)<br /> <br /> #---- Output result<br /> print '**********************************'<br /> #print rTable.toString()<br /> #outfn = os.path.join(outdir, 'stdata_test.csv')<br /> #rTable.saveAsCSVFile(outfn)<br /> #outfn = os.path.join(outdir, 'stdata_test.txt')<br /> #rTable.saveAsASCIIFile(outfn)<br /> <br /> print 'Run time: %s'%(datetime.datetime.now() - stime)<br /> print 'Finished!'

二爷名声在外

王老师太给力了

sysun

功能越来越强大了

风行浪尖

这个给力，我用SQL、C#和MacthCAD也写过

风行浪尖

这个给力，我用SQL、C#和MacthCAD也写过

elegisite

大神们在前面各种给力，渣渣们表示各种福利哇。。。

張顧.煒Zane

先收藏一下。。以后慢慢看

1032610020

先收藏le ...

梦纯

王老师，温度和露点来求相对湿度的思路是啥啊？能否告知一下

MeteoInfo

梦纯 发表于 2015-10-21 17:43
王老师，温度和露点来求相对湿度的思路是啥啊？能否告知一下

记不住了，自己找相关资料看吧