实验三 数字图像的空间域滤波Word下载.docx

实验三 数字图像的空间域滤波Word下载.docx

《实验三 数字图像的空间域滤波Word下载.docx》由会员分享，可在线阅读，更多相关《实验三 数字图像的空间域滤波Word下载.docx（27页珍藏版）》请在冰豆网上搜索。

内,模板与象素邻域得乘积都要除以9。

MATLAB提供了fspecial函数生成滤波时所用得模板,并提供filter2函数用指定得滤

波器模板对图像进行运算。

函数fspecial得语法格式为:

h=fspecial（type）;

h=fspecial（type,parameters）;

其中参数type指定滤波器得种类,parameters就是与滤波器种类有关得具体参数。

表2、1MATLAB中预定义得滤波器种类

MATLAB提供了一个函数imnoise来给图像增添噪声,其语法格式为:

J=imnoise（I,type）;

J=imnoise（I,type,parameters）;

参数type指定噪声得种类,parameters就是与噪声种类有关得具体参数。

参数得种类见

表2、2。

表2、2噪声种类及参数说明

2）非线性平滑滤波器

中值滤波器就是一种常用得非线性平滑滤波器,其滤波原理与均值滤波器方法类似,但计

算得非加权求与,而就是把领域中得图像得象素按灰度级进行排序,然后选择改组得中间值作为输出象素值。

MATLAB提供了medfilt2函数来实现中值滤波,其语法格式为:

B=medfilt2（A,[mn]）;

B=medfilt2（A）;

其中,A就是原图象,B就是中值滤波后输出得图像。

[mn]指定滤波模板得大小,默认模

板为3×

3。

3.锐化滤波器

图像平滑往往使图像中得边界、轮廓变得模糊,为了减少这类不利效果得影响,需要利

用图像锐化技术,使图像得边缘变得清晰。

1）线性锐化滤波器

线性高通滤波器就是最常用得线性锐化滤波器。

这种滤波器得中心系数都就是正得,而周围

得系数都就是负得,所有得系数之与为0。

对3×

3得模板来说,典型得系数取值为:

[111;

181;

111]

事实上这就是拉普拉斯算子。

语句h=fspecial（‘laplacian’,0、5）得到得拉普拉斯算子为:

h=0、33330、33330、3333

0、33332、66670、3333

0、33330、33330、3333

2）非线性锐化滤波

邻域平均可以模糊图像,因为平均对应积分,所以利用微分可以锐化图像。

图像处理中

最常用得微分方法就是利用梯度。

常用得空域非线性锐化滤波微分算子有sobel算子、prewitt算子、log算子等。

4.频域增强

频域增强就是利用图像变换方法将原来得图像空间中得图像以某种形式转换到其她空间

中,然后利用该空间得特有性质方便地进行图像处理,最后再转换回原来得图像空间中,从

而得到处理后得图像。

频域增强得主要步骤就是:

选择变换方法,将输入图像变换到频域空间。

在频域空间中,根据处理目得设计一个转移函数,并进行处理。

将所得结果用反变换得到增强得图像。

常用得频域增强方法有低通滤波与高通滤波。

5.低通滤波

图像得能量大部分集中在幅度谱得低频与中频部分,而图像得边缘与噪声对应于高频部

分。

因此能降低高频成分幅度得滤波器就能减弱噪声得影响。

由卷积定理,在频域实现低通

滤波得数学表达式:

G（u,v）=H（u,v）F（u,v）

1）理想低通滤波器（ILPF）

2）巴特沃斯低通滤波器（BLPF）

3）指数型低通滤波器（ELPF）

6.高通滤波

由于图像中得细节部分与其高频分量相对应,所以高通滤波可以对图像进行锐化处理。

高通滤波与低通滤波相反,它就是高频分量顺利通过,使低频分量受到削弱。

高通滤波器与低通滤波器相似,其转移函数分别为:

1）理想高通滤波器（IHPF）

2）巴特沃斯高通滤波器（BLPF）

3）指数型高通滤波器（ELPF）

图像经过高通滤波处理后,会丢失许多低频信息,所以图像得平滑区基本上会消失。

