一種改進的水下偏振圖像融合算法研究

楊威，王曉曼，趙海麗，張谷祥

（長春理工大學(xué) 電子信息工程學(xué)院，長春 130022）

我國海域遼闊，海洋中儲藏著各種資源，人們的研究重點不斷的被吸引到對海洋的探索。傳統(tǒng)觀念認(rèn)為，海洋環(huán)境相對于陸地環(huán)境，較為單一，海洋干擾因素較少，海洋目標(biāo)更容易被探測識別。但是海洋觀測容易受到氣候條件影響，有時會造成目標(biāo)與背景對比度反差較小，此時采用傳統(tǒng)的成像技術(shù)獲得的光強圖像，很難將目標(biāo)區(qū)分出來。采用偏振成像技術(shù)則可以很好的克服此缺陷，它可以有效地降低海洋環(huán)境對被測目標(biāo)的干擾，如偏振度圖像可以提高海洋目標(biāo)與背景環(huán)境的對比度，使被測目標(biāo)凸顯出來［1-2］。

圖像融合是指將包含同一場景的不同圖像進行綜合，使圖像能夠更好地還原被測目標(biāo)物信息的過程。自上世紀(jì)80年代以來，圖像融合技術(shù)在各個領(lǐng)域就已經(jīng)開始有所應(yīng)用，如偏振遙感圖像融合、多模態(tài)醫(yī)學(xué)圖像融合、海面霧天軍事目標(biāo)識別等。Ruxue Zhao等人提出一種基于邊緣信息增強的多尺度強度和偏振差的圖像融合算法，將偏振強度圖像與局部放電圖像融合，該方法能很好的實現(xiàn)渾濁環(huán)境下圖像的邊緣增強；Gemine Vivone等人將低空間分辨率的多光譜圖像與高空間分辨率的全色圖像進行融合，得到與全色圖像具有相同空間分辨率、與多光譜數(shù)據(jù)具有相同光譜分辨率的合成產(chǎn)物，提高了可見、近紅外和短波紅外波段的空間分辨率［3-5］。

目前，圖像融合技術(shù)雖然用途廣泛，涉及面廣，但是還不夠系統(tǒng)完善，一些方法也只是針對具體的應(yīng)用環(huán)境提出的。本文針對水下復(fù)雜環(huán)境，提出一種改進的基于區(qū)域和能量的圖像融合算法。對于提高水下運動目標(biāo)與背景對比度，將目標(biāo)物從復(fù)雜背景中分離出來，準(zhǔn)確識別被測目標(biāo)具有重要的意義。

1 光強圖像和偏振圖像的獲取

偏振是光的基本性質(zhì)之一，與普通光強相比，它包含了被測物的大量特征信息，一直是國內(nèi)外研究的熱點。偏振光的描述通常用Stokes矢量法來表示。即：

式中，S0或I代表光波的強度，即在不加檢偏裝置的情況下原始的光強信息；S1或Q代表光波的水平偏振分量與垂直偏振分量之差；S2或U代表光波沿45°方向的偏振分量與沿135°方向的偏振分量強度之差；S3或V代表光波的右旋圓偏振分量與左旋圓偏振分量的強度之差，近似為零［6］。

根據(jù)Stokes矢量計算線偏振度（Dolp）和偏振角（Aop）的公式分別為：

對于水下光強圖像和偏振圖像的獲取，本文擬采用普通光強度相機和分焦平面相機分別對水下運動目標(biāo)實時采集。為使實驗效果顯著，目標(biāo)物選擇玩具手槍射出的子彈，本實驗在靜水環(huán)境下進行。實驗平臺搭建完成，利用中型強光手電發(fā)出的自然光經(jīng)過偏振片，波片后照射在運動物體上，通過圖像采集裝置分別采集光強圖像和偏振圖像。圖1為水下光強圖像和偏振圖像采集原理圖，主要由輔助光源系統(tǒng)和圖像采集與處理系統(tǒng)兩部分組成。

圖1 水下光強圖像和偏振圖像采集原理圖

圖2為水下光強圖像和偏振圖像采集部分實物圖，圖（a）中左側(cè)為光強度相機，右側(cè)為分焦平面相機。

圖2 水下光強圖像和偏振圖像采集部分實物圖

圖3、圖4為偏振片旋轉(zhuǎn)90°時，利用光強度相機和分焦平面相機分別采集到的第21幀時的水下運動目標(biāo)的光強圖像和偏振圖像。圖4為采集到的不同角度的偏振分量圖像，與圖3相比，圖4更好的展現(xiàn)了目標(biāo)物信息，對于圖像的進一步處理，提供了較為準(zhǔn)確的數(shù)據(jù)。

