数字图像处理试验报告概况

1、数字图像处理实验报告实验一数字图像根本操作及灰度调整一、实验目的1)掌握读、写图像的根本方法.2)掌握MATLAB语言中图像数据与信息的读取方法.3)理解图像灰度变换处理在图像增强的作用.4)掌握绘制灰度直方图的方法,理解灰度直方图的灰度变换及均衡化的方法.二、实验内容与要求1,熟悉MATLAB言中对图像数据读取,显示等根本函数特别需要熟悉以下命令：熟悉imread()函数、imwrite()函数、size()函数、Subplot()函数、Figure()函数.1) 将MATLAB目录下work文件夹中的forest.tif图像文件读出,用到imread,imfinfo等文件,观察一以下图像数

2、据,了解一下数字图像在MATLAB中的处理就是处理一个矩阵.将这个图像显示出来(用imshow).尝试修改map颜色矩阵的值,再将图像显示出来,观察图像颜色的变化.2) 将MATLAB目录下work文件夹中的b747.jpg图像文件读出,用rgb2gray()将其转化为灰度图像,记为变量Bo2 .图像灰度变换处理在图像增强的作用读入不同情况的图像,请自己编程和调用Matlab函数用常用灰度变换函数对输入图像进行灰度变换,比拟相应的处理效果.3 .绘制图像灰度直方图的方法,对图像进行均衡化处理请自己编程和调用Matlab函数完成如下实验.1)显示B的图像及灰度直方图,可以发现其灰度值集中在一段区

3、域,用imadjust函数将它的灰度值调整到0,1之间,并观察调整后的图像与原图像的差异,调整后的灰度直方图与原灰度直方图的区别.2)对B进行直方图均衡化处理,试比拟与源图的异同.3对B进行如下图的分段线形变换处理,试比拟与直方图均衡化处理的异同.图1.1分段线性变换函数三、实验原理与算法分析1.灰度变换灰度变换是图像增强的一种重要手段,它常用于改变图象的灰度范围及分布,是图象数字化及图象显示的重要工具.1图像反转灰度级范围为0,L-1的图像反转可由下式获得s=L-1-r2对数运算：有时原图的动态范围太大,超出某些显示设备的允许动态范围,如直接使用原图,那么一局部细节可能丧失.解决的方法是对原

4、图进行灰度压缩,如对数变换：s=clog1+r,c为常数,r>03哥次变换：s=cr,c-0,:-04比照拉伸：在实际应用中,为了突出图像中感兴趣的研究对象,常常要求局部扩展拉伸某一范围的灰度值,或对不同范围的灰度值进行不同的拉伸处理,即分段线性拉伸:其对应的数学表达式为:ab'-a'b-aMa式")="f,y)-h+b',b<fry)<M2.直方图均衡化灰度直方图的横坐标是灰度级,纵坐标是该灰度级出现的频度,它是图像最根本的统计特征.依据定义,在离散形式下,用rk代表离散灰度级,用Pr(rk)代表Pr(r),并且有下式成立：Pr(

