1. 程式人生 > >數字影象處理matlab版第四章

數字影象處理matlab版第四章

第四章主要講述頻域處理,濾波都是通過傅立葉變換在頻域中實現的。

在頻域原點處變換的值稱為傅立葉變換的直流分量。

在實際應用中,DFT及其逆變換可以通過使用快速傅立葉變換(FFT)演算法來實現。

matlab中函式fft2實現對陣列M×N影象陣列f的傅立葉變換,F=fft2(f)。

傅立葉頻譜可以使用函式abs來獲得:s=abs(F) 該函式計算陣列的每一元素的幅度(實部和虛部平方和的平方根)。

函式ifft2可以計算傅立葉逆變換,該函式的基本語法為f=ifft2(F),ifft2(F,P,Q)。

空間域和頻域線性濾波的基礎都是卷積定理,該定理可以寫為:f(x,y)*h(h,y)=H(u,v)F(u,v),f(x,y)為影象,h(x,y)為濾波掩模。大寫後為傅立葉變換後結果。

兩個空間函式的卷積可以通過計算兩個傅立葉變換函式的乘積的逆變換得到。相反地,兩個空間函式的卷積的傅立葉變換恰好等於兩個函式的傅立葉變換的乘積。

頻域濾波的目的是選擇一個濾波器傳遞函式H(u,v),以便按照指定的方式修改F(u,v)。

衰減F(u,v)的高頻分量,而保持低頻分量相對不變,具有這種特性的濾波器稱為低通濾波器。

空間卷積通常使用較小的掩模來簡化,使用這種較小的掩模的目的是儘可頻域對應內容的顯著特性。

函式paddedisze用於計算滿足前述等式的P和Q的最小偶數值,通常處理偶數維陣列以加快FFT的計算速度。

DFT濾波的基本步驟,1使用函式paddedsize獲得填充引數2得到使用填充的傅立葉變換3生成一個大小為PQ(1)×PQ(2)的濾波函式H4將變換乘以濾波函式5獲得G的傅立葉逆變換的實部6將左上部的矩形修剪為原始大小。

有限衝擊響應(FIR)濾波器,由線性系統理論可知,在某種適度條件下,輸入到線性系統的一個衝擊完全可以表徵系統。

頻域濾波函式dftfilt:g=dftfilt(f,H)

空間濾波器h,頻域濾波器H,轉換的明顯方法是令H=fft2(h,PQ(1),PQ(2)),其中向量PQ的值取決於我們想要對其濾波的影象的大小。

函式freqz2用於計算FIR濾波器的頻率響應,如:H=freqz2(h,R,C),其中h是一個二維空間濾波器,H是相應的二維頻域濾波器,R為H的行數,C為H的列數。

通過建立一副閾值二值影象,我們可更清楚地看到邊緣。

M函式的核心是需要在頻率矩形中計算任意點到指定點的距離函式,FFT的計算假設變換的原點在頻率矩陣的左上角,資料通過函式fftshift重新排序,以實現視覺化(因而原點的值被轉換為頻率矩形的中心)。

低通頻域濾波器:n階巴特沃茲低通濾波器(BLPF),高斯低通濾波器(GLPF)

函式lpfilter用於生成本節中討論的所有低通濾波器的傳遞函式

三維線框圖及表面圖,這些圖形可用於視覺化二維濾波器的傳遞函式。

對於給定的二維函式H,繪製線框圖的最簡單方法是使用函式mesh,該函式的語法為mesh(H)

函式surf(H)繪製表面圖,代替線框圖。

低通濾波器使影象變模糊,高通濾波器使影象變清晰。

高通濾波器函式hpfilter。

若偏移量與將濾波器乘以一個大於1的常數結合起來,則這種方法就稱為高頻強調濾波,因為該常量乘數突出了高頻部分。

將高頻濾波與直方圖均衡化結合起來使用,所得到的結果要好於單獨使用任何一種方法所得到的結果。