s095.m

简单的数字图像处理程序

% 经典频域指数低通滤波器  s095
clear;
I95=imread('2.jpg');
I950=rgb2gray(I95);
i2=imnoise(I950,'salt & pepper');
f=double(i2);
g=fft2(f);
g=fftshift(g);
[N1,N2]=size(g);
n=2;
d0=20;
n1=fix(N1/2);
n2=fix(N2/2);
for i=1:N1
for j=1:N2
d=sqrt((i-n1)^2+(j-n2)^2);
h=exp(-(d/d0)^n/0.414);
result(i,j)=h*g(i,j);
end
end
result=ifftshift(result);
x2=ifft2(result);
x3=uint8(real(x2));
subplot(1,2,1),imshow(i2);
subplot(1,2,2),imshow(x3)