圖3 水下光強圖像

圖4 水下偏振圖像

圖5是根據(jù)公式（1）計算得到的三個偏振分量圖像。圖6是根據(jù)3個偏振分量圖像進一步解算出的線偏振度和偏振角圖像，Dolp反映了圖像的邊緣輪廓信息，Aop反映了圖像的表面特征信息，將二者進行融合有利于全面掌握圖像的偏振特征，對于提取水下偏振圖像的邊緣信息具有重要作用。本文采用Toet融合方法，有利于突出原圖像的特征部分。其融合公式為：

圖5 I、Q、U分量圖像

圖6 線偏振度圖像和偏振角圖像

其中，a，b，c為權(quán)重因子，仿真實驗表明，a，b，c分別取 0.6，0.2，0.2時，仿真效果最好；f(x,y)為Dolp與Aop融合得到的偏振特征圖像Dolp-Aop；

圖7 偏振特征圖像Dolp-Aop

2 水下偏振圖像融合

本文針對水下偏振圖像成像質(zhì)量差、易受水下環(huán)境干擾、邊緣細節(jié)模糊等問題，提出一種改進的基于二維離散小波變換的圖像融合算法，二維離散小波變換主要由離散小波分解與離散小波重構(gòu)兩個部分組成。二維離散小波分解主要通過采用選定的算法將源圖像分解為低頻圖像與對角細節(jié)圖像，通過查閱文獻發(fā)現(xiàn)，小波分解的程度并非越深越好，分解的程度過深時，最底層的圖像在融合的過程中，不可避免的會造成信息的丟失，通常這些丟失的信息是不可能被恢復(fù)的，因此本次實驗將進行三層小波分解。二維離散小波重構(gòu)是二維離散小波分解的逆過程，主要通過對二維離散小波分解的反演得到。

本文通過Mallat算法進行圖像的三層小波分解與重構(gòu)。Mallat算法主要是對圖像的每一行、每一列進行濾波處理，而濾波通常借助一組低通濾波器和高通濾波器實現(xiàn)。每一層的小波分解均可得到圖像的低頻分量、水平細節(jié)分量、垂直細節(jié)分量、對角細節(jié)分量，其中低頻分量涵蓋了圖像的主要特征，包含了源圖像的大量能量，其他三個分量則包含了源圖像的邊緣細節(jié)信息［7］。其一尺度分解如圖8所示。

圖8 一尺度小波分解算法

圖中An為低頻圖像，Hn為水平細節(jié)圖像，Vn為垂直細節(jié)圖像，Dn為對角細節(jié)圖像，且An-1是An和Hn、Vn、Dn之和。圖中H和G分別為低通和高通濾波器，當(dāng)小波函數(shù)確定后即可推導(dǎo)低通和高通濾波器。1↓2表示降采樣，水平方向二取一，2↓1表示降采樣，垂直方向二取一。通過圖8得到的是一尺度二維離散小波分解圖像，尺度為三的二維離散小波分解可通過對低頻圖像進行三次小波分解得到，如圖9（a）所示，圖9（b）為經(jīng)過三層小波分解后得到的一系列低頻圖像和細節(jié)圖像。

圖9 三尺度小波分解及分解系數(shù)顯示

二維離散小波重構(gòu)是分解的逆運算，其算法如圖10所示。

圖10 一尺度小波重構(gòu)算法

2.1 改進的二維離散小波變換算法

由離散小波變換算法的原理可知，分解得到的低頻分量和高頻分量所包含的信息側(cè)重點不同，針對低頻分量和高頻分量的特點，應(yīng)采取不同的處理方法，融合方法的選擇會有所不同，融合方法應(yīng)該主要針對它們所包含的信息特點進行選擇，這樣才能保證最終得到的融合圖像達到預(yù)期效果。

通過查閱文獻發(fā)現(xiàn)，傳統(tǒng)的圖像融合算法通常采用加權(quán)平均法，大多應(yīng)用在背景環(huán)境單一的情況下，優(yōu)點是算法簡單，易于實現(xiàn)，但是加權(quán)平均法沒有詳細分析原始圖像各個頻帶的信息，融合出來的圖像質(zhì)量較差，帶有很強的噪聲，并且容易出現(xiàn)塊狀現(xiàn)象［8］。本文在綜合各算法的優(yōu)缺點后，提出一種改進的二維離散小波變換算法，用于對光強圖像和偏振特征圖像Dolp-Aop進行融合處理。

