高斯白噪声滤波

1、西安电子科技大学课程论文数字图像处理高斯白噪声滤波班级：070821作者：董文凯学号： 07082001时间：2011-06-30高斯白噪声滤波实验要求对实际Lena图像分别加入噪声标准差 仃=15, 20, 25的高斯白噪声，用理想 低通滤波器、高斯低通滤波器、算术均值滤波器和中值滤波器对实际 Lena图像 进行去噪，比较其去噪效果。实验内容1 .对Lena图像加高斯白噪声1.1 原始图例:Lena图像见图1.采用经典Lena图像作为实验样例进行本实验的操作，原始图1原始Lena图1.2 加噪结果:图2。=0.15的噪声图图3。=0.20的噪声图图4仃=0.25的噪声图结论：经过对以上三图的

2、分析知。越大图像越不清晰1.3 源程序：X=imread( 'Lena.jpg' );J1=imnoise(X,'gaussian',0,0.15A2);imshow(J1) ；J2=imnoise(X,'gaussian',0,0.20A2);imshow(J2)；J2=imnoise(X,'gaussian',0,0.25A2); imshow(J2)2对高斯白噪声进行滤波2.1 理想低通滤波器：2.1.1 滤波原理理想低通滤波器：其传递函数为:"产« <eoc cA8c cH(js) = 0卜(。尸

3、at。理想低通滤波器的冲激响应为：c sin c t-t。h t 二一二JISa% t -to2.1.2滤波结果- c t-to图5不同d0(5,15,3。)值下图2的滤波结果图6不同d0(5,15,30)值下图3的滤波结果图7不同d0(5,15,30)值下图4的滤波结果2.1.3 源代码X=imread( 'Lena(25).jpg' );% 读取图像I=rgb2gray(X);%将图像变为灰度图figure;%创建图形图像对象imshow(I);% 显示灰度图像 title('原始图像');%加标题 %>灰度图像的二维不连续 Fourier变换的零频率

4、成分引导频谱的中心s=fftshift(fft2(I); figure;%创建图形图像对象imshow(log(abs(s),);%显示对s的绝对值取对数后的图像title('傅里叶频谱图');%加标题 M,N=size(s);%分别返回s的行数到Ml,列数到N中n1=floor(M/2);%XM/2 进行取整n2=floor(N/2);%XN/2 进行取整%ILPF 滤波，d0=5,15,30 d0=XX;%初始化d0 for i=1:Mfor j=1:Nd=sqrt(i-n1)A2+(j-n2)A2);%点(i,j)到傅里叶变换中心的距离if d<=d0 % 点(i,

5、j)在通带内的情况h=1;%通带变换函数else % 点(i,j) 在阻带内的情况 h=0;%阻带变换函数end s(i,j)=h*s(i,j);%ILPF滤波后的频域表示end end s=ifftshift(s);%X：' s进行反 FFI 移动%对进行二维反离散的Fourier变换后，取复数的实部转化为无符号8位整数s=uint8(real(ifft2(s); figure; imshow(s);% 显示ILPF滤波后的图像 title( 'ILPF 滤波(d=XX)');2.2 高斯低通滤波器：2.2.1 滤波原理由于高斯函数的傅里叶变换仍是高斯函数，因此高斯函

6、数能够成为在时域和 频域都具有平滑性能的低通滤波器。2.2.2 滤波结果M .金图8均值为1,方差依次为 0.01,0.02,0.03,0.04,0.05,0.06的高斯低通滤波对图2滤波的结果图9均值为1,方差依次为 0.01,0.02,0.03,0.04,0.05,0.06的高斯低通滤波对图3滤波的结果U iM图10均值为1,方差依次为 0.01,0.02,0.03,0.04,0.05,0.06的高斯彳氐通滤波对图4滤波的结果2.2.3 源代码function d=gaussfilt(k,n,s)%s是需要滤波的图像，n是均值，k是方差 Img=double(s);n1=floor(n+1

7、)/2);% 计算图像中心 b=zeros(n,n); for i=1:n for j=1:n b(i,j)=exp(-(i-n1)A2+(j-n1)A2)/(4*k)/(4*pi*k); end end %生成高斯序列 Img1=conv2(Img,b, 'same' );%用生成的高斯序列卷积运算，进行高斯滤波 d=uint8(Img1); 主函数：h=imread( 'LenaXX.jpg' ); c=rgb2gray(h); figure imshow(c); title( '加噪图像'); n=1;k;%n是均值，k是方差，n选定1,

