图像处理频域锐化（理想和巴特沃斯）

图像频域平滑（锐化）：选择一幅图像，例如rice.png，分别使用理想高通滤波器、巴特沃斯高通滤波器、指数高通滤波器和梯形高通滤波器（至少使用两种高通滤波器），显示滤波前后的频域能量分布图，空间图像。

%图像频域锐化    理想  和  巴特沃斯
%导入图像
I=imread('rice.png');
A=I;
AN = A;
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 = double(zeros(m,n));
H = double(zeros(m,n));
for i = 1:m
    for j = 1:n
        d=sqrt((i-round(m/2))^2+(j-round(n/2))^2);
        if d>D0
            H(i,j)=1;
        end
    end
end
AAA=uint8(255.*H);
subplot(334)
imshow(AAA);
title('理想高通滤波器')
 
AIDFT = ADFT.*H;
subplot(335);
imshow(AIDFT);
title('理想滤波后的频域能量分布')
AIDFT=ifftshift(AIDFT);%直流分