clear all;
clc
am=imread('1.tiff');
a=double(am);
c=size(a);
n=c(1)
D0=input('Enter cutoff of frequency');
for u=1:1:c(1)
for v=1:1:c(2)
D=((u-(n/2))^2+(v-(n/2))^2)^0.5;
if D<D0;
H(u,v)=1;
else H(u,v)=0;
end
end
end
b=fft2(a);
d=fftshift(b);
x=d.*H;
p=abs(ifft2(x));
figure(1),imshow(uint8(am))
%figure(2),mesh(H)
figure(2), imshow(uint8(p))
%figure(4), imagesc(H),colormap(gray)
评论0