5、rk)J0<rk<1k=0,1,2,l-1n式中：nk为图像中出现rk级灰度的像素数,n是图像像素总数,而n/n即为频数.直方图均衡化处理是以累积分布函数变换法为根底的直方图修正法.假定变换函数为r(a)Lena图像(b)Lena图像的直方图s=T(r)=0Pr()d图1.2Lena图像及直方图当灰度级是离散值时,可用频数近似代替概率值,即Pr(rk)0k<1k=0,1,l-1n式中：l是灰度级的总数目,Pr(rk)是取第k级灰度值白概率,nk是图像中出现第k级灰度的次数,n是图像中像素总数.所以积分可以表示为以下累计分布函数(cumulativedistributionfu

6、nction,CDF)knjkSk=T(rJ/Pr(rj)03三1k=0,1,l-1j-0nj=0四、实验步骤1 .熟悉MATLA皿言中对图像数据读取,显示等根本函数1)文件读取与信息显示:loadtrees;X,map=imread('forest.tif');subimage(X,map);I=imread('forest.tif');imshow(I);imfinfo('forest.tif');2 )map颜色矩阵的修改X,map=imread('forest.tif');map1=map+map;subimage(X,m

7、ap1);3)灰度图像的转化RGB=imread('b747.jpg');B=rgb2gray(RGB);2 .图像灰度变换处理在图像增强的作用g1=imadjust(I,01,10);g2=imcomplement(g1);g3=im2uint8(mat2gray(log(1+double(I);3 .绘制图像灰度直方图的方法,对图像进行均衡化处理1)图像灰度直方图的显示与灰度调整imhist(B);J=imadjust(B,01);imhist(J);subimage(J);2)对B进行直方图均衡化处理,试比拟与原图的异同.I=imread('pout.tif

8、9;);J,T=histeq(I);figure,plot(0:255)/255,T);3)对B进行如下图的分段线形变换处理,试比拟与直方图均衡化处理的异同.x1=0:0.01:0.125;x2=0.125:0.01:0.75;x3=0.75:0.01:1;y1=2*x1;y2=0.25+0.6*(x2-0.125);y3=0.625+1.5*(x3-0.75);x=x1,x2,x3;y=y1,y2,y3;subplot(2,2,4);plot(x,y);五、实验结果分析与讨论1,熟悉MATLAB言中对图像数据读取,显示等根本函数1图像文件的读出与图像数据的观察RGBImageGrayImag

9、e图1,3真彩色图像与灰度图像显示ans=Filename:'F:MATLABR2007atoolboximagesimdemosforest.tif'FileModDate:'04-Dec-200013:57:58'FileSize:124888Format:'tifFormatversion:Width:447Height:301BitDepth:8ColorType:'indexed'Formatsignature:7373420ByteOrder:'little-endian'NewSubFileType:0Bit

10、sPerSample:8Compression:'PackBits'Photometricinterpretation:'RGBPalette'StripOffsets:17x1doubleSamplesPerPixel:1RowsPerStrip:18StripByteCounts:17x1doubleXResolution:72YResolution:72ResolutionUnit:'Inch'Colormap:256x3doublePlanarConfiguration:'Chunky'TileWidth:TileLeng

11、th:TileOffsets:TileByteCounts:Orientation:1FillOrder:1GrayResponseUnit:0.0100MaxSampleValue:255MinSampleValue:0Thresholding:1ImageDescription:'CarmanahAncientForest,BritishColumbia,Canada'map颜色矩阵修改后图像颜色的变化图1.4原图像与map巨阵值增强一倍后的图像新的颜色矩阵值变成原文件的2倍,图像明显变亮,颜色的R、G、B值增强.2图像文件转化为灰度图像图1.5真彩色图像与转化成的灰度图像

12、2.图像灰度变换处理在图像增强的作用图1.6灰度变化增强图像反转、求补、对数变换s=L-1-r得到的图像反转,对反图中对图像文件进行了根本的灰度变换,包括用式转图像的求补,以及对数变换的采用.3.绘制图像灰度直方图的方法,对图像进行均衡化处理1图像灰度及灰度直方图的调整图1.7灰度范围的调整与直方图显示在原始图像中,直方图的组成成分集中在高灰度等级亮一侧,且图像灰度范围为0,1,故将灰度值调整到0,1间后直方图无明显变化.类似的,将灰度值调整到0,0.5时,整个图像变暗,直方图横向压缩1倍.2直方图均衡化OriginalGreyImageGreyRangeEqualizedImageGreyV

