数字图像去噪典型算法及matlab实现希望得到大家的指点和帮助图像去噪是数字图像处理中的重要环节和步骤。

去噪效果的好坏直接影响到后续的图像处理工作如图像分割、边缘检测等。

图像信号在产生、传输过程中都可能会受到噪声的污染，一般数字图像系统中的常见噪声主要有：高斯噪声（主要由阻性元器件内部产生）、椒盐噪声（主要是图像切割引起的黑图像上的白点噪声或光电转换过程中产生的泊松噪声）等；目前比较经典的图像去噪算法主要有以下三种：均值滤波算法：也称线性滤波，主要思想为邻域平均法，即用几个像素灰度的平均值来代替每个像素的灰度。

有效抑制加性噪声，但容易引起图像模糊，可以对其进行改进，主要避开对景物边缘的平滑处理。

中值滤波：基于排序统计理论的一种能有效抑制噪声的非线性平滑滤波信号处理技术。

中值滤波的特点即是首先确定一个以某个像素为中心点的邻域，一般为方形邻域，也可以为圆形、十字形等等，然后将邻域中各像素的灰度值排序，取其中间值作为中心像素灰度的新值，这里领域被称为窗口，当窗口移动时，利用中值滤波可以对图像进行平滑处理。

其算法简单，时间复杂度低，但其对点、线和尖顶多的图像不宜采用中值滤波。

很容易自适应化。

Wiener维纳滤波：使原始图像和其恢复图像之间的均方误差最小的复原方法，是一种自适应滤波器，根据局部方差来调整滤波器效果。

对于去除高斯噪声效果明显。

