实验三 数字图像的空间域滤波.docx
- 文档编号:7344267
- 上传时间:2023-01-23
- 格式:DOCX
- 页数:27
- 大小:1.43MB
实验三 数字图像的空间域滤波.docx
《实验三 数字图像的空间域滤波.docx》由会员分享,可在线阅读,更多相关《实验三 数字图像的空间域滤波.docx(27页珍藏版)》请在冰豆网上搜索。
实验三数字图像的空间域滤波
实验三、四数字图像得空间域滤波与频域滤波
1.实验目得
1.掌握图像滤波得基本定义及目得。
2.理解空间域滤波得基本原理及方法。
3.掌握进行图像得空域滤波得方法。
4.掌握傅立叶变换及逆变换得基本原理方法。
5.理解频域滤波得基本原理及方法。
6.掌握进行图像得频域滤波得方法。
2.实验基本原理
1.空间域增强
空间域滤波就是在图像空间中借助模板对图像进行领域操作,处理图像每一个像素得取值都就是根据模板对输入像素相应领域内得像素值进行计算得到得。
空域滤波基本上就是让图像在频域空间内某个范围得分量受到抑制,同时保证其她分量不变,从而改变输出图像得频率分布,达到增强图像得目得。
空域滤波一般分为线性滤波与非线性滤波两类。
线性滤波器得设计常基于对傅立叶变换
得分析,非线性空域滤波器则一般直接对领域进行操作。
各种空域滤波器根据功能主要分为
平滑滤波器与锐化滤波器。
平滑可用低通来实现,平滑得目得可分为两类:
一类就是模糊,目
得就是在提取较大得目标前去除太小得细节或将目标内得小肩端连接起来;另一类就是消除噪
声。
锐化可用高通滤波来实现,锐化得目得就是为了增强被模糊得细节。
结合这两种分类方法,
可将空间滤波增强分为四类:
线性平滑滤波器(低通)
非线性平滑滤波器(低通)
线性锐化滤波器(高通)
非线性锐化滤波器(高通)
空间滤波器都就是基于模板卷积,其主要工作步骤就是:
1)将模板在图中移动,并将模板中心与图中某个像素位置重合;
2)将模板上得系数与模板下对应得像素相乘;
3)将所有乘积相加;
4)将与(模板得输出响应)赋给图中对应模板中心位置得像素。
2.平滑滤波器
1)线性平滑滤波器
线性低通平滑滤波器也称为均值滤波器,这种滤波器得所有系数都就是正数,对3×3
得模板来说,最简单得就是取所有系数为1,为了保持输出图像任然在原来图像得灰度值范围
内,模板与象素邻域得乘积都要除以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.平滑空间滤波:
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('eight、tif');
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,'replicate');
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,'replicate');
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次均值滤波,查瞧其特点,显示均值处理后得图像。
J=imnoise(I,'salt&pepper',0、05);
h=fspecial('average');%AveragingFiltering
J1=imfilter(J,h);
fori=1:
10
J1=imfilter(J,h);
end
fori=1:
20
J2=imfilter(J,h);
图2、5椒盐噪声污染图像经10次、20次均值滤波图像
由图2、5可得,20次滤波后得效果明显好于10次滤波,但模糊程度也更强。
5)对加入椒盐噪声得图像分别采用均值滤波法,与中值滤波法对有噪声得图像做处理,要求在同一窗口中显示结果。
h=fspecial('average');%AveragingFiltering
J1=imfilter(J,h);
J2=medfilt2(J);%MedianFiltering
图2、6椒盐噪声污染图像及均值、中值滤波图像
从图2、6中可以瞧出,对于椒盐噪声污染得图像处理,中值滤波效果要明显好于均值滤波。
经均值滤波器处理后得图像比均值滤波器中结果图像更加模糊。
6)设计平滑空间滤波器,并将其对噪声图像进行处理,显示处理后得图像。
domain=[00800;
00800;
88888;
00800;
00800];
K1=ordfilt2(J,5,domain);
图2、7椒盐噪声污染图像及5*5平滑滤波器掩模
掩模值为w=1/25*[11111;11111;11111;11111;11111]
图2、8椒盐噪声污染图像及5*5平滑滤波器掩模
掩模值为w=[00800;00800;88888;00800;00800]
2.锐化空间滤波
1)采用3×3得拉普拉斯算子w=[1,1,1;1–81;1,1,1]滤波
I=imread('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初始图像与拉普拉斯算子锐化图像
2)编写函数w=genlaplacian(n),自动产生任一奇数尺寸n得拉普拉斯算子,如5×5得拉普拉斯算子:
w=[11111
11111
112411
11111
11111]
functionw=genlaplacian(n)
%putestheLaplacianoperator
w=ones(n);
x=ceil(n/2);
w(x,x)=1*(n*n1);
3)分别采用5×5,9×9,15×15与25×25大小得拉普拉斯算子对blurry_moon、tif进行锐化滤波,并利用式完成图像得锐化增强,观察其有何不同,要求在同一窗口中显示。
不同尺寸拉普拉斯算子滤波以及图像增强
w1=genlaplacian(5);
I=imread('blurry_moon、tif');
T=double(I);
K=conv2(T,w1,'same');
J=TK;
图2、10初始图像与不同拉普拉斯算子锐化图像
图像锐化得实质就是将原图像与梯度信息叠加,相当于对目标物得边缘进行了增强。
图2、11拉普拉斯算子锐化与锐化增强图像
4)采用不同得梯度算子对blurry_moon、tif进行锐化滤波,并比较其效果
[I,map]=imread('blurry_moon、tif');
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
K=find(G>=7);
J3(K)=255;
imshow(J3,map);
J4=I;%gradient4
K=find(G<=7);
J4(K)=255;
imshow(J4,map);
J5=I;%gradient5
K=find(G<=7);
J5(K)=0;
Q=find(G>=7);
J5(Q)=255;
imshow(J5,map);
图2、12原始图像与不同梯度子锐化图像
作为二阶微分算子,拉普拉斯变换在图像细节得增强处理上有明显得优点,但会产生更多得噪声。
梯度变换在灰度变化区域得响应更强,但对噪声与细节得响应比拉普拉斯变换弱。
5)自己设计锐化空间滤波器,并将其对噪声图像进行处理,显示处理后得图像;
图2、13原始图像与不同边缘锐化图像
3.傅立叶变换
1)读出woman、tif这幅图像,对其进行快速傅立叶变换,分别显示其幅度图像与相位图像。
domain=[88088;
88088;
00000;
88088;
88088];
K1=ordfilt2(J,5,domain);
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)仅对相位部分进行傅立叶反变换后查瞧结果图像。
F1=fft2(F);
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置为其共轭后进行反变换,比较新生成图像与原始图像得差异。
F1=fft2(F);
F2=log(1+abs(F1));%amplitudespectrum
F3=fftshift(F1);
F4=angle(F1);%phasespectrum
F5=F4
F6=double(F3*exp(F4));%theplexconjugateofthefouriertransformF7=ifft2(F6);%inversefouriertransform
imshow(real(F7),[]);
图2、18共轭傅立叶逆变换结果
傅立叶变换得相位谱为对称得,原变换与其共轭变换间仅频率谱互为相反。
4.平滑频域滤波
1)设计理想低通滤波器、巴特沃斯低通滤波器与高斯低通滤波器
理想低通滤波器
I=imread('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);
fori=1:
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);
end
end
result=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);
nn=2;%2gradeButterworthlowpssfilter
d0=100;
m=fix(M/2);n=fix(N/2);
fori=1:
M
forj=1:
N
d=sqrt((im)^2+(jn)^2);
h=1/(1+0、414*(d/d0)^(2*nn));%filtertransformfunction
%h=1、/(1+(d、/d0)、^(2*n))
%h=exp((d、^2)、/(2*(d0^2)));
result(i,j)=h*g(i,j);
end
end
result=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:
M
forj=1:
N
d=sqrt((im)^2+(jn)^2);
h=exp((d、^2)、/(2*(d0^2)));%gaussianfiltertransform
result(i,j)=h*g(i,j);
end
end
result=ifftshift(result);
J1=ifft2(result);
J2=uint8(real(J1));
图2、19理想低通滤波器透视图
图2、20巴特沃斯低通滤波器透视图
图2、21高斯低通滤波器透视图
由各类低通滤波器透视图可见,高斯滤波器剖面线最平滑,二阶巴特沃斯低通滤波器函数剖面线较为紧凑,而理想滤波器完全为圆筒状结构,未考虑选择范围内不同信息频率得差别化处理。
图2、22理想低通滤波器滤波效果(
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验三 数字图像的空间域滤波 实验 数字图像 空间 滤波
![提示](https://static.bdocx.com/images/bang_tan.gif)