8、k可以变化 A2=gaussfilt(k,n,c); figure imshow(A2); title( '高斯滤波)2.3 算术均值滤波器：2.3.1 滤波原理均值滤波也称为线性滤波，其采用的主要方法为领域平均法。其基本原理是用均值代替原图像中的各个像素值，即对待处理的当前像素点(x,y),选择一个模板，该模板由其近邻的若干像素组成， 求模板中所有像素的均值，再把该均值 赋予当前朦昊点(x,y),作为疝理后囱彳象在该点工的反低个 g(x,y),即：g(x,y) =1/m" f(x,y)m为该模板中包含当前像素在内的像素总个数。2.3.2 滤波结果图11模板为分别为5和3的均

9、值滤波器对图 2进行滤波的结果图12模板为分别为5和3的均值滤波器对图3进行滤波的结果图13模板为分别为5和3的均值滤波器对图4进行滤波的结果2.3.3 源代码function d=avefilt(x,n)a(1:n,1:n)=1;%a 即n*n模板，元素全为1p=size(x);%输入图像是 p*q 的，且 p>n,q>nx1=double(x);x2=x1;%A(a:b,c:d)表示矩阵A的第a到b行，第c到d列的所有元素for i=1:(p(1)-n+1)for j=1:(p(2)-n+1)c=x1(i:i+(n-1),j:j+(n-1).*a;%取出x1中从(i,j)开始的

10、n行n列元素与模板相乘s=sum(sum(c);%求c矩阵(即模板)中个元素之和x2(i+(n-1)/2,j+(n-1)/2)=s/(n*n);%>模板各元素的均值赋给模板中心位置的元素endend麻被赋值的元素取原值d=uint8(x2);主函数：h=imread( 'LenaXX.jpg' );c=rgb2gray(h);figureimshow(c);title( '加噪图像');A2=avefilt(c,n);%n表示模块，在本实验中分别取5和3figureimshow(A2);title( '均值滤波')2.4 中值滤波器:2.4

11、.1 滤波原理中值滤波是基于排序统计理论的一种能有效抑制噪声的非线性信号处理技 术，基本原理是把数字图像或数字序列中一点的信用该点的一个邻域中各点值的 中值代替，让周围的像素值接近真实值，从而消除孤立的噪声点。方法是在某种 结构的二维滑动模板下，将板内像素按照像素值的大小进行排序，生成单调上升(或下降)的二维数据序列。二维中值滤波输出为：g(x, y)=mid f(x -k, y -l),(k,l W)其中f(x,y),g(x,y)分别为原始图像和处理后图像。W为二维模板，通常为2*2, 3*3区域，也可以是不同的形状，如线状，圆形，十字形，圆环形等。2.4.2 滤波结果图14模板为分别为5和

12、3的中值滤波器对图2进行滤波的结果图15模板为分别为5和3的中值滤波器对图3进行滤波的结果图16模板为分别为5和3的中值滤波器对图4进行滤波的结果2.4.3 源代码 function d=midfilt(x,n) p=size(x); %输入图像是 p*q 的，且 p>n,q>nx1=double(x); x2=x1; for i=1:(p(1)-n+1) for j=1:(p(2)-n+1) c=x1(i:i+(n-1),j:j+(n-1); %取出 x1 中从(i,j) 开始的 n行n 列元素e=c(1,:);%c矩阵的第一行for k=2:n e=e,c(k,:);%将c矩阵

13、变为一个行矩阵end mm=median(e);%mm 是中值 x2(i+(n-1)/2,j+(n-1)/2)=mm;%将模板各元素的中值赋给模板中心位置元素end end %未被赋值的元素取原值 d=uint8(x2); 主函数：h=imread( 'LenaXX.jpg' ); c=rgb2gray(h); figure imshow(c); title( '加噪图像'); A2=midfilt(c,n);%n表示模块，在本实验中分别取 5和3figure imshow(A2); title( '中值滤波') 实验结果分析：通过以上四种滤波方法的实验及其结果，可以看出理想低通滤波器的滤波效果并不理想，特别是当d0