s099.m

来自「简单的数字图像处理程序」· M 代码 · 共 28 行

% 经典频域带通滤波器  s099
clear;
I99=imread('2.jpg');
I990=rgb2gray(I99);
f=im2double(I990);
g=fft2(f);
g=fftshift(g);
[N1,N2]=size(g);
n=2;
d1=5;
d2=40;
n1=fix(N1/2);
n2=fix(N2/2);
for i=1:N1
  for j=1:N2
     d=sqrt((i-n1)^2+(j-n2)^2);
     if d1<d<d2
         h(i,j)=1;
     else h(i,j)=0;
     end
     result(i,j)=h(i,j)*g(i,j);
  end
end
result=ifftshift(result);
x2=ifft2(result);
x3=im2uint8(real(x2));
subplot(1,2,1),imshow(I990);
subplot(1,2,2),imshow(x3)