區(qū)域方差和點銳度相結合的多聚焦圖像融合

趙迪迪，季軼群*

(1.蘇州大學 光電科學與工程學院 & 蘇州納米科技協(xié)同創(chuàng)新中心，江蘇 蘇州 215006；2．江蘇省先進光學制造技術重點實驗室 & 教育部現(xiàn)代光學技術重點實驗室，江蘇 蘇州 215006)

1 引 言

對于光學成像系統(tǒng)，由于其景深范圍有限，無法使場景內所有目標同時清晰成像，需對場景內不同區(qū)域分別聚焦成像，通過融合算法，將不同區(qū)域的聚焦圖像融合成一張全視場的清晰光學圖像，這種算法即多聚焦圖像融合。多聚焦圖像融合方法，對于實現(xiàn)光學系統(tǒng)在大視場范圍內獲得高分辨率的光學圖像具有重要意義。

多聚焦圖像的融合算法可分為基于空間域[1-3]和基于變換域[4-6]兩種?；诳臻g域，是在圖像灰度空間上直接對像素點進行融合處理，速度較快，但當某些區(qū)域噪聲較大或灰度變化平緩時，易產(chǎn)生塊效應；基于變換域，是先對源圖像進行變換，如金字塔變換、小波變換等，然后對變換系數(shù)按一定規(guī)則進行融合，得到融合圖像的變換系數(shù)，最后通過逆變換得到融合圖像。其中，基于變換域的圖像融合是目前研究的熱點，張曉菲等[7]提出基于區(qū)域能量特征的小波變換圖像融合方法, 彌補了傳統(tǒng)的小波變換融合算法中出現(xiàn)的信息熵減小、清晰度降低的現(xiàn)象,但此方法在處理邊緣細節(jié)時發(fā)現(xiàn)，對比度下降、易產(chǎn)生重影。葉明等[8]提出區(qū)域清晰度的小波變換圖像融合算法，可很好地消除傳統(tǒng)空域融合方法存在的塊效應，但圖像噪聲較大、像質下降明顯。

拉普拉斯金字塔變換是Burt和Adelson提出的一種多分辨率、多尺度的圖像融合方法[9]，本文基于拉普拉斯金字塔變換提出一種區(qū)域方差和點銳度相結合的融合算法。對金字塔頂層圖像，即低頻系數(shù)，采用區(qū)域方差作為融合度量，進行融合；其余各層為高頻系數(shù)，采用點銳度作為融合度量，進行融合。最后做拉普拉斯金字塔逆變換，得到完整的全視場清晰光學圖像。

2 拉普拉斯金字塔分解

對同一場景不同位置分別聚焦成像得到前景聚焦圖像圖1(a)和后景聚焦圖像圖1(b)，作為源圖像A和源圖像B。在前景聚焦圖像中，黑色測試卡刻度清晰，灰色測試卡模糊。后景聚焦圖像中，灰色測試卡刻度清晰，黑色測試卡模糊。

(1)

(2)

(3)

圖1 源圖像Fig.1 Source images

3 融合規(guī)則

拉普拉斯金字塔分解將兩幅源圖像分別分解到不同頻率層，在不同頻率層采用不同的融合規(guī)則，從而在融合圖像中保留各源圖像在不同頻率層的特征信息。

(4)

式中，LPN為金字塔頂層圖像，即低頻系數(shù)矩陣，μ表示鄰域像素均值，w為加權矩陣。

當0≤l

(5)

式中：m,n為所選局部區(qū)域的長和寬，df為灰度變化幅值，dx為像元間距離增量。公式(5)可表示為：對區(qū)域中的每點與其8鄰域像素點做差，對差值取絕對值后求加權和，最后將所有點所得值相加除以像素總數(shù)。

由于區(qū)域方差和點銳度的數(shù)量級不同，采用式(6)進行歸一化處理：

(6)

式中:X(i,j)為區(qū)域方差或點銳度值，Xmax和Xmin分別為矩陣X中的最大值和最小值，Y(i,j)為歸一化后的區(qū)域方差或點銳度值。

(7)

(8)

本文融合方法流程圖如圖2所示。

圖2 區(qū)域方差和點銳度相結合的圖像融合方法流程圖Fig.2 Image fusion framework combining regional variance and EAV

4 拉普拉斯金字塔逆變換

(9)

圖3 本文方法獲得的融合圖像Fig.3 Proposed method fusion image

5 仿真對比

對同一場景，分別采用基于絕對值取大的離散小波變換、基于區(qū)域能量特征的小波變換和區(qū)域清晰度的小波變換圖像融合算法進行圖像融合，得到如圖4(b)、(c)、(d)所示相應的融合圖像。

上述3種方法均實現(xiàn)了前景聚焦和后景聚焦圖像的融合，保留了源圖像的互補信息。但前兩種方法對比度降低，刻度周圍有重影。第三種方法細節(jié)信息模糊，圖像噪聲較大，像質下降明顯。通過對比，本文采用區(qū)域方差和點銳度相結合的算法，得到的融合圖像具有細節(jié)信息豐富，無重影，清晰度高，對比度高等優(yōu)點。

圖4 不同融合方法的性能比較Fig.4 Fusion performance comparison of different methods

表1給出了融合圖像的評價指標，包括標準差、互信息、邊緣保持度。標準差反映了像素值在其均值上下的波動程度；互信息反映了融合圖像從源圖像中獲得的信息量；邊緣保持度反映了源圖像轉移到融合圖像的邊緣信息量。3種評價指標的值越大，融合效果越好?？梢钥闯霰疚姆椒ǖ玫降娜诤蠄D像評價指標均高于其他3種。

表1 不同融合方法的性能指標比較Tab.1 Comparison of performance indexes of different methods

6 結 論

針對傳統(tǒng)光學成像系統(tǒng)景深范圍有限，難以在大視場范圍實現(xiàn)高分辨率成像的問題。本文提出一種基于區(qū)域方差和點銳度相結合的拉普拉斯金字塔變換融合算法。首先對多聚焦源圖像進行拉普拉斯金字塔變換；然后，對拉普拉斯金字塔的頂層圖像，采用區(qū)域方差考量鄰域像素間的關聯(lián)性,有效保留圖像的細節(jié)信息；對其余各層，采用點銳度對像素點周圍的灰度擴散程度做統(tǒng)計，提高融合圖像清晰度；最后，通過拉普拉斯金字塔逆變換，得到了噪聲低、邊緣信息豐富的多聚焦融合圖像。進一步從標準差、互信息和邊緣保持度等評價函數(shù)對其進行定量評價，結果表明本文提出的融合算法優(yōu)于現(xiàn)有的常用圖像融合方法，標準差、互信息和邊緣保持度分別提高了0.24%、6.8%和8.4%。對光學成像系統(tǒng)在大視場范圍獲得高分辨率的光學圖像具有重要意義。