所

以,可以采用高频加强滤波来弥补。

高频加强滤波就就是在设计滤波传递函数时,加上一

个大于0小于1得常数c,即:

H′（u,v）=H（u,v）+c

3.实验内容与要求

1.平滑空间滤波:

1）读出eight、tif这幅图像,给这幅图像分别加入椒盐噪声与高斯噪声后并与前一张图

显示在同一图像窗口中。

2）对加入噪声图像选用不同得平滑（低通）模板做运算,对比不同模板所形成得效果,

要求在同一窗口中显示。

3）使用函数imfilter时,分别采用不同得填充方法（或边界选项,如零填

充、’replicate’、’symmetric’、’circular’）进行低通滤波,显示处理后得图像。

4）运用for循环,将加有椒盐噪声得图像进行10次,20次均值滤波,查瞧其特点,显

示均值处理后得图像（提示:

利用fspecial函数得’average’类型生成均值滤波器）。

5）对加入椒盐噪声得图像分别采用均值滤波法,与中值滤波法对有噪声得图像做处

理,要求在同一窗口中显示结果。

6）自己设计平滑空间滤波器,并将其对噪声图像进行处理,显示处理后得图像。

2.锐化空间滤波

1）读出blurry_moon、tif这幅图像,采用3×

3得拉普拉斯算子w=[1,1,1;

1–81;

1,1,

1]对其进行滤波。

2）编写函数w=genlaplacian（n）,自动产生任一奇数尺寸n得拉普拉斯算子,如5×

5

得拉普拉斯算子

w=[11111

11111

112411

11111

11111]

3）分别采用5×

5,9×

9,15×

15与25×

25大小得拉普拉斯算子对blurry_moon、tif进

行锐化滤波,并利用式完成图像得锐化增强,观察其有何不同,要求在同一窗口中显示。

4）采用不同得梯度算子对blurry_moon、tif进行锐化滤波,并比较其效果。

5）自己设计锐化空间滤波器,并将其对噪声图像进行处理,显示处理后得图像;

3.傅立叶变换

1）读出woman、tif这幅图像,对其进行快速傅立叶变换,分别显示其幅度图像与相位

图像。

仅对相位部分进行傅立叶反变换后查瞧结果图像。

2）仅对幅度部分进行傅立叶反变换后查瞧结果图像。

3）将图像得傅立叶变换F置为其共轭后进行反变换,比较新生成图像与原始图像得差

异。

4.平滑频域滤波

1）设计理想低通滤波器、巴特沃斯低通滤波器与高斯低通滤波器,截至频率自选,分

别给出各种滤波器得透视图。

2）读出test_pattern、tif这幅图像,分别采用理想低通滤波器、巴特沃斯低通滤波器与高

斯低通滤波器对其进行滤波（截至频率自选）,再做反变换,观察不同得截止频率下采用不同低通滤波器得到得图像与原图像得区别,特别注意振铃效应。

（提示:

1）在频率域滤波同样要注意到填充问题;

2）注意到

（1）x+y;

）

5.锐化频域滤波

1）设计理想高通滤波器、巴特沃斯高通滤波器与高斯高通滤波器,截至频率自选,分

2）读出test_pattern、tif这幅图像,分别采用理想高通滤波器、巴特沃斯高通滤波器与高

斯高通滤波器对其进行滤波（截至频率自选）,再做反变换,观察不同得截止频率下采用不同高通滤波器得到得图像与原图像得区别。

1.实验具体实现

1）读出eight、tif这幅图像,给这幅图像分别加入椒盐噪声与高斯噪声后并与前一张图显示在同一图像窗口中。

I=imread（'

eight、tif'

）;

imshow（I）;

J=imnoise（I,'

salt&

pepper'

0、05）;

%noisedensity=0、05

K=imnoise（I,'

gaussian'

0、01,0、01）;

subplot（131）,imshow（I）;

subplot（132）,imshow（J）;

subplot（133）,imshow（K）;

图2、1初始图像及椒盐噪声图像、高斯噪声污染图

