天津理工大學(xué)中環(huán)信息學(xué)院 喬兆亮 蓋琦

針對低照度條件下彩色圖像的細(xì)節(jié)丟失與視覺效果不佳等問題，提出了一種基于四元數(shù)傅里葉變換的彩色圖像增強(qiáng)算法。把彩色圖像表示為一個四元數(shù)矩陣，整體做傅里葉變換，設(shè)計頻域濾波器，對圖像做濾波處理。實驗結(jié)果表明，該算法能夠明顯提升低照度彩色圖像的視覺效果。

隨著計算機(jī)和網(wǎng)絡(luò)技術(shù)的發(fā)展，圖像已經(jīng)成為信息表達(dá)的一種重要方式，圖像作為一種載體，是信息的最直觀的表示。但圖像在拍攝過程中會由于光照不足、噪聲等因素影響，使得圖像的視覺質(zhì)量無法令人滿意，需要做增強(qiáng)處理。圖像增強(qiáng)的方法按照作用域的不同，可分為空間域和頻率域，其中頻率域濾波需要對圖像做傅里葉變換，在變換域?qū)崿F(xiàn)濾波處理[1-4]。

彩色圖像由三個獨立分量組成，如在RGB空間，由R、G、B三個分量組成，一般對于彩色圖像的處理方法，是各個分量單獨進(jìn)行處理，如果在頻率域濾波，需要針對三個分量做三次傅里葉變換，這樣就忽略了彩色圖像各個分量之間的內(nèi)在聯(lián)系。應(yīng)用四元數(shù)理論，可以把彩色圖像的像素作為一個矢量整體進(jìn)行處理，只需要做一次四元數(shù)傅里葉變換即可，體現(xiàn)了圖像的色彩關(guān)聯(lián)[5-7]。

1 頻域濾波器設(shè)計

圖像頻域濾波原理框圖如圖1所示。

圖1 頻域濾波原理框圖Fig.1 Principle block diagram of frequency domain filtering

圖中f(x,y)表示原始圖像，g(x,y)表示處理后的圖像，DFT代表離散傅里葉變換，DFT-1代表離散逆傅里葉變換，F(xiàn)(u,v)表示f(x,y)的頻譜，H(u,v)代表濾波器。只要設(shè)計一個合適的濾波器H(u,v)，即可消除低頻或高頻成分，得到增強(qiáng)后的圖像。本文選擇的H(u,v)三維立體圖如圖2所示。

圖2 H(u,v)三維立體圖Fig.2 Three dimensional stereogram of H(u,v)

濾波函數(shù)可由公式（1）表示為[8]：

其中，D0代表截止頻率，rl和rh的取值決定了濾波器對低頻和高頻分量的影響，c用來控制濾波函數(shù)斜面的銳化程度，若選擇rl<1，rh>1，將同時實現(xiàn)抑制低頻和增強(qiáng)高頻的目的，可以使圖像在壓縮動態(tài)范圍的同時增強(qiáng)圖像的對比度[8]。

2 彩色圖像的四元數(shù)傅里葉變換

一幅彩色圖像可以表示為一個四元數(shù)矩陣，數(shù)學(xué)表示式如公式（2）所示[5]：

對于彩色圖像，fr(x,y)=0。

彩色圖像f(x,y)的離散四元數(shù)傅里葉變換定義如公式（3）所示[5]：

對應(yīng)的逆變換定義見如公式(4)所示[5]：

式中(x,y)和(u,v)分別為空間域和頻率域的坐標(biāo)。μ是一個單位純四元數(shù)，即μ的實部為0，模為1，且μ2=-1。F(u,v)也是一四元數(shù)，可表示如公式(5)所示：F(u,v)=Fr(u,v)+Fi(u,v)i+Fj(u,v)j+Fk(u,v)k（5）包含4個分量，即Fr(u,v)、Fi(u,v)、Fj(u,v)、Fk(u,v)。

3 算法實現(xiàn)步驟

本算法實現(xiàn)步驟如下：

(1)讀取原始彩色圖像f(x,y)，并做四元數(shù)傅里葉變換，得到F(u,v)的4個分量；

(2)選擇不同的參數(shù)，設(shè)計4個濾波器H1(u,v)、H2(u,v)、H3(u,v)、H4(u,v)；

(3)對F(u,v)的4個分量分別做濾波處理，得到濾波后的結(jié)果F'(u,v)；

(4)對F'(u,v)做四元數(shù)逆傅里葉變換，得到濾波后的圖像。

4 實驗結(jié)果

本實驗應(yīng)用MATLAB 2019a軟件編程實現(xiàn)，原始圖像選擇MATLAB圖像庫中自帶的低照度圖像，實驗結(jié)果如圖3所示。

圖3 實驗結(jié)果Fig.3 Experimental results

5 結(jié)語

本文利用四元數(shù)理論，把彩色圖像表示為一個四元數(shù)矩陣，整體做傅里葉變換，并設(shè)計頻域濾波器，對圖像做濾波處理。實驗結(jié)果表明，此算法原理簡單，實現(xiàn)方便，增強(qiáng)效果較好。