实验一：均值滤波对高斯噪声的效果l=imread（'C:\Documents and 桌面\1.gif'）;% 读取图像J=imnoise(l,'gaussian',0,0.005);% 加入均值为0 ,方差为 0.005 的高斯噪声subplot(2,3,1);imshow(l);title(' 原始图像');subplot(2,3,2); imshow(J);ti tle('加入高斯噪声之后的图像’)；%采用MATLAB 中的函数filter2 对受噪声干扰的图像进行均值滤波K1=filter2(fspecial('average',3),J)/255; % 模板尺寸为3K2=filter2(fspecial('average',5),J)/255;% 模板尺寸为5K3=filter2(fspecial('average',7),J)/255; % 模板尺寸为7K4= filter2(fspecial('average',9),J)/255; % 模板尺寸为9subplot(2,3,3);imshow(K1);ti tle(' 改进后的图像1');subplot(2,3,4); imshow(K2);title(' 改进后的图像2');subplot(2,3,5);imshow(K3);title(' 改进后的图像3');subplot(2,3,6);imshow(K4);title(' 改进后的图像4');J=imnoise(l,'gaussian',0,0.005);% 加入均值为0 ,方差为 0.005 的高斯噪声PS:filter2用法：filter2用法fspecial函数用于创建预定义的滤波算子，其语法格式为: h = fspecial(type)h = fspecial(type,parameters)参数type制定算子类型，parameters指定相应的参数，具体格式为：K2=wiener2(J,[7 7]); %对加噪图像进行二维自适应维纳滤波type='average'，为均值滤波，参数为n,代表模版尺寸，用向量表示，默认值为［3,3］。

type= 'gaussian'，为高斯低通滤波器，参数有两个，n表示模版尺寸，默认值为［3,3］，sigma表示滤波器的标准差，单位为像素，默认值为0.5。

type= 'laplacian'，为拉普拉斯算子，参数为alpha，用于控制拉普拉斯算子的形状，取值范围为［0,1］，默认值为0.2。

type= 'log'，为拉普拉斯高斯算子，参数有两个，n表示模版尺寸，默认值为［3,3］，sigma为滤波器的标准差，单位为像素，默认值为0.5type= 'prewitt'，为prewitt算子，用于边缘增强，无参数。

type= 'sobel'，为著名的sobel算子，用于边缘提取，无参数。

type= 'unsharp'，为对比度增强滤波器，参数alpha用于控制滤波器的形状，范围为［0,1］，默认值为0.2。

据我目测，使用均值滤波去噪(高斯噪声)效果选用的邻域半径越大效果越好，当然其代价也会更大，另外确切的去噪效果的好坏还需要用SNR等数据来度量。

实验二：二维自适应维纳滤波对高斯噪声的滤除效果l=imread('C:\Documents and 桌面\1.gif'); % 读取图像J=imnoise(l,'gaussian',0,0.005); % 加入均值为0，方差为0.005 的高斯噪声K2=wiener2(J,［3 3］); %对加噪图像进行二维自适应维纳滤波K2=wiener2(J,［5 5］); %对加噪图像进行二维自适应维纳滤波K2=wiener2(J,[9 9]); %对加噪图像进行二维自适应维纳滤波subplot(2,3,1);imshow(l);title(' 原始图像');subplot(2,3,2);imshow(J);title(' 加噪图像');subplot(2,3,3);imshow(K1);ti tle(' 恢复图像1');subplot(2,3,4);imshow(K2);title(' 恢复图像2');subplot(2,3,5);imshow(K3);title(' 恢复图像3');subplot(2,3,6);imshow(K4);title(' 恢复图像3');PS:维纳滤波的两个函数wiener2与deconvwnr都能够完成维纳滤波的功能，deconvwnr强调图象复原方面，wiener2强调图象空间域锐化的作用，其中J=wiener2(l,[m,n])返回有噪声图像I经过wierner(维纳)滤波后的图像，[m,n]指定滤波器窗口大小为m*n,默认值为3*3，J=wiener2(l,[m,n],noise)指定噪声的功率，[J,noise]=wiener2(l,[m,n])在图像滤波的同时，返回噪声功率的估计值noise。

:维纳滤波的两个函数wiener2与deconvwnr都能够完成维纳滤波的功能，deconvwnr强调图象复原方面，wiener2强调图象空间域锐化的作用，其中J=wiener2(l,[m,n])返回有噪声图像I经过wierner(维纳)滤波后的图像，[m,n]指定滤波器窗口大小为m*n,默认值为3*3，J=wiener2(l,[m,n],noise)指定噪声的功率，[J,noise]=wiener2(l,[m,n])在图像滤波的同时，返回噪声功率的估计值noise。

imnoise的语法格式为K2=wiener2(J,[7 7]); %对加噪图像进行二维自适应维纳滤波J = imno ise(l,type)J = imn oise(l,type,parameters)其中J = imnoise(l,type)返回对原始图像I添加典型噪声的有噪图像J。

参数type和parameters用于确定噪声的类型和相应的参数。

F面的命令是对图像1.gif分别加入高斯噪声、椒盐噪声和乘性噪声，其结果实验三：对加入椒盐噪声的图像分别作均值、中值和维纳滤波l=imread(1.gif');J1=im noise(l,'gaussia n',0,0.02);J2=imnoise(l,'salt & pepper',0.02);J3=im no ise(l,'speckle',0.02);运行效果见图2l=imread('C:\Documents and 桌面\1.gif');J=imnoise(l,'salt & pepper',0.02);%h=ones(3,3)/9% 产生3*3 的全1 数组%B=co nv2(J,h);% 卷积运算K2=filter2(fspecial('average',3),J)/255; % 均值滤波模板尺寸为3K= medfilt2(J);% 采用二维中值滤波函数medfilt2 对受椒盐噪声干扰的图像滤波K1=wiener2(J,[3 3]); %对加噪图像进行二维自适应维纳滤波subplot(2,3,1);imshow(l);title(' 原始图像');subplot(2,3,2);imshow(J);title(' 加噪图像');subplot(2,3,3);imshow(K2);ti tle('均值滤波后的图像'); subplot(2,3,4);imshow(K); ti tle('中值滤波后的图像');subplot(2,3,5);imshow(K1);ti tle(' 维纳滤波后的图像');PS:MATLAB中提供了卷积运算的函数命令conv2,其语法格式为：:MATLAB中提供了卷积运算的函数命令co nv2,其语法格式为：C = con v2(A,B)C = conv2(A,B)返回矩阵A和B的二维卷积C。

若A为ma x na的矩阵，B为mb x nb 的矩阵，贝U C 的大小为(ma+mb+1) x (na+nb+1)。

MATLAB图像处理工具箱提供了基于卷积的图象滤波函数filter2 , filter2的语法格式为：Y = filter2(h,X)其中丫 = filter2(h,X)返回图像X经算子h滤波后的结果，默认返回图像丫与输入图像X大小相同。

例如：其实filter2和conv2是等价的。

MATLAB在计算filter2时先将卷积核旋转180 度，再调用conv2函数进行计算。

Fspecial函数用于创建预定义的滤波算子，其语法格式为：h = fspecial(type)h = fspecial(type,parameters)参数type制定算子类型，parameters指定相应的参数，具体格式为前文已有叙述。