1 一維與二維離散傅里葉變換
以周期
對函數 f(t) 采樣可表示為
,
對采樣函數進行傅里葉變換得
,
整理得
。
由於對函數 f(t) 的采樣周期為
,采樣函數的傅里葉變換的一個完整周期為
,
同樣的,
也是采樣函數的傅里葉變換的一個完整周期,只是這個周期不是以原點對稱的。
在
區間中取 M 個點,則第 m 個點的頻率為
,
帶入公式得
,
其中,
為連續函數 f(t) 對應的 M 個離散值,
為取樣函數
的傅里葉變換對應的 M 個離散值,
整理公式得
(由於函數僅在 [0,M-1] 上有非零值,故真實求和區間為 [0,M-1])。
因此,一維離散傅里葉變換對為
,
。
類似的,二維離散傅里葉變換對為
,
。
2 傅里葉變換的性質
1)傅里葉變換平移特性
,
用指數項乘以 f(t) 使得傅里葉變換后原點移動到
處,
使用負指數乘以
使得反傅里葉變換后原點移動到
處,證明如下:
,
使用
替換
得
,
因此有
,類似推導可得
。
將平移特性擴展到二維離散變量上有
。
2)離散傅里葉變換一定具有周期特性,因為離散傅里葉變換的頻率取值在
區間內,有限頻率導致必然具有周期性,
連續傅里葉變換頻率取值為無窮大,所以連續傅里葉變換一般不具有周期性(但也有所有頻率都一樣的函數)。
離散傅里葉變換周期性可表示為
。
觀察公式
或
,
發現頻率取值在
之間,而一個完整的頻率應該在
之間,如下圖:

如果直接應用公式進行傅里葉變換,得到的頻率為 [0,M-1]區間,這是兩個半周期組成的一個周期。
在圖像中則表現為低頻信號分布在4個角落,這顯然不便於觀察頻率信息。
結合傅里葉變換的平移特性,可以將原函數乘以一個正指數項,使得平移后傅里葉變換再 [0,M-1]區間正好是一個完整的周期。
將原函數平移 M/2 可以實現該目標,具體分析如下:
原函數平移 M/2 得
,
由於 x 為非負整數,
,
最終得到
。
對於二維離散變量有相似結論
。
3)原函數(二維及以上)旋轉一定角度,其傅里葉變換也旋轉對應角度。
令
為原函數變量的列向量,
為傅里葉變換函數變量的列向量,對
的傅里葉變換可表示為
,
對
旋轉一定角度可表示為
,其中 R 為旋轉矩陣,
對
的傅里葉變換可表示為
,
由
得
,並將其帶入上式得
,
由於
,
因此
,使得傅里葉變換旋轉相應角度。
4)傅里葉變換具有對稱性,對於二維圖像來說,由於圖像值為實數,其傅里葉變換具有共軛對稱性。
,
由於 f(x,y) 為取值為實數,因此
。
在將傅里葉級數從三角函數轉換為指數函數過程中,通過分析指數函數系數組成部分,可以知道傅里葉級數的共軛對稱性。
而傅里葉變換是傅里葉級數在周期無限大情況下的極限表示,因此,傅里葉變換也應該滿足共軛對稱性。
3 圖像上簡單應用
1)由於離散傅里葉變換的平移與周期特性,在傅里葉變換乘以
可將頻譜中心化。
由於頻譜變化范圍很廣(在不同數量級),為了便於觀察頻譜信息,在可視化之前一般對頻譜進行取對數處理,如
。
2)由於傅里葉變換的共軛對稱性,在求解圖像傅里葉變換時,只需要求解四分之一頻譜信息,其他部分可通過共軛對稱推導而來。
3)圖像頻譜可表示為
,其中,
為傅里葉變換幅度,
為傅里葉變換頻譜。
將傅里葉變換信息拆分為幅度信息與頻譜信息,幅度信息表征了圖像亮度特征,頻譜信息表征了圖像形狀特征。
令
,然后進行反傅里葉變換得到只包含形狀的圖像。
令
,然后進行反傅里葉變換得到只包含亮度的圖像。
4)當圖像前景平移時,根據
可知,其傅里葉變換不發生改變。
5)當圖像前景旋轉時,根據
可知,其傅里葉變換進行相應旋轉。
參考資料 Digital Image Processing Rafael C. Gonzalez & Richard E. Woods