實驗通過光強度相機采集光強圖像，由分焦平面相機采集偏振圖像，進一步處理得到偏振特征圖像Dolp-Aop，對二者采用本文提出的一種改進的二維離散小波變換。具體算法步驟如下所示：（1）對光強圖像和偏振特征圖像Dolp-Aop進行配準(zhǔn)處理，使光強圖像和偏振特征圖像Dolp-Aop達到像素級匹配；（2）分別對光強圖像和偏振特征圖像Dolp-Aop圖像進行三層小波分解，根據(jù)低頻分量和高頻分量的特點，采用不同的融合規(guī)則。對第三層低頻分量采用基于區(qū)域方差加局部能量相結(jié)合的融合規(guī)則，其他兩層低頻分量采用加權(quán)平均法。高頻分量采用基于Sobel算子的融合方法；（3）將融合后的低頻分量和高頻分量進行二維離散小波重構(gòu)，最終得到融合后圖像。融合流程如圖11所示。

圖11 改進的二維離散小波變換流程

2.2 低頻融合策略

低頻分量與源圖像最為接近，源圖像的有用特征信息都集中在低頻部分。對于低頻分量常采用像素絕對值取大或取小法，這種方法算法簡單，易于操作，常用于圖像灰度差異不大的情況下，但是容易引入新的噪聲，對比度也會得到削弱。為了更好的抑制噪聲，避免產(chǎn)生新的噪聲干擾，提高融合性能，本文擬采用一種基于區(qū)域方差加局部能量相結(jié)合的融合策略。假設(shè)M、N為光強圖像和偏振特征圖像經(jīng)過小波分解后的低頻分量，F(xiàn)為融合后的低頻分量，則該方法的具體步驟為：

（1）分別計算M(x,y)和N(x,y)以點(x,y)為中心的3×3鄰域內(nèi)的區(qū)域方差，記為VarM(x,y)和

其中，Var為圖像的某一鄰域內(nèi)的區(qū)域方差；S(x,y)代表圖像在(x,y)點的像素值；(x,y)代表在鄰域內(nèi)的像素平均值；ω(a,b)代表窗口函數(shù)，且A×B的窗口大小為3×3鄰域，窗口函數(shù)ω取為：

（2）當(dāng)L

當(dāng)L≥T時，根據(jù)式（9）計算像素點(x,y)周圍3×3區(qū)域內(nèi)的局部能量，分別記為EM(x,y)、EN(x,y)。

式中，E(x,y)是以(x,y)為像素中心點的區(qū)域局部能量；Q(x,y)為低頻分量；w(x,y)表示以(m,n)為中心點時Q(x,y)對應(yīng)的權(quán)重系數(shù)。

進一步計算兩幅圖像在融合后圖像中的權(quán)重系數(shù)，記為WM和WN，如式（11）所示：

2.3 高頻融合策略

高頻分量與低頻分量不同，它包含的能量較少，所含信息較為分散，但是細節(jié)信息豐富，圖像的邊緣輪廓信息就包含在高頻分量中。本實驗選取的目標(biāo)物為玩具手槍射出的長方形子彈，為準(zhǔn)確識別目標(biāo)物，對于其邊緣輪廓的獲取至關(guān)重要。邊緣檢測的方法有很多，每一種方法都有優(yōu)缺點，需要視具體情況選擇合適的方法。例如常見的Robert算子、Prewitt算子。Rob?ert算子是一個二階矩陣，主要利用的是對兩個對角線上的像素值作差來檢測邊緣，但是處理后的圖像邊緣不平滑，通常應(yīng)用在噪聲較少且邊緣明顯的情況下；Prewitt算子是一個三階矩陣，它通常由水平邊沿檢測模板和垂直邊沿檢測模板組成，水平邊沿檢測模板通過對水平方向上的像素值作差來檢測邊緣；垂直邊沿檢測模板主要通過對垂直方向上的像素值作差來檢測邊緣。采用Prewitt算子檢測邊緣，對于噪聲具有很好的平滑作用，使圖像的邊緣連續(xù)性得到加強，但是容易出現(xiàn)拼接現(xiàn)象。結(jié)合目標(biāo)物的特征，本文對于高頻部分的處理擬采用一種基于Sobel算子的融合策略［9］。

Sobel算子是一種一階微分邊緣檢測算子，它的檢測模板由三階矩陣組成，與Robert、Prewitt算子不同的是，Sobel算子采用的是一種像素平均的方法，能夠很好地抑制環(huán)境噪聲的影響，更好的顯示圖像的邊緣輪廓。同時，它還加入了像素權(quán)重的概念，Sobel算子認(rèn)為與某一像素點距離不同的像素點，會對該像素點產(chǎn)生不同程度的影響，即距離越遠，影響越??；距離越近，影響越大。Sobel算子能夠使邊緣連續(xù)性得到增強，從而突出圖像的邊緣輪廓，獲得較好的邊緣檢測效果。

