1、数字图像处理第六次作业数字图像处理第六次作业摘要本次报告主要记录第六次作业中的各项任务完成情况。本次作业以Mat lab 2013为平 台，结合matlab函数编程实现对所给图像文件的相关处理：1.在测试图像上产生高斯噪声 lena图-需能指定均值和方差；并用滤波器(自选)恢复图像：2.实现下边要求：(a)实现 模糊滤波器如方程Eq. (5.6-11) ； (b)模糊lena图像：45度方向，1=1： (c)在模糊的 lena图像中增加高斯噪声，均值;0 ,方差=10 pixels以产生退化图像；(d)分别利用方 程Eq. (5.8-6)和(5.9-4),恢复图像。以上任务完成后均得到了预期的
2、结果。1.在测试图像上产生高斯噪声lena图-需能指定均值和方差:并用滤波器(自选)恢复图像： (1)实验原理与方法图像复原处理是建立在图像退化的数学模型基础上的，这个退化数学模型能够反映图像 退化的原因。图像的退化过程可以理解为施加于原图像上的运算和噪声两者联合作用的结 果，图像退化模型如图1所示，可以表示为：g(x,y) = Hf(x, y )+n( x,y) = f( x, y ) * h( x,y)+n( x,y )图1图像退化模型高斯噪声是指它的概率密度函数服从高斯分布(即正态分布)的一类噪声。一个高斯随 机变量z的PDF可表示为：P(Z)= /，2猿 2(?其中Z代表灰度，U是Z的
3、均值，是Z的标准差。高斯噪声的灰度值多集中在均值附近。本文采用5X5模板的中值滤波器和高斯滤波器(。=1.5)作业四中已经介绍过，中值滤波器是使用一个像素邻域中灰度级的中值来替代该像素值, 即，f(x,y)=median!。高斯滤波是一种根据高斯函数的形状来选择模板权 (SJ用5孙值的线性平滑滤波方法，具体操作是：用一个模板(或称卷积)扫描图像中的每一个像素， 用模板确定的邻域内像素的加权平均灰度值去替代模板中心像素点的值。利用mat lab中imnoise函数加入高斯噪声：g=imnoise(f, type,parameters)调用格式：g = imnoise(I, type)g = im
4、noise(I, type, parameters)参数Type对应的噪声类型如下：gaussian高斯白噪声localvar 0均值白噪声poisson，泊松噪声 salt & pepper1 盐椒噪声speckle乘性噪声滤波程序同作业四（2）处理结果原始图像lena.bmp加 入 gaussian噪声后的 lena.bmp(u=0.5 f s2=0.01)原始图像lena.bmp力口 入 gaussian噪 声后的 lena.bmp(u=O , s2 =0.01)原始图像lena.bmp加入 gaussian噪 声后的 lena.bmp(ii=0.5 , s2=0.1)lena力口 入
5、gauss ianO 声后(u=0.5 , s2 =0.01). bmp原始图像lena.bmp中值滤波(5x5)高斯滤波5x5(3)结果分析通过imnoise函数产生了被均值和方差可选的高斯噪声污染的图像。当高斯噪声均值 不变为0时，随着方差增加，图像噪声越严重：当高斯噪声方差不变时，均值会影响到整个 图像的灰度值，使整个图像变亮。与理论上均值和方差对图像的影响一致。分别使用高斯滤 波器和中值滤波器对加噪图像进行恢复。两种方法在一定程度上都可以降低噪声。高斯滤波 器降低噪声的同时保存的图像细行更丰富，亮度比原噪声图像和中值滤波后图像暗更接近原 始图像，中值滤波后图像亮度基本与原噪声图像相同。
6、2.实现下边要求：(a)实现模糊滤波器如方程Eq. (5. 6-11).(b)模糊模na图像:45度方向,T=l；(c)在模糊的lena图像中增加高斯噪声，均值二0 ,方差二10 pixels以产生退化图 像：(d)分别利用方程Eq. (5.8-6)和(5. 9-4),恢复图像：H(u9v) =-sin网ua + vb)产力)(5. 6-11)7i( ua + vb )利用上式模糊lena图像：45度方向，T=L即使a=b=0. 1, T=l,对原始用像的模像矩 阵做傅里叶变换并移至图像中心得到频域矩阵F,使H与F相乘后反傅里叶变换到空域得到 变换后图像。维纳滤波综合了退化函数和噪声统计特性两
7、个方面进行复原处理，其目标是寻找一个滤波器，使得复原后图像f( x, y)与原始图像，x, y )的均方误差最小：E_f(xfy)-f(x,y =min因此维纳滤波器又称为最小均方误差滤波器，在频率中用下式表达:F(uy) = 1 H( uy)H( u9v )|2 +Sn( ufv )/Sf( w,v)其中G(u,v)是退化图像的傅里叶变换，H(u,v)是退化函数 SQjXpV(2为噪声功率谱，Sf(uy) = F(uy为未退化图像的功率谱.式5. 8-6为人 1 H(u9v)2F(u9v) = -G(uy) H(u,v)H(uy) +K属于维纳滤波式5. 9-4为八 H ( U v)F(U9
8、V)= - 一一7G(uy)H(u9v)+yP(u9v)其中，丫是一个参数，必须对它进行调整以满足k-而=忖的条件，属于约束最 小二层方滤波利用以上两式恢复图像的流程与对模糊原始图像的流程相似:对原始图像的图像矩阵做 傅里叶变换并移至图像中心得到频域矩阵F,通过H得到户，使户与F相乘后反傅里叶变 换到空域得到变换后图像。在实现5. 9-4时借助matlab工具包以得到更好的效果。用 fspecial和imnoise函数得到45度方向，T=1的模糊lena图像，并在此图像上实现维纳滤 波和约束最小二乘方滤波。1)imfilter功能：对任意类型数组或多维图像进行滤波。用法:B = imfilte
9、r(A, H)B = imfilter (A, H, optionl, option2,.)或写做 g = imfilter (f, w, filtering_mode, boundary_options, size_options) 其中，f为输入图像，W为滤波掩模，g为滤波后图像。filtering_mode用于指定在波波过 程中是使用“相关”还是“卷积”。boundary.options用于处理边界充零问题，边界的大小 由滤波器的大小确定。2)fspecial功能：fspecial函数用于建立预定义的滤波算子。用法：h = fspecial (type) h = fspecial (ty
10、pe, para) 其中type指定算子的类型，para指定相应的参数。(2)处理结果lena.bm限始图像运动模楣化lena.bmp运动模糊化lena.bmp 模糊lena.bmp加入高斯噪声(u=0 , s2=0.01)lena运动模糊+高斯噪声维纳滤波的结果(K=0.06)lena运动模柳+高斯噪声约束最小二乘滤波的结果(3)结果分析1.按照书上公式编写的模糊函数图像是斜向下45度运动模糊，mat lab函数是斜向上45度 运动模糊的，公式的程序得到图像棱角比较分明边界比较明显。2.使用自己编写的函数进行维纳滤波，难点在于寻找令信噪比最大的K值，报告中显示了 K=0. 06时的滤波结果，
11、从结果看，视觉上的效果并不是很理想，噪声依然很大，要想达到更 好的效果可能需要寻找更加合适的K值或者直接使用matlab的deconvreg函数实现。3.最后采用MATL.UB提供的deconvreg函数进行约束最小二乘方滤波。从滤波后的结果看, 约束最小二乘方滤波得到了比维纳滤波更好的结果,噪声基本消除，图像变得模糊但是平滑。附录：参考文献：1 Rafael C. Gonzalez., et al.数字图像处理(第三版),电子工业出版社,2011.2周品.MATLAB数字图像处理北京，清华大学出版社，2012源代码：1. imgl. m产生高斯噪声并用高斯滤波器和中值滤波器滤波I=imrea
12、dC lena. bmp);figured);subplot (1, 2, 1)imshow(I);title(原始图像lena. bmp*);inwrite(I,原始图像lena. bmp*);I2=imnoise(I, , gaussian , 0. 5, 0. 01);subplot(1, 2,2)imshow(I2);title(,加入gaussian噪声后的lena. bmp(u=0. 5, s*2=0. 01)J;imwrite(I2,加入gaussian噪声后的lena. bmp(u=0. 5, s*2=0. 01). bmp*);figure (2);subplot 2, 1)
13、imshow(I);title (原始图像luna. bmp);subplot (2, 2, 2)imshow(I2);title( lena加入gaussian噪声后的(u=0. 5, s*2=0. 01). bmp);n=5;a=ones(n, n);p=size(I2);xl二double (12);x2=xl;for i=l:p(l)-n+lfor j=l:p(2)-n+lc=xl (i:i+(n-1), j:j+(n-l);e=c(l,:);for u=2:n e=e, c(u,:);endmm=median(e);x2(i+(n-l)/2, j+(n-l)/2)=mm;endend
14、I3=uint8(x2);subplot (2, 2, 3)imshow(I3);title(中值派波(5x5) ();imwrite(I3,中值滤波(5x5) . bmp*);12, map二imreadC I川入gaussian噪声后的lena. bmp(u=0. 5, s*2=0. 01). bmp,);k=1.5;Img=double(I2);n=5;nl=floor(n+l)/2); for i=l:nfor j=l:nb(i, j) =exp(-(i-nl)2+ (j-nl)*2)/(4*k)/(4*pi*k); endendImgl=conv2(Img, b, same);d=u
15、int8(Imgl);subplot (2,2,4);imshow(d, map);title(高斯滤波5x5)2.img2.m运动模糊+高斯噪声、维纳滤波、约束最小二乘滤波I=imread(* lena. bmp);figured);subplot (1, 2, 1)imshow(I);title (? lena. bmp原始图像);inwrite (I, lena原始图像.bmp);f=double(I);F=fft2(f);F=fftshift(F);M,N=size(F);a=0. 1;b=0. 1;T=l;for u=l:Mfor v=l:NH (u, v) = (T/(pi* (u
16、*a+v*b) *sin(pi* (u*a+v*b) *exp (-sqrt (-1) *pi* (u*a+v*b); G(u, v) =H (u, v) *F (u, v);endendG=ifftshift(G);g=ifft2(G);g=256. *g. /max (max(g);g=uint8(real(g);subplot(1,2,2);imshow(g);title。运动模糊化lena. bmp);imwrite(g, lena运动模糊的结果.bmp,);figure (2)subplot (1, 2, 1);imshow(g);title。运动模糊化lena. bmp);imwr
17、ite (g, lena运动模糊的结果.bmp);I2=imnoise(g, gaussian , 0, 0, 01);subplot(1, 2,2)imshow(I2);title(,模糊lena. b【np加入高斯噪声(u=0, s*2=0. 01);imwrite (12,模糊lena. bmp/jl 入高斯噪声(u=0, s*2=0. 01). bmp，);figure (3)I=imread(* lena. bmp);h=fspecial ( motion, 50, 45);Il=imfilter (I, h, circular, conv);I2=imnoise(Il, gauss
18、ian, 0, 0. 01);subplot (1,2, 1)imshow(I2);title( lena运动模糊+高斯噪声);imwrite (12, lena运动模糊+高斯噪声.bmp);gl=double(I2);Gl=fft2(gl);Gl=fftshift(Gl);M, N=size(Gl);a=0. 1;b=0. 1;T=l;K=0. 06;for u=l:Mfor v=l:NHl (u, v) = (T/(pi*(u*a+v*b)*sin(pi*(u*a+v*b)*exp(-sqrt (-l)*pi*(u*a+v*b); F(u, v)=l/Hl (u, v)*(abs (Hl
19、(u, v) *2/(abs (Hl (u, v) *2+K)*G1 (u, v);endendF=ifftshift；f=ifft2(F);f=256. *f. /max (max(f);f=uint8(real(f);subplot(1, 2, 2)imshow(f);titleC维纳滤波的结果(K=0. 06),);imwrite (f, 维纳滋波的结果(K=0. 06). bmp );I二imreadC lena. bmp);h=fspecialC motion, 50, 45);Il=imfilter(I, h, circular*, conv);I2=imnoise(Il, gaussian*, 0, 0. 01);figure (4);subplot (1,2,1)imshow(I2);title。lena运动模糊+高斯噪声);imwrite(I2, * len&m丁-高.斯噪声.bmp);V=0. 0001;No i s ePower=V*prod(s i z e (I);g, LAGRA=deconvreg(11, h, NoisePower);subplot(1, 2, 2)imshow(g);titleC约束最小二乘/波的姑果);imwrite(g,约束最小二乘滤波的结果.bmp);（注：可编辑下载，若有不当之处，请指正，谢谢!）
