数字图像处理 空域滤波 实验报告

1、数字图像处理实验三空域滤波实验报告一、实验目的了解空域滤波的方法；掌握几种模板的基本原理。二、实验内容使用函数fspecial()生成几种特定的模板；使用函数inifUter()配合模板对图象数据进行二维卷积:比较各种滤波器的效果。三、实验步骤及结果线性平滑(低通)滤波器1用h=fspecial(averagej得到的h为3X3的邻域平均模板，然后用h进行平滑处理。x=iniiead(cameraman,tif);h=fspecial(caverage);%h=fspecial(average7,7);y=inifilter(xji);figuie(l);iiiishow(x)figure；i

2、iiishow(y)原图：滤波效果:改变模板人小重试x=inuead(,cameianian.tif);%h=fspecial();h=fspecial(raverage7,刀);y=inifiltei(x,h);figuie(l);iiiishow(x)figure；imshow(y)原图：处理效果:Gaussian平滑改变fspecial()的参数为高斯函数，再进行同样的平滑处理，观察其结果。x=inuead();h=fspecial(gaussian,)/3X3的邻域平均模板(7X7的邻域平均模板h=fspecial(gaussian,7,7)y=inifilter(xji);实际程序：

3、x=inuead(,cameianian.tif);h=fspecial(gaussianj%3X3的邻域平均模板(7X7的邻域平均模板h=fspecial(aussiaif,7,7)y=inifilter(xji);subplot(121);iiiishow(x);subplot(122);iiiishow(y)；模板人小为3,左为原图，右为处理后效果模板人小为7,7,左为原图，右为处理后效呆模板人小为11,11，左为原图，右为处理后效果改变模板人小为7,7、9,9、11,11,观察实验得到的结构冇何变化？结果：高斯模板对结果影响较小，为什么？答：简单地说，这是由于高斯滤波中越靠近中心点的点

4、的权重越人，在很人程度上减弱了较远点的干扰，所以即使模版大小较人，其效果仍远远好过“不知亲疏远近变化”的邻域平均模板。具体地说，高斯滤波是一种线性平滑滤波，适用于消除高斯噪声，广泛应用于图像处理的减噪过程。通俗的讲，高斯滤波就是对整幅图像进行加权平均的过程，每一个像素点的值，都由其本身和邻域内的其他像素值经过加权平均后得到。高斯滤波的具体操作是：用一个模板（或称卷积、掩模）打描图像中的每一个像素，用模板确定的邻域内像素的加权平均灰度值去替代模板中心像素点的值。一般的模板为3X3或5X5人小，其权值分布如下图：1X16121242121127314741414261447264126741626

5、ltf414741若使用3X3模板，则计算公式如下：g(x,y)=f(x-l,y-l)+f(x-1,y+1)+f(x+1,y-l)+f(x+1,v+1)+f(x-l,y)+f(x,y-1)+f(x+1,y)+f(x,y+1)*2+f(x,y)*4/16;其中，f(x,y)为图像中(x,y)点的灰度值，g(x、y)为该点经过高斯滤波后的值。用m亡sh函数直观观察模板hfspeciaXgaussiaiiMOOJ);x刃=meshgiid(l:100);mesh(xyhl);h2=fspecial(lgaussiaii100,10);figuremesh(x,yji2)非线性平滑滤波器中值(ndia

6、n)滤波器是非线性滤波器的一种，它工作原理如下：将模板在图中漫游，并将模板中心与图中某个彖素位置重合；读取模板下各对应彖素的灰度值；将这些灰度值从小到人排成一列；找出这些值中排在中间的1个；将这个中间值赋给对应模板中心位置的彖素。非线性中值滤波步骤对图象添加胡椒盐噪声(salt&peppei),然后分别用前面的均值滤波器和中值滤波器medfilt2()进行去噪声处理，并给出比较结果。中值滤波和均值滤波代码：x=inuead(,eight.tif);h=imnoise(x/salt&pepper*,0.02);%改变噪声的大小试一试?y=medfilt2(h);subplot(2,2,l)iii

7、ishow(x)subplot(2,2,2)iiiishow(y)y1=fspecial(caverage?;subplot(2,2,3)iiiishow(x)subplot(2,2,4)iiiishow(yl)椒盐噪声人小是0.02,左侧两图是原图，右上是中值滤波，右下是均值滤波椒盐噪声人小是0.10,左侧两图是原图，右上是中值滤波，右下是均值滤波发现中值滤波对图像边缘的噪声无能为力，而均值滤波则可以较好地去除边缘噪声。对图像的清晰度，只能说差别不人，半斤八两，都勉强及格吧。线性锐化滤波器线性锐化滤波器的模板为h=-1-1-1；-18-1;-1-1-1；用该模板对图象进行锐化处理；代码：x=

8、inuead(,cameianian.tif);h=-1-1-1;-18-1;-1-1-1y=inifilter(xji);subplot(l,2,l)iiiishow(x)subplot(1,2,2)iiiishow(y)也可以fspecial(rlaplacian)函数得到锐化模板；代码：x=inuead(,cameranian.tif);h=fspeciaXlaplacian*)y=inifiltei(x,h);subplot(l,2,l)iiiishow(x)subplot(1,2,2)imshow(y)sobel算子索贝尔(sob亡1)算子hl=-l0l;-202;-l01;h2=1

9、2l;000;-l-2-1;分别使用以上算子对图彖进行处理可以得到图彖的边缘。x=inuead(,cameianian.tif);hl=-l01;-202;-101;y=inifilter(xji);subplot(l,2,l)imshow(x)subplot(1,2,2)iiiishow(y)x=inuead(,cameianian.tif);112=12l;000;-l-2-1;y=inifilter(xji2);subplot(l,2,l)imshow(x)subplot(1,2,2)iiiishow(y)Camiv算子Helpedge”canny、sob亡1、PrewittsLOG等Canny边缘检测算子是JolmF.Caimy于1986年开发出来的一个多级边缘检测算法。Sobel算子主要用作边缘检测。在技术上，它是一离散性差分算子，用来运算图像亮度函数的梯度之近似值。在图像的任何一点使用此算子，将会产生对应的梯度矢量或是其法矢量Prewitt算子是一种一阶微分算子的边缘检测，利用像素点上卞、左右邻点的灰度差，在边缘处达到极值检测边缘，