1、基于双变量模型的子带自适应滤波方法侯建华 熊承义（中南民族大学电子信息工程学院,武汉 430074）摘 要 利用贝叶斯最大后验估计理论研究了双变量层间模型、模型的子带自适应参数估计方法，推导了对应的萎缩函数；在此基础上设计了一种子带自适应图像滤波方法，并与两种基于模型的类似算法进行了比较。实验结果表明：利用小波系数层间模型可以有效地改善图像去噪质量。关键词 小波系数；最大后验估计；双变量模型；子带自适应中图分类号 TP391 文献标识码 A 文章编号 1672-4321(2005)04-0042-03A Sub-Band Adaptive Image Filtering MethodBased
2、 on Bivariate ModelHou Jianhua Xiong ChengyiAbstract The paper has investigated the bivariate model using Bayesian maximum a posteriori (MAP) estimation theory. The sub-band adaptive method for parameter estimation is discussed and the corresponding shrinkage function is derived. A sub-band adaptive
3、 image filtering method is designed and compared with other two similar algorithms which are also based on statistics modeling of wavelet coefficients. Experiment results show that image denoising can be improved by exploiting wavelet interscale models.Keywords wavelet coefficients; maximum a poster
4、iori estimation; bivariate model; sub-band adaptiveHou Jianhua Assoc Prof, College of Electronics and Information Engineering, SCUFN, Wuhan 430074,China大量研究表明，图像小波系数之间存在着一定的相关性，利用这种相关性可以显著改善有关算法的性能。具体地讲有 2 点：(1)小波系数具有空间聚集性(Clustering)，即如果某一系数较大( 或小)，则与其相邻的系数也往往较大(或小)；(2)小波系数具有尺度传递性(Persistence) ，即大(
5、或小)的系数有沿尺度传播的趋向 1。依据这些相关性和小波系数的非高斯特性，已提出了多种图像小波系数的先验分布模型，例如广义高斯分布 2、拉普拉斯分布 3这些层内模型利用了子带内的空间聚集性，在图像编码、去噪领域得到了成功的应用。另一类层间模型则根据小波系数的尺度传递性,对“父” “子”系数间的依赖性进行描述。例如 Shapiro 提出的小波零树编码算法是应用该特性的成功典范 4，Crouse 的小波域隐马尔科夫树模型在信号估计、图像分割、纹理分析等方面也表现出了出色的性能 1。近年来，Sendur 提出了一种双变量模型，用非高斯密度函数对“父”“子”系数分布进行建模，利用贝叶斯统计理论得到了其
6、最大后验估计的解析表达式 5。本文对这种双变量模型进行了深入研究，并设计了一种子带自适应图像滤波方法。设 w1 是当前的小波系数，w 2 是上一个粗糙尺度中空间位置与 w1 相同的系数,即 w1 的“父”系数。y1 = w1 + n1 ，y 2 = w2 + n2 (1)用矢量表示为： y = w + n (2)其中 y = ( y1 , y2 ) , w = ( w1 , w2 ) , n = (n1 , n2 )。为了更好地利用小波系数的相关性，Sendur 等统计了大量自然图像小波系数直方图中“父” “子”系数间的概率分布，提出了一种双变量混合模型来刻画这种关系5：, (3)212w3e
7、xp)( wp2 基于双变量模型的最大后验估计2.1 贝叶斯最大后验估计理论为简单起见，先考虑标量形式的 MAP 估计。设 y=w+n 分别代表观测到的含噪图像、真实图像以及独立同分布高斯噪声的小波系数。n 服从 N(0, )分布。2n我们的目的是从观测系数 y 中得到真实系数 w 的估计 。)(y设 pw(w)、p n(n)分别是真实系数 w、噪声系数 n 的概率密度函数 (pdf)。在贝叶斯统计理论中。最大后验估计(MAP)是一种常用的方法。MAP 估计即在给定y 条件下，使后验概率密度 pw|y(w|y)取最大的 w：, (4)|(maxrg)(|yw利用贝叶斯准则，得到。 (5)(ma
8、xrg)(|(r)(| pypy wnwwyw , (6)(w)lgp)-(ylgpmax r(y)wnw, (7)2e21pnn定义： 。 (8)()lgf()w将（7） （8）代人（6）得：， (9)(2(maxr)( fyywnw， (10)0)(2fn利用 pw(w)得到 ,进而求出方程(10)的解,即得对应的 MAP 估计 (y)。f w2.2 双变量萎缩函数利用噪声独立同分布假设：。 (11)2exp21)(1nnnp将矢量 y、w、n 以及(11)式代入(6) 式：， (12)(2()(maxrg)(21wfynn上式的等价方程为：(13)0)(121wfyn(14)(22fn其
9、中 分别代表 对 w1、w 2 的一阶导数：1wf2)(f(15)213)(f(16)2123)(wf将(15)、(16)代入(13)、(14)式,经简单推导，得：(17)121213yywn其中符号(g) +表示：(18)otherisg,0)(17)式就是在双变量模型下当前系数 w1 的 MAP 估计，它实质上也是一种软阈值萎缩函数，称为双变量萎缩函数；但与 Donoho 的经典软阈值萎缩规则不同，在做阈值处理时还考虑了“父”系数的影响。3 子带自适应滤波算法为了计算 MAP 估计表达式 (17)，要已知模型(3)式中的 2 个参数 和 。其中2n噪声方差可采用 Donoho 提出的鲁棒性
10、中值估计方法 :the finest subband HH (19)iinymeda,6745.0)(模型的边缘方差 2 则以子带 S 为单位，先估计含噪系数的标准差：, (20)syiyiN2其中 N2 是子带 S 中的系数个数。则：。 （21）)0,max(2ny基于双变量模型的子带自适应滤波算法描述如下：(1)对含噪图像进行小波变换；(2)对每个细节子带按(19) (21)式估计 和 按(17)式系数进行软阈值2n处理；(3)小波反变换。由于模型的边缘方差是以子带为单位进行估计的,故该算法具有子带自适应性。4 实验结果采用 3 幅 512512 Lena、Boat 和 Bridge 图像
11、为实验对象，噪声水平分别为 10、15、20、25,用 Daubechies8 小波进行 4 尺度分解。参与比较实验的去噪方法有 BayesShrink 算法 2,6和 MapShrink 算法 3,它们分别采用广义高斯分布 6、拉普拉斯分布，均为子带自适应方法；本文方法用 BiShrink 表示。去噪性能用峰值信噪比 PSNR 来评价，实验结果见表 1，其中黑体部分代表同一条件下的最好结果。从表 1 可看出,在这 3 种基于模型的子带自适应算法中，从去噪后图像的峰值信噪比改善程度看，本文的 BiShrink 最好，其次是BayesShrink、MapShrink。相对 ayesShrink
12、和 MapShrink，本文方法的 PSNR 增益最高可达 0.22 和 0.71dB。5 结语利用小波系数的统计模型可以改善图像处理的质量。双变量模型能够较好地刻画“父” “子”系数间的层间相关性，本文的实验结果证明了这一点。由于BiShrink 算法只考虑了 “父”系数对当前系数的影响，模型的边缘方差根据子带自适应调整，算法简洁；若进一步考虑子带内的空间聚集性，可得到基于混合模型的局部自适应算法,去噪性能将得到更大提高，但算法复杂度也会因此而增加。参 考 文 献1 Crouse M S, Nowak R D, Baraniuk R G. Wavelet-based statistical
13、signal processing using hiddenMarkov models J. IEEE Trans Signal Processing,1998, 46 (4) : 8869022 Chang S G, Yu B, Vetterli M. Adaptive wavelet thresholding for image denoising and compression J. IEEE Trans Image Processing, 2000, 9 (9):1 5321 5463 Moulin P, Liu J. Analysis of multiresolution image
14、 denoising schemes using generalized-gaussian and complexity priors J. IEEE Trans Information Theory,1999,45 (3):909 9194 Shapiro J M, “Embedded image coding using zerotrees of wavelet coefficients”, IEEE Trans Signal Processing, 1993, 41 (12) ,344534625 Sendur L, Selesnick I W. Bivariate shrinkage
15、func-tions for wavelet-based denoising exploiting interscale dependency J. IEEE Trans Signal Processing,2002, 50 (11) : 2 7442 7566 侯建华,熊承义 ,田晓梅.广义高斯分布及其在图像去噪中的应用J.中南民族大学学报(自然科学版),2005, 24 (3) : 4447Goosen B,Pizurica A, “Removal of correlated noise by modeling the signal of interest in the wavelet d
16、omain”,IEEE Trans, 2009,6(18),1153-1165Yu Yue, Xue Yonglin, Li Fengtinget, “Noise image coding using soft-threshold based denoising”, Journal of Image and Graphics(in Chinese), 2001, 6(1), 4650David L Donoho, Iain M Johnstone, “Ideal spatial adaptation by wavelet shrinkage”, Biometrika, 1994, 81(3), 425455
