上上下下的享受
close all;
clear all;
%加上两个close,clear 最后显⽰图像时候可以正常显⽰,否则图像会因为过⼤⽽被拉伸 img=imread('c.jpg');
img=rgb2gray(img);
劲松六中
figure,imshow(img);
g=fft2(double(img));
陕西理工大学魏乐%进⾏傅⾥叶变换,fft2是⼆维的傅⽴叶变换。
g=fftshift(g);
%中⼼化
imshow(abs(g),[])
[N1,N2]=size(g);
n=2;
d0=30;
%d0是终⽌频率
n1=fix(N1/2);
n2=fix(N2/2);
%n1,n2指中⼼点的坐标,fix()函数是往 0 取整 for i=1:N1
for j=1:N2
轻点CHINA学生CHINESE>中华
菊头蝠d=sqrt((i-n1)^2+(j-n2)^2);
h=1-exp(-d*d/(2*d0*d0));
龙年社区result(i,j)=h*g(i,j);
end
end
result=ifftshift(result);
X2=ifft2(result);
%傅⾥叶去中⼼化以及反变换
final=uint8(real(X2));
%最终的结果归⼀化到0-7之间
figure,imshow(final);