锐化.m

来自「在matlab中最简单的功能实现去噪、平滑和锐化」· M 代码 · 共 17 行

a=imread('F:\学习\数字图象\作业\第一次上机\去噪之后的图.jpg');             %读取
b=double(a);                                      %数据类型转换便于傅立叶变换   
F=fftshift(fft2(b));                              %傅立叶变换
A=sqrt(real(F).^2+imag(F).^2);                    %求频谱    
A=(A-min(min(A)))/(max(max(A))-min(min(A)))*255;  %频谱归一化
F=F*(1/(1+(100/sqrt(i^2+j^2))^50));
A=sqrt(real(F).^2+imag(F).^2);
A=(A-min(min(A)))/(max(max(A))-min(min(A)))*255; %频谱处理
figure;                                          %打开新窗口
imshow(A);                                      %显示处理后频谱

f=ifft2(ifftshift(F));                           %反变换
c=uint8(f);                                      %数据类型还原
figure;                                          %打开新窗口
imshow(c);                                       %去噪之后的图