实验四图像的傅立叶变换与频域滤波

实验四图像的傅立叶变换与频域滤波

ID:47644552

大小:451.40 KB

页数:9页

时间:2019-08-26

实验四图像的傅立叶变换与频域滤波_第1页
实验四图像的傅立叶变换与频域滤波_第2页
实验四图像的傅立叶变换与频域滤波_第3页
实验四图像的傅立叶变换与频域滤波_第4页
实验四图像的傅立叶变换与频域滤波_第5页
资源描述:

《实验四图像的傅立叶变换与频域滤波》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库

1、实验四图像的傅立叶变换与频域滤波一、实验目的1了解图像变换的意义和手段;2熟悉傅里叶变换的基本性质;3熟练掌握FFT方法的应用;4通过实验了解二维频谱的分布特点;5通过本实验掌握利用MATLAB编程实现数字图像的傅立叶变换。6、掌握怎样利用傅立叶变换进行频域滤波7、掌握频域滤波的概念及方法8、熟练掌握频域空间的各类滤波器9、利用MATLAB程序进行频域滤波二、实验原理1应用傅立叶变换进行图像处理傅里叶变换是线性系统分析的一个有力工具,它能够定量地分析诸如数字化系统、采样点、电子放大器、卷积滤波器、噪音和显示点等的作用。

2、通过实验培养这项技能,将有助于解决大多数图像处理问题。对任何想在工作屮有效应用数字图像处理技术的人来说,把吋间用在学习和掌握博里叶变换上是很有必要的。2傅立叶(Fourier)变换的定义对于二维信号,二维Fourier变换定义为:F{f(x,y)]=F@,v)=J£f(x9y)e~i2m}dxdy二维离散傅立叶变换为:M-N-.r=0y=0图像的傅立叶变换与一维信号的傅立叶变换变换一样,有快速算法,具体参见参考书目,有关傅立叶变换的快速算法的程序不难找到。实际上,现在有实现傅立叶变换的芯片,可以实时实现傅立叶变换。

3、3利用MATLAB软件实现数字图像傅立叶变换的程序:l=imread(z原图像名.gif);%读入原图像文件imshow(l);%显示原图像fftl=fft2(l);%二维离散傅立叶变换sfftl=fftshift(fftl);%直流分量移到频谱中心RR=real(sfftl);%取傅立叶变换的实部ll=imag(sfftl);%取傅立叶变换的虚部A=sqrt(RR.A2+ll.A2);%计算频谱幅值A=(A-min(min(A)J)/(max(max(A))-min(min(A)))*225;%归一化figure;%

4、设定窗口imshow(A);%显示原图像的频谱域滤波分为低通滤波和高通滤波两类,对应的滤波器分别为低通滤波器和高通滤波器。频域低通过滤的基本思想:G(u,v)=F(u,v)H(uzv)F(u,v)是需要钝化图像的傅立叶变换形式,H(u,v)是选取的一个低通过滤器变换函数,G(u,v)是通过H(u,v)减少F(u”v)的高频部分来得到的结果,运用傅立叶逆变换得到钝化后的图像。理想地通滤波器(ILPF)具有传递函数:z/D(w,v)Do其中,为指定的非负数,D(S)为(u,v)到滤波器的中心的距离

5、。£>(w,v)=D0的点的轨迹为一个圆。n阶巴特沃兹低通滤波器(BLPF)(在距离原点D.处出现截至频率)的传递函数为H(u,v)=1l+[D(S)/D.]2lt与理想地通滤波器不同的是,巴特沃兹率通滤波器的传递函数并不是在Q处突然不连续。高斯低通滤波器(GLPF)的传递函数为H(u,v)=eD2^v)/2a2其中,o•为标准差。相应的高通滤波器也包括:理想高通滤波器、n阶巴特沃兹高通滤波器、高斯高通滤波器。给定一个低通滤波器的传递函数H沁⑴,通过使用如下的简单关系,可以获得相应高通滤波器的传递函数:Hhp=-Hl

6、pM利用MATLAB实现频域滤波的程序f=imread(,room.tif);F=fft2(f);%对图像进行傅立叶变换S=fftshift(log(l+abs(F)));%对变换后图像进行队数变化,并对其坐标平移,使其中心化S=gscale(S);%将频谱图像标度在0-256的范围内imshow(S)%显示频谱图像h=special(,sobe「);%产生空间'sobeP模版freqz2(h)%查看相应频域滤波器的图像PQ=paddedsize(size(f));%产生滤波时所需大小的矩阵H=freqz2(h,PQ(

7、l),PQ(2));%产生频域中的'sobel'滤波器Hl=ifftshift(H);%重排数据序列,使得原点位于频率矩阵的左上角imshow(abs(H),[])%以图形形式显示滤波器figure,imshow(abs(Hl),[])gs=imfilter(double(f)?h);%用模版h进彳亍空域滤波gf=dftfilt(f,Hl);figure,imshow(gs,[])%用滤波器对图像进行频域滤波figure,imshow(gf,[])figure,imshow(abs(gs),[])figure,imsh

8、ow(abs(gf),[])f=imread(,number.tif');%读取图片PQ=paddedsize(size(f));%产生滤波时所需大小的矩阵D0=0.05*PQ(l);%设定高斯高通滤波器的阈值H=hpfilter(,gaussian,,PQ(l),PQ(2),D0);%产生高斯高通滤波器%对图像进行滤波figur

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。