13-11.m

matlab6.5图形图像处理源程序.很多的。很好用的。很有帮助的。

I = imread('saturn.tif');
J = imnoise(I,'salt & pepper',0.02);
%图像添加盐椒噪声
subplot(121),imshow(J)
title('含有噪声的原图像')
J=double(J);
f=fft2(J);
g=fftshift(f);
[M,N]=size(f);
n=3;d0=20;
n1=floor(M/2);n2=floor(N/2);
for i=1:M
    for j=1:N
        d=sqrt((i-n1)^2+(j-n2)^2);
        h=1/(1+0.414*(d/d0)^(2*n));
        g(i,j)=h*g(i,j);
    end
end
g=ifftshift(g);
g=uint8(real(ifft2(g)));
subplot(122),imshow(g)
title('三阶Butterworth滤波图像')