MeteoInfoLab脚本示例：从micaps 1数据中提取站点数据

MeteoInfo

micaps 1类数据是地面全要素观测数据，每个时次一个文件，每个文件包含很多站点。要提取某个站点的时间序列数据需要用循环读取多个时次的micaps 1数据，循环中用addfile_micaps函数读取micaps文件对象，然后用文件对象的read_dataframe函数将数据读入一个DataFrame对象，再从中提取一个或者多个站点的数据（也可以对提取那些要素进行限定），用DataFrame的append方法合并不同时次的数据。最终的DataFrame对象用to_csv方法保存到文件中。

1. #---- Set stations
   
2. stids = ['54511']
   
3. 
   
4. #---- Set variables
   
5. # Avialable variables in MICAPS 1 data:
   
6. #  "Altitude", "Grade", "CloudCover", "WindDirection", "WindSpeed", "Pressure",
   
7. #  "PressVar3h", "WeatherPast1", "WeatherPast2", "Precipitation6h", "LowCloudShape",
   
8. #  "LowCloudAmount", "LowCloudHeight", "DewPoint", "Visibility", "WeatherNow",
   
9. #  "Temperature", "MiddleCloudShape", "HighCloudShape"
   
10. variables = ['Visibility','WeatherNow','WindDirection','WindSpeed','Temperature','DewPoint']
    
11. 
    
12. #Read data
    
13. datadir = r'D:\Temp\micaps'
    
14. year = 2016
    
15. month = 1
    
16. day = 1
    
17. t = datetime.datetime(year, month, day)
    
18. rdf = None
    
19. for hour in range(2, 24, 3):
    
20.     t = t.replace(hour=hour)
    
21.     fn = os.path.join(datadir, t.strftime('%y%m%d%H') + '.000')
    
22.     print fn
    
23.     f = addfile_micaps(fn)
    
24.     df = f.read_dataframe()
    
25.     t = f.gettime(0)
    
26.     r = df.loc[stids,variables]
    
27.     r.insert(0, 'Time', t)
    
28.     if rdf is None:
    
29.         rdf = r
    
30.     else:
    
31.         rdf = rdf.append(r)
    
32. 
    
33. rdf.index = rdf['Time']
    
34. rdf = rdf.drop('Time')
    
35. print rdf
    
36. rdf.to_csv('D:/Temp/test/test.csv')

#---- Set stations<br /> stids = ['54511']<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 = ['Visibility','WeatherNow','WindDirection','WindSpeed','Temperature','DewPoint']<br /> <br /> #Read data<br /> datadir = r'D:\Temp\micaps'<br /> year = 2016<br /> month = 1<br /> day = 1<br /> t = datetime.datetime(year, month, day)<br /> rdf = None<br /> for hour in range(2, 24, 3):<br />     t = t.replace(hour=hour)<br />     fn = os.path.join(datadir, t.strftime('%y%m%d%H') + '.000')<br />     print fn<br />     f = addfile_micaps(fn)<br />     df = f.read_dataframe()<br />     t = f.gettime(0)<br />     r = df.loc[stids,variables]<br />     r.insert(0, 'Time', t)<br />     if rdf is None:<br />         rdf = r<br />     else:<br />         rdf = rdf.append(r)<br /> <br /> rdf.index = rdf['Time']<br /> rdf = rdf.drop('Time')<br /> print rdf<br /> rdf.to_csv('D:/Temp/test/test.csv')

输出的数据：              Visibility WeatherNow WindDirection WindSpeed Temperature DewPoint
2016-01-01 02       26.0          0         360.0       1.0        -3.4    -18.1
2016-01-01 05       30.0          0         340.0       2.0        -4.8    -18.6
2016-01-01 08        1.8         10          70.0       1.0        -4.9     -9.0
2016-01-01 11        3.8          5          20.0       1.0         1.4     -8.4
2016-01-01 14        2.8          5          90.0       1.0         4.9     -8.2
2016-01-01 17        1.7          5         180.0       1.0         2.8     -6.3
2016-01-01 20        1.1         10         180.0       1.0        -1.2     -6.5
2016-01-01 23        0.9         10           0.0       0.0        -2.9     -6.5

luohan552200

赞赞赞！

孤蓝et

老师，MeteoInfoLab脚本编写的py代码，单纯只做气象数据批量处理，不进行绘图操作，这类py代码可以脱离GUI界面，直接后台调用吗？

MeteoInfo

孤蓝et 发表于 2018-11-28 09:17
老师，MeteoInfoLab脚本编写的py代码，单纯只做气象数据批量处理，不进行绘图操作，这类py代码可以脱离GUI ...

包括绘图操作也可以后台调用

孤蓝et

谢谢解答，我再仔细看下文档，找一下具体怎么操作。

孤蓝et

MeteoInfo 发表于 2018-11-28 11:24
包括绘图操作也可以后台调用

抱歉，王老师，又打扰了。
之前问您的后台脚本，我看了论坛和meteothink，都没有找到MeteoInfoLab进行后台调用操作的文档唉，是不是我找的关键字不对？
还是，您说的后台调用，指的是C#版本 MeteoInfo下的Script脚本？

MeteoInfo

孤蓝et 发表于 2018-11-29 09:00
抱歉，王老师，又打扰了。
之前问您的后台脚本，我看了论坛和meteothink，都没有找到MeteoInfoLab进行后 ...

参考这里：http://www.meteothink.org/docs/i ... nfomap-meteoinfolab

孤蓝et

On Windows paltform, you can launch MeteoInfoMap by double-clicking “MeteoInfoMap.exe” file, or by run the batch command file “mimap.bat”.

The batch command file with parameters will run a Jython script (“mimap.bat test.py”).

看到了，谢谢老师指点。

lasamm

王老师，我在MeteoInfoLab的脚步里面运行这个代码，总是提示：
NameError: name 'addfile_micaps' is not defined

是还需要导入什么插件或者包么？首次使用，希望王老师指点一下

MeteoInfo

lasamm 发表于 2020-7-2 14:31
王老师，我在MeteoInfoLab的脚步里面运行这个代码，总是提示：
NameError: name 'addfile_micaps' is not  ...

应该不需要了，软件启动的时候已经导入了。