基于Contourlet變換和區域能量區域方差的圖像融合方法

王莉

摘 要： 針對圖像融合技術中小波變換方法不能很好的表征圖像邊緣、輪廓和紋理等高維幾何奇異性特征的問題，引入了Contourlet 變換方法。這種變換通過對源圖像的分解、重構等步驟來實現對圖像的融合；融合規則對高頻子帶采用基于區域能量匹配原則，低頻子帶采用區域方差匹配度原則來確定融合系數，更好的表示了圖像的幾何特性和邊緣特征，包含了更多的圖像信息。仿真結果表明，該方法是一種優于小波變換的新的圖像融合方法。

關鍵詞： 圖像融合； 小波變換； Contourlet變換； 融合規則； 融合評價

中圖文分類號：TP391.9 文獻標志碼：A 文章編號：1006-8228（2016）08-60-02

Abstract： Because the Wavelet transform method cannot be a good representation of the image edge， contour and texture and other high dimensional geometric singularity characteristics， so the Contourlet transform method is introduced. The method achieves image fusion by using decomposition， reconstruction and other steps to the source image， and the fusion coefficients are determined by the fusion rule that for the high frequency sub-band by using the regional energy matching principle and for the low frequency sub-band by using the regional variance matching principle， which the geometric features and edge features of the image are represented better， and more information of the image are contained. Simulation results show that proposed method is a new image fusion method that better than the Wavelet transform method.

Key words： image fusion； Wavelet transform； Contourlet transform； fusion rules； fusion evaluation

0 引言

小波變換由于其多分辨率和時-頻局部特性能夠很好的處理一維分段連續信號[1]。然而令人失望的是，在處理二維或更高維圖像信號時，小波變換在一維信號處理時所表現出來的優異特性卻不能得以體現。Contourlet變換作為一個新發展起來的離散圖像融合算法，不但延續了小波多尺度分析的精髓，而且還具有多方向分析的優勢，具有比小波變換更加稀疏的二維信號表示方法[1-2]。

1 Contourlet變換

作為一種新的圖像二維表示方法，多分辨率、局部定位、多方向性、近鄰界采樣和各向異性是Contourlet變換十分顯著的特征。自然圖像的一個非常重要的特性是邊緣輪廓，在圖像處理過程中往往希望能夠利用較少的系數就可以獲取到圖像的這些重要信息，具有多尺度和多方向屬性的Contourlet變換基函數恰恰能做到這一點。Contourlet變換的精髓在于，首先對圖像進行多尺度分解以便獲取相應的邊緣奇異點，然后再利用方向信息判斷這些奇異點的特性，最后利用位置相近的點形成其輪廓段[3]。

2 Contourlet變換圖像融合

假設A、B和F分別表示已經配準的源圖像和融合后的圖像，則Contourlet變換算法的步驟主要由分解、融合和重構三個步驟組成。

⑴ 分解：利用Contourlet變換規則對兩幅源圖像A和B分別進行相應分解，在此過程中將會獲取到低頻部分圖像和高頻部分圖像。

⑵ 融合：采用不同的融合規則，分別處理分解后的低頻和高頻系數，形成各尺度方向上的Contourlet融合系數。

⑶ 重構：利用Contourlet逆變換的相關規則對步驟。

3 融合規則

除多尺度多方向分解和重構對圖像融合結果起重要作用外，融合規則的選擇屬于圖像融合的又一個關鍵核心，它將直接影響圖像融合的效果。

高頻子帶反映的是圖像的細節信息，通常圖像清晰度越高的區域其分解后的高頻子帶能量也越大。故高頻子帶融合可選擇基于區域能量匹配的融合規則，其表達式如下。

匹配度表達式為：

一般來講，圖像的方差與圖像包含的信息量成正比關系，即方差越大，圖像包含的信息越豐富。所以可以采用區域方差匹配度來確定低頻子帶融合系數。需先定義區域方差匹配度為：

4 評價指標

融合后的圖像如何判斷其優劣，通常需要通過相應的評價指標和評價方法來判定。

⑴ 互信息MI

MI計算源圖像有多少信息轉移到了融合結果中，其值越高表明融合的效果越好[6]。

用MIFA和MIFB分別表示融合圖像F與待融合圖像A、B的互信息：

⑵ 邊緣信息保留測度QAB/F

QAB/F利用Sobel邊緣檢測來衡量有多少邊緣信息從源圖像轉移到融合圖像[7]。

邊緣信息保留測度可以用下式表示：

5 實驗結果

圖2（a）和圖2（b）為兩幅待融合的原始圖像。圖2（a）為CT圖像，圖2（b）為MRI圖，其中圖2（a）較清晰的顯示出骨骼組織部分，圖2（b）較清晰的顯示出軟組織部分，這兩幅圖從不同的角度反映了人體結構的不同信息。通過對兩幅圖像的融合，實現了信息的互補，使得融合后圖像的信息更全面。表1是三種不同融合算法即：直接融合算法、小波變換融合算法和Contourlet變化融合算法得到的融合圖像和評價指標，通過客觀指標中的值可以看出，Contourlet變化融合算法具有較好的效果和優勢。

6 結束語

Contourlet變化融合算法提供了一種更加稀疏的圖像表示方法，其優于以前的加權平均或小波變化等圖像融合算法，并針對高頻和低頻部分，分別采用不同的融合規則，使得圖像融合后的效果更佳，實驗結果數據證明，這種方法是一種更優的圖像融合新方法。本文在理論研究和實驗研究中由于研究時間的不足尚存在一些值得改進的地方，還可將研究的范圍進行拓展，本文主要是針對兩幅醫學圖像進行了融合再進行效果分析，今后可以擴展到對遙感圖像等進行融合處理，這方面的融合研究也具有一定的必要性。

參考文獻（References）：

[1] 蘇瑞文.基于小波分析方法的圖像邊緣檢測[J].硅谷，2011.15：168，184

[2] 劉楊，張毅，潘寧.基于小波變換的圖像壓縮技術[J].現代電子技術，2007.24：112-115

[3] Minh N. Do， Martin Vetterli. The contourlet transform： an efficient directional multiresolution image representation[J].IEEE Trans on Image Processing，2005.14（12）：2091-2106

[4] 郁梅，易文娟，蔣剛毅.基于Contourlet變換尺度間相關的圖像去噪[J].光電工程，2006.33（6）：73-78

[5] 張華，黃輝，劉邦明.結合Contourlet變換和區域特征的圖像融合算法[J].計算機工程與應用，2013.49（22）：160-163

[6] 胡俊峰，唐鶴云，錢建生.基于小波變換醫學圖像融合算法的對比分析[J].中國生物醫學工程學報，2011.30（2）：196-205

[7] 苗鳳美.基于NSCT和區域分割相結合的圖像融合新算法[D].西安電子科技大學，2009.