云导风数据的检验方法和装置
技术领域
本发明涉及气象技术领域,具体而言,特别涉及一种云导风数据的检验方法和装置。
背景技术
风是气象中最重要的要素之一,它与其他气象要素有着紧密的联系。风对人们的生产、生活产生着重大的影响。众所周知常规手段获取的风场信息主要来自于测站的测量,对于像高山、海洋、沙漠和极地等缺乏测站的地区,就很难获取真实的风场资料,卫星资料反演风场是重要的信息源。近几十年来,利用气象卫星资料反演的云导风已广泛应用于数值预报、台风分析等,并越来越显示其重要的作用。
然而,当云导风数据本身存在问题,例如准确性较低时,如果采用这样的数据进行数值预报、台风分析等,得到的结果必然也不准确。为此,本申请提出了一种云导风数据的检验方法和装置。
发明内容
本发明的主要目的在于提供一种云导风数据的检验方法和装置,以在使用云导风数据之前对其的准确性进行检验。
为解决上述技术问题,本发明公开了一种云导风数据的检验方法。该方法包括:从待检验云导风数据和检验数据中,匹配时间差在预设时间差范围内的一个云导风数据文件和一个检验数据文件;计算所述云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号,根据最临近法,获取整数型的网络索引号;根据所述检验数据文件和所述网络索引号获取所述检验数据文件的风速数据;计算所述检验数据文件的风速数据的平均值,以得到检验值;计算所述云导风数据文件的风速数据的平均值,以得到待检验值;计算所述检验值与所述待检验值两者之间差的统计量。
进一步地,所述检验数据为T639数据、欧洲中期天气预报中心数据或探空数据。
进一步地,在执行计算所述云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号之前,该方法还包括:去除所述云导风数据文件和所述检验数据文件中气压层小于预设气压值的数据。
进一步地,在执行计算所述云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号之前,该方法还包括:判断所述云导风数据文件的风速数据是否在预设云导风数据有效值范围内;其中,当判断结果为所述云导风数据文件的风速数据在所述预设云导风数据有效值范围内时,执行计算所述云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号的步骤。
进一步地,在执行计算所述云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号之前,该方法还包括:查找所述云导风数据文件中是否存在相邻气压层的数据,其中,当查找结果为所述云导风数据文件存在相邻气压层的数据时,执行计算所述云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号的步骤。
进一步地,在执行计算所述云导风数据文件的风速数据的平均值,以得到待检验值的步骤之前,该方法还包括:判断所述云导风数据文件的经纬度是否在预设的卫星观测范围内,其中,仅当判断结果为所述云导风数据文件的经纬度在所述预设的卫星观测范围内时,执行计算所述云导风数据文件的风速数据的平均值,以得到待检验值的步骤。
为解决上述技术问题,本发明还公开了一种云导风数据的检验方装置。该装置包括:匹配模块,用于从待检验云导风数据和检验数据中,匹配时间差在预设时间差范围内的一个云导风数据文件和一个检验数据文件;第一获取模块,用于计算所述云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号,根据最临近法,获取整数型的网络索引号;第二获取模块,用于根据所述检验数据文件和所述网络索引号获取所述检验数据文件的风速数据;第一计算模块,用于计算所述检验数据文件的风速数据的平均值,以得到检验值;第二计算模块,用于计算所述云导风数据文件的风速数据的平均值,以得到待检验值;第三计算模块,用于计算所述检验值与所述待检验值两者之间差的统计量。
进一步地,所述检验数据为T639数据、欧洲中期天气预报中心数据或探空数据。
进一步地,该装置还包括:处理模块,用于在所述第一获取模块执行计算所述云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号之前,去除所述云导风数据文件和所述检验数据文件中气压层小于预设气压值的数据。
进一步地,该装置还包括:判断模块,用于在所述第一获取模块执行计算所述云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号之前,判断所述云导风数据文件的风速数据是否在预设云导风数据有效值范围内,其中,当所述判断模块的判断结果为所述云导风数据文件的风速数据在所述预设云导风数据有效值范围内时,所述第一获取模块执行计算所述云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号的步骤。
通过本发明,能够对云导风数据的质量进行检验,从而保证采用高质量的云导风数据进行数值预报、台风分析,进而能够得到准确的数值预报和台风分析结果。
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,而可依照说明书的内容予以实施,并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂,以下特举本发明的具体实施方式。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1是根据本发明实施例一的云导风数据的检验方法的流程图;
图2是根据本发明实施例二的云导风数据的检验方法的流程图;
图3是根据本发明实施例六的云导风数据的检验装置的框图。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步说明。需要指出的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。
实施例一
首先,本发明实施例一提供了一种云导风数据的检验方法,该方法主要描述采用检验数据对云导风数据进行检验的过程,参见图1,该方法可以包括以下步骤:
步骤S102:从待检验云导风数据和检验数据中,匹配时间差在预设时间差范围内的一个云导风数据文件和一个检验数据文件。
其中,待检验云导风数据包括多个云导风数据文件,各个云导风数据文件之间相差约几个小时,检验数据包括多个检验数据文件,各个检验数据文件之间相差约几个小时,例如将预设时间差范围设置为30分钟,则对待检验云导风数据和检验数据中的数据文件进行时间匹配时,可匹配到时间差在30分钟内的一个云导风数据文件和一个检验数据文件。
云导风数据文件和检验数据文件的文件名中通常带有文件的时间,因而,在对待检验云导风数据和检验数据中的数据文件进行时间匹配时,只需提取云导风数据文件和检验数据文件的文件名中的时间,根据提取的时间进行匹配,找到两个文件的时间差在预设时间差范围,也即30分钟之内即可。
步骤S104:计算云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号,根据最临近法,获取整数型的网络索引号。
步骤S106:根据检验数据文件和网络索引号获取检验数据文件的风速数据。
具体采用网络索引号获取检验数据文件的风速数据的方法可采用现有技术中任一种方法,此处不再赘述。
步骤S108:计算检验数据文件的风速数据的平均值,以得到检验值。
步骤S110:计算云导风数据文件的风速数据的平均值,以得到待检验值。
步骤S112:计算检验值与待检验值两者之间差的统计量。
具体包括计算两者之间的误差值、标准方差以及差的标准差,通过统计量来反映云导风数据的质量,其中,误差值越小、标准方差越小、差的标准差越小,云导风数据的质量越好。
采用该实施例提供的云导风数据的检验方法,能够对云导风数据的质量进行检验,从而保证采用高质量的云导风数据进行数值预报、台风分析,进而能够得到准确的数值预报和台风分析结果。
实施例二
该实施例为在实施例一的基础上进一步优选的云导风数据的检验方法,参见图2,该方法可以包括以下步骤:
步骤S202:获取待检验云导风数据对应的卫星观测范围、云导风数据有效值范围和相邻两个气压层阈值。
其中,卫星观测范围为卫星参数,数据有效值范围和相邻两个气压层阈值为经验值。例如,该方法通过计算机程序实现时,将卫星观测范围、数据有效值范围和相邻两个气压层阈值这些参数写入程序的配置文件,从而在程序初始,解析配置文件即可获取到待检验云导风数据对应的卫星观测范围、数据有效值范围和相邻两个气压层阈值。
步骤S204:从待检验云导风数据和检验数据中,匹配时间差在预设时间差范围内的一个云导风数据文件和一个检验数据文件。
该步骤的具体实现方式与实施例一中步骤S102的实现方式相同,此处不再赘述。
步骤S206:去除云导风数据文件和检验数据文件中气压层小于预设气压值的数据。
在匹配到一个云导风数据文件和一个检验数据文件之后,根据文件格式对两个文件分别进行解析,解析后得到数据文件中的数据本身,进而将解析得到的数据中气压层小于预设气压值的数据去除,例如将预设气压值设为30。
步骤S208:判断云导风数据文件的风速数据是否在预设云导风数据有效值范围内。
具体地,将云导风数据文件的风速数据与步骤S202中获取到的云导风数据有效值范围进行比较,当云导风数据文件的风速数据在云导风数据有效值范围内时,执行步骤S210。
步骤S210:查找云导风数据文件中是否存在相邻气压层的数据。
具体地,针对在云导风数据有效值范围内的风速数据,按照步骤S202中相邻两个气压层阈值进行查找,当查找到有相邻气压层的数据时,执行步骤S212。
步骤S212:计算云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号,根据最临近法,获取整数型的网络索引号。
步骤S214:判断云导风数据文件的经纬度是否在预设的卫星观测范围内。
具体地,将云导风数据文件的经纬度与步骤S202中的卫星观测范围进行比较,当云导风数据文件的经纬度是在预设的卫星观测范围内,执行步骤S216。
步骤S216:根据检验数据文件和网络索引号获取检验数据文件的风速数据。
步骤S218:计算检验数据文件的风速数据的平均值,以得到检验值,计算云导风数据文件的风速数据的平均值,以得到待检验值,计算检验值与待检验值两者之间差的统计量。
实施例三
该实施例为在实施例一和实施例二的基础上进一步优选的云导风数据(也即AMV数据)的检验方法,在该方法中,采用T639数据作为检验数据,对云导风数据进行检验,具体检验流程如下:
步骤S302:解析配置文件;
配置文件中主要记录卫星观测范围,AMV数据有效值范围及相邻两个气压层阈值设置等参数。
步骤S304:时间匹配;
提取AMV数据和T639数据的数据文件名中的时间,根据所提取的时间进行匹配,两个文件的时间相差在30分钟以内。
步骤S306:数据解析;
根据相应的格式,对匹配到的AMV数据文件和T639数据文件进行解析,并对气压层小于30的数据进行过滤。
步骤S308:空间匹配;
由于AMV数据是散点数据,所以通过遍历AMV数据做空间匹配。
首先判断AMV风速数据是否为有效数据,如果为无效值跳过不做处理;再次通过查找是否有相邻气压层的数据,如果没有,该数据跳过不做处理;最后计算AMV数据投影到全球等经纬度二维平面中的浮点型网格索引号,根据最临近法,获取整数型的网格索引号。
步骤S310:数值计算;
判断AMV数据经纬度是否在观测范围内,如果不在观测范围,直接跳过。根据空间匹配所得到的网格索引号,获取T639数据.U和.V文件分别对应的值,并合成T639数据风速数据。
分别计算,AMV数据的风速平均值,T639数据的风速平均值;两者之间的误差值、标准方差和差的标准差。
步骤S312:数据统计。
具体可通过图像的方式展示出两者之前的误差值、标准方差和差的标准差。例如,可采用如下的展示方法:
a)时序误差统计:是根据检验数据的类型(也即数据类型为T639),通道类型(包括红外、水汽),开始时间、结束时间、气压层分布(高、中、低)、纬度范围等参数,以时间为X轴,统计值(均值、误差、标准差、样本点)为Y轴进行数据统计展示的方法。
b)时间平均地域平均场:是根据检验数据类型(T639),通道类型(红外、水汽)、数据类型(均值、误差、标准差)、气压范围(0-400hPa、400-700hPa、700-1000hPa),以经纬度为XY轴进行数据统计展示的方法。
c)纬向平均场(时间):是根据检验数据类型(T639),通道类型(红外、水汽)、数据类型(均值、误差、标准差)、气压范围(0-400hPa、400-700hPa、700-1000hPa),以时间为X轴,纬度为Y轴进行数据统计展示的方法。
d)纬向平均场(气压层):是根据检验数据类型(T639),通道类型(红外、水汽)、数据类型(均值、误差、标准差)、以纬度为X轴,气压层为Y轴进行数据统计展示的方法。
实施例四
该实施例为在实施例一和实施例二的基础上进一步优选的云导风数据的检验方法,在该方法中,采用欧洲中期天气预报中心数据(ECMWF数据)作为检验数据,对云导风数据进行检验,具体检验流程如下:
步骤S402:解析配置文件;
配置文件中主要记录卫星观测范围,AMV数据有效值范围及相邻两个气压层阈值设置等参数。
步骤S404:时间匹配;
提取AMV数据和ECMWF数据文件名中的时间,根据所提取的时间进行匹配,两个文件的时间相差在30分钟以内。
步骤S406:数据解析;
根据相应的格式,对AMV数据和ECMWF数据文件进行解析,并对气压层小于30的数据进行过滤。
步骤S408:空间匹配;
由于AMV数据是散点数据,所以通过遍历AMV数据做空间匹配。
首先判断AMV风速数据是否为有效数据,如果为无效值跳过不做处理;再次通过查找是否有相邻气压层的数据,如果没有,该数据跳过不做处理;最后计算AMV数据投影到全球等经纬度二维平面中的浮点型网格索引号,根据最临近法,获取整数型的网格索引号。
步骤S410:数值计算;
判断AMV数据经纬度是否在观测范围内,如果不在观测范围,直接跳过。根据空间匹配所得到的网格索引号,获取ECMWF数据文件数据风速数据。
分别计算,AMV数据的风速平均值,ECMWF数据的风速平均值;两者之间的误差值;标准方差;差的标准差;
步骤S412:数据统计;
a)时序误差统计:是根据检验数据类型(ECMWF),通道类型(红外、水汽),开始时间、结束时间、气压层分布(高、中、低)、纬度范围等参数,以时间为X轴,统计值(均值、误差、标准差、样本点)为Y轴进行数据统计展示的方法。
b)时间平均地域平均场:是根据检验数据类型(ECMWF),通道类型(红外、水汽)、数据类型(均值、误差、标准差)、气压范围(0-400hPa、400-700hPa、700-1000hPa),以经纬度为XY轴进行数据统计展示的方法。
c)纬向平均场(时间):是根据检验数据类型(ECMWF),通道类型(红外、水汽)、数据类型(均值、误差、标准差)、气压范围(0-400hPa、400-700hPa、700-1000hPa),以时间为X轴,纬度为Y轴进行数据统计展示的方法。
d)纬向平均场(气压层):是根据检验数据类型(ECMWF),通道类型(红外、水汽)、数据类型(均值、误差、标准差)、以纬度为X轴,气压层为Y轴进行数据统计展示的方法。
实施例五
该实施例为在实施例一和实施例二的基础上进一步优选的云导风数据的检验方法,在该方法中,采用探空数据作为检验数据,对云导风数据进行检验。
步骤S502:解析配置文件;
配置文件中主要记录卫星观测范围,AMV数据有效值范围及相邻两个气压层阈值设置等参数。
步骤S504:时间匹配;
提取AMV数据和探空数据文件名中的时间,根据所提取的时间进行匹配,两个文件的时间相差在30分钟以内。
步骤S506:数据解析;
根据相应的格式,对AMV数据和探空数据文件进行解析,并对气压层小于30的数据进行过滤。
步骤S508:空间匹配;
由于AMV数据是散点数据,所以通过遍历AMV数据做空间匹配。
首先判断AMV风速数据是否为有效数据,如果为无效值跳过不做处理;再次通过查找是否有相邻气压层的数据,如果没有,该数据跳过不做处理;最后计算AMV数据投影到等经纬度二维平面中的浮点型网格索引号,根据最临近法,获取整数型的网格索引号。
步骤S510:数值计算;
判断AMV数据经纬度是否在观测范围内,如果不在观测范围,直接跳过。根据空间匹配所得到的网格索引号,获取探空数据风速数据。
分别计算,AMV数据的风速平均值,探空数据的风速平均值;两者之间的误差值;标准方差;差的标准差。
步骤S512:数据统计;
时序误差统计:是根据检验数据类型(探空),通道类型(红外、水汽),开始时间、结束时间、气压层分布(高、中、低)、纬度范围等参数,以时间为X轴,统计值(均值、误差、标准差、样本点)为Y轴进行数据统计展示的方法。
时间平均地域平均场:是根据检验数据类型(探空),通道类型(红外、水汽)、数据类型(均值、误差、标准差)、气压范围(0-400hPa、400-700hPa、700-1000hPa),以经纬度为XY轴进行数据统计展示的方法。
纬向平均场(时间):是根据检验数据类型(探空),通道类型(红外、水汽)、数据类型(均值、误差、标准差)、气压范围(0-400hPa、400-700hPa、700-1000hPa),以时间为X轴,纬度为Y轴进行数据统计展示的方法。
纬向平均场(气压层):是根据检验数据类型(探空),通道类型(红外、水汽)、数据类型(均值、误差、标准差)、以纬度为X轴,气压层为Y轴进行数据统计展示的方法。
以上是对本发明所提供的云导风数据的检验方法进行的描述。下面将对本发明提供的云导风数据的检验装置进行描述,需要说明的是,该装置可用于执行上述任意一种云导风数据的检验方法。
实施例六
与本发明实施例一提供的云导风数据的检验方法相对应,本发明实施例还提供了一种云导风数据的检验装置,参见图3,该装置可以包括匹配模块610、第一获取模块620、第二获取模块630、第一计算模块640、第二计算模块650和第三计算模块660。
匹配模块610用于从待检验云导风数据和检验数据中,匹配时间差在预设时间差范围内的一个云导风数据文件和一个检验数据文件。第一获取模块620用于计算云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号,根据最临近法,获取整数型的网络索引号;第二获取模块630用于根据检验数据文件和网络索引号获取检验数据文件的风速数据;第一计算模块640用于计算检验数据文件的风速数据的平均值,以得到检验值;第二计算模块650用于计算云导风数据文件的风速数据的平均值,以得到待检验值;第三计算模块660用于计算检验值与待检验值两者之间差的统计量。
优选地,检验数据为T639数据、欧洲中期天气预报中心数据或探空数据。
优选地,该检验装置还包括:处理模块,用于在第一获取模块执行计算云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号之前,去除云导风数据文件和检验数据文件中气压层小于预设气压值的数据。
优选地,该检验装置还包括:第一判断模块,用于在第一获取模块执行计算云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号之前,判断云导风数据文件的风速数据是否在预设云导风数据有效值范围内。
其中,当第一判断模块的判断结果为云导风数据文件的风速数据在预设云导风数据有效值范围内时,第一获取模块执行计算云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号的步骤。
优选地,该检验装置还包括:查找模块,用于在第一获取模块执行计算云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号之前,查找模块查找云导风数据文件中是否存在相邻气压层的数据。
其中,当查找结果为云导风数据文件存在相邻气压层的数据时,第一获取模块执行计算云导风数据文件中的数据投影到全球等经纬度二维平面中的浮点型网络索引号的步骤。
优选地,该检验装置还包括:第二判断模块,用于第二计算模块执行计算云导风数据文件的风速数据的平均值,以得到待检验值之前,判断云导风数据文件的经纬度是否在预设的卫星观测范围内。
其中,仅当判断结果为云导风数据文件的经纬度在预设的卫星观测范围内时,第二计算模块执行计算云导风数据文件的风速数据的平均值,以得到待检验值的步骤。
需要说明的是,上述装置实施例属于优选实施例,所涉及的模块并不一定是本申请所必须的。
本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。对于本申请的装置实施例而言,由于其与方法实施例基本相似,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。