MATLAB图片处理(五.完)
图像的频域滤波
f=imread('b.jpg');
imshow(f)
g=rgb2gray(f);
figure,imshow(g)
g=imresize(g,[600 600]);
p=imnoise(g,'gaussian',0.01);
figure,imshow(p)
F=fft2(p);
fc=fftshift(F);
figure,imshow(abs(fc),[])
s2=log(1+ abs(fc));
figure,imshow(s2,[])
[M,N]=size(p)
sig=30
h=lpfilter('gaussian',M,N,sig);
G=h.*F;
b=real(ifft2(G));
figure,imshow(b,[])
sig=40
h=lpfilter('gaussian',M,N,sig);
G=h.*F;
b=real(ifft2(G));
figure,imshow(b,[])