orgImage=imread('tire.tif');
figure(1); imshow(orgImage);
fftI1=fft2(orgImage); %二维离散傅立叶变换
sfftI1=fftshift(fftI1); %直流分量移到频谱中心
RR1=real(sfftI1); %取傅立叶变换的实部
II1=imag(sfftI1); %取傅立叶变换的虚部
A1=sqrt(RR1.^2+II1.^2); %计算频谱幅值
a1=(A1-min(min(A1)))/(max(max(A1))-min(min(A1)))*225;%归一化
figure(2); imshow(a1); %显示原图像的频谱
image=ifftshift(a1);
image3=real(image);
figure(3); imshow(image3);
我按你的意思改了程序,可还是得不到原图啊???
追答这个我忽略了一点,你把最后一行的那个figure, imshow(g),改成figure,imshow(g, [ ])就行了。