my_fft2_2.m

两个滤波器的代码

clear all;
close all;
I=imread('lena.tif');
koko=size(I);
intRows=koko(1);
intCols=koko(2);
I2=imnoise(I,'gaussian',0.10,0.01);
I2fft=fft2(fftshift(I2));
suodin=fspecial('gaussian',[intRows intCols],0.1);
fftsuodin=fft2(suodin);
I3fft=I2fft.*fftsuodin;
IM=ifft2(I3fft);
figure,subplot(221),
imshow(I),title('Original','color','r'),
subplot(222),imshow(I2),title('Gaussian noise figure','color','r'),
subplot(223),imshow(log10(I3fft)),title('Filter','color','r'),
subplot(224),imshow(unit8(IM)),title('Filter image','color','r');