圖12為Sobel算子的兩組檢測模板，分別為水平邊緣檢測模板、垂直邊緣檢測模板。計算時，分別將光強圖像和偏振特征圖像Dolp-Aop與兩個矩陣作水平方向卷積和垂直方向卷積，即可得到圖像的邊緣輪廓，最終獲得包含邊緣輪廓信息的高頻融合圖像。

圖12 Sobel算子邊緣檢測模板

3 實驗結(jié)果與分析

3.1 融合結(jié)果分析

將已經(jīng)配準(zhǔn)好的光強圖像和偏振特征圖像在matlab R2010b軟件平臺上進行仿真實驗。為了更有效的驗證本文提出的算法對圖像的處理效果，分別采用基于加權(quán)平均法、像素絕對值取小法、參考論文［10］中提出的基于Prewitt算子的圖像融合算法對圖像進行融合處理，將得到的結(jié)果與本文提出的算法進行對比，實驗結(jié)果如圖13所示。

圖13 實驗仿真結(jié)果

由實驗結(jié)果發(fā)現(xiàn)，圖13（a）采用簡單的加權(quán)平均法雖然消除了部分水體噪聲的干擾，但是子彈的輪廓變得模糊，子彈身上出現(xiàn)塊狀現(xiàn)象，難以識別被測目標(biāo)物；圖 13（b）和圖 13（c）是分別采用像素絕對值取小法、基于Prewitt算子的圖像融合算法處理后的結(jié)果，可以看出融合圖像亮度得到了一定的提升，顯示了更多的邊緣輪廓信息，但是還是存在一些環(huán)境噪聲的干擾；圖13（d）為本文算法得到的結(jié)果，觀察發(fā)現(xiàn)，圖像整體亮度得到明顯提升，子彈與背景的對比度也得到加強，細節(jié)信息更加豐富，子彈上尾部凹槽已初步顯現(xiàn)，子彈后端圓形尾部也可明顯觀察到。

3.2 融合圖像的評價指標(biāo)

對于融合圖像的結(jié)果如果僅僅通過直觀的觀察來評價融合圖像的質(zhì)量，有時容易受到主觀因素的干擾，不能客觀反映事實。因此常引入一些客觀評價指標(biāo)，將主觀判斷與客觀評價相結(jié)合，綜合判斷圖像融合結(jié)果。本文擬采用信息熵、平均梯度、空間頻率、標(biāo)準(zhǔn)差來評價融合結(jié)果。

四幅采用不同算法得出的圖像評價指標(biāo)值如表1所示。由表1數(shù)據(jù)發(fā)現(xiàn)，采用加權(quán)平均法得到的圖像，反映其包含的信息量，圖像的清晰度的指標(biāo)都不理想，融合效果最差，相較于后三種算法，差距明顯，這與之前直觀觀察融合圖像的結(jié)論一致。基于Prewitt算子的圖像融合得到的標(biāo)準(zhǔn)差最大，但是通過對比兩幅融合圖像發(fā)現(xiàn)，本文算法得到的圖像包含的信息量更加豐富。而在信息熵、平均梯度、空間頻率等指標(biāo)方面，本文的算法都要優(yōu)于其他三種算法。說明通過本文算法得到的融合圖像對比度高，包含的信息量豐富，可以顯現(xiàn)出更多的細節(jié)信息。

表1 圖像融合評價指標(biāo)

4 結(jié)論

圖像融合技術(shù)是近些年研究的熱點，尤其是在醫(yī)療、軍事等領(lǐng)域，例如水下軍事目標(biāo)的探測識別，未來將有廣闊的前景。本文針對水下運動目標(biāo)提出一種改進的二維離散小波變換，對于低頻分量采用一種區(qū)域方差加局部能量相結(jié)合的圖像融合算法，較好的保留了源圖像的大部分信息；而對于高頻部分，針對其含有的少量邊緣輪廓信息，提出一種基于Sobel算子的圖像融合方法。通過實驗仿真，本文提出的算法相對于其他三種算法有一定的優(yōu)勢，融合圖像顯現(xiàn)出了更多的細節(jié)信息，這對于準(zhǔn)確識別被測物有重要的意義。但是在實驗過程中也發(fā)現(xiàn)，本文提出的算法會引入一些新的噪聲，造成一些新的干擾，所以本文的算法還有待改善。