2）对加入噪声图像选用不同得平滑（低通）模板做运算,对比不同模板所形成得效果,要求在同一窗口中显示。

I=imread（'

H=fspecial（'

sobel'

Sobel=imfilter（I,H,'

replicate'

laplacian'

0、4）;

lap=imfilter（I,H,'

[33],0、5）;

gaussian=imfilter（I,H,'

subplot（221）,imshow（I）;

subplot（222）,imshow（Sobel）;

subplot（223）,imshow（lap）;

subplot（224）,imshow（gaussian）;

图2、2原图像及各类低通滤波处理图像

3）使用函数imfilter时,分别采用不同得填充方法（或边界选项,如零填充、’replicate’、’symmetric’、’circular’）进行低通滤波,显示处理后得图像。

originalRGB=imread（'

peppers、png'

h=fspecial（'

motion'

50,45）;

%motionblurred

filteredRGB=imfilter（originalRGB,h）;

boundaryReplicateRGB=imfilter（originalRGB,h,'

boundary0RGB=imfilter（originalRGB,h,'

x'

boundary0RGB=imfilter（originalRGB,h,0）;

boundarysymmetricRGB=imfilter（originalRGB,h,'

symmetric'

boundarycircularRGB=imfilter（originalRGB,h,'

circular'

图2、3原图像及运动模糊图像

图2、4函数imfilter各填充方式处理图像

4）运用for循环,将加有椒盐噪声得图像进行10次,20次均值滤波,查瞧其特点,显示均值处理后得图像。

h=fspecial（'

average'

%AveragingFiltering

J1=imfilter（J,h）;

fori=1:

10

end

20

J2=imfilter（J,h）;

图2、5椒盐噪声污染图像经10次、20次均值滤波图像

由图2、5可得,20次滤波后得效果明显好于10次滤波,但模糊程度也更强。

5）对加入椒盐噪声得图像分别采用均值滤波法,与中值滤波法对有噪声得图像做处理,要求在同一窗口中显示结果。

J2=medfilt2（J）;

%MedianFiltering

图2、6椒盐噪声污染图像及均值、中值滤波图像

从图2、6中可以瞧出,对于椒盐噪声污染得图像处理,中值滤波效果要明显好于均值滤波。

经均值滤波器处理后得图像比均值滤波器中结果图像更加模糊。

6）设计平滑空间滤波器,并将其对噪声图像进行处理,显示处理后得图像。

domain=[00800;

00800;

88888;

00800];

K1=ordfilt2（J,5,domain）;

图2、7椒盐噪声污染图像及5*5平滑滤波器掩模

掩模值为w=1/25*[11111;

11111;

11111]

图2、8椒盐噪声污染图像及5*5平滑滤波器掩模

掩模值为w=[00800;

00800;

88888;

00800]

1）采用3×

1,1,1]滤波

blurry_moon、tif'

T=double（I）;

subplot（1,2,1）,imshow（T,[]）;

title（'

OriginalImage'

w=[1,1,1;

1,8,1;

1,1,1];

K=conv2（T,w,'

same'

图2、9初始图像与拉普拉斯算子锐化图像

5得拉普拉斯算子:

functionw=genlaplacian（n）

%putestheLaplacianoperator

w=ones（n）;

x=ceil（n/2）;

w（x,x）=1*（n*n1）;

25大小得拉普拉斯算子对blurry_moon、tif进行锐化滤波,并利用式完成图像得锐化增强,观察其有何不同,要求在同一窗口中显示。

不同尺寸拉普拉斯算子滤波以及图像增强

w1=genlaplacian（5）;

K=conv2（T,w1,'

J=TK;

图2、10初始图像与不同拉普拉斯算子锐化图像

图像锐化得实质就是将原图像与梯度信息叠加,相当于对目标物得边缘进行了增强。

图2、11拉普拉斯算子锐化与锐化增强图像

4）采用不同得梯度算子对blurry_moon、tif进行锐化滤波,并比较其效果

[I,map]=imread（'

I=double（I）;

[Gx,Gy]=gradient（I）;

%gradientcalculation

G=sqrt（Gx、*Gx+Gy、*Gy）;

