[0004] 本发明针对现有技术的不足,提出了一种基于图滤波器的两通道近似正交图滤波器组频率选择性误差交替优化设计方法。该方法具体包含以下步骤:
[0005] 步骤一、根据设计要求,分别确定图滤波器组中低通滤波器h0(x)的通带截止频率ωp,阻带截止频率ωs,提升滤波器的系数sn、tn的最大长度Ls、Lt。
[0006] 步骤二、由提升滤波器的系数sn、tn求出低通滤波器h0(x)与低通滤波器g0(x)的表达式,进而求出此时滤波器组的频率选择性误差以及重构误差。
[0007] 2‑1.由提升滤波器的系数sn、tn求出低通滤波器h0(x)与低通滤波器g0(x)的表达式低通滤波器h0(x)与低通滤波器g0(x)的表达式分别为:
[0008]
[0009]
[0010] 其中s(x)和t(x)分别为:
[0011]
[0012]
[0013] 2‑2.确定近似正交图滤波器组的完美重构条件
[0014] 图滤波器组的理想重构条件为:在x∈[0,2]的频域上满足:当x=0时,当x=1时,h0(x)=1。同时,对于近似正交图滤波器组结构,低通滤波器h0(x)、g0(x)与高通滤波器h1(x)、g1(x)存在着严格的等式关系:
[0015] h1(x)=g0(2‑x),g1(x)=h0(2‑x)
[0016] 对于图滤波器组结构,完美重构条件为:
[0017] h0(x)g0(x)+h1(x)g1(x)=c2
[0018] 令 并将低通滤波器与高通滤波器之间的等式关系代入到上式,得到对于近似正交图滤波器组的完美重构条件公式:
[0019] h0(x)g0(x)+h0(2‑x)g0(2‑x)=2
[0020] 2‑3.求近似正交图滤波器组的实际重构误差
[0021] 实际设计过程中往往并不能达到完美重构条件,因此两者之间存在一个差值。定义两者之间的差值E(x)为实际重构误差:
[0022] E(x)=|h0(x)g0(x)+h0(2‑x)g0(2‑x)‑2|
[0023] 在以重构误差数值为核心比较图滤波器组性能优劣的情况下,只需要取重构误差中最大的值进行比较,即最大重构误差Emax:
[0024] Emax=||E(x)||∞
[0025] 2‑4.求近似正交图滤波器组的实际频率选择性误差
[0026] 由于设计出的滤波器通带与阻带应平坦才能达到很好的滤波效果,因此在考虑滤波器组指标时频率选择性也是一个重要衡量指标。在给定低通滤波器h0(x)通带截止频率ωp,阻带截止频率ωs后,低通滤波器h0(x)的通带误差ε_h0_pass和阻带误差ε_h0_stop,g0(x)的通带误差ε_g0_pass和阻带误差ε_g0_stop分别为:
[0027]
[0028] ε_h0_stop=‖[h0(xs=ws) h0(xs+1) … h0(2)]T‖∞
[0029]
[0030] ε_g0_stop=‖[g0(xs=ws) g0(xs+1) … g0(2)]T‖∞
[0031] 步骤三、设计交替优化算法
[0032] 3‑1.设计优化目标函数为:
[0033]
[0034] 其中s、t分别为提升滤波器初始系数 该目标函数是一个最小最大化问题,同时也是一个线性问题,可以使用CVX函数进行优化求解。
[0035] 3‑2.固定s系数,优化t系数
[0036] 先固定初始系数 对系数 进行优化,优化目标函数可以写为:
[0037]
[0038] 其中 为经过一次优化以后新的t系数,使用CVX函数进行优化后,令新的系数替代原来的系数。
[0039] 3‑3.固定t系数,优化s系数
[0040] 再固定优化好的系数 对系数 进行优化,优化目标函数可以写为:
[0041]
[0042] 其中 为经过一次优化以后新的s系数,使用CVX函数进行优化后,令新的系数替代原来的系数。
[0043] 步骤四、将最终输出的最大重构误差值以及最大频率选择性误差值作为最终优化结果,同时输出交替优化结束时得到的提升滤波器系数s、t,优化结束。
[0044] 本发明具有以下有益效果:
[0045] 相较于以往通过牺牲频率选择性误差值来获取较好重构误差值的滤波器组设计方法,本发明在已经具备极佳重构误差值的基础上,优化图滤波器组的频率选择性误差值,同时保证重构误差值小的指标能够非常好的保留,在两个指标中找到了一个极佳的平衡点,使得优化出的图滤波器组不仅在滤波性能上同时也在图像美观程度上有了进一步的改观。