图像去噪算法综述

图像降噪算法总结

分析各种算法的优点和缺点

1、BM3D 降噪
2、DCT 降噪
3、PCA 降噪
4、K-SVD 降噪
5、非局部均值降噪
6、WNNM 降噪
7、基于主成分分析和双边滤波的图像降噪算法
8、小波变换
9、小波阈值降噪
10、Contourlet 变换
11、基于平移不变 Contourlet 变换的 SAR 图像降噪**

1、BM3D 降噪
BM3D 是一种降噪方法提高了图像在变换域的稀疏表示。BM3D 降噪方法的优点是更好的保留图像中的一些细节，BM3D采用了不同的去噪策略。通过搜索相似块并在变换域进行滤波，得到块评估值，最后对图像中每个点进行加权得到最终去噪效果。 
原理：首先将一幅图像分割成尺寸较小的小像素片，选定参考片后，寻找与参考片相似的小片组成 3D 块。此过程过后将得到 3D 块。然后将所有相似块进行 3D 变换。将变换后的 3D 块进行阈值收缩，这也是除去噪声的过程。然后进行 3D 逆变换。最后将所有的 3D 块通过加权平均后还原到图像中。
BM3D算法的大致流程：
第1步. 初始估计
（1）逐块估计。对含噪图像中的每一块 
 （i）分组。找到它的相似块然后把它们聚集到一个三维数组。 
（ii）联合硬阈值。对形成的三维数组进行三维变换，通过对变换域的系数进行硬阈值处理减弱噪声，然后逆变换得到组中所有图像块的估计值，然后把这些估计值返回到他们的原始位置。
（2）聚集。对得到的有重叠的块估计，通过对他们进行加权平均得到真实图像的基础估计。 
第2步. 最终估计 
（1）逐块估计。对基础估计图像中的每一块 
  （i） 分组。通过块匹配找到与它相似的相似块在基础估计图像中的位置，通过这些位置得到两个三维数组，一个是从含噪图像中得到的，一个是从基础估计图像中得到的。 
   （ii）联合维纳滤波。对形成的两个三维数组均进行三维变换，以基础估计图像中的能量谱作为能量谱对含噪三维数组进行维纳滤波，然后逆变换得到组中所有图像块的估计，然后把这些估计值返回到他们的原始位置。
（2）聚集。对得到的有重叠的局部块估计，通过对他们进行加权平均得到真实图像的最终估计。

2、DCT 降噪
DCT 变换是正交变换。这个方法是对滑动的窗口内的图像做 DCT 变换，在局部变换域做系数收缩。由滑动窗口产的重叠部分满足图像的过完备表示。如果采用的变换不能将普通图像的细节过完备表示，那么这种变换的过完备性就不能抵消无效的收缩。因此只使用二维变换获得稀疏表示只适于特定的图像模式。
B-DCT图像降噪方法一般是处理逐像素滑动的窗口得到的像素矩阵，对矩阵做二维DCT变换，（这以过程可以看作是在矩阵左右两侧同时乘正交矩阵，左乘的矩阵相当于对像素矩阵做列DCT变换，右乘的正交矩阵相当于对像素矩阵做行DCT变换）在变换域中对系数进行硬阈值收缩，然后对收缩后的矩阵做二维DCT逆变换。与B-DCT变换不同的是 SA-DCT变换可以对任意形状的像素片做DCT变换，它的方法是分别对每行和每列做一维DCT变换做完后进行系数收缩，然后再相应一维 DCT 变换的顺序做逆变换。
形状自适应变换域系数收缩降噪方法面临两点问题：1、变换能否适应像素片形状，2、像素片形状是否适应图像的纹理特点。PSA-DCT 方法针对第二点采用 LPA-ICI 方法获得与图像纹理特点对应的像素小片。实验虽然获得了较好的降噪效果，但是 DCT 变换降噪这种方法本质上的问题：1、由于是逐像素滑动变换动窗口带来的数据存储和 DCT 变换计算量较大。2、还有DCT变换基底并没有将自然图像信息与噪声做出良好的分割导致了这类方法的局限性。

3、PCA 降噪
PCA应用于图像降噪的方法的提出克服了标准正交变换带来的缺点，这种方法对于高结构性的图像细节都有很好的结果。但是由于噪声的存在势必影响主成分分析的准确性。

4、K-SVD 降噪
K-SVD降噪方法解决了固定变换矩阵基底不能自适应图像纹理信息的缺点，，这种算法的缺点是更新字典的计算量较大。
5、非局部均值降噪
早期的降噪方法一般为局部平滑滤波方法，例如，高斯滤波降噪，局部均值滤波降噪。这类方法很难处理图像的非平滑部分，像添加的噪声和图像的纹理信息，去噪的同时不能有效保留纹理信息。非局部均值降噪方法就为解决这一问题。

6.WNNM 降噪
设y=x+n，y是含噪声的图像，x是清晰图像，n是均值为零标准差为σ_n的高斯白噪声。y中局部小片yj，在y寻找yj的相似片，形成矩阵Yj，我们就 Yj=Xj+Nj,Xj和Nj是原图像和噪声形成的矩阵，Xj是一个低阶矩阵，通过这个性质可以解析下式达到降噪的目的：

