【图像去噪】小波变换和多种阙值改进的小波变换彩色图像去噪（含PSNR）【含Matlab源码 3700期】

⛄一、图像去噪及滤波简介

1 图像去噪
1.1 图像噪声定义
噪声是干扰图像视觉效果的重要因素，图像去噪是指减少图像中噪声的过程。噪声分类有三种：加性噪声，乘性噪声和量化噪声。我们用f(x,y）表示图像，g(x,y）表示图像信号，n(x,y）表示噪声。
图像去噪是指减少数字图像中噪声的过程。现实中的数字图像在数字化和传输过程中常受到成像设备与外部环境噪声干扰等影响，称为含噪图像或噪声图像。去噪是图像处理研究中的一个重点内容。在图像的获取、传输、发送、接收、复制、输出等过程中，往往都会产生噪声，其中的椒盐噪声是比较常见的一种噪声，它属于加性噪声。

1.2 图像噪声来源
（1）图像获取过程中
图像传感器CCD和CMOS采集图像过程中受传感器材料属性、工作环境、电子元器件和电路结构等影响，会引入各种噪声。
（2）图像信号传输过程中
传输介质和记录设备等的不完善，数字图像在其传输记录过程中往往会受到多种噪声的污染。

1.3 噪声分类
噪声按照不同的分类标准可以有不同的分类形式：
基于产生原因：内部噪声，外部噪声。
基于噪声与信号的关系：
加性噪声：加性噪声和图像信号强度是不相关的，这类带有噪声的图像g可看成为理想无噪声图像f与噪声n之和：
g = f + n；
乘性嗓声：乘性噪声和图像信号是相关的，往往随图像信号的变化而变化，载送每一个象素信息的载体的变化而产生的噪声受信息本身调制。在某些情况下，如信号变化很小，噪声也不大。为了分析处理方便，常常将乘性噪声近似认为是加性噪声，而且总是假定信号和噪声是互相统计独立。
g = f + f*n
按照基于统计后的概率密度函数：
是比较重要的，主要因为引入数学模型这就有助于运用数学手段去除噪声。在不同场景下噪声的施加方式都不同，由于在外界的某种条件下，噪声下图像-原图像（没有噪声时）的概率密度函数（统计结果）服从某种分布函数，那么就把它归类为相应的噪声。下面将具体说明基于统计后的概率密度函数的噪声分类及其消除方式。

2 小波滤波
随着小波理论的日益完善，其以自身良好的时频特性在图像去噪领域受到越来越多的关注，开辟了用非线性方法去噪的先河。具体来说，小波能够去噪主要得益于小波变换有如下特点：
（1）低熵性。小波系数的稀疏分布，使图像变换后的熵降低。 意思是对信号（即图像）进行分解后，有更多小波基系数趋于0（噪声），而信号主要部分多集中于某些小波基，采用阈值去噪可以更好的保留原始信号。
（2）多分辨率特性。由于采用了多分辨方法，所以可以非常好地刻画信号的非平稳性，如突变和断点等（例如0-1突变是傅里叶变化无法合理表示的），可以在不同分辨率下根据信号和噪声的分布来消除噪声。
（3）去相关性。小波变换可对信号去相关，且噪声在变换后有白化趋势，所以小波域比时域更利于去噪。
（4）基函数选择灵活。小波变换可灵活选择基函数，也可根据信号特点和去噪要求选择多带小波和小波包等（小波包对高频信号再次分解，可提高时频分辨率），对不同场合，选择不同小波基函数。

根据基于小波系数处理方式的不同，常见去噪方法可分为三类：
（1）基于小波变换模极大值去噪（信号与噪声模极大值在小波变换下会呈现不同变化趋势）
（2）基于相邻尺度小波系数相关性去噪（噪声在小波变换的各尺度间无明显相关性，信号则相反）
（3）基于小波变换阈值去噪

小波去噪实现步骤：
（1）二维信号的小波分解。选择一个小波和小波分解的层次N，然后计算信号s到第N层的分解。
（2）对高频系数进行阈值量化。对于从1~N的每一层，选择一个阈值，并对这一层的高频系数进行软阈值量化处理。
（3）二维小波重构。根据小波分解的第N层的低频系数和经过修改的从第一层到第N的各层高频系数，计算二维信号的小波重构。

⛄二、部分源代码

% 对彩色图像进行去噪
I = imread(‘lena_color.png’, ‘png’); % 读入图像
X = im2double(I); % 转换成双精度类型

x_noise = noise(X, ‘gaussian’, 0.01); % 加入高斯噪声

% 提取三个通道信息
xr = x_noise(:, :, 1); % R通道
xg = x_noise(:, :, 2); % G通道
xb = x_noise(:, :, 3); % B通道

% 估计三个通道的阈值
[Cr, Sr] = wavedec2(xr, 2, ‘sym4’);
[Cg, Sg] = wavedec2(xg, 2, ‘sym4’);
[Cb, Sb] = wavedec2(xb, 2, ‘sym4’);
thr_r = Donoho(xr); % R通道全局阈值
thr_g = Donoho(xg); % G通道全局阈值
thr_b = Donoho(xb); % B通道全局阈值
thr_lvd_r = Birge_Massart(Cr, Sr); % R通道局部阈值
thr_lvd_g = Birge_Massart(Cg, Sg); % G通道局部阈值
thr_lvd_b = Birge_Massart(Cb, Sb); % B通道局部阈值

% 对三个通道分别进行去噪===================================================
% Donoho全局阈值 软阈值公式----------------------------------------------
x_soft_r = wdenoise(xr, ‘gbl’, ‘s’, thr_r, ‘sym4’, 2);
x_soft_g = wdenoise(xg, ‘gbl’, ‘s’, thr_g, ‘sym4’, 2);
x_soft_b = wdenoise(xb, ‘gbl’, ‘s’, thr_b, ‘sym4’, 2);
% -----------------------------------------------------------------------

% Donoho全局阈值 硬阈值公式----------------------------------------------
x_hard_r = wdenoise(xr, ‘gbl’, ‘h’, thr_r, ‘sym4’, 2);
x_hard_g = wdenoise(xg, ‘gbl’, ‘h’, thr_g, ‘sym4’, 2);
x_hard_b = wdenoise(xb, ‘gbl’, ‘h’, thr_b, ‘sym4’, 2);
% -----------------------------------------------------------------------

% Bige-Massa策略 软阈值公式----------------------------------------------
x_soft_lvd_r = wdenoise(xr, ‘lvd’, ‘s’, thr_lvd_r, ‘sym4’, 2);
x_soft_lvd_g = wdenoise(xg, ‘lvd’, ‘s’, thr_lvd_g, ‘sym4’, 2);
x_soft_lvd_b = wdenoise(xb, ‘lvd’, ‘s’, thr_lvd_b, ‘sym4’, 2);
% -----------------------------------------------------------------------

% Bige-Massa策略 硬阈值公式----------------------------------------------
x_hard_lvd_r = wdenoise(xr, ‘lvd’, ‘h’, thr_lvd_r, ‘sym4’, 2);
x_hard_lvd_g = wdenoise(xg, ‘lvd’, ‘h’, thr_lvd_g, ‘sym4’, 2);
x_hard_lvd_b = wdenoise(xb, ‘lvd’, ‘h’, thr_lvd_b, ‘sym4’, 2);
% -----------------------------------------------------------------------

% 半软阈值---------------------------------------------------------------
x1_r = den1(xr, ‘sym4’, 2, thr_r);
x1_g = den1(xg, ‘sym4’, 2, thr_g);
x1_b = den1(xb, ‘sym4’, 2, thr_b);
% -----------------------------------------------------------------------

% 半软阈值 + 均值滤波----------------------------------------------------
x1_5_r = den1_5_1(xr, ‘sym4’, 2, thr_r, 0.5thr_r);
x1_5_g = den1_5_1(xg, ‘sym4’, 2, thr_g, 0.5thr_g);
x1_5_b = den1_5_1(xb, ‘sym4’, 2, thr_b, 0.5*thr_b);
% -----------------------------------------------------------------------

% 自适应阈值-------------------------------------------------------------
x4_r = den4(xr, ‘sym4’, 2);
x4_g = den4(xg, ‘sym4’, 2);
x4_b = den4(xb, ‘sym4’, 2);

⛄三、运行结果

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 蔡利梅.MATLAB图像处理——理论、算法与实例分析[M].清华大学出版社，2020.
[2]杨丹,赵海滨,龙哲.MATLAB图像处理实例详解[M].清华大学出版社，2013.
[3]周品.MATLAB图像处理与图形用户界面设计[M].清华大学出版社，2013.
[4]刘成龙.精通MATLAB图像处理[M].清华大学出版社，2015.
[5]基于matlab的传统算法图像去噪的实现原理

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除