1. 程式人生 > >系統學習數字影象處理之頻域濾波

系統學習數字影象處理之頻域濾波

最近在看模板匹配,雖然很簡單,但還是想認真過下基礎,因此把訊號處理頻域相關的內容,接著影象處理再過一遍。

理論上,對連續變數t的連續函式f(t)的傅立葉變換為F(u),利用f(t)取樣後的函式重建f(t),則必須滿足取樣定理,取樣函式的傅立葉變換為F'(U),它是連續週期的,因此F‘(u)*H(u),就是F(U),即可重建。要注意的是,f(t)必須是連續的帶限函式。

關於混淆

一維情況下,在取樣過的函式中,混淆總是存在的,儘管原始取樣過的函式是帶限的,但實踐中,我們要做的是限制時間的函式,在我們限制函式的持續時間時,總會引進無限的頻率分量(例如對f(t)*h(t),則實現限制函式的持續時間,但由卷積定理知道它們的結果會產生無限的頻率分類,因為H(u)包含無限的頻率分類,卷積後也是)。因此,沒有有限持續時間的函式是帶限的。反之,一個帶限函式,一定從-∞擴充套件到﹢∞。所以,用有限長度的取樣記錄,混淆不可避免。因為我們不可能得到帶限函式。

正如一維情況下一樣,影象中(二維情況)混淆也是不可避免的。包括空間混淆和時間混淆。主要研究空間混淆,即X,Y方向的有害頻率分量,而時間的混淆,主要體現在車輪看上去與實際轉速相反的例子上,因為人眼取樣率低於車輪轉速,導致高頻被混淆成低頻了。

影象的空間混淆主要是認為引入的缺陷,如線狀特徵中的鋸齒,偽高光及原影象中不存在的模式。

通過稍微散焦被數字化的場景來削弱高頻可以降低混淆的影響,但一定是在取樣前。

空間混淆,主要是欠取樣引起的,因此,採用畫素複製的方式,進行影象重取樣時,混淆的影響會更壞。這在很強邊緣處看起來就是鋸齒。對於內插,採用雙線性內插,則鋸齒現象會減小。同時,為減小混淆,在縮小影象之前,稍微模糊一下影象。

在數字影象中,當掃描介質印刷物(如報紙)時,或者在具有周期分量的影象中,如果它的間隔與取樣間隔可比的時候,就會出現莫爾模式。也就是波紋效應。

對於影象來說,傅立葉譜對影象的平移不敏感,但隨著影象的旋轉而旋轉。但相角很敏感。

影象的灰度資訊由傅立葉譜攜帶,對於影象特性,相位起支配作用。

因為DFT表示式中是周期函式,因此通過DFT的乘積的IDFT計算卷積會出現因為週期的靠近而互相干擾的纏繞問題,需要0填充P>= A + B - 1.一般來說DFT演算法對偶數尺寸的陣列執行較快,因此最好選擇P為滿足條件的最小偶整數。

頻域濾波

給定M*N的影象f(x,y),填充為P*Q,平移變換中心,計算F(u,v),設計濾波器H(u,v),大小與處理後的f(x,y)相同,計算乘積的IDFT,取實部,再移回變換中心,裁剪為M*N。

頻域平滑

理想低通(尖銳),巴特沃斯(高階尖銳,低階平滑),高斯(平滑)。

實際中,由於理想低通的截止頻率不平滑,會有振鈴現象,巴特沃斯濾波器在1階處沒有振鈴,2階有但看不出來,隨著階數升高,振鈴越明顯。高斯則無振鈴,因為其DFT也是高斯。但高斯與2階巴特沃斯相比,效果差點,因此,在需要嚴格控制低頻和高頻截止頻率的過渡的情況下,巴特沃斯低通是最好的選擇。除此之外,高斯最佳。

頻域銳化

即高通 = 1-低通。與低通(理想,巴特沃斯,高斯)類似,但要注意的是巴特沃斯D0/Du。

另外還有,拉普拉斯(要注意歸一化標定問題),與空域類似的鈍化模板,高提升濾波,高頻強調濾波。

同態濾波

同時增強對比度和壓縮動態範圍,利用對影象取自然對數後做頻域濾波在取指數,通過降低支配照射分量的影響,允許低灰度變得可見,同時由於高頻分量被增強,故反射分量(邊緣資訊)被明顯銳化了。也就是影象的照射分量通常由慢的空間變化來表徵,反射分量則往往引起突變這一特性。