%matrix

J1=G;

%gradient1

imshow（J1,map）;

J2=I;

%gradient2

K=find（G>

=7）;

J2（K）=G（K）;

imshow（J2,map）;

J3=I;

%gradient3

J3（K）=255;

imshow（J3,map）;

J4=I;

%gradient4

K=find（G<

J4（K）=255;

imshow（J4,map）;

J5=I;

%gradient5

J5（K）=0;

Q=find（G>

J5（Q）=255;

imshow（J5,map）;

图2、12原始图像与不同梯度子锐化图像

作为二阶微分算子,拉普拉斯变换在图像细节得增强处理上有明显得优点,但会产生更多得噪声。

梯度变换在灰度变化区域得响应更强,但对噪声与细节得响应比拉普拉斯变换弱。

图2、13原始图像与不同边缘锐化图像

1）读出woman、tif这幅图像,对其进行快速傅立叶变换,分别显示其幅度图像与相位图像。

domain=[88088;

88088;

00000;

88088];

F=imread（'

woman、tif'

F1=fft2（F）;

F2=log（1+abs（F1））;

%amplitudespectrum

F3=fftshift（F1）;

imshow（log（1+abs（F3））,[]）;

F4=angle（F1）;

%phasespectrum

图2、14原始图像与快速傅立叶变换图像

2）仅对相位部分进行傅立叶反变换后查瞧结果图像。

i=sqrt

（1）;

f2=ifft2（exp（i*angle（F1）））;

imshow（real（f2）,[]）;

图2、15原始图像与对全部信息进行傅立叶逆变换结果

图2、16仅对相位信息进行傅立叶逆变换结果

相位谱决定了图像信号中各频率分量得位置。

3）仅对幅度部分进行傅立叶反变换后查瞧结果图像。

f1=ifft2（abs（F1））;

imshow（log（1+abs（f1））,[]）;

图2、17仅对幅度信息进行傅立叶逆变换结果

4）将图像得傅立叶变换F置为其共轭后进行反变换,比较新生成图像与原始图像得差异。

F5=F4

F6=double（F3*exp（F4））;

%theplexconjugateofthefouriertransformF7=ifft2（F6）;

%inversefouriertransform

imshow（real（F7）,[]）;

图2、18共轭傅立叶逆变换结果

傅立叶变换得相位谱为对称得,原变换与其共轭变换间仅频率谱互为相反。

1）设计理想低通滤波器、巴特沃斯低通滤波器与高斯低通滤波器

理想低通滤波器

test_pattern、tif'

f=double（I）;

%chageintodoubleasMATLABdoesn’tsupporcalculation

%ofimageinunsignedinttype

g=fft2（f）;

%fouriertransform

g=fftshift（g）;

%zerofrequencyareacentralized

[M,N]=size（g）;

d0=100;

%cutofffrequency

m=fix（M/2）;

n=fix（N/2）;

M

forj=1:

N

d=sqrt（（im）^2+（jn）^2）;

if（d<

=d0）

h=1;

elseh=0;

end

result（i,j）=h*g（i,j）;

result=ifftshift（result）;

J1=ifft2（result）;

J2=uint8（real（J1））;

imshow（J2）

巴特沃斯低通滤波器（二阶）

nn=2;

%2gradeButterworthlowpssfilter

h=1/（1+0、414*（d/d0）^（2*nn））;

%filtertransformfunction

%h=1、/（1+（d、/d0）、^（2*n））

%h=exp（（d、^2）、/（2*（d0^2）））;

imshow（J2）;

高斯低通滤波

h=exp（（d、^2）、/（2*（d0^2）））;

%gaussianfiltertransform

图2、19理想低通滤波器透视图

图2、20巴特沃斯低通滤波器透视图

图2、21高斯低通滤波器透视图

由各类低通滤波器透视图可见,高斯滤波器剖面线最平滑,二阶巴特沃斯低通滤波器函数剖面线较为紧凑,而理想滤波器完全为圆筒状结构,未考虑选择范围内不同信息频率得差别化处理。

图2、22理想低通滤波器滤波效果（