二維傅立葉變換

傅立葉轉換（英语：Fourier transform）是一種幫助我們分析訊號頻域成分的積分變換，詳細內容詳見傅立葉轉換一文。一般教科書所教的通常是一維的傅立葉轉換，但我們也可以將傅立葉轉換推廣到多維的空間。而二維傅立葉轉換即是由一維傅立葉轉換推廣而來，近幾十年來常被運用在影像處理上。其他相關的數學工具，例如二維餘弦轉換、二維濾波器……等等，均是建立在二維傅立葉轉換的概念上而得到的。 考慮一個信號二維的信號， s ( t 1 , t 2 ) {\displaystyle s(t_{1},t_{2})} ，其中 t 1 , t 2 {\displaystyle t_{1},t_{2}} 為兩個獨立變數，且 s ( t 1 , t 2 ) {\displaystyle s(t_{1},t_{2})} 滿足下列方程式: s ( t 1 , t 2 ) = s ( t 1 + k T 1 , t 2 + l T 2 ) , {\displaystyle s(t_{1},t_{2})=s(t_{1}+kT_{1},t_{2}+lT_{2}),} for all t 1 , t 2 {\displaystyle t_{1},t_{2}} , 其中k,l為整數 也就是說 s ( t 1 , t 2 ) {\displaystyle s(t_{1},t_{2})} 在 t 1 − t 2 {\displaystyle t_{1}-t_{2}} 平面為週期函數，在 t 1 {\displaystyle t_{1}} 方向的週期為 T 1 {\displaystyle T_{1}} ，在 t 2 {\displaystyle t_{2}} 方向的週期為 T 2 {\displaystyle T_{2}} ，而信號的傅立葉級數為 s ( t 1 , t 2 ) = ∑ n = − ∞ ∞ ∑ m = − ∞ ∞ a m n e j ( m ω 1 t 1 + n ω 2 t 2 ) ,   ω 1 = 2 π / T 1 , ω 2 = 2 π / T 2 {\displaystyle s(t_{1},t_{2})=\sum _{n=-\infty }^{\infty }\sum _{m=-\infty }^{\infty }a_{mn}e^{j(m\omega _{1}t_{1}+n\omega _{2}t_{2})},\ \omega _{1}=2\pi /T_{1},\omega 2=2\pi /T_{2}} 而 a m n {\displaystyle a_{mn}} 可藉由積分求得 a m n = 1 T 1 T 2 ∫ 0 T 2 ∫ 0 T 1 s ( t 1 , t 2 ) e − j m ω 1 t 1 e − j n ω 2 t 2 d t 1 d t 2 {\displaystyle a_{mn}={\frac {1}{T_{1}T_{2}}}\int _{0}^{T_{2}}\int _{0}^{T_{1}}s(t_{1},t_{2})e^{-jm\omega _{1}t_{1}}e^{-jn\omega _{2}t_{2}}\,dt_{1}\,dt_{2}} 對於平面上的連續週期信號，我們可以使用二維傅立葉級數來分析，但對於平面上的連續非週期信號 s ( t 1 , t 2 ) {\displaystyle s(t_{1},t_{2})} ，我們則需使用二維連續傅立葉變換。 二維連續傅立葉變換定義為 S ( j ω 1 , j ω 2 ) = ∫ − ∞ ∞ ∫ − ∞ ∞ s ( t 1 , t 2 ) e − j ( ω 1 t 1 + ω 2 t 2 ) d t 1 d t 2 {\displaystyle S(j\omega _{1},j\omega _{2})=\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }s(t_{1},t_{2})e^{-j(\omega _{1}t_{1}+\omega _{2}t_{2})}\,dt_{1}\,dt_{2}} 二維連續傅立葉逆變換定義為 s ( t 1 , t 2 ) = 1 4 π 2 ∫ − ∞ ∞ ∫ − ∞ ∞ S ( j ω 1 , j ω 2 ) e j ( ω 1 t 1 + ω 2 t 2 ) d ω 1 d ω 2 {\displaystyle s(t_{1},t_{2})={\frac {1}{4\pi ^{2}}}\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }S(j\omega _{1},j\omega _{2})e^{j(\omega _{1}t_{1}+\omega _{2}t_{2})}\,d\omega _{1}\,d\omega _{2}} 為了方便，我們會將一些變數以向量的形式表示 t → = ( t 1 , t 2 ) t , ω → = ( ω 1 , ω 2 ) t , ω → t t → = ω 1 t 1 + ω 2 t 2 {\displaystyle {\overrightarrow {t}}=(t_{1},t_{2})^{t},{\overrightarrow {\omega }}=(\omega _{1},\omega _{2})^{t},{\overrightarrow {\omega }}^{t}{\overrightarrow {t}}=\omega _{1}t_{1}+\omega _{2}t_{2}} 則上述兩式則可表示為 S ( j ω → ) = ∫ − ∞ ∞ ∫ − ∞ ∞ s ( t → ) e j ω → t t → d t → ,     {\displaystyle S(j{\overrightarrow {\omega }})=\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }s({\overrightarrow {t}})e^{j{\overrightarrow {\omega }}^{t}{\overrightarrow {t}}}\,d{\overrightarrow {t}},\ \ } s ( t → ) = 1 4 π 2 ∫ − ∞ ∞ ∫ − ∞ ∞ S ( j ω → ) e j ω → t t → d ω → {\displaystyle s({\overrightarrow {t}})={\frac {1}{4\pi ^{2}}}\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }S(j{\overrightarrow {\omega }})e^{j{\overrightarrow {\omega }}^{t}{\overrightarrow {t}}}\,d{\overrightarrow {\omega }}} 一般在作影像處理的影像大多不是連續信號，而對於平面上的不連續信號，我們則需使用二維離散傅立葉變換。 假設輸入的影像s[n,m]水平方向長度是N，垂直方向長度是M 二維離散傅立葉變換定義為 S ( j ω 1 , j ω 2 ) = ∑ m = − ∞ ∞ ∑ n = − ∞ ∞ s [ n , m ] e − j ( ω 1 n + ω 2 m ) {\displaystyle S(j\omega _{1},j\omega _{2})=\sum _{m=-\infty }^{\infty }\sum _{n=-\infty }^{\infty }s[n,m]e^{-j(\omega _{1}n+\omega _{2}m)}} 二維離散傅立葉逆變換定義為 s [ n , m ] = 1 4 π 2 ∫ − ∞ ∞ ∫ − ∞ ∞ S ( j ω 1 , j ω 2 ) e j ( ω 1 n + ω 2 m ) d ω 1 d ω 2 {\displaystyle s[n,m]={\frac {1}{4\pi ^{2}}}\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }S(j\omega _{1},j\omega _{2})e^{j(\omega _{1}n+\omega _{2}m)}\,d\omega _{1}\,d\omega _{2}} 同樣為了方便，我們可將上述兩式改為向量形式 S ( j ω → ) = ∫ − ∞ ∞ ∫ − ∞ ∞ s ( n → ) e j ω → t n → d n → ,     {\displaystyle S(j{\overrightarrow {\omega }})=\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }s({\overrightarrow {n}})e^{j{\overrightarrow {\omega }}^{t}{\overrightarrow {n}}}\,d{\overrightarrow {n}},\ \ } s ( n → ) = 1 4 π 2 ∫ − ∞ ∞ ∫ − ∞ ∞ S ( j ω → ) e j ω → t n → d ω → {\displaystyle s({\overrightarrow {n}})={\frac {1}{4\pi ^{2}}}\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }S(j{\overrightarrow {\omega }})e^{j{\overrightarrow {\omega }}^{t}{\overrightarrow {n}}}\,d{\overrightarrow {\omega }}} 其中, n → = ( n , m ) t {\displaystyle {\overrightarrow {n}}=(n,m)^{t}} 。 下列為二維連續傅立葉變換的性質，大致與一維傅立葉變換的性質相似 假如 s ( t 1 , t 2 ) {\displaystyle s(t_{1},t_{2})} 可分解為 f ( t 1 ) g ( t 2 ) {\displaystyle f(t_{1})g(t_{2})} ，因為自然指數部分也可分解，則 S ( j ω 1 , j ω 2 ) = ∫ − ∞ ∞ ∫ − ∞ ∞ s ( t 1 , t 2 ) e j ( ω 1 t 1 + ω 2 t 2 ) d t 1 d t 2 = ∫ − ∞ ∞ f ( t 1 ) e j ω 1 t 1 d t 1 ∫ − ∞ ∞ g ( t 2 ) e j ω 2 t 2 d t 2 = F ( j ω 1 ) G ( j ω 2 ) = F { f ( t 1 ) } F { g ( t 2 ) } {\displaystyle S(j\omega _{1},j\omega _{2})=\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }s(t_{1},t_{2})e^{j(\omega _{1}t_{1}+\omega _{2}t_{2})}\,dt_{1}\,dt_{2}=\int _{-\infty }^{\infty }f(t_{1})e^{j\omega _{1}t_{1}}\,dt_{1}\int _{-\infty }^{\infty }g(t_{2})e^{j\omega _{2}t_{2}}\,dt_{2}=F(j\omega _{1})G(j\omega _{2})={\mathcal {F}}\{f(t_{1})\}{\mathcal {F}}\{g(t_{2})\}} 稱為分離性。 假設 s ( t → ) {\displaystyle s({\overrightarrow {t}})} 為實數，則 S ( j ω → ) = S ∗ ( − j ω → ) {\displaystyle S(j{\overrightarrow {\omega }})=S^{*}(-j{\overrightarrow {\omega }})} 假設 s ( t → ) {\displaystyle s({\overrightarrow {t}})} 為實數且偶對稱，則 S ( j ω → ) {\displaystyle S(j{\overrightarrow {\omega }})} 為實數且偶對稱 假設 s ( t → ) {\displaystyle s({\overrightarrow {t}})} 為實數且奇對稱，則 j S ( j ω → ) {\displaystyle jS(j{\overrightarrow {\omega }})} 為實數且奇對稱 F { s ( t → − t 0 → ) } = e − j ω → t t 0 → S ( ω → ) {\displaystyle {\mathcal {F}}\{s({\overrightarrow {t}}-{\overrightarrow {t_{0}}})\}=e^{-j{\overrightarrow {\omega }}^{t}{\overrightarrow {t_{0}}}}S({\overrightarrow {\omega }})} 一維傅立葉變換的平移概念是讓信號在時間軸上延遲、提前，相對的，二維傅立葉變換的平移概念是讓影像在空間域作位移，但兩者對應在頻域的概念卻是相同的，也就是將頻域作相位位移。 F { ∂ n s ( t → ) ∂ t i n } = ( j ω i ) n S ( j ω → ) {\displaystyle {\mathcal {F}}\{{\frac {\partial ^{n}s({\overrightarrow {t}})}{\partial t_{i}^{n}}}\}=(j\omega _{i})^{n}S(j{\overrightarrow {\omega }})} 原函數n階偏微分的傅立葉變換，為原函數的傅立葉變換乘上 ( j ω i ) n {\displaystyle (j\omega _{i})^{n}} 。 由於二維傅立葉變換與一維傅立葉變換同為積分運算，所以仍有線性特性。 F { a f ( t → ) + b g ( t → ) } = ∫ − ∞ ∞ ∫ − ∞ ∞ [ a f ( t → ) + b g ( t → ) ] e j ω → t t → d t → = a ∫ − ∞ ∞ ∫ − ∞ ∞ f ( t → ) e j ω → t t → d t → + b ∫ − ∞ ∞ ∫ − ∞ ∞ g ( t → ) e j ω → t t → d t → {\displaystyle {\mathcal {F}}\{af({\overrightarrow {t}})+bg({\overrightarrow {t}})\}=\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }[af({\overrightarrow {t}})+bg({\overrightarrow {t}})]e^{j{\overrightarrow {\omega }}^{t}{\overrightarrow {t}}}\,d{\overrightarrow {t}}=a\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }f({\overrightarrow {t}})e^{j{\overrightarrow {\omega }}^{t}{\overrightarrow {t}}}\,d{\overrightarrow {t}}+b\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }g({\overrightarrow {t}})e^{j{\overrightarrow {\omega }}^{t}{\overrightarrow {t}}}\,d{\overrightarrow {t}}} = a F { f ( t → ) } + b F { g ( t → ) } {\displaystyle =a{\mathcal {F}}\{f({\overrightarrow {t}})\}+b{\mathcal {F}}\{g({\overrightarrow {t}})\}} S ( j ω → ) = ∫ − ∞ ∞ ∫ − ∞ ∞ s ( t → ) e j ω → t t → d t → ,     {\displaystyle S(j{\overrightarrow {\omega }})=\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }s({\overrightarrow {t}})e^{j{\overrightarrow {\omega }}^{t}{\overrightarrow {t}}}\,d{\overrightarrow {t}},\ \ } F { s ( a t → ) } = 1 | a | S ( j ω → a ) {\displaystyle {\mathcal {F}}\{s(a{\overrightarrow {t}})\}={\frac {1}{|a|}}S({\frac {j{\overrightarrow {\omega }}}{a}})} 一維傅立葉變換的線性比例調整是讓信號在時間軸上加速或減速，相對的，二維傅立葉變換的線性比例調整則是讓影像在空間域以某個點放大或縮小，但不論是在時域或空間域作線性比例調整，均也會在頻域產生線性比例調整。 F { f ( t → ) ∗ g ( t → ) } = F { f ( t → ) } F { g ( t → ) } {\displaystyle {\mathcal {F}}\{f({\overrightarrow {t}})*g({\overrightarrow {t}})\}={\mathcal {F}}\{f({\overrightarrow {t}})\}{\mathcal {F}}\{g({\overrightarrow {t}})\}} 卷積特性大至與一維傅立葉變換相同。 聲音是只有時間一個維度的信號，而靜態的平面影像則是具有兩個維度的信號。一維傅立葉轉換常被用來分析隨時間變化的信號，將信號變換至頻域並得到信號的頻域成分強度與相位。然而，對一張的靜態的平面影像而言，雖然沒有時間的概念在裡面，取而代之的是二維空間域的觀念。二維傅立葉變換的目的是將信號由空間域轉到頻域，除此之外，我們也引入了空間頻率這個新的概念來取代傳統一維空間中的頻率。 對傳統隨時間變化的信號而言，以聲音為例，低頻成分的物理意義即為聲音的低音，而高頻成分的物理意義則是聲音的高音。但對空間頻率而言，空間頻率中的低頻成分指的是圖片中顏色緩慢變化的部分。相對的，空間頻率中的高頻成分則是指圖片中顏色迅速變化的部分，比方說物體的邊線。 我們可以用二維的濾波器分別將圖片的低頻與高頻成分濾掉以幫助我們了解空間頻率的概念。我們可以發現，當圖片通過低通濾波器後，被濾出來的圖片是一個模糊的影像，這就是圖片的低頻成分。而當圖片通過高通濾波器後，被濾出來的圖片僅剩下邊緣，這就是圖片的高頻成分，一般而言，圖片在頻域的能量大多集中在低頻，離散傅立葉餘弦轉換變是利用這個特性來進行有損壓縮。 二維傅立葉變換的作用是將影像由空間域變換到頻域，並對於影像不同頻段的成分進行分析與處理，所以二維傅立葉變換在影像處理領域有著舉足輕重的地位，其主要的應用為 影像分析 影像濾波 影像重建 影像壓縮 Oppenheim & Willsky, “Signals & Systems”, 2nd Ed. 1997, Prentice-Hall http://fourier.eng.hmc.edu/e101/lectures/image_processing/node6.html（页面存档备份，存于互联网档案馆） http://fourier.eng.hmc.edu/e101/lectures/image_processing/node7.html（页面存档备份，存于互联网档案馆） http://homepages.inf.ed.ac.uk/rbf/HIPR2/fourier.htm（页面存档备份，存于互联网档案馆） http://www.cs.toronto.edu/~jepson/csc320/notes/linearFilters2.pdf（页面存档备份，存于互联网档案馆）