13、alueGreyValue图1.8原图像与灰度直方图均衡化结果原图像中目标物的灰度主要集中于高亮度局部,而且象素总数较多,所占的灰度等级较少.经过直方图均衡后,目标物的所占的灰度等级得到扩展,比照度增强,使整个图像得到增强.数字图像均衡化后,其直方图并非完全均匀分布,这是由于图像的象素个数和灰度等级均为离散值,而且均衡化使灰度级并归.因此,均衡化后,其直方图并非完全均匀分布.3对B进行如下图的分段线性变换处理OriginalGreyImageLinearTransformedImage图1.9原图像与灰度直方图均衡化结果通过在所关心范围内为所有灰度值指定一个较高值,而为其他灰度指定一个较低值,

14、或将所需范围变亮,分段线性变换可提升图像中特定灰度范围的亮度,常用于图像特征值的提取.这里将原始图像位于0.125,0.75间的灰度值调低,放大其余的灰度值,突出显示图像低频域和高频域的局部.六、参考文献美RafaelC.Gonzalez.数字图像处理第二版M.阮秋琦阮宇智,译.北京：电子工业出版社,2003.3.实验二数字图像的空间域滤波和频域滤波一、实验目的1 .掌握图像滤波的根本定义及目的.2 .理解空间域滤波的根本原理及方法.3 .掌握进行图像的空域滤波的方法.4 .掌握傅立叶变换及逆变换的根本原理方法.5 .理解频域滤波的根本原理及方法.6 .掌握进行图像的频域滤波的方法.二、实验内

15、容与要求1.平滑空间滤波：1读出eight.tif这幅图像,给这幅图像分别参加椒盐噪声和高斯噪声后并与前一张图显示在同一图像窗口中.2对参加噪声图像选用不同的平滑低通模板做运算,比照不同模板所形成的效果,要求在同一窗口中显示.3使用函数imfilter时,分别采用不同的填充方法或边界选项,如零填充、replicate'、'symmetric'、'circular'进行低通滤波,显示处理后的图像.10次,20次均值滤波,查看其特点,显averages型生成均值滤波器.和中值滤波法对有噪声的图像做处4运用for循环,将加有椒盐噪声的图像进行示均值处理后的图像

16、提示利用fspecial函数的5对参加椒盐噪声的图像分别采用均值滤波法,理,要求在同一窗口中显示结果.6)2.1)自己设计平滑空间滤波器,并将其对噪声图像进行处理,显示处理后的图像.锐化空间滤波读出blurry_moon.tif这幅图像,采用3X3的拉普拉斯算子w=1,1,1;1-81;1,1,1对其进行滤波.2编写函数的拉普拉斯算子w=genlaplacian(n),自动产生任一奇数尺寸n的拉普拉斯算子,如5X53分别采用111115X5,111119X9,11-241115X111111111115和25X25大小的拉普拉斯算子对blurry_moon.tif进行锐化滤波,并利用式gx,y

17、=fx,y$2fx,y完成图像的锐化增强,观察其有何不同,要求在同一窗口中显示.4采用不同的梯度算子对blurry_moon.tif进行锐化滤波,并比拟其效果.5自己设计锐化空间滤波器,并将其对噪声图像进行处理,显示处理后的图像;3 .傅立叶变换1读出woman.tif这幅图像,对其进行快速傅立叶变换,分别显示其幅度图像和相位图像.仅对相位局部进行傅立叶反变换后查看结果图像.2仅对幅度局部进行傅立叶反变换后查看结果图像.3将图像的傅立叶变换F置为其共轲后进行反变换,比拟新生成图像与原始图像的差异.4 .平滑频域滤波1设计理想低通滤波器、巴特沃斯低通滤波器和高斯低通滤波器,截至频率自选,分别给出

18、各种滤波器的透视图.2读出test_pattern.tif这幅图像,分别采用理想低通滤波器、巴特沃斯低通滤波器和高斯低通滤波器对其进行滤波截至频率自选,再做反变换,观察不同的截止频率下采用不同低通滤波器得到的图像与原图像的区别,特别注意振铃效应.提示:1在频率域滤波同样要注意到填充问题；2注意到-1x+y；5.锐化频域滤波1设计理想高通滤波器、巴特沃斯高通滤波器和高斯高通滤波器,截至频率自选,分别给出各种滤波器的透视图.2读出test_pattern.tif这幅图像,分别采用理想高通滤波器、巴特沃斯高通滤波器和高斯高通滤波器对其进行滤波截至频率自选,再做反变换,观察不同的截止频率下采用不同高通