選擇濾波

即處理指定頻段或者頻率矩形的小區域。第一類是帶阻和帶通,第二類是陷波濾波器。其中陷波更有用,例如消除波紋效應。陷波直接處理實際的DFT,無須填充平移,卻不產生任何纏繞問題。它的主要應用之一是選擇性的修改DFT的區域性區域。

以上,對頻域濾波做個整體補充個總結,實際中頻域濾波是原型化一個實際問題的濾波器的最快速,直接,有效的方法。

相關推薦

系統學習數字影象處理濾波

最近在看模板匹配,雖然很簡單,但還是想認真過下基礎,因此把訊號處理頻域相關的內容,接著影象處理再過一遍。 理論上,對連續變數t的連續函式f(t)的傅立葉變換為F(u),利用f(t)取樣後的函式重建f(t),則必須滿足取樣定理,取樣函式的傅立葉變換為F'(U),它是連續週期的

系統學習數字影象處理基礎

常用數字影象處理例項有: 伽馬射線成像,即由同位素衰變放射出伽馬射線,由伽馬射線收集儀收集的伽馬射線形成影象。 X射線:穿過物體時吸收量調製X射線的強度,最後能量落在膠片上,形成影象。 紫外線:能量轉換髮光 可見光: 微波:即雷達,反射到雷達的微波能量 無線電波:穿透物體,

數字影象處理空間濾波和銳化(Octave實現)

濾波這一概念可以結合數字訊號處理這一領域中的濾波。而在數字影象處理中濾波可以分為空間域濾波和頻率域濾波。這篇博文主要來學習下空間域濾波。 空間域濾波機理 *空間濾波器由一個鄰域(典型的是一個較小的矩形)構成,對該鄰域所包圍的畫素按照一定的操作計算出目標畫素的值,這一過程就是空

數字影象處理筆記——濾波、取樣和頻譜混疊( Frequency domain filtering; sampling and aliasing)

頻域濾波 頻域濾波就是將訊號先做傅立葉變換再與濾波器頻域響應相乘,最後再做傅立葉反變換得到 低通濾波器 讓我們先來看看理想低通濾波器,理想低通濾波器的頻率響應是一箇中間是1,周圍是0的正方形或圓形,而在時域上是sinc函式,,我們看經過低通濾波器後的影象變得模糊了,但是會發現圖中多

數字影象處理直方圖均衡化(Octave)

直方圖的均衡化是什麼呢? 舉個簡單的例子:在一個圓中有很多石頭,都集中在圓心附近,對其均衡化就是讓這些石頭儘可能的均勻分佈在圓這個區域內。 並且還有一個原則:如果石頭A在原來狀態下距離圓心的距離在所有石頭是第4位,那麼均衡化後仍然是第4位,相對順序不能變。(大概就是一個拉的更寬了)

影象處理頻率數學基礎

複數 複數C的定義如下: 其中R和I是實數,j是虛數,即。 C的共軛複數C*: 極座標下表示覆數: 其中,是該向量和實軸(x軸)的夾角。 根據尤拉公式: 有: ‘ 另外,複函式F(u)可表述為: 其中,R(u)和I(u)分別表示實分量函

數字影象處理(三)】用影象增強談灰度變換

調整索引色影象的調色盤map。如果low_in, high_in, low_out, high_out 和 gamma 都是標量,那麼對 r,g,b 分量同時都做此對映。對於每個顏色分量都有唯一的對映,當 low_in 和 high_in 同時為1*3向量或者 low_out 和 high_out 同時為1

數字影象處理伽馬(gamma)變換

輸入影象灰度值 r,輸出影象灰度值s,變換公式 s = T(r),如下 : 原影象“arial.jpg“,取0.7,3.0,變換後圖像如下,並附三通道分佈直方圖(B,G,R). Arial                                        

數字影象處理二維碼影象提取演算法(二)

二維碼標準: 國外:pdf417,Datamatrix,QR 國內:龍貝,漢信。 影象輸入-> 影象預處理-> 二維碼影象切割(定位)-> 二維碼影象校正-> 譯碼 譯碼:二維碼碼字提取-> 糾錯譯碼-> 資訊譯碼 糾錯譯碼:求解伴隨因

