数字图像处理，相位 频谱

本资料主要对lena图像进行处理，用MATLAB处理相位频谱，并给出结果。

clear all;
close all;
I=imread（‘C:\Documents and Settings\Administrator\桌面\lena.bmp‘）;
[heightwidth]=size（I）;
figureimshow（I）title（‘原始图像‘）;

fftI=fft2（double（I））;
sfftI=fftshift（fftI）;
A=abs（sfftI）;
E=（A-min（min（A）））/（max（max（A））-min（min（A）））*255 ;
figureimshow（E）title（‘图像幅频‘）;
B=angle（sfftI）;
F=（B-min（min（B）））/（max（max（B））-min（min（B）））*255;
figureimshow（F）title（‘图像相频‘）;

C=abs（ifft2（A））;
C=（C-min（min（C）））/（max（max（C））-min（min（C）））*255 ;
figureimshow（C）title（‘幅频反变换的图像‘）;
for i=1:height                              %B的范围由（-pipi）调整至（02pi）
     for j=1:width
         if B（ij）<0
             B（ij）=B（ij）+2*pi;
         end
    end
end                   
D=ifft2（ifftshift（exp（-j*B）））;

%D=ifft2（exp（-j*B））;
D=（D-min（min（D）））/（max（max（D））-min（min（D）））*255;

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

     文件     231322  2010-10-26 08:34  图像处理\lena.bmp

     文件        878  2010-11-19 20:32  图像处理\第一题\fourier.m

     文件     366258  2010-11-19 20:31  图像处理\第一题\图像幅频.bmp

     文件     366258  2010-11-19 20:31  图像处理\第一题\图像相频.bmp

     文件     366258  2010-11-19 20:30  图像处理\第一题\幅频反变换的图像.bmp

     文件     366258  2010-11-19 20:32  图像处理\第一题\相频反变换的图像.bmp

     文件     366258  2010-11-19 21:16  图像处理\第二题\GLPF变换后的图像.bmp

     文件       1137  2010-11-19 21:15  图像处理\第二题\GPLF.m

     目录          0  2010-12-03 09:26  图像处理\第一题

     目录          0  2010-12-03 09:26  图像处理\第二题

     目录          0  2010-12-03 09:26  图像处理

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

              2064627                    11