学生实验报告开课学院及实验室：2013年5月24日学院机械与电气工程学院年级、专业、班10通信姓名王敏华学号1007400030实验课程名称多媒体通信技术成绩实验项目名称实验四：频域滤波增强指导老师高星辉一、实验目的图像的边缘和噪声都对应于傅里叶变换中的高频部分，而灰度变换平缓的区域对应傅里叶变换中的低频部分。利用不同的滤波器对图像进行滤波，就可以达到频域增强的效果。通过本实验，理解一幅图像在空域及其傅里叶变换之间的关系，学习构造不同的滤波器，掌握频域滤波实现图像增强的方法。二、实验原理在MATLAB中构造频域滤波器，需先建立一个与欲增强图像相同大小、用于构建频域滤波器的网格数组；利用该网格数组，可求出每点距离原点的欧氏距离，进而构造各种频域滤波器。方法说明：使用函数dftuv(M,N)，产生网格数组。求出距离函数，根据不同滤波器的定义，生成滤波器。频域滤波器的显示需先使用fftshift（）函数移动变换原点。对理想频域滤波器使用ifft2（）函数变换到空域，会产生振铃现象。函数说明：dftuv(M,N)：产生M×N的网格数组。例如：[U,V]=dftuv(500,500);构造理想低通滤波器和巴特沃斯低通滤波器并分别和图像的傅里叶变换相乘，即可实现频域低通滤波。方法说明：理想低通滤波器：H=double(D<=5);截止频率为5的理想低通滤波器巴特沃斯低通滤波器：H=1./(1+(D./5).^(2*1));截止频率为5，阶为1的BLPF以上两式中，D为每点与原点欧氏距离网格。构造理想高通滤波器和巴特沃斯高通滤波器并分别和图像的傅里叶变换相乘，即可实现频域高通滤波。方法说明：理想高通滤波器：H=double(D>5);截止频率为5的理想高通滤波器巴特沃斯高通滤波器：H=1-1./(1+(D./5).^(2*1));截止频率为5，阶为1的BHPF以上两式中，D为每点与原点欧氏距离网格。三、使用仪器、材料计算机，MATLAB6.5四、实验步骤构造频域滤波器使用dftuv（）函数构造一个100×100的网格，并根据其求距离数组D。代码参考如下：[U,V]=dftuv(100,100);D=sqrt(U.^2+V.^2);根据a步骤产生的距离函数，构造截止频率为5的理想低通滤波器H，并参考如下语句，以2维和3维方式显示该滤波器：figure,imshow(fftshift(H))figure,mesh(fftshift(H(1:100,1:100)))分别保存以上两图为1b-1.tif，1b-2.tif；对b步骤中的理想低通滤波器H使用ifft2（）函数进行二维傅里叶反变换，对变换结果使用fftshift（）改变变换原点，再使用real（）函数取实部，显示该结果图像，观察理解理想低通滤波器在空域产生的振铃现象。将图像保存为1c.tif；根据a步骤产生的距离函数，构造截止频率为5、阶为1的巴特沃斯低通滤波器H，并分别以2维和3维方式显示该滤波器，保存两图为1d-1.tif，1d-2.tif；对d步骤中的巴特沃斯低通滤波器H使用ifft2（）函数进行二维傅里叶反变换，对变换结果使用fftshift（）改变变换原点，再使用real（）函数取实部，显示该结果图像，观察理解巴特沃斯低通滤波器的空域图像。将图像保存为1e.tif。低通滤波读入图像test.tif并显示图像。使用fft2（）求其傅里叶变换F，使用如下代码显示其频谱：figure,imshow(log(1+abs(fftshift(F))),[])并使用Export菜单保存为2a.tif；构造一个与图像大小相同、截止频率为5的理想低通滤波器，将其与a步骤中图像傅里叶变换F进行点乘运算，对结果再使用ifft2（）函数进行反变换，使用real（）函数取实部得到滤波后图像并显示，并使用Export菜单保存为2b.tif；如b步骤，分别构造截止频率为20、50的理想低通滤波器并对图像test.tif进行滤波，滤波结果使用Export菜单保存为2c-20.tif，2c-50.tif；分别构造阶为1、截止频率为5、20、50的巴特沃斯低通滤波器并对图像test.tif进行滤波，滤波结果使用Export菜单保存为2d-5.tif，2d-20.tif，2d-50.tif。高通滤波读入图像test.tif并显示图像，使用fft2（）求其傅里叶变换F；分别构造截止频率为5、20、50的理想高通滤波器并对图像test.tif进行滤波，滤波结果使用Export菜单保存为3b-5.tif，3b-20.tif，3b-50.tif；分别构造阶为1、截止频率为5、20、50的巴特沃斯高通滤波器并对图像test.tif进行滤波，滤波结