步骤:
(1)计算原图像f(x,y)的DFT,
(2) 讲频谱的零频点移动到频谱图的中心位置;
(3)计算滤波器函数H(U,V)与F(U,V)的乘积G(U,V);
(4)讲频谱G(U,V)的零频点移回到频谱图的坐上角。
(5)计算(4)的结果的傅立叶反变换g(x,y);
(6)取g(x,y)的实部作为最终的滤波后的结果图像。
- clc;
- clear all;
- close all;
- if(ndims==3 && size(I,3)==3)
- I=rgb2gray(I);
- end
- f=fft2(I);
- s=fftshift(f);
- out=s.*ff
- out=ifftshift(out);
- out=ifft2(out);
- out=abs(out);
- out=out/max(out(:));