




版權說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權,請進行舉報或認領
文檔簡介
1、第第3 3章章頻域處理頻域處理2022-3-6YUAN Guo-wu , YunNan University23.1 3.1 傅里葉變換基礎傅里葉變換基礎 法國數(shù)學家傅立葉生于1768年,他對世人的最大貢獻在于他指出任何周期函數(shù)都可以表示為不同頻率的正弦和(或余弦和)的形式,每個正弦(或余弦)乘以不同的系數(shù)。 甚至非周期的函數(shù)(曲線是有限的情況下)也可以用正弦(或余弦)乘以加權函數(shù)的積分來表示。 一個恰當?shù)谋扔魇菍⒏盗⑷~變換比作一個玻璃棱鏡,棱鏡是可以將光分成不同顏色成分的物理儀器,每個顏色由光的頻率(波長)決定;而傅立葉變換可看作“數(shù)學的棱鏡”,將函數(shù)基于頻率分成不同的成分。2022-3-6
2、YUAN Guo-wu , YunNan University32022-3-6YUAN Guo-wu , YunNan University4 傅立葉于1807年在法國科學學會上發(fā)表了相關的論文。當時審查這個論文的人,其中有兩位是歷史上著名的數(shù)學家拉格朗日(Joseph Louis Lagrange, 1736-1813)和拉普拉斯(Pierre Simon de Laplace, 1749-1827),當拉普拉斯和其它審查者投票通過并要發(fā)表這個論文時,拉格朗日堅決反對,在近50年的時間里,拉格朗日堅持認為傅立葉的方法無法表示帶有棱角的信號,如在方波中出現(xiàn)非連續(xù)變化斜率。法國科學學會屈服于拉
3、格朗日的威望,拒絕了傅立葉的工作。直到拉格朗日死后15年這個論文才被發(fā)表出來。 誰是對的呢?拉格朗日是對的:正弦曲線無法組合成一個帶有棱角的信號。但是,可以用正弦曲線來非常逼近地表示它,逼近到兩種表示方法不存在能量差別,基于此,傅立葉是對的。2022-3-6YUAN Guo-wu , YunNan University53.1 3.1 傅里葉變換基礎傅里葉變換基礎 分解信號的方法是無窮的,例如也還可以用方波或三角波,但分解信號的目的是為了更加簡單地處理原來的信號。 為什么選正弦曲線? 一個正弦曲線信號輸入后,輸出的仍是正弦曲線,只有幅度和相位可能發(fā)生變化,但是頻率和波的形狀仍是一樣的。且只有正
4、弦曲線才擁有這樣的性質(zhì),正因如此我們才不用方波或三角波來表示。 2022-3-6YUAN Guo-wu , YunNan University63.1 3.1 傅里葉變換基礎傅里葉變換基礎 傅立葉變換:針對信號的的長度是無窮大的。因為正余弦波被定義成從負無窮小到正無窮大 計算機:只能處理有限數(shù)據(jù)。無法把一個長度無限的信號組合成長度有限的信號 解決方法:把長度有限的信號表示成長度無限的信號(1)可以把信號無限地從左右進行延伸,延伸的部分用零來表示;(2)可以把信號用復制的方法進行延伸,這樣信號就變成了周期性離解信號。2022-3-6YUAN Guo-wu , YunNan University7
5、3.1 3.1 傅里葉變換基礎傅里葉變換基礎 傅里葉變換在數(shù)學中的定義是嚴格的。設f(x)為x的函數(shù),如果滿足下面的狄里赫萊條件: ()具有有限個間斷點; ()具有有限個極值點; ()絕對可積。 則有下列二式成立:dxexfuFuxj2)()(dueuFxfuxj2)()(2022-3-6YUAN Guo-wu , YunNan University8 式中x是時域(空間域)變量,u為頻率變量。 如令 , 則有 2 udxexfuFxj)()(deFxfxj)(21)(通常把以上公式稱為傅里葉變換對。通常把以上公式稱為傅里葉變換對。2022-3-6YUAN Guo-wu , YunNan Un
6、iversity9函數(shù)函數(shù)f(x)的傅里葉變換一般是一個復量,它可以由下式表示:的傅里葉變換一般是一個復量,它可以由下式表示:FRjI( )( )( )或?qū)懗芍笖?shù)形式或?qū)懗芍笖?shù)形式FRI( )( )( )22FFej( )( )() ( )( )( ) arctgIR把 叫做 的傅里葉譜,而 叫相位譜。F( )f x( )(sincosjej附 歐拉公式2022-3-6YUAN Guo-wu , YunNan University103.2 3.2 二維離散傅里葉變換(二維離散傅里葉變換(DFTDFT)3.2.1 二維連續(xù)傅里葉變換二維連續(xù)傅里葉變換l 二維連續(xù)函數(shù) f (x, y)的傅里葉變
7、換定義如下:l設 是獨立變量 的函數(shù),且在 上絕對可積,則定義積分 為二維連續(xù)函數(shù) 的傅里葉變換,并定義 為 的反變換。 和 為傅里葉變換對。),(yxfyx, dxdyeyxfvuFvyuxj)(2),(),(),(yxf dudvevuFyxfvyuxj)(2),(),(),(vuF),(yxf),(vuF2022-3-6YUAN Guo-wu , YunNan University113.2.2 3.2.2 二維離散傅里葉變換二維離散傅里葉變換l尺寸為MN的離散圖像函數(shù)的DFT l反變換可以通過對F(u,v) 求IDFT獲得 1010)/(2),(1),(MxNyNvyMuxjeyxfM
8、NvuF1010)/(2),(),(MuNvNvyMuxjevuFyxf2022-3-6YUAN Guo-wu , YunNan University12 lDFT變換進行圖像處理時有如下特點:l(1)直流成分為F(0,0)。l(2)幅度譜|F(u,v)|對稱于原點。l(3)圖像f (x, y)平移后,幅度譜不發(fā)生變化,僅有相位發(fā)生了變化。 ),(),(),(vujIvuRvuF),(),(arctan),(vuRvuIvu2022-3-6YUAN Guo-wu , YunNan University13先來看一個變換實例: 一個原始信號的長度是16,于是可以把這個信號分解9個余弦波和9個正弦
9、波(一個長度為N的信號可以分解成N/2+1個正余弦信號,這是為什么呢?結合下面的18個正余弦圖,我想從計算機處理精度上就不難理解,一個長度為N的信號,最多只能有N/2+1個不同頻率,再多的頻率就超過了計算機所能所處理的精度范圍),如下圖:2022-3-6YUAN Guo-wu , YunNan University149個正弦信號:2022-3-6YUAN Guo-wu , YunNan University159個余弦信號:2022-3-6YUAN Guo-wu , YunNan University16把以上所有信號相加即可得到原始信號。對于以上的變換結果,在程序中又是該怎么表示的,可以看
10、看下面這個示例圖:上圖中左邊表示時域中的信號,右邊是頻域信號表示方法,從左向右表示正向轉換(Forward DFT),從右向左表示逆向轉換(Inverse DFT),用小寫x表示信號在每個時間點上的幅度值數(shù)組, 用大寫X表示每種頻率的幅度值數(shù)組, 因為有N/2+1種頻率,所以該數(shù)組長度為N/2+1,X數(shù)組又分兩種,一種是表示余弦波的不同頻率幅度值:Re X,另一種是表示正弦波的不同頻率幅度值:Im X。Re是實數(shù)(Real)的意思,Im是虛數(shù)(Imagine)的意思,采用復數(shù)的表示方法把正余弦波組合起來進行表示,目的是為了便于表達(在后面我們會知道,復數(shù)形式的傅立葉變換長度是N,而不是N/2+
11、1)。2022-3-6YUAN Guo-wu , YunNan University173.2.3 3.2.3 二維離散傅里葉變換的性質(zhì)二維離散傅里葉變換的性質(zhì)1周期性和共軛對稱性l 周期性和共軛對稱性來了許多方便。l 我們首先來看一維的情況。l設有一矩形函數(shù)為,求出它的傅里葉變換: 其他00)(XxAxfuXjXuxjuxjeuXuXAXdxeAdxexfuF022sin)()(2022-3-6YUAN Guo-wu , YunNan University18幅度譜: uXuXAXuFsin)( (a)幅度譜 (b)原點平移后的幅度譜頻譜圖 2022-3-6YUAN Guo-wu , Yun
12、Nan University19n DFT取的區(qū)間是0,N-1,在這個區(qū)間內(nèi)頻譜是由兩個背靠背的半周期組成的 ,要顯示一個完整的周期,必須將變換的原點移至u=N/2點。n 根據(jù)定義,有 n 在進行DFT之前用(-1)x 乘以輸入的信號 f (x) ,可以在一個周期的變換中(u0,1,2,N1),求得一個完整的頻譜。10102)2/(2)() 1(1)(1)2/(NxNxxuNjxNuxNjexfNexfNNuF2022-3-6YUAN Guo-wu , YunNan University20l 推廣到二維情況。在進行傅里葉變換之前用(-1)x+y 乘以輸入的圖像函數(shù),則有: l DFT的原點,
13、即F(0,0)被設置在u=M/2和v=N/2上。l (0,0)點的變換值為: 即 f (x,y) 的平均值。l 如果是一幅圖像,在原點的傅里葉變換F(0,0)等于圖像的平均灰度級,也稱作頻率譜的直流成分。 )2/, 2/() 1)(,(NvMuFyxfDFTyx1010),(1)0 , 0(MxNyyxfMNF2022-3-6YUAN Guo-wu , YunNan University21(a)原始圖像 (b) 中心化前的頻譜圖 (c) 中心化后的頻譜圖圖4.5 圖像頻譜的中心化 2022-3-6YUAN Guo-wu , YunNan University222可分性可分性l 離散傅里葉變
14、換可以用可分離的形式表示 這里l 對于每個x值,當v0,1,2,N1時,該等式是完整的一維傅里葉變換。 1010/2/2),(11),(MxNyNvyjMuxjeyxfNeMvuF10/2),(1MxMuxjevxFM10/2),(1),(NyNvyjeyxfNvxF2022-3-6YUAN Guo-wu , YunNan University23n二維變換可以通過兩次一維變換來實現(xiàn)。n同樣可以通過先求列變換再求行變換得到2D DFT。 圖 二維DFT變換方法2022-3-6YUAN Guo-wu , YunNan University243離散卷積定理離散卷積定理l 設f(x,y)和g(x,
15、y) 是大小分別為AB和CD的兩個數(shù)組,則它們的離散卷積定義為l 卷積定理卷積定理 1010),(),(),(*),(MmNnnymxgnmfyxgyxf),(),(),(*),(vuGvuFyxgyxfDFT2022-3-6YUAN Guo-wu , YunNan University25【例】【例】用MATLAB實現(xiàn)圖像的傅里葉變換。 解:解:MATLAB程序如下: A=imread(pout.tif); %讀入圖像 imshow(A); %顯示圖像 A2=fft2(A); %計算二維傅里葉變換 A2=fftshift(A2); %將直流分量移到頻譜圖的中心figure, imshow(l
16、og(abs(A2)+1),0 10); %顯示變換后的頻譜圖2022-3-6YUAN Guo-wu , YunNan University26 (a)原始圖像 (b)圖像頻譜傅里葉變換2022-3-6YUAN Guo-wu , YunNan University273.3 3.3 頻域濾波頻域濾波 頻域處理法的基礎是卷積定理。它采用修改圖像傅里葉變換的方法實現(xiàn)對圖像的增強處理。由卷積定理可知,如果原始圖像是 ,處理后的圖像是 ,而 是處理系統(tǒng)的沖激響應,那么,處理過程可由下式表示 其中代表卷積。 f x y( , )g x y( , )h x y(,)g x yh x yf x y( , )
17、( , )( , )2022-3-6YUAN Guo-wu , YunNan University28 如果 , , 分別是 的傅立葉變換,那么,上面的卷積關系可表示為變換域的乘積關系,即 式中, 為傳遞函數(shù)。 ),(vuG),(vuH),(vuFg x y h x yf x y( , ), ( , ),( , ) ,() ,() ,(vuFvuHvuG),(vuH2022-3-6YUAN Guo-wu , YunNan University29 在增強問題中, 是給定的原始數(shù)據(jù),經(jīng)傅立葉變換可得到 。選擇合適的 ,使得由式 得到的 比 在某些特性方面更加鮮明、突出,因而更加易于識別、解譯。
18、例如,可以強調(diào)圖像中的低頻分量使圖像得到平滑,也可以強調(diào)圖像中的高頻分量使圖像的邊緣得到增強等等。以上就是頻域處理法的基本原理。f x y( , ),(vuF),(vuH) ,() ,(F ) ,(1vuFvuHyxgg x y( , )f x y( , )2022-3-6YUAN Guo-wu , YunNan University302022-3-6YUAN Guo-wu , YunNan University313.3.1 低通濾波法 這種方法是一種頻域處理法。在分析圖像信號的頻率特性時,一幅圖像的邊緣、跳躍部分以及顆粒噪聲代表圖像信號的高頻分量,而大面積的背景區(qū)則代表圖像信號的低頻分量
19、。用濾波的方法濾除其高頻部分就能去掉噪聲,使圖像得到平滑。 濾除高頻成分,保留低頻成分2022-3-6YUAN Guo-wu , YunNan University322022-3-6YUAN Guo-wu , YunNan University33 常用的幾種低通濾波器 (1)理想低通濾波器 一個理想的二維低通濾波器的傳遞函數(shù)由下式表示:001 ( , )( , )0 ( , )D u vDH u vD u vD式中 是一個規(guī)定的非負的量,叫做理想低通濾波器的截止頻率。 是從頻率頰的原點到(u,v)點的距離,即 2122u=) ,(vvuD D0) ,(vuD2022-3-6YUAN Guo
20、-wu , YunNan University34理想低通濾波器傳遞函數(shù)徑向剖面圖2022-3-6YUAN Guo-wu , YunNan University35 所謂理想低通濾波器是指以截頻 為半徑的圓內(nèi)的所有頻率都能無損地通過,而在截頻之外的頻率分量完全被衰減。理想低通濾波器可以用計算機模擬實現(xiàn),但是卻不能用電子元器件來實現(xiàn)。 理想低通濾波器在處理過程中會產(chǎn)生較嚴重的模糊和振鈴現(xiàn)象。這種現(xiàn)象正是由于傅立葉變換的性質(zhì)決定的。原因:略D02022-3-6YUAN Guo-wu , YunNan University362022-3-6YUAN Guo-wu , YunNan Universi
21、ty37 (2)巴特沃思(Butterworth)低通濾波器 一個n階布特沃斯低通濾波器的傳遞函數(shù)由下式表示 式中 為截止頻率, 的值由下式?jīng)Q定nDvuDvuH20),(+11=) ,(2122=),(vuvuD0D),(vuD2022-3-6YUAN Guo-wu , YunNan University38 布特沃斯低通濾波器又稱最大平坦濾波器。它與理想低通濾波器不同,它的通帶與阻帶之間沒有明顯的不連續(xù)性。也就是說,在通帶和阻帶之間有一個平滑的過渡帶。通常把 下降到某一值的 那一點定為截止頻率 。上式中是把 下降到原來值的 時的 定為截頻點 。) ,(vuHD0) ,(vuH12v),(uD
22、D02022-3-6YUAN Guo-wu , YunNan University39布特沃斯低通濾波器剖面圖2022-3-6YUAN Guo-wu , YunNan University40 布特沃斯低通濾波器的特點: 1)、由于有平緩的過渡帶,圖像將不會有振鈴現(xiàn)象。 2)、模糊程度大大減小。2022-3-6YUAN Guo-wu , YunNan University41(3)指數(shù)低通濾波器(4)梯形低通濾波器 用低通濾波器進行平滑處理可以使噪聲減低到不顯眼的程度,但是由于低通濾波器對噪聲等成分濾除的同時,對有用高頻成分也濾除,因此,這種去噪的美化處理是以犧牲清晰度為代價而換取的。2022-3-6YUAN Guo-wu , YunNan University423.3.2 高通濾波法 常用的幾種高通濾波器: (1)理想高通濾波器 一個理想的二維高通濾波器的傳遞函數(shù)由下式表示
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經(jīng)權益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
- 6. 下載文件中如有侵權或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 2025至2030中國百貨行業(yè)發(fā)展研究與產(chǎn)業(yè)戰(zhàn)略規(guī)劃分析評估報告
- 2025至2030中國生物貂行業(yè)產(chǎn)業(yè)運行態(tài)勢及投資規(guī)劃深度研究報告
- 2025至2030中國玉米淀粉行業(yè)發(fā)展研究與產(chǎn)業(yè)戰(zhàn)略規(guī)劃分析評估報告
- 影樓團隊培訓課件
- 技術賦能教育實現(xiàn)個性化教學的突破
- 財務報銷流程培訓
- 教育數(shù)據(jù)挖掘的潛力學生在多元評價體系中的應用
- 年度培訓計劃編寫課件
- 技術創(chuàng)新助力教育混和教學模式新發(fā)展
- 智慧城市服務中智能公共服務設施的可持續(xù)發(fā)展融資策略
- 2025至2030全球及中國企業(yè)文件共享和同步(EFSS)行業(yè)產(chǎn)業(yè)運行態(tài)勢及投資規(guī)劃深度研究報告
- 上海金山區(qū)屬國有企業(yè)招聘筆試真題2024
- 2025至2030中國碳化硅陶瓷膜行業(yè)發(fā)展趨勢分析與未來投資戰(zhàn)略咨詢研究報告
- 2025至2030中國生石灰行業(yè)市場深度調(diào)研及發(fā)展趨勢與投資方向報告
- 一通三防管理課件
- 2025秋二年級上冊語文上課課件 2 我是什么
- 胖東來總值班管理制度
- 口腔診室終末消毒流程
- 2024年廣州市荔灣區(qū)社區(qū)專職招聘考試真題
- 切口感染案例分析
- 2025-2030年中國管道運輸行業(yè)市場深度分析及發(fā)展前景與投資研究報告
評論
0/150
提交評論