下边为三个滤波器的使用,截至频率半径值分别设置30,90,160,,先对灰度图像进行加高斯噪声,然后进行滤波处理
F=fftshift(F);%将变换的原点移到频率矩形的中心
[M,N]=size(f);%截至半径30理想低通滤波
h1=zeros(M,N);
for i=1:M
for j=i:N
if(sqrt(((i-M/2)^2+(j-N/2)^2))<30)
h1(i,j)=1;
end
end
end
G1=F.*h1;
G1=ifftshift(G1);
g1=real(ifft2(G1)); %截至半径90理想低通滤波
h2=zeros(M,N);
for i=1:M
for j=i:N
if(sqrt(((i-M/2)^2+(j-N/2)^2))<90)
h2(i,j)=1;