數字影象處理四】影象的放大縮小旋轉

縮放放大: 修改紅框裡面的引數,2.0表示放大兩倍,0.5表示縮小兩倍。 旋轉: 修改紅框裡面的引數,表示旋轉多少度。逆時針。 旋轉編寫程式碼的時候可能會出現這樣的情況: 因為旋轉的時候,

數字影象處理(一)】數字影象處理與相關領域概述

機器視覺就是用機器代替人眼來做測量和判斷。機器視覺系統是指通過機器視覺產品(即影象攝取裝置,分 CMOS 和CCD 兩種)將被攝取目標轉換成影象訊號,傳送給專用的影象處理系統,根據畫素分佈和亮度、顏色等資訊,轉變成數字化訊號;影象系統對這些訊號進行各種運算來抽取目標的特徵,進而根據判別的結果來實現高效的機器控

數字影象處理二維影象的傅立葉變換(1)

1.1、訊號變化越快,說明頻率越大,訊號變化越慢,說明頻率越小。這裡的頻率不一定是通常意義上的頻率,通常的頻率是指週期的倒數,我們把通常意義上的頻率叫時間頻率。廣義上的頻率是指變化的快慢,比如圖片來說,從這個畫素到另外一個畫素的灰度值差距比較大,那麼頻率就比較高

數字影象處理歸一化方法

1.把數變為(0,1)之間的小數 為了資料處理方便提出來的,把資料對映到0~1範圍之內處理,更加便捷快速 例1:{2.5 3.5 0.5 1.5}歸一化後變成了{0.3125 0.4375 0.0625 0.1875} 解法:2.5+3.5+0.5+1.5=8;   2.5

數字影象處理點陣圖在計算機中的儲存結構

點陣圖是windows中廣泛應用的一種影象格式,其後綴名為.bmp.點陣圖也稱為位對映圖片. 一張點陣圖包含了許許多多的畫素點,每個畫素點有不同的顏色。由此構成了五彩斑斕的點陣圖影象, 然而正是因為由畫素點這一基本元素構成,所以在放大圖片的時候,人眼便能夠區分一個一

數字影象處理特徵提取及常用方法

形狀特徵   (一)特點:各種基於形狀特徵的檢索方法都可以比較有效地利用影象中感興趣的目標來進行檢索,但它們也有一些共同的問題,包括:①目前基於形狀的檢索方法還缺乏比較完善的數學模型;②如果目標有變形時檢索結果往往不太可靠;③許多形狀特徵僅描述了目標區域性的性質,要全面描述目標常對計算時間和儲存量有較高的要

Matlab濾波下三大LPF

目錄 一、初始知識 二、理想LPF 三、高斯LPF 四、巴特沃斯(Butterworth)LPF 五、濾波器和頻譜的卷積 六、完整程式碼   Python中OpenCV2下實現參見這裡 一、初始知識 關於 distmatrix

數字影象處理中的空間濾波

影象的空間濾波 空間卷積:                                                     卷積是處理線性系統的基礎,以下介紹的濾波計算方式,都與卷積計算有密切關係,其計算形式與卷積的計算方式有高度的相似性。如果下面介紹的濾波運算

數字影象處理- 3.4 空間濾波 and 3.5 平滑空間濾波器

3.4 空間濾波基礎 • Images are often corrupted by random variations in intensity, illumination, or have poor contrast and can’t be used directly. • Filte

影象增強空間濾波

1、為什麼進行影象增強 影象增強是指增強影象中某些特徵,同時削弱或去除某些不需要的資訊,即為了某種應用而去改善影象的質量,消除噪聲,顯現那些被模糊了的細節或簡單的突出一副影象中感興趣的特徵。所以影象增強並不是增強影象的原始資訊,而是隻針對某一特徵以提高某種資訊的辨別能力,影象增強需要根據需求採用特定的方法,當

頻率濾波基礎二(讀數字影象處理學習halcon)

二維離散傅立葉(DFT)變化及其反傅立葉變換(IDFT) 傅立葉變換通用形式 對通用形式來講,c=1,s=-1即為傅立葉變化(影象空間域轉頻域);c=1/MN,s=1即為逆變換(頻域轉空間域) halcon運算元fft_generic(Image : ImageFFT