（全新）海南大学数字图像处理实验四.docVIP

实验四、频域滤波 实验目的 了解频域滤波的方法； 掌握频域滤波的基本步骤。 二、实验内容 使用二维快速傅立叶变换函数fft2( )及其反变换函数ifft2( )对图象进行变换； 自己编写函数生成各种频域滤波器； 比较各种滤波器的特点。 三、实验步骤 图象的傅立叶变换 对图象 1.bmp 做傅立叶变换。 x=imread(‘1.bmp’); f=fft2(x); imshow(real(f)) %显示变换后的实部图象 figure f1=fftshift(f); imshow(real(f1)) x=imread(1.bmp); f=fft2(x); subplot(1,2,1) imshow(real(f)) %显示变换后的实部图象 title(实部图象); subplot(1,2,2) f1=fftshift(f); imshow(real(f1)) title(中心平移图象); 对图象 cameraman.tif 进行傅立叶变换，分别显示变换后的实部和虚部图象。Help real x=imread(cameraman.tif); f=fft2(x); subplot(1,2,1) imshow(real(f)) %显示变换后的实部图象 title(实部图象); subplot(1,2,2) imshow(imag(f)) title(虚部图象); 思考： 对图象 cameraman.tif 进行傅立叶变换，并显示其幅度谱|F(U,V)|。结果类似下图。Abs() 显示结果命令imshow（uint8（y./256）,[]） x=imread(cameraman.tif); f=fft2(x); ff=abs(f); subplot(1,3,1) imshow(x) subplot(1,3,2) imshow(uint8(ff./256),[]) title(|F(u,v)|); subplot(1,3,3) f1=fftshift(ff); imshow(uint8(f1./256),[]) title(中心平移图象); 频域滤波的步骤 求图象的傅立叶变换得 F=fft2(x) 用函数F=fftshit(F) 进行zhongxin移位 生成一个和F 一样大小的滤波矩阵 H . 用F和H相乘得到 G , G=F.*H 求G的反傅立叶变换得到 g 就是我们经过处理的图象。 这其中的关键就是如何得到H 。 理想低通滤波器 函数dftuv( )在文件夹中，它用生成二维变量空间 如： [U V]=dftuv(11,11) 生成理想低通滤波器 [U V]=dftuv(51,51); D=sqrt(U.^2+V.^2); H=double(D=15); Mesh(U,V,H) 应用以上方法，对图象cameraman.tif进行低通滤波; close all Q=0.7 F=imread(cameraman.tif) [U V]=dftuv(size(F,1),size(F,2)); D=sqrt(U.^2+V.^2); H=double(D= size(F,1)/2*Q); %修改系数Q为0.5,0.3,0.2 FF=fft2(F); G=FF.*H; imshow(real(fftshift(FF))) figure imshow(real(fftshift(G))) g=real(ifft2(G)); figure imshow(uint8(g)) close all Q=0.7 F=imread(cameraman.tif) [U V]=dftuv(size(F,1),size(F,2)); D=sqrt(U.^2+V.^2); H=double(D= size(F,1)/2*Q); %修改系数Q为0.5,0.3,0.2 FF=fft2(F); G=FF.*H; subplot(1,3,1) imshow(real(fftshift(FF))) subplot(1,3,2) imshow(real(fftshift(G))) g=real(ifft2(G)); subplot(1,3,3) imshow(uint8(g)) 思考： 观察理想低通滤波器不同滤波半径对滤波结果的影响。即系数Q分别为0.7、0.5、0.3、0.2时，结果如何？并加以解