19、滤波器得到的图像与原图像的区别.三、实验原理与算法分析1 .空间域增强处理图像每一个像素的取值空域滤波根本上是让图像在从而改变输出图像的频率分空间域滤波是在图像空间中借助模板对图像进行领域操作,都是根据模板对输入像素相应领域内的像素值进行计算得到的.频域空间内某个范围的分量受到抑制,同时保证其他分量不变,布,到达增强图像的目的.空域滤波一般分为线性滤波和非线性滤波两类.线性滤波器的设计常基于对傅立叶变换的分析,非线性空域滤波器那么一般直接对领域进行操作.各种空域滤波器根据功能主要分为平滑滤波器和锐化滤波器.平滑可用低通来实现,平滑的目的可分为两类：一类是模糊,目的是在提取较大的目标前去除太小的

20、细节或将目标内的小肩端连接起来；另一类是消除噪声.锐化可用高通滤波来实现,锐化的目的是为了增强被模糊的细节.结合这两种分类方法,可将空间滤波增强分为四类：线性平滑滤波器低通非线性平滑滤波器低通线性锐化滤波器高通非线性锐化滤波器高通空间滤波器都是基于模板卷积,其主要工作步骤是：1将模板在图中移动,并将模板中央与图中某个像素位置重合;2将模板上的系数与模板下对应的像素相乘；3将所有乘积相加；4将和模板的输出响应赋给图中对应模板中央位置的像素.2 .平滑滤波器1)线性平滑滤波器线性低通平滑滤波器也称为均值滤波器,这种滤波器的所有系数都是正数,对3X3的模板来说,最简单的是取所有系数为1,为了保持输出

21、图像任然在原来图像的灰度值范围内,模板与象素邻域的乘积都要除以9.MATLAB提供了fspecial函数生成滤波时所用的模板,并提供filter2函数用指定的滤波器模板对图像进行运算.函数fspecial的语法格式为：h=fspecial(type);h=fspecial(type,parameters);其中参数type指定滤波器的种类,parameters是与滤波器种类有关的具体参数.Type睨明k三球均值溺ft,如果锚域为方阵.可hsia为标量.否剧由两元素向量Min指定邻域的打效和列数,di:kradium由3加济空+1)个边的阑形均值魂波器标准偏差为力吐H大小为的高斯/通滤波需lap

22、laciiDilpha系数由“口小】0)决定的二维拉普拉斯博灌Iceh力ze.二ig皿标准偏差为立大小为的离斯被波旋转对林拉氏算子motion加口-theta按角度位口移动len个教宗的运动沌选器pieivin无近似计重垂直梯度的水平功修强词异子sobel无近似计算垂直梯度光滑效应的水平边除强调算子alpha根据ilplu决定的拉氏箱子创立的推桢滤波器表2.1MATLAB中预定义的滤波器种类MATLAB提供了一个函数imnoise来给图像增添噪声,其语法格式为：J=imnoise(I,type);J=imnoise(I,type,parameters);参数type指定噪声的种类,parame

23、ters是与噪声种类有关的具体参数.参数的种类见表2.2.押类参数说明gaussianv均色为n】方茏为t的高斯嶂声localvarV均值为.+方差为中的高斯白噪声passion无泊松噪声saltpeppet椒盐啜声speckleV均值为0,方差为中的均匀分布随机啜声表2.2噪声种类及参数说明2)非线性平滑滤波器中值滤波器是一种常用的非线性平滑滤波器,其滤波原理与均值滤波器方法类似,但计算的非加权求和,而是把领域中的图像的象素按灰度级进行排序,然后选择改组的中间值作为输出象素值.MATLAB提供了medfilt2函数来实现中值滤波,其语法格式为：B=medfilt2(A,mn);B=medfi

