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 img woman')
subplot(2,4,7),imshow(phase_woman,),title('女人的相位')
ifft_womon_phase=ifft2(exp(1j.*phase_woman));
subplot(2,4,5),imshow(log(real(ifft_womon_phase))+1,);title('只用女人的相位重建的影象')
ifft_img=ifftshift(ifft2(woman_fre));
subplot(2,4,6),imshow(log(real(ifft_img))+1,),title('只用女人的振幅重建的影象')
rec_fre_woman_phase=ifft2(rec_fre.*exp(1j*phase_woman));
rec_phase_woman_fre=ifft2(woman_fre.*exp(1j*phase_rec));
subplot(2,4,8),imshow(log(real(rec_phase_woman_fre))+1,);title('用矩形相位和女人的振幅重建的影象')
figure,imshow(real(rec_fre_woman_phase),);title('用woman相位和矩形振幅重建的影象')
figure,imshow(ori_rec,);title('矩形原圖')
ps： 物體的相位對成像有著關鍵性的影響，相位決定了物體的形狀，換句話說相位資訊決定了這個是什麼物體，尤其是在晶體學解析中相位資訊非常重要，但是相位現在無法測量，只能測量振幅，不過現在有一種相位迭代恢復演算法，可以恢復相位。迭代演算法非常奇妙，如在分形學中簡單的公式加上限制條件，就可以產生非常奇妙的影象·····這個世界是如此的簡單（簡單的公式描述），又是如此的複雜（反覆的迭代），令人驚嘆不已。
傅利葉變換 2 傅利葉幅度譜與相位譜
快速傅利葉變換 void fft2image inputarray src,outputarray dst 求傅利葉變換的幅度譜 amplitude spectrum void amplitudespectrum inputarray srcfft,outputarray dstspectrum 傅...
傅利葉級數（一）
概述 傅利葉級數，起初是用來解決如何求解微分方程的問題的。學過高數的應該都知道，對於非齊次微分方程的求解，如果輸入是指數冪，sin或者cos系列，就可以先求其其次微分方程的通解，然後通過公式求出特解。如果是其他的函式，就很難解決。傅利葉級數就是將所有的週期函式轉化為sin和cos的乘積。這篇部落格就...
傅利葉級數與傅利葉變換
法國數學家傅利葉發現，任何週期函式都可以用正弦函式和余弦函式構成的無窮級數來表示 選擇正弦函式與余弦函式作為基函式是因為它們是正交的 後世稱傅利葉級數為一種特殊的三角級數，根據尤拉公式，三角函式又能化成指數形式，也稱傅利葉級數為一種指數級數。傅利葉變換，表示能將滿足一定條件的某個函式表示成三角函式 ...