CN108830444A - 一种探空观测数据的评估和修正方法及装置 - Google Patents
一种探空观测数据的评估和修正方法及装置 Download PDFInfo
- Publication number
- CN108830444A CN108830444A CN201810359842.4A CN201810359842A CN108830444A CN 108830444 A CN108830444 A CN 108830444A CN 201810359842 A CN201810359842 A CN 201810359842A CN 108830444 A CN108830444 A CN 108830444A
- Authority
- CN
- China
- Prior art keywords
- data
- raob
- barosphere
- coefficient
- suspicious
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Development Economics (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Operations Research (AREA)
- Marketing (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Complex Calculations (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明的实施例公开了一种探空观测数据的评估和修正方法及装置,该方法针对欲进行评估的目标探空观测数据中涉及到的每一气压层,分别通过在该气压层观测到的第二探空观测数据和相应的模式背景数据,计算出用于评价该气压层的探空观测数据的质量的评价参数,通过该评价参数对目标探空观测数据的质量进行评估和修正,以改善探空观测数据的质量。该方法针对同一气压层的探空观测数据,通过上述评价参数共同评价该层的探空观测数据的质量,评估结果更为准确。同时,依据这些评价参数对探空观测数据进行修正,使得探空观测数据更为合理,便于将这些数据应用到实验或者气象分析中。
Description
技术领域
本发明实施例涉及探空数据处理技术领域,尤其是涉及一种探空观测数据的评估和修正方法及装置。
背景技术
现有的对探空数据的处理都是基于数值预报同化系统提供的背景场,计算观测资料与背景场的差别(简称观测余差);对观测资料进行质量控制排除异常值;统计分析各类观测报告的接收率、观测报告时空分布、异常值比率、观测余差的平均偏差和标准差时空分布等。比如,全球观测资料质量监视评估系统以资料同化的6小时预报场为背景场,计算观测资料与背景场的偏差;运用质量控制技术排除观测资料中的可能是错误的资料;实时和非实时地统计分析各类观测报告的状况和观测资料的质量。该系统应用于日常观测资料质量监视评估及中国区域探空资料问题的详细分析、探空站综合评分、探空观测资料归档、探空观测资料偏差订正效果检验、L波段探空资料质量分析等,取得了显著效果。
除此之外,还有一些通过模式验证法来评估探空观测数据的质量的方法,例如,采用欧洲数值预报中心(ECMWF)全球数值天气预报模式6小时预报场(初估场)的100hPa位势高度和100-30hPa厚度资料作为参考,对已经运行L波段探测系统的高空气象台站的资料从平均偏差和标准偏差的角度进行了分析研究;以GRAPES、EC等不同模式分析场为背景场,评估中国L波段探空湿度观测资料的质量状况;利用初估场来剔除由天气过程引起的气象要素场的地区差异,从而进行探空仪的间接对比。
然而,模式验证法在上述的研究中心都只有从偏差,平均偏差和标准差三方面评估探空资料,没有更深入的研究探讨。而且这些研究都只是给出了一些定性的评估结论,没有科学地统计确定一个定量的评估标准。
在实现本发明实施例的过程中,发明人发现现有的对探空数据进行评估和修正的方法,只是定性的进行评估和修正,评估结果不准确,依据该评估对探空数据进行修正得到气象数据的质量较低。
发明内容
本发明所要解决的技术问题是现有的对探空数据进行评估和修正的方法,只是定性的进行评估和修正,评估结果不准确,依据该评估对探空数据进行修正得到气象数据的质量较低的问题。
针对以上技术问题,本发明的实施例提供了一种探空观测数据的评估和修正方法,包括:
获取与欲进行评估的目标探空观测数据有关的所有气压层;
针对每一气压层,从所述目标探空观测数据中获取在所述气压层下观测得到的探空观测数据,作为第一探空观测数据,并获取观测所述第一探空观测数据的第一观测时间;
获取在包括所述第一观测时间的第二观测时间观测到的对应于所述气压层的探空观测数据,作为第二探空观测数据,以及与每一第二探空观测数据对应的模式背景数据,作为目标模式背景数据;
根据所述第二探空观测数据和所述目标模式背景数据计算在所述气压层下观测的探空观测数据的至少一个评价参数;
根据计算出的所述气压层对应的所有评价参数,对所述第一探空观测数据的质量进行评估,并对所述第一探空观测数据进行修正;
其中,所述评价参数包括平均偏差、标准偏差、均方根误差、相关系数、偏度系数和峰度系数。
第二方面,本实施例提供了一种探空观测数据的评估和修正装置,包括:
第一获取模块,用于获取与欲进行评估的目标探空观测数据有关的所有气压层;
第二获取模块,用于针对每一气压层,从所述目标探空观测数据中获取在所述气压层下观测得到的探空观测数据,作为第一探空观测数据,并获取观测所述第一探空观测数据的第一观测时间;
第三获取模块,用于获取在包括所述第一观测时间的第二观测时间观测到的对应于所述气压层的探空观测数据,作为第二探空观测数据,以及与每一第二探空观测数据对应的模式背景数据,作为目标模式背景数据;
计算模块,用于根据所述第二探空观测数据和所述目标模式背景数据计算在所述气压层下观测的探空观测数据的至少一个评价参数;
处理模块,用于根据计算出的所述气压层对应的所有评价参数,对所述第一探空观测数据的质量进行评估,并对所述第一探空观测数据进行修正;
其中,所述评价参数包括平均偏差、标准偏差、均方根误差、相关系数、偏度系数和峰度系数。
第三方面,本发明的实施例还提供了一种电子设备,包括:
至少一个处理器;以及
与所述处理器通信连接的至少一个存储器,其中:
所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行以上所述的方法。
第四方面,本发明的实施例提供了一种非暂态计算机可读存储介质,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令使所述计算机执行以上所述的方法。
本发明的实施例提供了一种探空观测数据的评估和修正方法及装置,该方法针对欲进行评估的目标探空观测数据中涉及到的每一气压层,分别通过在该气压层观测到的第二探空观测数据和相应的模式背景数据,计算出用于评价该气压层的探空观测数据的质量的评价参数,通过该评价参数对目标探空观测数据的质量进行评估。同时,针对目标探空观测数据中处于该气压层的探空观测数据,根据上述评价参数和根据探空观测数据和相应的模式背景数据的差值对处于该气压层的探空观测数据进行修正,以改善探空观测数据的质量。该方法针对同一气压层的探空观测数据,通过上述评价参数共同评价该层的探空观测数据的质量,评估结果更为准确。同时,依据这些评价参数对探空观测数据进行修正,使得探空观测数据更为合理,便于将这些数据应用到实验或者气象分析中。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明一个实施例提供的探空观测数据的评估和修正方法的流程示意图;
图2是本发明另一个实施例提供的目标探空观测数据为2017年11月7日8时的温度数据和第一差值的示意图;
图3是本发明另一个实施例提供的2017年11月每天的8时在850hPa观测到的探空观测数据的示意图;
图4是本发明另一个实施例提供的对图3中的探空观测数据进行修正的示意图;
图5是本发明另一个实施例提供的2017年9月1日00时的温度数据和第一差值的示意图;
图6是本发明另一个实施例提供的各气压层对应的平均偏差的示意图;
图7是本发明另一个实施例提供的各气压层对应的标准偏差的示意图;
图8是本发明另一个实施例提供的每一层偏差的概率密度分布(PDF)示意图;
图9是本发明另一个实施例提供的对修正后的数据进行验证所计算的平均偏差示意图;
图10是本发明另一个实施例提供的对修正后的数据进行验证所计算的标准偏差示意图;
图11是本发明另一个实施例提供的探空观测数据的评估和修正装置的结构框图;
图12是本发明另一个实施例提供的电子设备的结构框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1是本实施例提供的一种探空观测数据的评估和修正方法的流程示意图,参见图1,该方法包括:
101:获取与欲进行评估的目标探空观测数据有关的所有气压层;
102:针对每一气压层,从所述目标探空观测数据中获取在所述气压层下观测得到的探空观测数据,作为第一探空观测数据,并获取观测所述第一探空观测数据的第一观测时间;
103:获取在包括所述第一观测时间的第二观测时间观测到的对应于所述气压层的探空观测数据,作为第二探空观测数据,以及与每一第二探空观测数据对应的模式背景数据,作为目标模式背景数据;
104:根据所述第二探空观测数据和所述目标模式背景数据计算在所述气压层下观测的探空观测数据的至少一个评价参数;
105:根据计算出的所述气压层对应的所有评价参数,对所述第一探空观测数据的质量进行评估,并对所述第一探空观测数据进行修正;
其中,所述评价参数包括平均偏差、标准偏差、均方根误差、相关系数、偏度系数和峰度系数。
需要说明的是,本实施例提供的方法通常由计算机或者服务器执行。目标探空观测数据可以是某个时间点观测到的一组探空数据,例如,如图2所示,目标探空观测数据为某年某月某日某时在某个地点观测到的不同气压层的温度数据,或者,如图3所示,目标探空观测数据为在某个地点观测到的某年某月的每一天某个时刻的温度数据。本实施例中的所说的探空观测数据可以是诸如温度、湿度这样的气象数据,本实施例对此不做具体限制。由于不同的气压层对应的探空观测数据差别较大,因此计算用于评价探空观测数据的数据质量的评价参数时,选取在该气压层下测量的到探空观测数据进行计算。
与目标探空观测数据有关的所有气压层指的是目标探空观测数据涉及到的所有气压层,例如,图2中的各点对应的气压层,图3中仅涉及850hPa一个气压层。针对每一气压层,第一探空观测数据是目标探空观测数据中,在该气压层下观测得到的探空观测数据,例如,针对图2中的800hPa这一气压层,第一探空观测数据仅为800hPa对应的数据点。该数据点对应的第一观测时间为2017年11月7日8时。本市实施例中选取获取2017年11月每一天8时在800hPa这一气压层观测到的一组温度数据作为第二探空观测数据。
评价参数是以模式背景数据作为参考,计算的用于评价某一气压层的数据质量的参数。模式背景数据用于对探空观测数据进行评估的参考数据。本实施例中计算的评价参数至少包括3个,实际中也可以加入更多的参数作为对探空数据质量进行评价的参数。
本实施例提供了一种探空观测数据的评估和修正方法,该方法针对欲进行评估的目标探空观测数据中涉及到的每一气压层,分别通过在该气压层观测到的第二探空观测数据和相应的模式背景数据,计算出用于评价该气压层的探空观测数据的质量的评价参数,通过该评价参数对目标探空观测数据的质量进行评估。同时,针对目标探空观测数据中处于该气压层的探空观测数据,根据上述评价参数和根据探空观测数据和相应的模式背景数据的差值对处于该气压层的探空观测数据进行修正,以改善探空观测数据的质量。该方法针对同一气压层的探空观测数据,通过上述评价参数共同评价该层的探空观测数据的质量,评估结果更为准确。同时,依据这些评价参数对探空观测数据进行修正,使得探空观测数据更为合理,便于将这些数据应用到实验或者气象分析中。
更进一步地,在上述实施例的基础上,所述根据计算出的所述气压层对应的所有评价参数对所述第一探空观测数据进行修正,包括:
针对每一所述第一探空观测数据,计算所述第一探空观测数据和与其相同观测条件对应的目标模式背景数据的第一差值,若所述第一差值大于第一阈值,则判定所述第一探空观测数据为错误数据,若所述第一差值小于或等于所述第一阈值,大于第二阈值,则判定所述第一探空观测数据为可疑数据;
分别根据每一所述评价参数判断在所述气压层下观测得到的探空观测数据是否属于可疑层数据,若根据每一所述评价参数均判断在所述气压层下观测得到的探空观测数据属于可疑层数据,则对所述可疑数据和所述错误数据进行修正,否则,对所述错误数据进行修正;
其中,所述第一阈值为所述气压层对应的平均偏差与3倍的标准偏差的差或者和,且所述第二阈值为所述气压层对应的平均偏差与2倍的标准偏差的差;或者,所述第一阈值为所述气压层对应的平均偏差与3倍的标准偏差的差或者和,且所述第二阈值为所述气压层对应的平均偏差与2倍的标准偏差的和。
相同观测条件通常指的是观测地点、观测时间、观测的气压层均相同。如图2左边的图所示,目标探空观测数据为2017年11月7日8时的温度数据,观测点为A地点。模式背景数据为GRAPES模式温度。对于每一探空观测数据,计算其与相同观测条件对应的目标模式背景数据的第一差值。O为实际观测的探空观测数据,B为模式背景数据,图2右边的图示出了每一第一探空观测数据对应的第一差值(即O-B的偏差分布图)。本实施例规定第一差值大于第一阈值(通常第一阈值为平均偏差±3倍标准偏差)为错误数据,大于第二阈值(通常第二阈值为平均偏差±2倍标准偏差)小于第一阈值为可疑数据。每一气压层对应的平均偏差和标准偏差均不相同,因此每一气压层对应的第一阈值和第二阈值也不相同。由图2可以看出,高空观测温度和模式温度之间的趋势差异很小,只有在0-100hPa存在少量的错误数据。
在对探空观测数据进行修正时,对错误数据直接进行修正。对可疑数据需要结合该可疑数据对应的气压层的评价参数确定是否对该可疑数据进行修正。若所有评价参数均判断该气压层对应的探空观测数据为可疑数据,则需要对该可疑数据进行修正,否则不对该可疑数据进行处理。
本实施例提供了一种探空观测数据的评估和修正方法,该方法限定了如何对探空观测数据进行修正,通过对错误数据的直接修正使得该目标探空观测数据不包括与模式背景数据相差较大的数据。通过气压层整体的观测数据的质量来确定是否修正可疑点,使得可疑点的修正充分考虑了探空观测数据呈现的整体趋势,保证了修正后数据的合理性。
更进一步地,在上述各实施例的基础上,对所述错误数据或者所述可疑数据进行修正,包括:
获取欲进行修正的所述错误数据或者所述可疑数据对应的观测条件,从所述目标模式背景数据中获取对应于所述观测条件下的模式背景数据,作为替换数据,用所述替换数据替换所述错误数据或者所述可疑数据。
例如,图3为了对2017年11月每天的8时,在850hPa观测到的探空观测数据的示意图,图4为通过本实施例提供的方法对图3中的探空观测数据进行修正后的示意图。由图3可以看出,11月整个月08时850hPa的温度质量很好,整个月只有19日误差很大,为错误点。通过图4质控后的图可以看出,质控后探空观测数据的平均偏差、标准偏差和均方根误差都变小了,相关系数增大,峰度系数更趋于3,偏度系数更趋于0,这些参数均表明质控后对探空观测数据的数据质量带来了明显的改善。
本实施例提供了一种探空观测数据的评估和修正方法,该方法针对需要修正的可疑数据或者错误数据,通过对应的模式背景数据进行替换,使得修正后的数据和模式背景数据能够更好的吻合。
更进一步地,在上述各实施例的基础上,所述分别根据每一所述评价参数判断在所述气压层下观测得到的探空观测数据是否属于可疑层数据,包括:
若所述平均偏差大于第三阈值,则判定在所述气压层下观测得到的探空观测数据属于可疑层数据;
若所述标准偏差大于第四阈值,则判定在所述气压层下观测得到的探空观测数据属于可疑层数据;
若所述均方根误差大于第五阈值,则判定在所述气压层下观测得到的探空观测数据属于可疑层数据;
若所述相关系数小于第六阈值,则判定在所述气压层下观测得到的探空观测数据属于可疑层数据;
若所述偏度系数大于第七阈值,则判定在所述气压层下观测得到的探空观测数据属于可疑层数据;
若所述峰度系数大于第八阈值,则判定在所述气压层下观测得到的探空观测数据属于可疑层数据。
更进一步地,在上述各实施例的基础上,所述根据所述第二探空观测数据和所述目标模式背景数据计算在所述气压层下观测的探空观测数据的至少一个评价参数,包括:
计算每一所述第二探空观测数据和与其相同观测条件下所对应的目标模式背景数据的第二差值,将所有所述第二探空观测数据对应的第二差值的平均值作为在所述气压层下观测的探空观测数据的平均偏差;
根据公式计算在所述气压层下观测的探空观测数据的标准偏差;其中,S1为所述标准偏差,dn为第n个第二探空观测数据对应的所述第二差值,d0为所述平均偏差,N为第二探空观测数据的数据点总个数;
根据公式计算在所述气压层下观测的探空观测数据的均方根误差;其中,S2为所述均方根误差,On为第n个第二探空观测数据,Bn为对应于On的模式背景数据;
对所述第二探空观测数据和所述目标模式背景数据进行拟合,得到所述相关系数;
根据所述第二探空观测数据的数据分布对称性得到所述偏度系数;
根据所述第二探空观测数据的数据聚集在中心数据的聚集程度得到所述峰度系数。
更进一步地,在上述各实施例的基础上,所述根据计算出的所述气压层对应的所有评价参数,对所述第一探空观测数据的质量进行评估,包括:
根据所述平均偏差评价所述第一探空观测数据偏离所述目标模式背景数据的偏差平均值;
根据所述标准偏差评价所述第一探空观测数据偏离所述目标模式背景数据的偏差相对于所述平均偏差的离散程度;
根据所述均方根误差评价所述第一探空观测数据偏离所述目标模式背景数据的偏离程度;
根据所述相关系数评价所述第一探空观测数据与所述目标模式背景数据的相关程度;
根据所述偏度系数评价所述第一探空观测数据的数据分布偏离对称性的程度;
根据所述峰度系数评价所述第一探空观测数据的数据聚集在中心数据的聚集程度。
例如,收集全国高空报文数据、探空秒数据、中国GRAPES模式背景场和美国GFS模式背景场等资料,对高空报文数据进行解码、数据补算、与探空秒数据相互验证等一系列预处理,统计到报情况,并基于中国GRAPES模式背景场、美国GFS模式背景场等对高空报文数据进行质量评估(其中,模式背景数据为B,探空观测数据为O)。具体评估方法包括:
1、计算第一差值
计算每个时刻的观测余差(O-B)。
2、针对涉及到的所有气压层中每一气压层计算评价参数
计算一段时间(月、季、年)内的参数:平均偏差、标准偏差、相关系数、均方根误差、概率密度分布的峰度系数和偏度系数,利用这些参数对高空观测数据进行评估,经过使用一年数据进行试验,统计选取最佳评估标准。
(1)平均偏差:查看一段时间内探空数据的平均质量。通过统计计算选择一个最佳评估标准。
其中,气温界限值:1℃<Tdvavg则说明该气压层的探空观测数据为可疑数据。湿度界限值:第一个对流层顶及其以下:5%<RHdvavg则明该气压层的探空观测数据为可疑数据;第一个对流层顶以上:10%<RHdvavg明该气压层的探空观测数据为可疑数据。
(2)标准偏差:标准偏差代表每个时刻的误差相对平均偏差的离散程度,标准偏差越小,说明每个时刻的误差越靠近平均偏差。
在平均偏差的基础上观察标准偏差廓线,针对平均偏差中找到的可疑点和错误点观察对应的标准偏差。如果标准偏差特别大,则说明这段时间内的一部分点或者个别点是错误的或者可疑的,需要每个时刻都检查,将对应的可疑点进行订正或者剔除;如果标准偏差特别小,则说明这段时间内的对应层的点都有问题,考虑需要全部剔除或者订正。其中温度的标准偏差的分割线为1.5℃,湿度的标准偏差的分割线为15%。
在平均偏差和标准偏差计算完毕后,对第一差值进行可疑值和错误值的标识。在数学统计中,误差分布一般遵循正态分布。根据正态分布的特征规律,统计学中规定与平均值的偏差超过两倍标准差的测定值称为异常值,与平均值的偏差超过三倍标准差的测定值,称为高度异常的异常值。在统计检验时,指定为检出异常值的显著性水平α=0.05,称为检出水平;指定为检出高度异常的异常值的显著性水平α=0.01,称为舍弃水平。因此,在处理数据时,应剔除高度异常的异常值。异常值是否剔除,视具体情况而定。由此,试验中将选择dvavg±2S1作为单个时刻温度的可疑值阈值,dvavg±3S1作为单个时刻温度的错误值阈值。
(3)均方根误差:探空观测值与模式值偏差的平方与观测次数n比值的平方根,即多个时刻计算,每层等压面有一个均方根误差。用来衡量每层等压面中高空观测数据偏离模式数据的程度。均方根误差越小,说明探空数据越靠近模式数据,说明探空数据质量越高。经过统计计算后给均方根误差定一个阈值,温度为1.5,S<=1.5认为数据质量较好,没有问题,S>1.5超过这个值认为该层数据存在问题,需要对该层的各个时刻的数据进行检查剔除或者订正可疑点;湿度为20,S<=20认为数据质量较好,没有问题,S>20超过这个值认为该层数据存在问题,需要对该层的各个时刻的数据进行检查剔除或者订正可疑点。
(4)相关系数:多个时刻计算,每层等压面有一个相关系数。相关系数|r|≤1,当越大,说明该层等压面的探空数据和模式数据相关程度越大,说明探空数据的质量越好。经过统计计算后给相关系数的绝对值定一个最小值,温度为0.9,湿度为0.7,小于最小值认为该层数据有问题,需要对该层的各个时刻的数据进行检查剔除或者订正可疑点。
绘制这段时间内每一层等压面偏差的概率密度分布图(PDF),结合对应的峰度系数和偏度系数进行评估分析。
(5)偏度系数(skew):偏度系数是描述分布偏离对称性程度的一个特征数。当分布左右对称时,偏度系数为0。当偏度系数大于0时,即重尾在右侧时,该分布为右偏。当偏度系数小于0时,即重尾在左侧时,该分布左偏。在探空数据的评估中,当skew>0时,说明该层等压面的探空数据比模式数据大的时间比较多;当skew<0时,说明该层等压面的探空数据比模式数据小的时间比较多。确定一个阈值为±1,如果超过这个范围,认为数据存在问题。
(6)峰度系数(kurtosis):峰度系数用来度量数据在中心聚集程度。在正态分布情况下,峰度系数值是3。由于计算峰度系数的函数(一年北京探空站数据试验)已经减去3,所以kurtosis>0的峰度系数说明观察量更集中,有比正态分布更短的尾部;kurtosis<0的峰度系数说明观测量不那么集中,有比正态分布更长的尾部,类似于矩形的均匀分布。在探空数据的评估质控中,温度为0±3.5,湿度为峰度系数0±4,如果超过这个范围,认为数据存在问题。
通过计算上述中的各个参数,对探空数据质量有一个定性定量的评估,在评估的基础上,根据评估结果进行质控:首先将每个时刻的错误数据进行剔除和订正,然后考虑可疑数据的处理方法,将所有的参数综合考虑。如果所有的参数都认为某一层的数据可疑,则将该层数据中的可疑数据剔除或者订正,如果只是其中几个参数认为可疑,则不剔除或订正。将质控之后的探空观测数据再用上述方法评估来判断质控效果。本实施例提供的方法为气象局评估探空业务质量提供科学依据,利用这个技术方法开发高空观测数据质控和评估的国家级软件,该软件具有基于模式初估场对高空观测数据质量控制与评估,以及定量化分析全网观测数据误差分布特征的能力,能够初步建立国家级准实时评估业务。
作为一个更为具体的示例,本实施例提供了A地点在2016.7.1-2017.6.30一年的高空观测资料,以FNL再分析资料和GRAPES预报场为背景场,使用上述方法进行评估指控再评估。将这一年的数据按照季节分为春夏秋冬四个数据块,再根据08时和20时一共分为八个数据块进行试验,通过试验确定了各个参数指标。若是扩大指标,则质控效果不好,若是缩小指标,则错误点太多,剔除之后不能体现探空数据本身的性质,没有意义。以秋季(9月-11月)08时的温度评估质控(FNL为背景场)为例说明。
图5是2017年9月1日00时的温度数据和第一差值的示意图,参见图5,图中示出了可疑数据(即dvavg+2S1<ΔT<dvavg+3S1和dvavg-3S1<ΔT<dvavg-2S1的数据)和错误数据(ΔT>dvavg+3S1和ΔT<dvavg-3S1的数据)。由图可知,整体上该时刻的一致性很好,探空观测数据和模式背景数据基本吻合,且很好的体现了大气温度廓线的性质。其中没有错误点,待定的一共有三个可疑点(这3个可疑点对应的气压层为950hPa,900hPa和700hPa的点)。对于其它任一时刻均可按照此图所示的方法处理。对可疑数据的处理方式包括:
图6为各气压层对应的平均偏差的示意图,由图可知,整个秋季来看,探空数据的质量很不错,误差基本在0℃-0.8℃之间,只有四个可疑点,且都很接近1℃,分别是950hPa,925hPa,900hPa,7hPa,其中950hPa和之前单个时刻9月1日00时(UTC)重合,由此判断9月1日00时的六个可疑点只有950hPa的点存在问题,需要继续考虑,其余五个可疑点保留不变。平均偏差中其他几层虽然没有和9月1日的可疑数据重合,但是都需要考虑,因为秋季其他时间的数据在这几层存在可疑点。从平均偏差判断,所有层次中,除了这几层(950hPa,925hPa,900hPa,7hPa)以外,其它层次的可疑点直接保留,不再考虑。这几层可疑点暂时留着,根据后面的参数来确定怎样处理。
图7为各气压层对应的标准偏差的示意图,关注平均偏差中标出可疑点的位置,950hPa的标准偏差为1.36<1.5,考虑认为950hPa这一气压层整个秋季的90个00时(UTC)数据中可疑点非常可疑,结合后面的几个参数来确定是否剔除或者订正,即9月1日00时(UTC)950hPa的温度数据列入重点考虑范围。发现7hPa的标准偏差为1.6>1.5,观测余差离散度较大,可能只是错误点造成的,考虑认为整个秋季的90个00时(UTC)数据中只有少部分有问题,这一气压层90个温度数据中的可疑数据继续保留,不再考虑这一层的可疑点。925hPa,900hPa,300hPa的标准偏差都小于1.5,和950hPa的情况相同,需要在后面的参数中重点考虑这几层的温度可疑数据,确定否将整个秋季这几层的可疑点剔除或者订正。即现在需要考虑订正或者剔除的是950hPa,925hPa,900hPa和300hPa上各层90个温度数据中的可疑点。
图8为每一层偏差的概率密度分布(PDF)示意图,参见图8,950hPaPDF的偏度系数:0.85>0,结合上面的PDF图,得出重尾在右侧时,该分布为右偏。说明该层等压面的探空数据比模式数据大的时间比较多,该层的平均偏差也验证了这一点。
950hPaPDF的峰度系数:0<1.90<3.5,结合上面的PDF图,发现PDF的形状符合“尖峰”形态。说明该层数据十分集中,非常靠近平均偏差(平均偏差可疑)。PDF图和其峰度系数、偏度系数都在表明这一层的数据存在问题。同样的方法去分析925hPa,900hPa和300hPa的PDF和峰度,偏度系数,判断这几层的可疑点是否需要继续关注。
950hPa、925hPa,900hPa和300hPa的均方根误差分别是1.70、1.46、1.60和0.88。其中前三层大于1,300hPa的小于1,说明除了300hPa以外三层的可疑点还需要继续考虑。
950hPa、925hPa和900hPa的相关系数都是0.99,大于0.9,说明这几层的数据相关性都非常好,认为这几层的可疑点不需要剔除。
将所有的参数综合考虑,如果所有的参数都认为某一层的数据可疑,则将该层数据中的可疑点剔除或者订正,如果只是其中几个参数认为可疑,则不剔除或订正。
在此次试验中,将错误点和筛选出来认为需要修改的可疑点进行了订正,没有选择剔除。订正方法是直接用那些点的模式资料代替探空资料
将进行过质控数据再次计算一次方案二中的参数。如图9中的平均偏差和图10中的标准偏差所示,950hPa PDF的偏度系数为-0.03,更接近0,峰度系数为-0.59,更接近0,两者都说明质控后的误差分布更符合正态分布,说明质量变得更好。
总体来说,质控前后的评估进行对比,发现质控后这些参数变得更加合理,更贴近理论。这些都表明探空数据的质量得到了提升。说明这种评估质控是有效,有意义的。
本实施例提供的方法和现有技术相比,从更多的角度评估了探空观测数据质量。根据均方根误差和相关系数,可以更好地衡量探空数据和背景场的相关关系,均方根误差越小,说明探空数据越靠近模式数据,说明探空数据质量越高,相关系数越大,该层等压面的探空数据和模式数据相关程度越大,说明探空数据的质量越好。根据偏差的概率密度分布以及对应的偏度系数和峰度系数能够分析误差分布的特征和数据特点。总体来说,对于探空数据质量的了解不再是单薄的平均偏差和偏离偏差的程度,能够从更多的角度,更全面的认识数据质量,误差分布特征和数据特点。
另一方面,相对于现有技术,这个评估方法不再是定性地评估探空质量,而是根据一年的数据试验统计结果科学地给定了确定的评估标准。为气象局评估探空业务质量提供科学依据。并在这个基础上开发了高空观测数据质控和评估的国家级软件,该软件具有基于模式初估场对高空观测数据质量控制与评估,以及定量化分析全网观测数据误差分布特征的能力,初步建立国家级准实时评估业务。
图11示出了本发明的实施例提供的一种探空观测数据的评估和修正装置的结构框图,参见图11,本实施例提供的探空观测数据的评估和修正装置,包括第一获取模块1101、第二获取模块1102、第三获取模块1103、计算模块1104和处理模块1104,其中,
第一获取模块1101,用于获取与欲进行评估的目标探空观测数据有关的所有气压层;
第二获取模块1102,用于针对每一气压层,从所述目标探空观测数据中获取在所述气压层下观测得到的探空观测数据,作为第一探空观测数据,并获取观测所述第一探空观测数据的第一观测时间;
第三获取模块1103,用于获取在包括所述第一观测时间的第二观测时间观测到的对应于所述气压层的探空观测数据,作为第二探空观测数据,以及与每一第二探空观测数据对应的模式背景数据,作为目标模式背景数据;
计算模块1104,用于根据所述第二探空观测数据和所述目标模式背景数据计算在所述气压层下观测的探空观测数据的至少一个评价参数;
处理模块1105,用于根据计算出的所述气压层对应的所有评价参数,对所述第一探空观测数据的质量进行评估,并对所述第一探空观测数据进行修正;
其中,所述评价参数包括平均偏差、标准偏差、均方根误差、相关系数、偏度系数和峰度系数。
本实施例提供的探空观测数据的评估和修正装置适用于上述实施例中提供的探空观测数据的评估和修正方法,在此不再赘述。
本发明的实施例提供了一种探空观测数据的评估和修正装置,该装置针对欲进行评估的目标探空观测数据中涉及到的每一气压层,分别通过在该气压层观测到的第二探空观测数据和相应的模式背景数据,计算出用于评价该气压层的探空观测数据的质量的评价参数,通过该评价参数对目标探空观测数据的质量进行评估。同时,针对目标探空观测数据中处于该气压层的探空观测数据,根据上述评价参数和根据探空观测数据和相应的模式背景数据的差值对处于该气压层的探空观测数据进行修正,以改善探空观测数据的质量。该装置针对同一气压层的探空观测数据,通过上述评价参数共同评价该层的探空观测数据的质量,评估结果更为准确。同时,依据这些评价参数对探空观测数据进行修正,使得探空观测数据更为合理,便于将这些数据应用到实验或者气象分析中。
第三方面,图12是示出本实施例提供的电子设备的结构框图。
参照图12,所述电子设备包括:包括:处理器(processor)1201、存储器(memory)1202和总线1203;
其中,
所述处理器1201、存储器1202通过所述总线1203完成相互间的通信;
所述处理器1201用于调用所述存储器1202中的程序指令,以执行上述各方法实施例所提供的方法,例如包括:获取与欲进行评估的目标探空观测数据有关的所有气压层;针对每一气压层,从所述目标探空观测数据中获取在所述气压层下观测得到的探空观测数据,作为第一探空观测数据,并获取观测所述第一探空观测数据的第一观测时间;获取在包括所述第一观测时间的第二观测时间观测到的对应于所述气压层的探空观测数据,作为第二探空观测数据,以及与每一第二探空观测数据对应的模式背景数据,作为目标模式背景数据;根据所述第二探空观测数据和所述目标模式背景数据计算在所述气压层下观测的探空观测数据的至少一个评价参数;根据计算出的所述气压层对应的所有评价参数,对所述第一探空观测数据的质量进行评估,并对所述第一探空观测数据进行修正;其中,所述评价参数包括平均偏差、标准偏差、均方根误差、相关系数、偏度系数和峰度系数。
第四方面,本实施例提供一种非暂态计算机可读存储介质,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令使所述计算机执行上述各方法实施例所提供的方法,例如包括:获取与欲进行评估的目标探空观测数据有关的所有气压层;针对每一气压层,从所述目标探空观测数据中获取在所述气压层下观测得到的探空观测数据,作为第一探空观测数据,并获取观测所述第一探空观测数据的第一观测时间;获取在包括所述第一观测时间的第二观测时间观测到的对应于所述气压层的探空观测数据,作为第二探空观测数据,以及与每一第二探空观测数据对应的模式背景数据,作为目标模式背景数据;根据所述第二探空观测数据和所述目标模式背景数据计算在所述气压层下观测的探空观测数据的至少一个评价参数;根据计算出的所述气压层对应的所有评价参数,对所述第一探空观测数据的质量进行评估,并对所述第一探空观测数据进行修正;其中,所述评价参数包括平均偏差、标准偏差、均方根误差、相关系数、偏度系数和峰度系数。
本实施例公开一种计算机程序产品,所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序,所述计算机程序包括程序指令,当所述程序指令被计算机执行时,计算机能够执行上述各方法实施例所提供的方法,例如,包括:获取与欲进行评估的目标探空观测数据有关的所有气压层;针对每一气压层,从所述目标探空观测数据中获取在所述气压层下观测得到的探空观测数据,作为第一探空观测数据,并获取观测所述第一探空观测数据的第一观测时间;获取在包括所述第一观测时间的第二观测时间观测到的对应于所述气压层的探空观测数据,作为第二探空观测数据,以及与每一第二探空观测数据对应的模式背景数据,作为目标模式背景数据;根据所述第二探空观测数据和所述目标模式背景数据计算在所述气压层下观测的探空观测数据的至少一个评价参数;根据计算出的所述气压层对应的所有评价参数,对所述第一探空观测数据的质量进行评估,并对所述第一探空观测数据进行修正;其中,所述评价参数包括平均偏差、标准偏差、均方根误差、相关系数、偏度系数和峰度系数。
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于一计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
以上所描述的电子设备等实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是:以上各实施例仅用以说明本发明的实施例的技术方案,而非对其限制;尽管参照前述各实施例对本发明的实施例进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明的实施例各实施例技术方案的范围。
Claims (9)
1.一种探空观测数据的评估和修正方法,其特征在于,包括:
获取与欲进行评估的目标探空观测数据有关的所有气压层;
针对每一气压层,从所述目标探空观测数据中获取在所述气压层下观测得到的探空观测数据,作为第一探空观测数据,并获取观测所述第一探空观测数据的第一观测时间;
获取在包括所述第一观测时间的第二观测时间观测到的对应于所述气压层的探空观测数据,作为第二探空观测数据,以及与每一第二探空观测数据对应的模式背景数据,作为目标模式背景数据;
根据所述第二探空观测数据和所述目标模式背景数据计算在所述气压层下观测的探空观测数据的至少一个评价参数;
根据计算出的所述气压层对应的所有评价参数,对所述第一探空观测数据的质量进行评估,并对所述第一探空观测数据进行修正;
其中,所述评价参数包括平均偏差、标准偏差、均方根误差、相关系数、偏度系数和峰度系数。
2.根据权利要求1所述的方法,其特征在于,所述根据计算出的所述气压层对应的所有评价参数对所述第一探空观测数据进行修正,包括:
针对每一所述第一探空观测数据,计算所述第一探空观测数据和与其相同观测条件对应的目标模式背景数据的第一差值,若所述第一差值大于第一阈值,则判定所述第一探空观测数据为错误数据,若所述第一差值小于或等于所述第一阈值,大于第二阈值,则判定所述第一探空观测数据为可疑数据;
分别根据每一所述评价参数判断在所述气压层下观测得到的探空观测数据是否属于可疑层数据,若根据每一所述评价参数均判断在所述气压层下观测得到的探空观测数据属于可疑层数据,则对所述可疑数据和所述错误数据进行修正,否则,对所述错误数据进行修正;
其中,所述第一阈值为所述气压层对应的平均偏差与3倍的标准偏差的差或者和,且所述第二阈值为所述气压层对应的平均偏差与2倍的标准偏差的差;或者,所述第一阈值为所述气压层对应的平均偏差与3倍的标准偏差的差或者和,且所述第二阈值为所述气压层对应的平均偏差与2倍的标准偏差的和。
3.根据权利要求2所述的方法,其特征在于,对所述错误数据或者所述可疑数据进行修正,包括:
获取欲进行修正的所述错误数据或者所述可疑数据对应的观测条件,从所述目标模式背景数据中获取对应于所述观测条件下的模式背景数据,作为替换数据,用所述替换数据替换所述错误数据或者所述可疑数据。
4.根据权利要求2所述的方法,其特征在于,所述分别根据每一所述评价参数判断在所述气压层下观测得到的探空观测数据是否属于可疑层数据,包括:
若所述平均偏差大于第三阈值,则判定在所述气压层下观测得到的探空观测数据属于可疑层数据;
若所述标准偏差大于第四阈值,则判定在所述气压层下观测得到的探空观测数据属于可疑层数据;
若所述均方根误差大于第五阈值,则判定在所述气压层下观测得到的探空观测数据属于可疑层数据;
若所述相关系数小于第六阈值,则判定在所述气压层下观测得到的探空观测数据属于可疑层数据;
若所述偏度系数大于第七阈值,则判定在所述气压层下观测得到的探空观测数据属于可疑层数据;
若所述峰度系数大于第八阈值,则判定在所述气压层下观测得到的探空观测数据属于可疑层数据。
5.根据权利要求1所述的方法,其特征在于,所述根据所述第二探空观测数据和所述目标模式背景数据计算在所述气压层下观测的探空观测数据的至少一个评价参数,包括:
计算每一所述第二探空观测数据和与其相同观测条件下所对应的目标模式背景数据的第二差值,将所有所述第二探空观测数据对应的第二差值的平均值作为在所述气压层下观测的探空观测数据的平均偏差;
根据公式计算在所述气压层下观测的探空观测数据的标准偏差;其中,S1为所述标准偏差,dn为第n个第二探空观测数据对应的所述第二差值,d0为所述平均偏差,N为第二探空观测数据的数据点总个数;
根据公式计算在所述气压层下观测的探空观测数据的均方根误差;其中,S2为所述均方根误差,On为第n个第二探空观测数据,Bn为对应于On的模式背景数据;
对所述第二探空观测数据和所述目标模式背景数据进行拟合,得到所述相关系数;
根据所述第二探空观测数据的数据分布对称性得到所述偏度系数;
根据所述第二探空观测数据的数据聚集在中心数据的聚集程度得到所述峰度系数。
6.根据权利要求1所述的方法,其特征在于,所述根据计算出的所述气压层对应的所有评价参数,对所述第一探空观测数据的质量进行评估,包括:
根据所述平均偏差评价所述第一探空观测数据偏离所述目标模式背景数据的偏差平均值;
根据所述标准偏差评价所述第一探空观测数据偏离所述目标模式背景数据的偏差相对于所述平均偏差的离散程度;
根据所述均方根误差评价所述第一探空观测数据偏离所述目标模式背景数据的偏离程度;
根据所述相关系数评价所述第一探空观测数据与所述目标模式背景数据的相关程度;
根据所述偏度系数评价所述第一探空观测数据的数据分布偏离对称性的程度;
根据所述峰度系数评价所述第一探空观测数据的数据聚集在中心数据的聚集程度。
7.一种探空观测数据的评估和修正装置,其特征在于,包括:
第一获取模块,用于获取与欲进行评估的目标探空观测数据有关的所有气压层;
第二获取模块,用于针对每一气压层,从所述目标探空观测数据中获取在所述气压层下观测得到的探空观测数据,作为第一探空观测数据,并获取观测所述第一探空观测数据的第一观测时间;
第三获取模块,用于获取在包括所述第一观测时间的第二观测时间观测到的对应于所述气压层的探空观测数据,作为第二探空观测数据,以及与每一第二探空观测数据对应的模式背景数据,作为目标模式背景数据;
计算模块,用于根据所述第二探空观测数据和所述目标模式背景数据计算在所述气压层下观测的探空观测数据的至少一个评价参数;
处理模块,用于根据计算出的所述气压层对应的所有评价参数,对所述第一探空观测数据的质量进行评估,并对所述第一探空观测数据进行修正;
其中,所述评价参数包括平均偏差、标准偏差、均方根误差、相关系数、偏度系数和峰度系数。
8.一种电子设备,其特征在于,包括:
至少一个处理器;以及
与所述处理器通信连接的至少一个存储器,其中:
所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行如权利要求1至6中任一项所述的方法。
9.一种非暂态计算机可读存储介质,其特征在于,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令使所述计算机执行如权利要求1至6任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810359842.4A CN108830444B (zh) | 2018-04-20 | 2018-04-20 | 一种探空观测数据的评估和修正方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810359842.4A CN108830444B (zh) | 2018-04-20 | 2018-04-20 | 一种探空观测数据的评估和修正方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108830444A true CN108830444A (zh) | 2018-11-16 |
CN108830444B CN108830444B (zh) | 2020-10-09 |
Family
ID=64154889
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810359842.4A Active CN108830444B (zh) | 2018-04-20 | 2018-04-20 | 一种探空观测数据的评估和修正方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108830444B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110989046A (zh) * | 2019-12-25 | 2020-04-10 | 中国气象局气象探测中心 | 一种锚碇浮标站数据质量控制方法系统 |
CN111830451A (zh) * | 2019-04-19 | 2020-10-27 | 中国科学院光电研究院 | 一种非成像传感器的检验方法 |
CN113592127A (zh) * | 2020-04-30 | 2021-11-02 | 国家气象中心 | 一种往返式探空系统轨迹预测方法及装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020107638A1 (en) * | 2000-04-18 | 2002-08-08 | Intriligator Devrie S. | Space weather prediction system and method |
CN101853290A (zh) * | 2010-05-25 | 2010-10-06 | 南京信息工程大学 | 基于地理信息系统的气象服务效益评估方法 |
CN102183800A (zh) * | 2011-02-28 | 2011-09-14 | 中南大学 | 基于大气混合层高度参数的城市通风状况的预测方法 |
CN102521403A (zh) * | 2011-12-26 | 2012-06-27 | 南京成风大气信息技术有限公司 | 精细化气象信息服务系统 |
CN104777488A (zh) * | 2015-03-13 | 2015-07-15 | 中国科学院上海天文台 | 天顶对流层延迟建模方法、装置及测量方法、装置 |
WO2015150747A1 (en) * | 2014-03-31 | 2015-10-08 | Imperial Innovations Limted | A computer implemented method of deriving performance from a local model |
CN105426970A (zh) * | 2015-11-17 | 2016-03-23 | 武汉理工大学 | 一种基于离散动态贝叶斯网络的气象威胁评估方法 |
CN106547840A (zh) * | 2016-10-13 | 2017-03-29 | 国家卫星气象中心 | 一种全球三维大气数据的解析及管理方法 |
-
2018
- 2018-04-20 CN CN201810359842.4A patent/CN108830444B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020107638A1 (en) * | 2000-04-18 | 2002-08-08 | Intriligator Devrie S. | Space weather prediction system and method |
CN101853290A (zh) * | 2010-05-25 | 2010-10-06 | 南京信息工程大学 | 基于地理信息系统的气象服务效益评估方法 |
CN102183800A (zh) * | 2011-02-28 | 2011-09-14 | 中南大学 | 基于大气混合层高度参数的城市通风状况的预测方法 |
CN102521403A (zh) * | 2011-12-26 | 2012-06-27 | 南京成风大气信息技术有限公司 | 精细化气象信息服务系统 |
WO2015150747A1 (en) * | 2014-03-31 | 2015-10-08 | Imperial Innovations Limted | A computer implemented method of deriving performance from a local model |
CN104777488A (zh) * | 2015-03-13 | 2015-07-15 | 中国科学院上海天文台 | 天顶对流层延迟建模方法、装置及测量方法、装置 |
CN105426970A (zh) * | 2015-11-17 | 2016-03-23 | 武汉理工大学 | 一种基于离散动态贝叶斯网络的气象威胁评估方法 |
CN106547840A (zh) * | 2016-10-13 | 2017-03-29 | 国家卫星气象中心 | 一种全球三维大气数据的解析及管理方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111830451A (zh) * | 2019-04-19 | 2020-10-27 | 中国科学院光电研究院 | 一种非成像传感器的检验方法 |
CN111830451B (zh) * | 2019-04-19 | 2023-04-18 | 中国科学院光电研究院 | 一种非成像传感器的检验方法 |
CN110989046A (zh) * | 2019-12-25 | 2020-04-10 | 中国气象局气象探测中心 | 一种锚碇浮标站数据质量控制方法系统 |
CN113592127A (zh) * | 2020-04-30 | 2021-11-02 | 国家气象中心 | 一种往返式探空系统轨迹预测方法及装置 |
CN113592127B (zh) * | 2020-04-30 | 2024-02-02 | 国家气象中心 | 一种往返式探空系统轨迹预测方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN108830444B (zh) | 2020-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110634080B (zh) | 异常用电检测方法、装置、设备及计算机可读存储介质 | |
CN112819207B (zh) | 基于相似性度量的地质灾害空间预测方法、系统及存储介质 | |
CN108830444A (zh) | 一种探空观测数据的评估和修正方法及装置 | |
CN104903262A (zh) | 玻璃基板生产管理系统以及玻璃基板生产管理方法 | |
CN114077970B (zh) | 一种基于城市形态的碳排放相关因子的确定方法及装置 | |
CN115830030B (zh) | 一种石英晶片的外观质量评估方法及系统 | |
CN103902798B (zh) | 数据预处理方法 | |
CN111612371A (zh) | 一种基于层次分析法的智能电表质量评估方法 | |
CN111914424B (zh) | 一种基于短期测风资料的设计风速取值方法及系统 | |
JP6192432B2 (ja) | リスク計量システム | |
CN108446213A (zh) | 一种静态代码质量分析方法和装置 | |
CN115796585A (zh) | 一种企业经营风险评估方法及系统 | |
CN115906669A (zh) | 一种顾及负样本选取策略的密集残差网络滑坡易发性评价方法 | |
CN113689153B (zh) | 一种基于灰靶决策的复杂装备交付中典型问题的筛选方法 | |
CN115204501A (zh) | 企业评估方法、装置、计算机设备和存储介质 | |
CN115292887A (zh) | 基于病斑面积比率的植物病害严重度评估方法及装置 | |
CN105069467B (zh) | 车辆油位数据的筛选方法和装置 | |
CN109284478B (zh) | 一种估计对数正态型单元可靠性参数的方法 | |
CN107578168A (zh) | 一种用于缺陷库移植的方法、装置及电子设备 | |
CN106301880A (zh) | 一种确定网络关系稳定度、互联网业务推荐方法和设备 | |
CN110363374A (zh) | 一种不合格产品影响因素的定量分析方法 | |
Siok et al. | Empirical study of embedded software quality and productivity | |
CN112070415B (zh) | 基于软件定量管理过程构建的场景化基线与模型实现方法 | |
CN106570625A (zh) | 服务缺陷管理方法 | |
US20240193460A1 (en) | Data processing method and data processing apparatus |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |