专利名称：基于双树离散小波包的图像去噪方法
技术领域：
本发明涉及图像处理技术领域，特别涉及一种基于双树离散小波包的图像去噪方法。
背景技术：
近年来，由于小波变换具有良好的时-频局部化特性，在信号和图像去噪领域得到了广泛的应用。传统的小波域去噪方法是对小波系数进行萎缩处理，如Donoho提出的硬阈值和软阈值去噪法。现有的方法存在的缺点是一方面，硬阈值函数具有不连续性，重构所得的信号会产生伪吉布斯效应，而软阈值方法估计后的小波系数和分解得到的小波系数总存在恒定的偏差，直接影响重构信号与真实信号的逼近程度；另一方面，在某些动态环境例如无人机自主飞行中，动态获取的图像不仅包含大量的噪声，同时自然场景下的图像富含方向性特征， 更多得高频细节和高频噪声很难区分，给小波去噪带来了较大困难。
发明内容
本发明的目的旨在至少解决上述技术缺陷之一。为达到上述目的，本发明提出一种基于双树离散小波包的图像去噪方法，包括以下步骤构造双树离散小波；根据所述双树离散小波对图像进行分解，获得多个子频带以及所述子频带的小波系数；计算各个子频带的小波系数的自相关系数，并选取所述小波系数的自相关系数小于预定的第一阈值的子频带进行小波包分解，获得第二层双树离散小波包；以及基于所述第二层双树离散小波包将所述图像变换到小波域中，并计算所述小波系数的邻域相关性，以及根据所述小波系数的邻域相关性对所述图像进行去噪和图像增强。根据本发明实施例的基于双树离散小波包的图像去噪方法，通过利用小波系数之间的相似性关系分析图像的噪声分布特征特性，设计出的小波包构造方案在实现自然场景去噪的情况下，能够同时实现图像增强。而且，本发明的方法只对源图像进行处理，不需要任何先验知识，通用性强。本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解，其中图1为本发明实施例的基于双树离散小波包的图像去噪方法的流程图；图2为本发明一个实施例的三级双树离散小波变换结构示意图；以及图3为本发明一个实施例的满树小波包分解结构示意图。
具体实施例方式下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能解释为对本发明的限制。如图1所示，根据本发明实施例的基于双树离散小波包的图像去噪方法，包括以下步骤步骤S101，构造双树离散小波。在本发明的一个实施例中，使用q-shift设计方案构造双树离散小波。如图2所示为一个示例的三级双树离散小波变换结构，其中，hob(n) = h0a(N-l-n)，N为滤波器长度为偶数，即hQb(n)是Iltla(Ii)的反转。自然地，hob(n)和Iltla(Ii)的傅里叶变换长度相等，但是两者不能达到半样本延迟， 即|Η0“ω) I = |H0a(co)Z hob (η) = Z Htla (ω) - (N-I) ω。为了实现半样本延迟，令 Zhtlb (η) =-ZHcia(G) )-0.5 ω，则有Z Ηθ3(ω)-0. 5ω = - Z Ηθ3(ω)-(Ν"1) ωZ Ηθ3(ω) =-0.5 (Ν_1) ω +0. 25 ω也就是说，Htla (ω)是近似线性相位滤波器，它的对称中心是η = 0. 5 (N-I) -0. 25， 其比一般线性相位滤波器的对称点偏移了四分之一个样本。双树离散小波近似地实现了两路滤波器的半样本点延迟，获得了近似解析的复小波变换，有诸多优点，是应用最为广泛的解析复小波变换，并且解析精度随分解级数的增加而增加。但是，双树离散小波分解与普通的离散小波变换一样，分解的时候采用二进式分解，对高频子带不做进一步分解，通过对低频子带进行迭代分解，信号被分解为一个低频子带和一系列高频子带。分片线性光滑信号主要包含低频成分，所以双树离散小波能非常紧致地表示这一类信号除了低频子带之外，高频子带仅在信号的奇异点附近产生较大的系数。但是，许多信号如自然图像、遥感图像、生物医学图像等，不仅包含显著的低频分量，同时也包含了大量的高频信号，双树离散小波难以高效地表达这类信号。为了克服双树离散小波这一缺点，考虑将小波包方法和双树离散小波方法结合起来，对高频子带和低频子带都进行分解。与双树离散小波相比，双树离散小波包能够提供更精细的频域表示，继承了双树离散小波的方向选择性，并增加了小波方向的数目，同时能够根据输入信号的特性自适应地优选分解结构。步骤S102，根据双树离散小波对图像进行分解，获得多个子频带以及所述子频带的小波系数。具体地，首先采用q-shift滤波器将图像分解为方向子带后，再采用普通滤波器进行各向异性分解，从而保证所得的小波具有方向选择性。由于在各向异性分解中，DDWT(Distributed Discrete Wavelet Transform 分散式离散小波变换)子带作为复子带进行分解，因此虚部和实部具有相同的分解结构。由于各DDWT子带独立地进行自适应各向异性分解，因此ADDWP (adaptive dual-tree discrete wavelet packet自适应双树离散小波包)在基函数优选中需要搜索的分解结构数目为各DDWT子带分解结构数目之和为0029]
权利要求
1.一种基于双树离散小波包的图像去噪方法，其特征在于，包括以下步骤51构造双树离散小波；52根据所述双树离散小波对图像进行分解，获得多个子频带以及所述子频带的小波系数；53计算各个子频带的小波系数的自相关系数，并选取所述小波系数的自相关系数大于预定的阈值的子频带进行小波包分解，获得第二层双树离散小波包，以此类推获得多层双树离散小波包；S4:基于所述多层双树离散小波包将所述图像变换到小波域中，并计算所述小波系数的邻域相关性，以及根据所述小波系数的邻域相关性对所述图像进行去噪和图像增强。
2.根据权利要求1所述的方法，其特征在于，所述步骤Sl进一步包括 使用q-shift方案构造所述双树离散小波。
3.根据权利要求2所述的方法，其特征在于，所述步骤S2进一步包括 使用q-shift滤波器将所述图像分解为方向子带；以及使用滤波器对所述方向子带进行各向异性分解，获得所述多个子频带以及所述子频带的小波系数。
4.根据权利要求1所述的方法，其特征在于，所述步骤S3进一步包括m η m ηΣΣΣΣ ( - )( -%)^^ /=1 /=1 /=1 /2=1 q -----1 m η m ηm ηΣΣΣΣΣΣΚ/=1 j=\ /=1 h=\ /=1 /=1根据所述子频带的小波系数，通过以下的公式，计算所述小波系数的自相关系数，其中，Wu^lh为所述小波系数，m是 和Wlh的行数，η是、和Wlh的列数，Wy是 为中心的一个窗口的小波系数的均值，q是自相关系数，μ是空间邻接矩阵且if Wij 和 Wlh邻接 Mya =i O if %和％不邻接 O if Wlj 就是wlh
5.根据权利要求1所述的方法，其特征在于，所述步骤S4进一步包括541根据所述小波系数，通过以下的公式，计算所述小波系数的邻域相关性，W2 n2 m2 n2f _ γ _、Σ Σ Σ Σ - ^V^wIh -TnlTn2UlH2 m2 η2 m2 η2m2 η2ΣΣΣΣ ^ijih ΣΣ(^-%·)4/=W1 j=nx I=Tnl Ii=Yil/=W1 ]=ηλ其中，Hipnpn^n2为以Wij为中心的窗口的顶点坐标值；542根据所述小波系数的邻域相关性，通过以下的公式对所述小波系数进行修正，ιw = w amlnlm2n2 mlnlm2n2 mlnlm2n全文摘要
本发明提出一种基于双树离散小波包的图像去噪方法，包括以下步骤构造双树离散小波；根据双树离散小波对图像进行分解，获得多个子频带及其小波系数；计算各个子频带的小波系数的自相关系数，并选取自相关系数小于预定的阈值的子频带进行小波包分解，获得第二层双树离散小波包，依此类推获得多层双树离散小波包；基于多层双树离散小波包将图像变换到小波域中，并计算小波系数的邻域相关性，以及根据邻域相关性对图像进行去噪和图像增强。本发明通过利用小波系数之间的相似性关系分析图像的噪声分布特征特性，设计出的小波包构造方案在实现自然场景去噪的情况下，能够同时实现图像边缘和细节的增强。
文档编号G06T5/00GK102236888SQ20111020743
公开日2011年11月9日 申请日期2011年7月22日 优先权日2011年7月22日
发明者刘芳, 戴琼海 申请人:清华大学