24、lt2(A);其中,A是原图象,B是中值滤波后输出的图像.mn指定滤波模板的大小,默认模板为3X3.3 .锐化滤波器图像平滑往往使图像中的边界、轮廓变得模糊,为了减少这类不利效果的影响,需要利用图像锐化技术,使图像的边缘变得清楚.1)线性锐化滤波器线性高通滤波器是最常用的线性锐化滤波器.这种滤波器的中央系数都是正的,而周围的系数都是负的,所有的系数之和为0.4 3X3的模板来说,典型的系数取值为：-1-1-1;-18-1;-1-1-1事实上这是拉普拉斯算子.语句h=-fspecial(laplacian',0.5)得到的拉普拉斯算子为：h=-0.3333-0.3333-0.3333-0

25、.33332.6667-0.3333-0.3333-0.3333-0.33332)非线性锐化滤波邻域平均可以模糊图像,由于平均对应积分,所以利用微分可以锐化图像.图像处理中最常用的微分方法是利用梯度.常用的空域非线性锐化滤波微分算子有sobel算子、prewitt算子、10g算子等.4 .频域增强频域增强是利用图像变换方法将原来的图像空间中的图像以某种形式转换到其他空间中,然后利用该空间的特有性质方便地进行图像处理,最后再转换回原来的图像空间中,从而得到处理后的图像.频域增强的主要步骤是:选择变换方法,将输入图像变换到频域空间.在频域空间中,根据处理目的设计一个转移函数,并进行处理.将所得结果

26、用反变换得到增强的图像.常用的频域增强方法有低通滤波和高通滤波.5 .低通滤波图像的能量大局部集中在幅度谱的低频和中频局部,而图像的边缘和噪声对应于高频部分.因此能降低高频成分幅度的滤波器就能减弱噪声的影响.由卷积定理,在频域实现低通滤波的数学表达式：1)2)G(u,v)=H(u,v)F(u,v)理想低通滤波器(ILPF)flon(n.v)<L>c>Dn3)指数型低通滤波器ELPF巴特沃斯低通滤波器BLPF6 .高通滤波使低频分量受到削弱.高通滤波器和低由于图像中的细节局部与其高频分量相对应,所以高通滤波可以对图像进行锐化处理.高通滤波与低通滤波相反,它是高频分量用通滤波器相

27、似,其转移函数分别为：D(ilt)<I).Du.v)>Dc1） 理想高通滤波器IHPF=2巴特沃斯高通滤波器BLPFAD(jkv)39=1 +3指数型高通滤波器ELPF=e忑图像经过高通滤波处理后,会丧失许多低频信息,所以图像的平滑区根本上会消失.所以,可以采用高频增强滤波来弥补.高频增强滤波就是在设计滤波传递函数时,加上一个大于0小于1的常数c,即：H'u,v=Hu,v+c四、实验步骤1 .平滑空间滤波：1) 椒盐噪声和高斯噪声的参加I=imread('eight.tif');imshow(I);J=imnoise(I,'salt&pepp

28、er',0.05);%noisedensity=0.05K=imnoise(I,'gaussian',0.01,0.01);2)不同的平滑(低通)模板的处理H=fspecial('sobel');Sobel=imfilter(I,H,'replicate');H=fspecial('laplacian',0.4);lap=imfilter(I,H,'replicate');H=fspecial('gaussian',33,0.5);gaussian=imfilter(I,H,'rep

29、licate');3)不同填充方法的滤波originalRGB=imread('peppers.png');h=fspecial('motion',50,45);%motionblurredfilteredRGB=imfilter(originalRGB,h);boundaryReplicateRGB=imfilter(originalRGB,h,'replicate');boundary0RGB=imfilter(originalRGB,h,'x');boundary0RGB=imfilter(originalRGB,h

