文档介绍:.
实验五二维DFT和频域滤波器实验报告
1. 一、实验目的理解二维离散傅立叶变换的概念,掌握Matlab屮进行二维离散傅立叶变换并显示频谱图的方法。
2. 理解频域滤波的基本概念,掌握频域滤波的基本方法。
3. 理解理想低通: .
实验五二维DFT和频域滤波器实验报告
1. 一、实验目的理解二维离散傅立叶变换的概念,掌握Matlab屮进行二维离散傅立叶变换并显示频谱图的方法。
2. 理解频域滤波的基本概念,掌握频域滤波的基本方法。
3. 理解理想低通滤波器、巴特沃斯低通滤波器的概念,并会利用其传递函数对已知图像进行滤波。
1. 二、实验要求在Matlab中对给定图像进行二维离散傅立叶变换,并显示其居中的频谱图。
2. 分别生成尺寸为MXN(给定图像的尺寸),截止频率DO二50的理想低通滤波器和1阶巴特沃斯低通滤波器。并以图像的形式显示。
3. 对给定图像分别利用上述滤波器进行滤波,显示滤波后图像的频谱以及经滤波后的图像,并与原始图像及其频谱进行对比。
三、实验原理根据卷积定理,两个空间函数的卷积可以通过计算两个傅立叶变换函数的乘积的逆变换得到,如果f(xy)和h(x,y)分别代表图像与空间滤波器,F(u,v)和H(u,v)分别为响应的傅立叶变换(H(u,v)又称为传递函数),那么我们可以利用卷积定理来进行频域滤波。
理想低通滤波器具有传递函数1若D(u,v)DoH(u,v)0若D(u,v)Do其屮D。为制定的非负数,D(u,v)为点(u,v)到滤波器中心的距离。n阶巴特沃斯低通滤波器的传递函数为
H(U,V)11「D(ii・v)/D°l2n四、实验结果图像
理弊假询龙巧器茁舖谱图匚贞询理匸恨诵总涯后的频谓
五、实验代码%原始图片及其居中后的频谱:
J=imread('d:\');J_gray=rgb2gray(J);figure;imshow(J_gray);%显示原始灰度图像titlef原始灰度图像,);J2=fft2(J_gray);%二维离散傅立叶变换J3二fftsh辻t(J2);%居中J4=log(abs(J3));%对图像做二维DFT变换并居中figure;imshow(J4,[]);%显不居中后的频谱图titlef离散傅立叶变换居中频谱图,);%理想低通滤波器:
[m,n]=size(J_gray);dO二50;x0=round(256);yO=round(256);%定义一个半径为256的圆形区域%截止频率dono的理想低通滤波器fori=l:mforj=l:nd=sqrt((i-xO)A2+(j-yO)A2);%到原点中心的距离if(d<dO)h(i,j)=l;elseh(i,j)=0;figure;imshow(h);%画出理想低通滤波器的频谱图title('理想低通滤波器的频谱图’);J5二J3.*h;J6=log(abs(J5));figure;imshow(J6,