MeteoInfo脚本示例：从MICAPS第一类数据中提取数据

MeteoInfo

MICAPS第一类数据是地面全要素观测数据，每个文件一个时次，包含很多站点和要素。这个脚本程序是要解决从多个MICAPS第一类数据文件中提取某些站点、某些要素的数据，并保存为文本文件。应该对业务科研工作比较有用。

示例脚本程序如下（需要MeteoInfo Java最新版本1.1.9R3的支持）：

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

#--------------------------------------------------------        <br /> # Author: Yaqiang Wang                                           <br /> # Date: 2015-1-14                                               <br /> # Purpose: Extract station data from micaps1 data <br /> # Note: Sample                                                   <br /> #-----------------------------------------------------------<br /> from org.meteoinfo.data.meteodata import MeteoDataInfo<br /> from org.meteoinfo.table import DataTable<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 = ['WindDirection','WindSpeed']<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, 10, 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 />         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!'

elegisite

对于编程一窍不通的，这就是捷径哇。。感谢大神。。。

王楠_20120014

顶起来顶起来！！王老师辛苦咯

香草拿铁

今天刚接到任务需要这个，太感谢了！

香草拿铁

试了一下，搞定！
非常感谢分享！！！！

无公害胖子

你好
不知道有没有能提取分钟级别micaps第一类数据的脚本呢？

悠然一梦123

香草拿铁 发表于 2015-4-16 17:47
今天刚接到任务需要这个，太感谢了！

我的没有成功，出现了这种错误，是什么原因>>> run script...
  File "<iostream>", line 39
    &Acirc;&nbsp; &Acirc;&nbsp; mfn = os.path.join(mdir, sdate.strftime('%y%m%d%H') + '.000')
    ^
SyntaxError: no viable alternative at character '&nbsp;'
>>>

悠然一梦123

老师，我试着用这个脚本提取数据，出现了这样的提示，是为什么>>> run script...
Start script...
D:\MeteoInfo\sample\MICAPS\plot\14100102.000
Traceback (most recent call last):
  File "<iostream>", line 46, in <module>
        at java.text.SimpleDateFormat.<init>(Unknown Source)
        at java.text.SimpleDateFormat.<init>(Unknown Source)
        at org.meteoinfo.table.DataTable.toString(DataTable.java:749)
        at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
        at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)
        at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
        at java.lang.reflect.Method.invoke(Unknown Source)

java.lang.NullPointerException: java.lang.NullPointerException
>>>

MeteoInfo

悠然一梦123 发表于 2015-12-19 18:03
老师，我试着用这个脚本提取数据，出现了这样的提示，是为什么>>> run script...
Start script...
D:\Met ...

这是软件的bug，你重新下载最新版本的MeteoInfo 1.2.9R1试试：http://www.meteothinker.com/

ayzqs

今天用上了，很给力