MeteoInfo脚本示例：格点数据平均

MeteoInfo · 发表于 2012-10-22 11:47:33

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 MeteoInfo 于 2012-10-22 12:38 编辑

MeteoInfo中计算多个格点数据的平均值有两种方法：一种是将所有的格点数据累加，然后再除以格点数据个数；另一种方法是将所有格点数据添加到一个列表中，然后用DataMath.Average(List(GridData), Boolean)方法对该列表进行平均计算。如果格点数据有缺测值，两种方法计算的结果可能不同。第二种方法的第二个参数表示是否忽略缺测值，如果为True，则和第一种方法是一样的。

第二种方法的第二个参数为False时，某个格点只有在所有格点场均为缺测值时其平均值才为缺测值，否则会按有效值个数来计算平均，这样可以获得尽可能多的有效数据。比如计算三个月的格点场的平均，三个月的数据中格点（2, 2）都是有效值，那么该格点的平均值为(m1(2, 2) + m2(2, 2) + m3(2, 2)) / 3。如果其中一个月该格点的数据是缺测值（比如m2为缺测值），那么该格点的平均值为(m1(2, 2) + m3(2, 2) ) / 2，缺测数据不参与平均计算。但在第一种平均方法或第二种平均方法中第二个参数为True时该格点的平均值为缺测值。

讲起来有点绕，仔细体会一下应该能理解。这里给出一个脚本示例，打开model.ctl文件，将所有时次的U变量格点数据进行平均，使用两种方法来看结果的不同。

第一种方法的平均结果，有效值较少：

第二种方法的平均结果，有效值较多：

完整的脚本程序如下：

#--------------------------------------------------------
# Author: Yaqiang Wang
# Date: 2012-10-22
# Purpose: Calculate average grid data - tow methods
# Note: Sample
#-----------------------------------------------------------
import clr
from System.Collections.Generic import *
clr.AddReference("MeteoInfoC.dll")
from MeteoInfoC import *
from MeteoInfoC.Data import *
from MeteoInfoC.Data.MeteoData import *
#---- Set directories
BaseDir = "C:\\Program Files (x86)\\MeteoInfo\"
DataDir = BaseDir + "Sample\\GrADS\"
#---- Average grid data
mydata = MeteoDataInfo()
mydata.OpenGrADSData(DataDir + "model.ctl")
tnum = mydata.GetTimeNumber()
datagrid = GridData()
datalist = List[GridData]()
for i in range(0, tnum):
mydata.TimeIndex = i
gdata = mydata.GetGridData("U")
if i == 0:
datagrid = gdata
else:
datagrid = datagrid + gdata
datalist.Add(gdata)
avedata1 = datagrid / tnum
avedata2 = DataMath.Average(datalist, False)
mipy.MeteoDataset.SetDrawType("shaded")
avelayer1 = mipy.MeteoDataset.Display(avedata1)
avelayer1.LayerName = "Ave_Method1"
avelayer2 = mipy.MeteoDataset.Display(avedata2)
avelayer2.LayerName = "Ave_Method2"
mipy.MapDocument.ActiveMapFrame.UpdateLayerNode(avelayer1)
mipy.MapDocument.ActiveMapFrame.UpdateLayerNode(avelayer2)
print 'Finished!'

复制代码

wbriver · 发表于 2012-10-22 21:05:29

端板凳过来学习~

开始起飞1 · 发表于 2013-4-25 22:03:43

错误提示

请问楼主，我把脚本文件中的model.ctl 改成我要用的1978data.ctl，然后把U改成了a，把路径（X86）去掉了改成我需要的路径，但为什么就出现这个错误呢？不知道怎么改，麻烦指点下吧

MeteoInfo · 发表于 2013-4-25 23:18:19

开始起飞1 发表于 2013-4-25 22:03

登录/注册后可看大图

请问楼主，我把脚本文件中的model.ctl 改成我要用的1978data.ctl，然后把U改成了a，把路径（X86）去掉了改成 ...

错误提示被打上了水印，看不清楚。

开始起飞1 · 发表于 2013-4-26 08:06:41

MeteoInfo 发表于 2013-4-25 23:18

登录/注册后可看大图

错误提示被打上了水印，看不清楚。

file:///C:/Documents%20and%20Settings/Administrator/Application%20Data/Tencent/Users/617039055/QQ/WinTemp/RichOle/ETY~D6B(B_62B5]A~TMWJCL.jpg 这样行了吗？

开始起飞1 · 发表于 2013-4-26 08:07:23

本帖最后由开始起飞1 于 2013-4-26 08:09 编辑

开始起飞1 发表于 2013-4-26 08:06

登录/注册后可看大图

这样行了吗？

Traceback (most recent call last)
File"<string>",line 38,in <module>
ValueError: 列“” 不属于表

开始起飞1 · 发表于 2013-4-26 08:10:21

MeteoInfo 发表于 2013-4-25 23:18

登录/注册后可看大图

错误提示被打上了水印，看不清楚。

Traceback (most recent call last)
File"<string>",line 38,in <module>
ValueError: 列“” 不属于表

MeteoInfo · 发表于 2013-4-26 09:05:19

开始起飞1 发表于 2013-4-26 08:10

登录/注册后可看大图

Traceback (most recent call last)
File"",line 38,in
ValueError: 列“” 不属于表

你更新一下MeteoInfo的最新文件试试，见置顶帖子：MeteoInfo最新文件
http://bbs.06climate.com/forum.p ... 126&fromuid=106

血嗜子木 · 发表于 2015-7-19 18:04:06

谢谢您，这个方法可以不用做庞大的编程工作

		自动登录	找回密码
密码			立即注册

MeteoInfo脚本示例：格点数据平均

登录后查看更多精彩内容~

评分