30、,0);boundarysymmetricRGB=imfilter(originalRGB,h,'symmetric');boundarycircularRGB=imfilter(originalRGB,h,'circular');4)屡次均值滤波J=imnoise(I,'salt&pepper',0.05);h=fspecial('average');%AveragingFilteringJ1=imfilter(J,h);fori=1:10J1=imfilter(J,h);endfori=1:20J2=imfilter(

31、J,h);5)均值、中值滤波h=fspecial('average');%AveragingFilteringJ1=imfilter(J,h);J2=medfilt2(J);%MedianFiltering6)自行设计平滑空间滤波器domain=00800;00800;88888;00800;K1=ordfilt2(J,5,domain);2.平滑空间滤波:1)3X3的拉普拉斯算子滤波I=imread('blurry_moon.tif');T=double(I);subplot(1,2,1),imshow(T,口);title('OriginalImag

32、e');w=1,1,1;1,-8,1;1,1,1;K=conv2(T,w,'same');2)奇数尺寸拉普拉斯算子随机产生函数functionw=genlaplacian(n)%ComputestheLaplacianoperatorw=ones(n);x=ceil(n/2);w(x,x)=-1*(n*n-1);3)不同尺寸拉普拉斯算子滤波以及图像增强w1=genlaplacian(5);I=imread('blurry_moon.tif');T=double(I);K=conv2(T,w1,'same');J=T-K;4)不同尺寸梯度算

33、子的锐化滤波I,map=imread('blurry_moon.tif');I=double(I);%gradientcalculation%matrixGx,Gy=gradient(I);G=sqrt(Gx.*Gx+Gy.*Gy);J1=G;%gradient1imshow(J1,map);J2=I;%gradient2K=find(G>=7);J2(K)=G(K);imshow(J2,map);J3=I;%gradient3K=find(G>=7);J3(K)=255;imshow(J3,map);J4=I;%gradient4K=find(G<=7);J

34、4(K)=255;imshow(J4,map);J5=I;%gradient5K=find(G<=7);J5(K)=0;Q=find(G>=7);J5(Q)=255;imshow(J5,map);5)自行设计锐化空间滤波器domain=88088;88088;00000;88088;88088;K1=ordfilt2(J,5,domain);3 .傅立叶变换1)图像的快速傅立叶变换,分别显示其幅度图像和相位图像.F=imread('woman.tif');F1=fft2(F);F2=log(1+abs(F1);%amplitudespectrumF3=fftshif

35、t(F1);imshow(log(1+abs(F3),);F4=angle(F1);%phasespectrum2)相位局部进行傅立叶反变换.F1=fft2(F);i=sqrt(-1);f2=ifft2(exp(i*angle(F1);imshow(real(f2),);3)幅度局部进行傅立叶反变换.f1=ifft2(abs(F1);imshow(log(1+abs(f1),口);4)将图像的傅立叶变换F置为其共轲后进行反变换F1=fft2(F);F2=log(1+abs(F1);F3=fftshift(F1);F4=angle(F1);F5=-F4%amplitudespectrum%pha

36、sespectrumF6=double(F3*exp(F4);%thecomplexconjugateofthefouriertransform%inversefouriertransformF7=ifft2(F6);imshow(real(F7),口);4 .平滑频域滤波理想低通滤波器、巴特沃斯低通滤波器和高斯低通滤波器的设计与滤波理想低通滤波器I=imread('test_pattern.tif');f=double(I);%chageintodoubleasMATLABdoesn'tsupporcalculation%ofimageinunsignedinttyp

37、eg=fft2(f);%fouriertransformg=fftshift(g);%zero-frequencyareacentralizedM,N=size(g);d0=100;%cutofffrequencym=fix(M/2);n=fix(N/2);fori=1:Mforj=1:Nd=sqrt(i-m)A2+(j-n)A2);if(d<=d0)h=1;elseh=0;endresult(i,j)=h*g(i,j);endendresult=ifftshift(result);J1=ifft2(result);J2=uint8(real(J1);imshow(J2)巴特沃斯低通滤波

38、器(二阶)I=imread('test_pattern.tif');f=double(I);g=fft2(f);g=fftshift(g);M,N=size(g);nn=2;%2-gradeButterworthlowpssfilterd0=100;m=fix(M/2);n=fix(N/2);fori=1:Mforj=1:Nd=sqrt(i-m)A2+(j-n)A2);h=1/(1+0.414*(d/d0)A(2*nn);%filtertransformfunction%h=1./(1+(d./d0).A(2*n)%h=exp(-(d.A2)./(2*(d0A2);result

39、(i,j)=h*g(i,j);endendresult=ifftshift(result);J1=ifft2(result);J2=uint8(real(J1);imshow(J2);高斯低通滤波I=imread('test_pattern.tif');f=double(I);g=fft2(f);g=fftshift(g);M,N=size(g);d0=100;m=fix(M/2);n=fix(N/2);fori=1:Mforj=1:Nd=sqrt(i-m)A2+(j-n)A2);h=exp(-(d.A2)./(2*(d0A2);%gaussianfiltertransform

40、result(i,j)=h*g(i,j);endendresult=ifftshift(result);J1=ifft2(result);J2=uint8(real(J1);5 .锐化频域滤波理想高通滤波器、巴特沃斯高通滤波器和高斯高通滤波器的设计与滤波理想高通滤波器I=imread('test_pattern.tif');f=double(I);g=fft2(f);g=fftshift(g);M,N=size(g);d0=80;%d0=15,25,80m=fix(M/2);n=fix(N/2);fori=1:Mforj=1:Nd=sqrt(i-m)A2+(j-n)A2);if

41、(d>=d0)h=1;elseh=0;endresult(i,j)=h*g(i,j);endendresult=ifftshift(result);J1=ifft2(result);J2=uint8(real(J1);巴特沃斯高通滤波器I=imread('test_pattern.tif');f=double(I);g=fft2(f);g=fftshift(g);M,N=size(g);nn=2;%2-gradenutterworthhighpassfilterd0=80;m=fix(M/2);n=fix(N/2);fori=1:Mforj=1:Nd=sqrt(i-m)A

42、2+(j-n)A2);if(d=0)h=0;elseh=1/(1+0.414*(d0/d)A(2*nn);%transformfuctioncalculationendresult(i,j)=h*g(i,j);endendresult=ifftshift(result);J1=ifft2(result);J2=uint8(real(J1);高斯高通滤波器I=imread('test_pattern.tif');f=double(I);g=fft2(f);g=fftshift(g);M,N=size(g);d0=80;m=fix(M/2);n=fix(N/2);fori=1:Mf

43、orj=1:Nd=sqrt(i-m)A2+(j-n)A2);h=1-exp(-(d.A2)./(2*(d0A2);%gaussianfiltertransformresult(i,j)=h*g(i,j);endendresult=ifftshift(result);J1=ifft2(result);J2=uint8(real(J1);五、实验结果分析与讨论1.平滑空间滤波：1)读出eight.tif这幅图像,给这幅图像分别参加椒盐噪声和高斯噪声后并与前一张图显示在同一图像窗口中.OriginalImageSalt&PepperNoiseGaussianNoise图2.1初始图像及椒盐噪

44、声图像、高斯噪声污染图2对参加噪声图像选用不同的平滑低通模板做运算,比照不同模板所形成的效果,要求在同一窗口中显示.图2.2原图像及各类低通滤波处理图像3） 使用函数imfilter时,分别采用不同的填充方法或边界选项,如零填充、'replicate'、'symmetric'、'circular'进行低通滤波,显示处理后的图像.图2.3原图像及运动模糊图像图2.4函数imfilter各填充方式处理图像,显4） 运用for循环,将加有椒盐噪声的图像进行10次,20次均值滤波,查看其特点图2.5椒盐噪声污染图像经10次、20次均值滤波图像示均值处理后

45、的图像.由图2.5可得,20次滤波后的效果明显好于10次滤波,但模糊程度也更强.5对参加椒盐噪声的图像分别采用均值滤波法,和中值滤波法对有噪声的图像做处理,要求在同一窗口中显示结果.图2.6椒盐噪声污染图像及均值、中值滤波图像从图2.6中可以看出,对于椒盐噪声污染的图像处理,中值滤波效果要明显好于均值滤波.经均值滤波器处理后的图像比均值滤波器中结果图像更加模糊.6设计平滑空间滤波器,并将其对噪声图像进行处理,显示处理后的图像.Salt&PepperNoiseImage5*5SmoothingSpatialFiteredImage图2.7椒盐噪声污染图像及5*5平滑滤波器掩模掩模值为w=

46、1/25*11111;11111;11111;11111;11111Salt3tPepperNoisedImageLov/passFilteredImage图2.8椒盐噪声污染图像及5*5平滑滤波器掩模掩模值为w=00800;00800;88888;00800;008002.锐化空间滤波1采用3X3的拉普拉斯算子w=1,1,1;1-81;1,1,1滤波2)图2.9初始图像与拉普拉斯算子锐化图像编写函数w=genlaplacian(n),自动产生任一奇数尺寸普拉斯算子：n的拉普拉斯算子,如5X5的拉1111111111-24111111111111113) 分别采用5X5,9X9,15X15和2

47、5X25大小的拉普拉斯算子对blurry_moon.tif进行锐2化滤波,并利用式g(x,y)=f(x,V)Nf(x,y)完成图像的锐化增强,观察其有何不同,要求在同一窗口中显示.图2.10初始图像与不同拉普拉斯算子锐化图像图像锐化的实质是将原图像与梯度信息叠加,相当于对目标物的边缘进行了增强.LaplacianImageEnhancedLaplacianImage图2.11拉普拉斯算子锐化与锐化增强图像4) 采用不同的梯度算子对blurrymoon.tif进行锐化滤波,并比拟其效果Operator!ImageOperator2ImageOperatorsImageOperator4Image

48、OperatoiSImage图2.12原始图像与不同梯度子锐化图像作为二阶微分算子,拉普拉斯变换在图像细节的增强处理上有明显的优点,但会产生更多的噪声.梯度变换在灰度变化区域的响应更强,但对噪声和细节的响应比拉普拉斯变换弱.5自己设计锐化空间滤波器,并将其对噪声图像进行处理,显示处理后的图像;图2.13原始图像与不同边缘锐化图像3.傅立叶变换1读出woman.tif这幅图像,对其进行快速傅立叶变换,分别显示其幅度图像和相位图像.OriginalImageAmplitudeSpectrumZero-FrequencyCentralizedSpectrumPhaseSpectrum图2.14原始图

49、像与快速傅立叶变换图像2仅对相位局部进行傅立叶反变换后查看结果图像.OrgmalImageRestorationImage图2.15原始图像与对全部信息进行傅立叶逆变换结果3)Amplitude-basedRestorationImagePhase-basedRestorationImage图2.16仅对相位信息进行傅立叶逆变换结果相位谱决定了图像信号中各频率分量的位置.仅对幅度局部进行傅立叶反变换后查看结果图像.图2.17仅对幅度信息进行傅立叶逆变换结果4)将图像的傅立叶变换F置为其共轲后进行反变换,比拟新生成图像与原始图像的差异.图2.18共轲傅立叶逆变换结果傅立叶变换的相位谱为对称的,原变换与其共轲变换间仅频率谱互为相反.4.平滑频域滤波1设计理想低通滤波器、巴特沃斯低通滤波器和高斯低通滤波器图2.19理想低通滤波器透视图图2.20巴特沃斯低通滤波器透视图图2.21高斯低通滤波器透视图IdealLowpassFilterd0=15IdealLowpassFilterd0=30IdealLowpassFilterd0=100al