图像频域平滑(去噪):使用自生成图像(包含白色区域,黑色区域,并且部分区域添加椒盐噪声),然后进行傅里叶变换,并且分别使用理想低通滤波器、巴特沃斯低通滤波器、指数低通滤波器和梯形低通滤波器(至少使用两种低通滤波器),显示滤波前后的频域能量分布图,空间图像。
注意: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