这是一个WNNM问题，可以通过收缩Yj的奇异值方法解得。首先对Yj做SVD分解，得到特征矩阵U和V还有奇异值矩阵Σ. 对每个奇异值Σij减去权值??，
即Sw (Σ)ij=max(Σ_ii-wi,0)，权值w与Xj的奇异值有关，记Xj的奇异值为：
大的奇异值为Xj的主成分，收缩的少。
小的奇异值收缩的

7、基于主成分分析和双边滤波的图像降噪算法
方法主要分为两步进行处理。第一步首先采用局部像素块匹配算法选出 PCA 的训练样本集，然后应用主成分分析法滤除绝大部分的高斯噪声，第二步对系数收缩重构图像的残余噪声方差进行估计，然后将一步去噪后的图像输入自适应双边滤波器中，最后输出降噪后的图像。
8、小波变换
小波变换作为一种新的时频分析方法，具有多尺度、多分辨率分析的特点，为信号处理提供了一种新的强有力手段。小波变换在图像降噪领域的成功应用主要得益于其具有低熵性、多分辨率特性、去相关性和选基灵活性的优点。
小波降噪本质上是一个信号的滤波问题，实际上是特征提取和低通滤波的综合。波降噪的处理流程是：首先对含有噪声的信号进行多尺度小波变换，在各尺度下尽可能提取出小波系数，最后利用逆小波变换重构信号。

9、小波阈值降噪
信号通常是低频或者平稳信号的形式，而噪声一般都是高频信号，所以降噪的过程也就是滤除高频信号的一个过程，小波阈值降噪就是利用的这一原理。
小波阈值降噪方法就是通过设置某种阈值，将小波系数与阈值进行比较，将小于阈值的系数设置为 0，而对大于阈值的小波系数，通过某种阈值函数(也称阈值规则)进行修正得到其估计系数。阈值的确定是最为关键的问题。阈值的选择对降噪效果有着很大的影响，如果阈值太小，降噪后的信号仍然有噪声存在，反之，阈值太大则重要的图像特征又将被滤除，造成图像模糊，阈值的选择始终是在抑制噪声和保护图像细节之间折中。目前应用较多的阈值方法包括针对多维正态变量联合分布的 Visushrink 阈值方法、利用均方差准则的无偏估计的 Sureshrink 阈值方法、利用贝叶斯最大后验概率估计理论的 Mapsbrink 阈值方法、利用邻域信息的分块阈值方法等。
10、Contourlet 变换
Contourlet 变换是一种多分辨率、多方向的图像稀疏表示方法。它能够用少量的系数有效地表示图像中的轮廓等重要特征。它主要利用了拉普拉斯塔形分解和方向滤波器组来实现图像的稀疏表示。
①改进阈值法：改进阈值法主要有两种改进方法，第一种就是对阈值的选取进行改进，改进后仍然是一个硬阈值，这种方法与全局阈值法有所不同，它是根据系数及其领域系数的关系来计算出阈值，这种方法在一定程度上改善了降噪效果，但是该方法的阈值选取具有一定的随机性，降噪效果不稳定。第二种就是自适应的阈值方法，如引入高频能量比来不断更新阈值，但是这种改进的阈值方法也有一定的缺陷，因为它只是在全局阈值的基础上加了一个调节因子，得到的阈值时不够准确的，降噪效果不是很理想，也会损失图像的细节信息。
②Wiener 滤波：在 Contourlet 分解的基础上加入了 Wiener 滤波这一步骤，它的原理是对 Contourlet 分解的高频系数进行椭圆窗口的 Wiener 滤波；利用自适应窗口的 Wiener 滤波对 Contourlet 分解得到的高频子带进行滤波。
③建模法：据高频系数及其邻域系数的相关性，对这些系数进行 HMT 建模，然后应用到图像中对图像进行降噪处理；对信号方差进行最优估计，然后对图像进行阈值去噪。此类方法的优点是利用了系数间的相关性，缺点是并没有充分利用系数间的相关性。
④改进 Contourlet 降噪算法由于 Contourlet 变换的下采样操作导致 Contourlet没有平移不变性，LP 滤波器组缺乏局部化能力，这使 DFB 频谱中出现了频谱混叠。抗混叠 Contourlet 变换对 LP 的滤波器组进行了一些改进，采用可分离低通滤波器以及满足 Nyquist 抽样定律的高通滤波器来代替 Contourlet 中的 LP 滤波器组，图像通过高通滤波器后级联 DFB，这样使各方向的子带中无法出现由多尺度分解造成的频谱混叠，抑制了频谱混叠问题。
11、基于平移不变 Contourlet 变换的 SAR 图像降噪
Contourlet 变换会在重构图像中会出现严重的划痕原因是Contourlet 变换在 LP 和 DFB 分解中采用了下采样操作，导致了局部性不够理想，不能够平移不变 。 鉴于这个原因，本文采用了 Eslami 提出的平移不变 Contourlet变换（TICT)，该变换克服了 Contourlet 不能够平移不变的缺点，不仅能很好的抑制噪声，还能达到更强的鲁棒性。