dft 图像的傅里叶变换。在代码中实现了将两幅图像的幅度谱和相位谱进行交换

图像的傅里叶变换。在代码中实现了将两幅图像的幅度谱和相位谱进行交换，进而理解两者在图像显示中的作用

clear;
i=imread（‘Tulips.jpg‘）;
grayim=rgb2gray（i）;%转为灰度图%
[width height]=size（grayim）;
figure（1）;imshow（i）;title（‘原图‘）;%原图%
figure（2）;imshow（grayim）;title（‘原灰度图‘）;%原图%
f=fft2（grayim）;%dft%
f_shift=log（abs（fftshift（f）））;%把直流分量移到中间%
figure（3）;imshow（f_shift[]）;title（‘DFT结果‘）;colormap（jet（64）） colorbar;%彩色图直流分量移至中间dft结果%
f1=f;f1（log（abs（f1））>11）=0;%f1为f去低频分量%
f2=f;count=0;pixels=width*height;%去高频分量%
for j=1:width;
    for k=1:height;
        if log（abs（f2（jk）））<8
            f2（jk）=0;
            count=count+1;
        end;
    end;
end;
f3=f;f3（11）=0;%去直流分量%

iff=ifft2（f）;%idft%
iff1=ifft2（f1）;%idft%
iff2=ifft2（f2）;%idft%
iff3=ifft2（f3）;%idft%
 iff4=abs（f）;%求幅值%
pyy=angle（f）;%求相位%
p=180/pi*pyy;
figure（4）;imshow（iff[]）;title（‘没有压缩重构的图像‘）;%图像复原%
figure（5）;i

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     595284  2011-04-26 21:39  dft\Hydrangeas.jpg

     文件      72660  2011-05-01 20:01  dft\H傅里叶变换幅度.jpg

     文件     117190  2011-05-01 20:01  dft\H傅里叶变换相位.jpg

     文件      54885  2011-05-01 20:01  dft\H幅度+T相位.jpg

     文件     620888  2009-07-14 12:52  dft\Tulips.jpg

     文件      68860  2011-05-01 20:02  dft\T傅里叶变换幅度.jpg

     文件     107379  2011-05-01 20:02  dft\T傅里叶变换相位.jpg

     文件      51593  2011-05-01 20:01  dft\T幅度+H相位.jpg

     文件      90737  2011-03-20 17:06  dft\傅里叶变换.jpg

     文件       1255  2011-04-26 21:27  dft\傅里叶变换.m

     文件      47600  2011-03-20 17:05  dft\原图灰度图.jpg

     文件      53453  2011-03-20 17:06  dft\去低频轮廓.jpg

     文件      49399  2011-03-20 17:07  dft\去直流分量.jpg

     文件        964  2011-05-01 20:00  dft\幅相变换1.m

     文件      48598  2011-03-20 17:06  dft\无压缩重构.jpg

     文件      49922  2011-03-20 17:07  dft\进行图像压缩.jpg

     目录          0  2011-05-01 20:02  dft

----------- ---------  ---------- -----  ----

              2030667                    17