基于互素多项式的两通道非下采样图滤波器组的设计方法技术

本发明专利技术公开一种基于互素多项式的两通道非下采样图滤波器组的设计方法，其主要从滤波器组的完全重构条件考虑，基于互素多项式进行非下采样图滤波器组的设计，设计所得的两通道图滤波器组能够应用于任意图。实验表明，本发明专利技术设计所得的非下采样图滤波器组的重构性能更好，且具备较好的去噪性能。

【技术实现步骤摘要】
基于互素多项式的两通道非下采样图滤波器组的设计方法
本专利技术涉及图信号处理
，具体涉及一种基于互素多项式的两通道非下采样图滤波器组的设计方法。
技术介绍
图信号处理是分析非规则域数据的有效工具。它弥补了传统信号处理在处理非规则域数据中的不足。如何将传统信号处理中滤波器组的概念应用到图信号处理中引起了广大研究者的兴趣。Narang和Ortega首次提出了临界采样图小波滤波器组。随后他们后来又提出了双正交图小波滤波器组，其中分析滤波器和综合滤波器可以具有频域紧支撑。这两种设计都是为二分图定义的。对于任意图而言，它需要被分解成几个二分图，并且将滤波器组分别应用于这些二分图，其中滤波器组的数量取决于从输入图分解的二分图的数量。为了分析单个滤波器组中的图信号，Saikiyama和Tanaka提出了一种关于图信号过采样的技术。有学者针对循环图提出了一种样条图小波滤波器组的设计方法。为了提供图多分辨率分析，H.Q.Nguyen和M.N.Do提出了最大生成树方法来重新连接定义下采样顶点。N.Tremblay和P.Borgnat提出了类似Haar小波滤波器组的滤波器组设计，该滤波器组的设计基于将输入图划分为连接的子图。除了图小波滤波器之外，一些研究还以其他方式设计图滤波器组。在一些研究中，图的频谱被分割成几个子带，在分析部分，图信号在每个子带上进行滤波，然后进行下采样。D.I.Shuman等人提出了拉普拉斯金字塔变换并应用于图信号，这种变换可以提供图和图信号的多尺度分析。A.G.Marques等人提出了一种称为聚合采样的独特采样方法，这种方法修复了一个节点，并考虑了该节点上信号的采样，因为递归运算符是递归应用的。最近，Teke和Vaidyanathan在M块循环图上提出了图信号的M通道最大抽取滤波器组，通过利用相似变换，他们进一步使他们的设计适用于任意图。然而，上述提到的所有滤波器组的设计都包含图信号的采样运算，它们是很难考虑图结构的重构，其主要问题是分析后下采样图的结构信息不能用于重构输入图。输入图的拓扑结构必须在整个过程中进行记录和维护，并且对于现有的图滤波器组设计会导致巨大的存储价格。
技术实现思路
本专利技术针对目前图信号处理中难以准确定义一般图的图信号下采样运算的问题，提供一种基于互素多项式的两通道非下采样图滤波器组的设计方法。为解决上述问题，本专利技术是通过以下技术方案实现的：基于互素多项式的两通道非下采样图滤波器组的设计方法，具体包括步骤如下：步骤1，根据给定的图，计算该图的归一化拉普拉斯矩阵该归一化拉普拉斯矩阵的行数和列数均为N；步骤2，根据步骤1中的归一化拉普拉斯矩阵设计两通道非下采样图滤波器组的分析滤波器组，该分析滤波器组包括低通分析滤波器H0和高通分析滤波器H1，其中步骤3，根据两通道非下采样图滤波器组的完全重构条件和贝祖恒等式，基于互素多项式设计两通道非下采样图滤波器组的综合滤波器组，该综合滤波器组包括低通综合滤波器G0和高通综合滤波器G1，其中式中，IN表示N×N的单位阵，n表示滤波器的阶数。上述方案中，低通分析滤波器H0和高通分析滤波器H1均为线性滤波器。上述方案中，低通综合滤波器G0和高通综合滤波器G1均为线性滤波器。上述方案中，滤波器的阶数n≥1。本专利技术主要从滤波器组的完全重构条件考虑，基于互素多项式进行非下采样图滤波器组的设计，设计所得的两通道图滤波器组能够应用于任意图。实验表明，本专利技术设计所得的非下采样图滤波器组的重构性能更好，且具备较好的去噪性能。附图说明图1为两通道非下采样图滤波器组的结构。图2为图滤波器的幅度响应。其中(a)为本专利技术中n＝1时图滤波器的幅度响应；(b)为本专利技术中n＝2时图滤波器的幅度响应。图3为本专利技术中Minnesota交通图去噪仿真结果。其中(a)为输入信号；(b)为噪声信号(σ＝1/2)；(c)为基于现有方法1(临界采样双正交图滤波器组)设计的图滤波器组去噪后的输出信号；(d)为基于现有方法2(M通道过采样图滤波器组)设计的图滤波器组去噪后的输出信号；(e)为本专利技术中n＝1时设计所得的两通道非下采样图滤波器组去噪后的输出信号；(f)为本专利技术n＝2时设计所得的两通道非下采样图滤波器组去噪后的输出信号。图4为本专利技术中实测的美国温度数据去噪仿真结果。其中(a)为输入信号；(b)为噪声信号(σ＝10)；(c)为基于现有方法1(临界采样双正交图滤波器组)设计的图滤波器组去噪后的输出信号；(d)为本专利技术n＝2时设计所得的两通道非下采样图滤波器组去噪后的输出信号。具体实施方式为使本专利技术的目的、技术方案和优点更加清楚明白，以下结合具体实例，并参照附图，对本专利技术进一步详细说明。在一些实际的应用中，一个图滤波器A可以是关于归一化拉普拉斯矩阵的多项式根据图1中表示的两通道非下采样图滤波器组，其中两个线性滤波器H0,H1组成分析滤波器组，综合滤波器组由两个线性滤波器G0,G1组成。图信号f为图滤波器组的输入信号，输出信号为f0,f1分别表示低频子带系数和高频子带系数。记两通道非下采样图滤波器组的输入输出关系可表示为T是图滤波器组的传递函数。两通道非下采样图滤波器组的完全重构条件可写为G0H0+G1H1＝IN(4)根据样条图小波的定义，非下采样图滤波器组的分析滤波器组可表示为H0,H1分别为分析滤波器组的低通和高通滤波器，n表示滤波器的阶数。则有根据上述的设计的非下采样图滤波器组的分析滤波器组，记综合滤波器组为根据Bezout(贝祖)恒等式当且仅当多项式P0,P1没有共同的根时，Bezout(贝祖)恒等式(10)有解，利用基于互素多项式设计综合滤波器组根据上述设计出的分析和综合滤波器组，以下证明设计的非下采样图滤波器组满足重构条件。对于n≥1时，完全重构条件满足基于以上原理，本专利技术所设计的一种基于互素多项式的两通道非下采样图滤波器组的设计方法，其具体包括步骤如下：步骤1，根据给定的图，计算该图的归一化拉普拉斯矩阵该归一化拉普拉斯矩阵的行数和列数均为N，即图的归一化拉普拉斯矩阵为一个N×N矩阵；步骤2，根据步骤1中的归一化拉普拉斯矩阵设计两通道非下采样图滤波器组的分析滤波器组，该分析滤波器组包括低通分析滤波器H0和高通分析滤波器H1，其中步骤3，根据两通道非下采样图滤波器组的完全重构条件和贝祖恒等式，基于互素多项式设计两通道非下采样图滤波器组的综合滤波器组，该综合滤波器组包括低通综合滤波器G0和高通综合滤波器G1，其中式中，IN表示N×N的单位阵，n表示滤波器的阶数，n≥1。下面通过具体仿真实例，对本专利技术的性能进行说明。仿真实例1：本专利技术中，根据基于互素多项式设计两通道非下采样图滤波器组。并以常用的Minnesota图信号作为输入信号，当n＝1时，其非下采样滤波器组的幅度响应如图2(a)所示，重构信噪比SNR＝331.70dB。当n＝2时，其非下采样滤波器组的幅度响应如图2(b)所示，重构信噪比SNR＝310.87dB。实验结果表明基于互素多项式设计所得的两通道非下采样图滤波器组具备完全重构特性，当n＝1时，滤波器组具备更好的重构特性。表1给出了本专利技术与现有方法设计的带采样的图滤波器组的重构性能，与带采样的滤波器组现有方法1(临界采样双正交图滤波器组)和现有方法2(M通道过采样图本文档来自技高网...

【技术保护点】
1.基于互素多项式的两通道非下采样图滤波器组的设计方法，其特征是，具体包括步骤如下：步骤1，根据给定的图，计算该图的归一化拉普拉斯矩阵

【技术特征摘要】
1.基于互素多项式的两通道非下采样图滤波器组的设计方法，其特征是，具体包括步骤如下：步骤1，根据给定的图，计算该图的归一化拉普拉斯矩阵该归一化拉普拉斯矩阵的行数和列数均为N；步骤2，根据步骤1中的归一化拉普拉斯矩阵设计两通道非下采样图滤波器组的分析滤波器组，该分析滤波器组包括低通分析滤波器H0和高通分析滤波器H1，其中步骤3，根据两通道非下采样图滤波器组的完全重构条件和贝祖恒等式，基于互素多项式设计两通道非下采样图滤波器组的综合滤波器组，该综合滤波器组包括低通...

【专利技术属性】
技术研发人员：蒋俊正，杨圣，欧阳缮，杨杰，赵海兵，李杨剑，李龙斌，
申请(专利权)人：桂林电子科技大学，
类型：发明
国别省市：广西,45