专利类型 | 发明授权 | 法律事件 | 公开; 实质审查; 授权; |
专利有效性 | 有效专利 | 当前状态 | 授权 |
申请号 | CN202210830796.8 | 申请日 | 2022-07-15 |
公开(公告)号 | CN114911789B | 公开(公告)日 | 2022-10-11 |
申请人 | 武汉大学; | 申请人类型 | 学校 |
发明人 | 黄煜; 陈华; 武泽恒; | 第一发明人 | 黄煜 |
权利人 | 武汉大学 | 权利人类型 | 学校 |
当前权利人 | 武汉大学 | 当前权利人类型 | 学校 |
省份 | 当前专利权人所在省份:湖北省 | 城市 | 当前专利权人所在城市:湖北省武汉市 |
具体地址 | 当前专利权人所在详细地址:湖北省武汉市武昌区八一路299号 | 邮编 | 当前专利权人邮编:430072 |
主IPC国际分类 | G06F16/215 | 所有IPC国际分类 | G06F16/215 ; G06F16/2458 ; G01F23/28 ; G01F23/292 ; G01S13/86 |
专利引用数量 | 1 | 专利被引用数量 | 0 |
专利权利要求数量 | 6 | 专利文献类型 | B |
专利代理机构 | 武汉科皓知识产权代理事务所 | 专利代理人 | 俞琳娟; |
摘要 | 本 发明 提供基于雷视多 传感器 数据的 水 位数据融合方法及装置,既能对传感器测量值进行有效性验证,剔除异常观测数据,同时又能充分考虑各数据源之间的离散性和相关性,实现全自动、自适应、高准确水位监测。水位数据融合方法包括:步骤1.采集待研究区域基于视频与雷达的多传感器水位原始数据集;步骤2.构建原始数据集理论函数和其约束条件,并对数据集样本进行可疑值检验和剔除,得到优选数据集S;步骤3.对S进行 频率 统计,计算熵值,并构建相对熵矩阵H;步骤4.对S进行相关计算,并构建相关性矩阵R;步骤5.基于H和R构建权重系数矩阵ωT;步骤6.基于ωT构建归一化系数矩阵W,将W与S子集矩阵做乘积处理,得到最终融合结果。 | ||
权利要求 | 1.基于雷视多传感器数据的水位数据融合方法,其特征在于,包括以下步骤: |
||
说明书全文 | 基于雷视多传感器数据的水位数据融合方法及装置技术领域背景技术[0002] 水位是水动力因素中最基本、最重要的观测要素。传统的水位测量仪器主要包括水尺,压力式、气泡式、浮子式水位计等。随着科学技术的发展,人们对新时期流域动态化、精细化和智慧化管理需求不断提高,更加先进的非接触式测量方法,如视频、雷达等开始应用于水位监测中。视频依据计算机视觉和数字图像处理技术对水位进行智能识别,对实时水位监控和构建河流水位监测预警系统有着重大意义,但摄像头传感器在采集数据时受天气和光照的影响较大,难以适应复杂多变的环境条件,鲁棒性较差。雷达根据观测回波来获取与目标的距离,进而推算出水位,其主要优势是全天候观测,不仅可以在灯光昏暗的夜间进行工作,而且可以在雨雪雾等复杂天气情况下进行工作,但测量精度易受波浪和水面漂浮物等影响。可见,单一传感器虽然各具优势,但受限于自身构造,在可探测范围、可探测对象及所获得的数据类型上存在一定局限性,且易受监测环境噪声干扰,测量结果稳定性较低,误差不可避免。此外,当传感器受到瞬时干扰、发生缓变故障、永久性故障等情况时,可能出现误检和漏跟等情况,无法胜任越来越复杂的应用场景。 [0003] 为了解决视频与雷达单一传感器测水位时存在的信息量少、抗干扰能力弱、鲁棒性差等问题,人们开始尝试将二者进行传感器融合。多传感器融合技术,通过对各个传感器的信息进行排列整合、冗余互补和优化组合,可得到所监测目标全方位、多层次、高准确率的观测信息,能扩展测量感知范围,增强系统自适应能力。目前,在雷达和视频水位数据融合算法方面的研究较少,最常用的方法为算术平均法,即简单地将各传感器数据进行累加取平均值。此种方法虽然简单,但难以得到可靠结果。这是因为一方面,其未考虑各传感器偶尔输出的异常、明显偏离真实值的数据,若将异常观测和正常数据进行融合通常会得到不准确的估计,甚至使系统崩溃。另一方面,由于各类传感器提供的多源信息与数据可能具有各不相同的特性或特征,这些特征可能是多样的,相互支持或互补的,确定或者模糊的,也可能是冲突矛盾的,具有差异性,不同传感器所被赋予的权重也应该有所不同,不能同一而论。 发明内容[0004] 本发明是为了解决上述问题而进行的,目的在于提供基于雷视多传感器数据的水位数据融合方法及装置,既能对传感器测量值进行有效性验证,剔除异常观测数据,同时又能充分考虑各数据源之间的离散性和相关性,实现全自动、自适应、高准确水位监测。 [0005] 本发明为了实现上述目的,采用了以下方案: [0006] 本发明提供基于雷视多传感器数据的水位数据融合方法,其特征在于,包括以下步骤: [0007] 步骤1.采集待研究区域基于视频与雷达的多传感器水位原始数据集o,o={o1(x)、o2(x)、o3(x)、…、om(x)},oi(x)为第i个传感器监测的数据子集,oi(x)=[x1、x2、x3、…、xn],xi为oi(x)的样本,m为数据子集个数,n为子集中数据样本个数; [0008] 步骤2.构建原始数据集理论函数和其约束条件,并对数据子集样本进行可疑值检验和剔除,得到优选数据集S; [0009] 对于任一子集oi(x),根据其所包含的样本,将其理论函数表示为: [0010] (2‑1) [0011] 式中,y表示数据子集的理论函数;αnt=[α1、α2、…、αnt]T为待定系数,βnt(x)=[β1T(x)、β2(x)、…、βnt(x)]为一组线性无关的函数;nt为一组中线性无关函数的个数; [0012] 为保证理论函数最优,其目标函数应满足: [0013] (2‑2) [0014] 目标函数的约束条件为: [0015] (2‑3) [0016] 式中,μ是待求解的常数;xi为oi(x)的样本; [0017] 可疑值检验和剔除过程为,首先建立目标函数的拉格朗日函数: [0018] (2‑4) [0019] 式中,φ为约束函数; [0020] 然后对拉格朗日函数进行求解: [0021] (2‑5) [0022] 验证原始数据集o各子集中各个样本是否满足上式,将不满足的剔除,将满足的保留记为xi’,构成优选子集Si(x’);即原始数据子集oi(x)=[x1、x2、x3、…、xn]经过可疑值检验和剔除后将变为优选子集Si(x’)=[x1’、x2’、x3’、…、xN’],xi’为优选子集Si(x’)的样本,N为优选子集中样本个数;所有优选子集形成优选数据集S,S={S1(x’)、S2(x’)、S3(x’)、…、Sm(x’)}; [0023] 步骤3.对优选数据集S进行频率统计,计算熵值,并构建相对熵矩阵H; [0024] 步骤4.对优选数据集S进行相关计算,并构建相关性矩阵R; [0025] 步骤5.基于相对熵矩阵H和相关性矩阵R构建权重系数矩阵ωT; [0026] 步骤6.基于权重系数矩阵ωT构建归一化系数矩阵W,并将W与S子集矩阵做乘积处理,得到最终融合结果。 [0027] 优选地,本发明提供的基于雷视多传感器数据的水位数据融合方法,还可以具有这样的特征:在步骤3中,频率统计为:记Si(x’)中任意x’i出现的次数为k,则其频率P(x’i)=k/N, [0028] 数据集的熵值为: [0029] (3‑1) [0030] 对于任意两个数据子集Si(x’)和Sj(x’),i≠j,其相对熵为: [0031] (3‑2) [0032] 对所有优选子集进行相对熵计算,并构建相对熵矩阵H: [0033] (3‑3)。 [0034] 优选地,本发明提供的基于雷视多传感器数据的水位数据融合方法,还可以具有这样的特征:在步骤4中,对于任意两个数据子集Si(x’)和Sj(x’)其相关计算为: [0035] (4‑1) [0036] 式中, ; [0037] 所有优选子集的相关性矩阵R为: [0038] (4‑2)。 [0039] 优选地,本发明提供的基于雷视多传感器数据的水位数据融合方法,还可以具有T这样的特征:在步骤5中,权重系数矩阵ω 为: [0040] (5‑1) [0041] 其归一化处理为: [0042] (5‑2)。 [0043] 优选地,本发明提供的基于雷视多传感器数据的水位数据融合方法,还可以具有这样的特征:在步骤6中,归一化系数矩阵W为: [0044] (6‑1) [0045] 融合结果表达式为: [0046] (6‑2)。 [0047] 进一步,本发明还提供基于雷视多传感器数据的水位数据融合装置,其特征在于,包括: [0048] 采集部,采集待研究区域基于视频与雷达的多传感器水位原始数据集o,o={o1(x)、o2(x)、o3(x)、…、om(x)},oi(x)为第i个传感器监测到的数据子集,oi(x)=[x1、x2、x3、…、xn],xi为oi(x)的样本,m为数据子集个数,n为子集中数据样本个数; [0049] 筛选部,构建原始数据集理论函数和其约束条件,并对数据子集样本进行可疑值检验和剔除,得到优选数据集S; [0050] 对于任一子集oi(x),根据其所包含的样本,将其数据集理论函数表示为: [0051] (2‑1) [0052] 式中,y表示数据子集的理论函数;αnt=[α1,α2,…,αnt]T为待定系数,βnt(x)=[β1T(x)、β2(x)、…、βnt(x)]为一组线性无关的函数;nt为一组中线性无关函数的个数; [0053] 为保证理论函数最优,其目标函数应满足: [0054] (2‑2) [0055] 目标函数的约束条件为: [0056] (2‑3) [0057] 式中,μ是待求解的常数;xi为oi(x)的样本; [0058] 可疑值检验和剔除过程为,首先建立目标函数的拉格朗日函数: [0059] (2‑4) [0060] 式中,φ为约束函数; [0061] 然后对拉格朗日函数进行求解: [0062] (2‑5) [0063] 验证原始数据集o各子集中各个样本是否满足上式,将不满足的剔除,将满足的保留记为xi’,构成优选子集Si(x’);即原始数据子集oi(x)=[x1、x2、x3、…、xn]经过可疑值检验和剔除后将变为优选子集Si(x’)=[x1’、x2’、x3’、…、xN’],xi’为优选子集Si(x’)的样本,N为优选子集中样本个数;所有优选子集形成优选数据集S,S={S1(x’)、S2(x’)、S3(x’)、…、Sm(x’)}; [0064] 相对熵矩阵构建部,对优选数据集S进行频率统计,计算熵值,并构建相对熵矩阵H; [0065] 相关性矩阵构建部,对优选数据集S进行相关计算,并构建相关性矩阵R; [0066] 权重系数矩阵构建部,基于相对熵矩阵H和相关性矩阵R构建权重系数矩阵ωT; [0067] 融合部,基于权重系数矩阵ωT构建归一化系数矩阵W,并将W与S子集矩阵做乘积处理,得到最终融合结果; [0068] 控制部,与采集部、筛选部、相对熵矩阵构建部、相关性矩阵构建部、权重系数矩阵构建部和融合部均通信相连,控制它们的运行。 [0069] 优选地,本发明所提供的基于雷视多传感器数据的水位数据融合装置,还可以具有这样的特征:在相对熵矩阵构建部中,记Si(x’)中任意x’i出现的次数为k,则其频率P(x’i)=k/N, [0070] 数据集的熵值为: [0071] (3‑1) [0072] 对于任意两个数据集Si(x’)和Sj(x’),i≠j,其相对熵为: [0073] (3‑2) [0074] 对所有优选子集进行相对熵计算,并构建相对熵矩阵H: [0075] (3‑3)。 [0076] 优选地,本发明所提供的基于雷视多传感器数据的水位数据融合装置,还可以具有这样的特征:在相关性矩阵构建部中,对于任意两个数据子集Si(x’)和Sj(x’)其相关计算为: [0077] (4‑1) [0078] 式中, ; [0079] 所有优选子集的相关性矩阵R为: [0080] (4‑2)。 [0081] 优选地,本发明所提供的基于雷视多传感器数据的水位数据融合装置,还可以具T有这样的特征:在权重系数矩阵构建部中,权重系数矩阵ω 为: [0082] (5‑1) [0083] 其归一化处理为: [0084] (5‑2)。 [0085] 优选地,本发明所提供的基于雷视多传感器数据的水位数据融合装置,还可以包括:输入显示部,与采集部、筛选部、相对熵矩阵构建部、相关性矩阵构建部、权重系数矩阵构建部、融合部和控制部均通信相连,用于让用户输入操作指令,并进行相应显示。 [0086] 本发明所提供的基于雷视多传感器数据的水位数据融合方法及装置,能充分利用和挖掘视频水位数据与雷达水位数据之间的深层次联系,考虑各数据源之间的离散性和相关性,自动对异常数据进行检测和剔除,自动优选原始数据集,根据优选数据集分别独立地计算相对熵值和相关值,然后将优选数据集组成一个相对熵矩阵、一个相关矩阵,最后基于相对熵矩阵和相关矩阵构建权重系数矩阵并完成各类数据的权重赋值,得到融合结果,达到水位结果的精确输出,对对象先验知识无要求,具有较强的自学能力和适应性。附图说明 [0088] 图2为本发明实施例涉及的装置的部分结构示意图; [0089] 图3为本发明实施例涉及的可疑值检验和剔除过程示意图。 具体实施方式[0090] 以下结合附图对本发明涉及的基于雷视多传感器数据的水位数据融合方法及装置的具体实施方案进行详细地说明。 [0091] 如图1所示,本实施例提供的基于雷视多传感器数据的水位数据融合方法包括如下步骤: [0092] 步骤1.采集待研究区域基于视频与雷达的多传感器水位原始数据集o,o={o1(x),o2(x)、o3(x)、…、om(x)},oi(x)为第i个传感器监测的数据子集,oi(x)=[x1、x2、x3、…、xn],xi为oi(x)的样本,m为数据子集个数,n为子集中数据样本个数; [0093] 步骤2.构建原始数据集理论函数和其约束条件,并对数据集样本进行可疑值检验和剔除,得到优选数据集S。 [0094] 对于任一子集oi(x),根据其所包含的所有样本,将其理论函数表示为: [0095] (2‑1) [0096] 式中,y表示数据子集的理论函数;αnt=[α1、α2、…、αnt]T为待定系数,βnt(x)=[β1T(x)、β2(x)、…、βnt(x)]为一组线性无关的函数;nt为一组中线性无关函数的个数。 [0097] 为保证理论函数最优,其目标函数应满足: [0098] (2‑2) [0099] 目标函数的约束条件为: [0100] (2‑3) [0101] 式中,μ是待求解的常数;xi为oi(x)的样本。 [0102] 可疑值检验和剔除过程为,首先建立目标函数的拉格朗日函数: [0103] (2‑4) [0104] 式中,φ为约束函数。 [0105] 然后对拉格朗日函数进行求解: [0106] (2‑5) [0107] 验证原始数据集o各子集中各个样本是否满足上式,将不满足的剔除,将满足的保留在相应子集中记为xi’,构成优选子集Si(x’);即原始数据子集oi(x)=[x1、x2、x3、…、xn]经过可疑值检验和剔除后将变为优选子集Si(x’)=[x1’、x2’、x3’、…、xN’],xi’为优选子集Si(x’)的样本,N为优选子集中样本个数;所有优选子集形成优选数据集S,S={S1(x’)、S2(x’)、S3(x’)、…、Sm(x’)}。 [0108] 步骤3.对优选数据集S进行频率统计,计算熵值,并构建相对熵矩阵H。 [0109] 频率统计为:记Si(x’)中任意x’i出现的次数为k,则其频率P(x’i)=k/N,[0110] 数据集的熵值为: [0111] (3‑1) [0112] 对于任意两个数据子集Si(x’)和Sj(x’),i≠j,其相对熵为: [0113] (3‑2) [0114] 对所有优选子集进行相对熵计算,并构建相对熵矩阵H: [0115] (3‑3)。 [0116] 当i≠j时,相对熵矩阵中的元素由相对熵值的倒数组成。 [0117] 步骤4.对优选数据集S进行相关计算,并构建相关性矩阵R。 [0118] 对于任意两个数据子集Si(x’)和Sj(x’)其相关计算为: [0119] (4‑1) [0120] 式中, 。 [0121] 所有优选子集的相关性矩阵R为: [0122] (4‑2)。 [0123] 步骤5.基于相对熵矩阵H和相关性矩阵R构建权重系数矩阵ωT。 [0124] 具体地,权重系数矩阵ωT为: [0125] (5‑1) [0126] 其归一化处理为: [0127] (5‑2)。 [0128] 步骤6.基于权重系数矩阵ωT构建归一化系数矩阵W,并将W与S子集矩阵做乘积处理,得到最终融合结果。 [0129] 归一化系数矩阵W为: [0130] (6‑1) [0131] 融合结果表达式为: [0132] (6‑2)。 [0133] 为了更好的解释及验证本发明的可靠性及有效性,引用中国发明专利公开号为CN110389971A中的试验数据作为原始数据集,并采用本实施例算法与CN110389971A中涉及算法进行对比。原始数据集,如下表1所示: [0134] 表1多传感器观测数据(真值900) [0135] [0136] 对表1中的数据,分别使用CN110389971A中最优融合集法、均值法、可靠性算法与本文算法进行数据融合,结果见表2: [0137] 表2本文算法与其他算法融合结果对比 [0138] [0139] 由上表2可以看出本实施例方法、最优融合集法、均值法和可靠性算法与真实值900的绝对误差分别为1.0957、1.9083、2.7389和2.7182。很明显,本实施例方法的绝对误差相对其他三种算法均更小,证明了本发明方法的优越性和可靠性。 [0140] 进一步,本实施例还提供能够自动实现上述方法的基于雷视多传感器数据的水位数据融合装置,该装置包括采集部、筛选部、相对熵矩阵构建部、相关性矩阵构建部、权重系数矩阵构建部、融合部、输入显示部、控制部。 [0141] 采集部执行上文步骤1所描述的内容,采集待研究区域基于视频与雷达的多传感器水位原始数据集o,o={o1(x)、o2(x)、o3(x)、…、om(x)},oi(x)为第i个传感器监测到的数据子集,oi(x)=[x1、x2、x3、…、xn],xi为oi(x)的样本,m为数据子集个数,n为子集中数据样本个数。 [0142] 筛选部执行上文步骤2所描述的内容,构建数据集理论函数和其约束条件,并对数据集样本进行可疑值检验和剔除,得到优选数据集S。 [0143] 相对熵矩阵构建部执行上文步骤3所描述的内容,对优选数据集S进行频率统计,计算熵值,并构建相对熵矩阵H。 [0144] 相关性矩阵构建部执行上文步骤4所描述的内容,对优选数据集S进行相关计算,并构建相关性矩阵R。 [0145] 权重系数矩阵构建部执行上文步骤5所描述的内容,基于相对熵矩阵H和相关性矩T阵R构建权重系数矩阵ω 。 [0146] 融合部执行上文步骤6所描述的内容,基于权重系数矩阵ωT构建归一化系数矩阵W,并将W与S子集矩阵做乘积处理,得到最终融合结果。 [0147] 输入显示部与采集部、筛选部、相对熵矩阵构建部、相关性矩阵构建部、权重系数矩阵构建部、融合部均通信相连,用于让用户输入操作指令,并对相应部获取和处理得到的数据以文字、列表、静态图或者动态趋势图等方式进行相应显示。 [0148] 控制部与采集部、筛选部、相对熵矩阵构建部、相关性矩阵构建部、权重系数矩阵构建部、融合部、输入显示部均通信相连,控制它们的运行。 |