频域滤波（一） 傅里叶谱与相位

ori_rec=imread('C:\Program Files\MATLAB\R2013a\bin\Original_Images\DIP3E_Original_Images_CH04\Fig0424(a)(rectangle).tif');
ori_rec=mat2gray(ori_rec(257:768,257:768)); % 归一化 转化为double 
fft_rec=fftshift(fft2(ori_rec)); % 矩形的傅里叶变换2D
phase_rec=angle(fft_rec);% 矩形的相位
rec_fre=abs(fft_rec);% 矩形的振幅
subplot(2,4,1),imshow(phase_rec,[]);title('ori rec phase 矩形的相位');
% figure,imshow(abs(fft_rec),[]),title('frequency img no log')
subplot(2,4,2),imshow(log(rec_fre+1),[]),title('frequency img  log 矩形的频谱') 
ori_woman=imread('C:\Program Files\MATLAB\R2013a\bin\Original_Images\DIP3E_Original_Images_CH04\woman.tif');
ori_woman=mat2gray(ori_woman);
subplot(2,4,3),imshow(ori_woman,[]);title('女人原图')
fft2_woman=fftshift(fft2(ori_woman));
phase_woman=angle(fft2_woman);  % 女人的相位
woman_fre=abs(fft2_woman);      %女人的振幅
subplot(2,4,4),imshow(log(woman_fre)+1,[]),title('frequency im