图像处理频域平滑（理想和巴特沃斯）

图像频域平滑（去噪）：使用自生成图像（包含白色区域，黑色区域，并且部分区域添加椒盐噪声），然后进行傅里叶变换，并且分别使用理想低通滤波器、巴特沃斯低通滤波器、指数低通滤波器和梯形低通滤波器（至少使用两种低通滤波器），显示滤波前后的频域能量分布图，空间图像。

注意：uint8类型的，会认为是0到255之间    double类型的  会认为是0到1之间！！   复数类型时，显示时一般先取幅值real，再用   uint8（255.*矩阵）！！

%图像频域平滑    理想  和  巴特沃斯
%自制黑白灰度图像
for i=1:300
    for j =1:300
        I(i,j)=0;
    end
end
for i=1:30
    for j =1:50
        I(100+i,100+j)=255;
    end
end
A=I;
%椒盐噪声处理
AN = imnoise(A,'salt & pepper',0.04);
D0 = 20;%截止频率
subplot(331);
imshow(A);
title('原图');
subplot(332);
imshow(AN);
title('椒盐噪声处理');
AN= double(AN);
ADFT  =fft2(AN);%傅里叶变换后的矩阵
ADFT=fftshift(ADFT);%直流分量移到频谱中心
subplot(333);
imshow(ADFT);
title('滤波前频率能量分布');
[m,n] = size(ADFT);
A1