2019年整理高斯白噪声滤波

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

2、以上三图的分析知 b越大图像越不清晰。1.3源程序：X=imread( 'Lena.j pg' );J1=imnoise(X,'gaussian',0,0.15A2);imshow(J1) ；J2=imnoise(X,'gaussian',0,0.20A2); imshow(J2)；J2=imnoise(X,'gaussian',0,0.252); imshow(J2)2对高斯白噪声进行滤波2.1理想低通滤波器：2.1.1滤波原理理想低通滤波器：其传递函数为:>c肿(问=0戶-at。理想低通滤波器的冲激响应为：sinc(tt

3、o )2.1.2滤波结果兀兀不同d0(5,15,30)值下对图2的滤波结果 VV J7J不同d0(5,15,30)值下对图3的滤波结果不同d0(5,15,30)值下对图4的滤波结果2.1.3 源代码X=imread( 'Lena(25).jpg');% 读取图像l=rgb2gray(X);% 将图像变为灰度图figure;% 创建图形图像对象imshow(I);%显示灰度图像title('原始图像);加标题%将灰度图像的二维不连续 Fourier变换的零频率成分引导频谱的中心figure;%创建图形图像对象显示对s的绝对值取对数后的图像imshow(log(abs(s)

4、,);%'傅里叶频谱图);加标题title(M,N=size(s);%n仁floor(M/2);%对M/2进行取整分别返回s的行数到M中，列数到N中对N/2进行取整n2=floor(N/2);%ILPF 滤波，d0=5,15,30 dO=XX;%初始化d0 for i=1:M点(i,j)到傅里叶变换中心的距离for j=1:Nd=sqrt(i-n1)A2+(j-n2)A2);%s(i,j)=h*s(i,j);%IL PFend滤波后的频域表示if d<=d0 % 点(i,j)在通带内的情况h=1;%通带变换函数else%点(i,j)在阻带内的情况h=0;%阻带变换函数endend

5、s=ifftshift(s);%对 s进行反 FFI 移动%对进行二维反离散的Fourier变换后，取复数的实部转化为无符号8位整数 s=uint8(real(ifft2(s);figure;imshow(s);% 显示ILPF滤波后的图像title('ILPF 滤波(d=XX)');2.2高斯低通滤波器：2.2.1滤波原理由于高斯函数的傅里叶变换仍是高斯函数，因此高斯函数能够成为在时域和 频域都具有平滑性能的低通滤波器。2.2.2滤波结果 *- F «i>z图8均值为1，方差依次为0.01,0.02,0.03,0.04,0.05,0.06的高斯低通滤波对图2滤

6、波的结果lO¥图9均值为1，方差依次为 0.01,0.02,0.03,0.04,0.05,0.06的高斯低通滤波对图3滤波的结果i.3鴛菽，图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是方差lmg=double(s);n1=floor(n+1)/2);% b=zeros(n,n);计算图像中心for i=1:nfor j=1:nb(i,j)=ex p(-(i-n1)A2+(j-n1)A2)/(4*k)/(4*

7、 pi *k);endend%生成高斯序列Img1=conv2(lmg,b, d=uint8(lmg1); 主函数：h=imread( c=rgb2gray(h);same' );%用生成的高斯序列卷积运算，进行高斯滤波'LenaXX.j pg' );figure imshow(c);title( '加噪图像'n=1;k;%n 是均值，A2=gaussfilt(k,n,c););k是方差，n选定1，k可以变化figureimshow(A2);title( '高斯滤波2.3算术均值滤波器：2.3.1滤波原理均值滤波也称为线性滤波，其采用的主要方法为

8、领域平均法。其基本原理是用均值代替原图像中的各个像素值，即对待处理的当前像素点(X, y)，选择一个模板，该模板由其近邻的若干像素组成， 求模板中所有像素的均值，再把该均值 赋予当前像素点(x,y)，作为处理后图像在该点上的灰度个g(x, y)，即：g(x,y)=1/沱 f(x, y)m为该模板中包含当前像素在内的像素总个数。232滤波结果11模板为分别为5和3的均值滤波器对图2进行滤波的结果12模板为分别为5和3的均值滤波器对图3进行滤波的结果13模板为分别为5和3的均值滤波器对图4进行滤波的结果2.3.3 源代码function d=avefilt(x,n)a(1:n,1:n)=1;%a

9、即n*n模板，元素全为1 p=size(x);%输入图像是 p*q 的，且 p>n,q>n x1=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) 开始的n行n列元素与模板相乘s=sum(sum(c);%求c矩阵(即模板)中个元素之和x2(i+(n-1)/2,j+(n-1)/2)=s/(n*n);%将模板各元素的均值赋给模板中心位置的元素endend %未被赋值的元素取原值 d

10、=uint8(x2);主函数：h=imread('LenaXX.j pg');c=rgb2gray(h);figureimshow(c);title( '加噪图像)；A2=avefilt(c,n);%n表示模块，在本实验中分别取5和3figureimshow(A2);title( '均值滤波')2.4中值滤波器：2.4.1滤波原理中值滤波是基于排序统计理论的一种能有效抑制噪声的非线性信号处理技 术，基本原理是把数字图像或数字序列中一点的值用该点的一个邻域中各点值的 中值代替，让周围的像素值接近真实值，从而消除孤立的噪声点。方法是在某种 结构的二维滑动模板

11、下，将板内像素按照像素值的大小进行排序，生成单调上升(或下降)的二维数据序列。二维中值滤波输出为：g(x,y) =mid f(xk, y l),(k"W)其中f (x,y), g(x,y)分别为原始图像和处理后图像。W为二维模板，通常为2*2，3*3区域，也可以是不同的形状，如线状，圆形，十字形，圆环形等。2进行滤波的结果2.4.2滤波结果If! I 图15模板为分别为5和3的中值滤波器对图3进行滤波的结果图14模板为分别为5和3的中值滤波器对图2图16模板为分别为5和3的中值滤波器对图4进行滤波的结果2.4.3 源代码function d=midfilt(x,n) p=size(x

12、); % 输入图像是 p*q 的，且 p>n,q>n x1=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); %矩阵的第一行e=c(1,:);%cfor k=2:n取出x1中从(i,j) 开始的n行n列元素e=e,c(k,:);%将c矩阵变为一个行矩阵endmm=median(e);%mm是中值x2(i+(n-1)/2,j+(n-1)/2)=mm;%end将模板各元素的中值赋给模板中心位置元素end %未被赋值的元素取原值 d=uint8(x2););主函数：h=imread( 'LenaXX.j pg'c=rgb2gray(h);figure imshow(c);A2=midfilt(c,n);%n figuretitle( '加噪图像)；表示模块，在本实验中分别取5和3 imshow(A2);title( '中值滤波')实验结果分析:通过以上四种滤波方法的实验及其结果， 可以看出理想低通滤波器的滤波效 果并不理想，特别是当dO的水平很低时，图像非常模糊，即使