基于小波變換的圖像融合方法研究

摘要:介紹了圖像融合的基本原理、結(jié)構(gòu)模型及應(yīng)用領(lǐng)域，簡述圖像小波分解與重構(gòu)的Mallat算法。在Matlab7.0環(huán)境下，采用Daubechies小波(dB4)濾波器組對待融合圖像進(jìn)行了3級正交小波分解，低頻系數(shù)采用邊緣保持，高頻系數(shù)采用基于區(qū)域能量最大的規(guī)則進(jìn)行融合，并對融合圖像進(jìn)行質(zhì)量評價。

關(guān)鍵字:圖像融合;小波變換;融合質(zhì)量評價

中圖分類號:TP391 文獻(xiàn)標(biāo)識碼:A 文章編號:1009-3044(2009)35-10075-02

Image Fusion Based on Wavelet Transform Method

WANG Chen，QI De-ning，CHU Bin-bin，PANG Lu-lu

(Artillery Academy of PLA Information Engineering Teaching and Research，Hefei 230031，China)

Abstract: This paper introduces the basic principles of image fusion， the structure models and applications， the brief image of the Mallat wavelet decomposition and reconstruction algorithms. The Matlab7.0 environment， the use of Daubechies wavelet (dB4) filter treatment of fusion images of the three orthogonal wavelet decomposition， using edge-preserving low-frequency coefficients， high-frequency coefficients of the largest energy， based on the rules of regional integration， and integration of image quality evaluation.

Key words: image fusion; wavelet transform; integration of quality assessment

20世紀(jì)90年代以來，隨著圖像傳感器技術(shù)的迅猛發(fā)展，多傳感器圖像融合技術(shù)引起了人們越來越多的關(guān)注。特別是近年來，數(shù)字圖像處理技術(shù)、數(shù)據(jù)融合技術(shù)及小波變換等理論的發(fā)展，使圖像融合技術(shù)研究更加成為眾多研究者的研究熱點。圖像融合在遙感、自動目標(biāo)識別、計算機(jī)視覺、機(jī)器人智能、網(wǎng)絡(luò)安全、工業(yè)檢測、等領(lǐng)域都起著重要的作用，尤其在軍事指揮領(lǐng)域，以多傳感器圖像融合為核心內(nèi)容的戰(zhàn)場感知技術(shù)已成為現(xiàn)代戰(zhàn)爭中最具影響力的軍事高科技技術(shù)[1-2]。

圖像融合就是根據(jù)某一算法，將來自不同傳感器(或同一傳感器在不同時間或不同觀測角度)對同一目標(biāo)或場景觀測得到的多幅圖像進(jìn)行處理，從而得到一幅新的、達(dá)到某種要求的、對目標(biāo)或場景的描述更加準(zhǔn)確、更加全面、更加可靠的圖像。數(shù)字圖像融合充分利用了多個被融合圖像包含的互補(bǔ)信息，大大增加了融合圖像包含的信息量，同時也將多幅被融合圖像中的冗余信息去除掉，提高了系統(tǒng)的可靠性，從而高效利用由多傳感器獲取的圖像信息。數(shù)字圖像融合系統(tǒng)結(jié)構(gòu)模型如圖1所示。該模型將圖像融合評價的信息加入到融合規(guī)則的選取和參數(shù)的選擇過程中，可以更充分地利用信息源提供的信息[1-3]。

1 圖像的小波分解與重構(gòu)[1，4]

對二維圖像信號進(jìn)行小波分解與重構(gòu)，即在空間L2(R2)對信號進(jìn)行二維可分離正交多分辨率分析，Mallat算法的實現(xiàn)使得小波在圖像處理領(lǐng)域的應(yīng)用成為可能。圖像Mallat算法一層小波分解即對圖像矩陣先進(jìn)行行小波變換，再進(jìn)行列小波變換得到圖像的四個頻帶——低頻近似部分子圖LL，高頻細(xì)節(jié)部分水平方向子圖HL、垂直方向子圖LH和對角線方向子圖HH，如圖2所示，下一層分解僅在LL子圖上進(jìn)行。

設(shè){Vj}j∈Z 是L2(R)的一個多分辨率分析，?準(zhǔn)為尺度函數(shù)，?鬃為小波函數(shù)，{hk}k∈Z 為對應(yīng)尺度函數(shù)的低通濾波器系數(shù)，{gk}k∈Z 為對應(yīng)小波函數(shù)的高通濾波器系數(shù)，圖像小波分解的Mallat算法如下:

圖像小波重構(gòu)的Mallat算法如下:

2 基于小波變換的圖像融合

2.1 融合規(guī)則的確定[1，5]

圖像融合過程中，融合規(guī)則的選擇對于最終的融合圖像的質(zhì)量是至關(guān)重要，設(shè)計圖像融合的融合規(guī)則的理論基礎(chǔ)是小波變換后低頻子帶表征的是圖像近似部分，而高頻子帶表征的是圖像的細(xì)節(jié)信息。高頻子帶的系數(shù)在零值左右波動，絕對值越大的系數(shù)表示該處灰度變化越劇烈，即包含圖像的重要信息，如圖像的邊緣、線條以及區(qū)域的邊界。另外，同一場景經(jīng)過不同的傳感器得到的圖像，其低頻近似部分的系數(shù)值差別不大，而高頻細(xì)節(jié)部分卻存在顯著差異。因此，本文采用小波域低頻系數(shù)采用邊緣保持，高頻系數(shù)采用基于區(qū)域能量最大的規(guī)則。

2.2 圖像融合步驟[3，6]

對二維圖像進(jìn)行N層小波分解，最終有(3N+1)個不同頻帶，其中包含3N個高頻帶和一個低頻帶?；谛〔ǘ喑叨确纸鈭D像融合方案如圖3所示，圖像融合的基本步驟為:

1)將待融合圖像進(jìn)行小波塔式分解;

2)按融合規(guī)則對各頻帶層分別進(jìn)行融合，得到融合后各高頻帶和低頻帶小波系數(shù);

3)將融合后各高頻帶和低頻帶小波系數(shù)進(jìn)行圖像重構(gòu)，得到融合后的圖像。

2.3 圖像融合實驗

待融合圖像如圖4(a)(b)所示，在Matlab7.0環(huán)境下，采用Daubechies小波(dB4)濾波器組對待融合圖像進(jìn)行了3級正交小波分解，對兩個分解后的圖像以小波域低頻系數(shù)取平均、高頻系數(shù)模值取大的融合規(guī)則進(jìn)行融合，經(jīng)過小波逆變換，重建圖像，實驗結(jié)果如圖4(d)所示。圖4(c)為對應(yīng)像素灰度值取大的空域直接融合結(jié)果，從視覺感受而言，小波域融合能夠取得較好的效果。

2.4 融合圖像質(zhì)量評價

融合后圖像效果的質(zhì)量評價研究，對于在實際應(yīng)用中選擇適當(dāng)?shù)娜诤纤惴?，以及對現(xiàn)有融合算法的改進(jìn)和研究新的融合算法都具有十分重要的意義[2，4]，一般有主觀評價和客觀評價兩類。主觀評價是由專家進(jìn)行諸如好、較好、一般、不好等定性的等級評分，主觀因素造成評價標(biāo)準(zhǔn)的不同一，具有很大局限性;從信息理論與圖像處理的角度出發(fā)，依據(jù)評定方法所需條件的不同，圖像融合效果的客觀質(zhì)量評價分為基于標(biāo)準(zhǔn)參考圖像的質(zhì)量評價和無參考質(zhì)量評價[2，3]。通常情況下，不存在標(biāo)準(zhǔn)參考圖像，因此無參考質(zhì)量評價更具有實際意義，本文在無標(biāo)準(zhǔn)參考圖像的情況下，根據(jù)圖像自身統(tǒng)計特性對其進(jìn)行評價，選取以下四個常用評價指標(biāo)(設(shè)融合圖像用F表示，L表示圖像F的總灰度級數(shù)，行數(shù)M、列數(shù)為N)。本文實驗選取各圖像的客觀評價指標(biāo)見表1。

1)灰度均值:指圖像中所有像素灰度值的算術(shù)平均，對人眼反映為平均亮度，其定義為:

2)標(biāo)準(zhǔn)差(Standard Error):反映灰度相對于灰度均值的離散情況，即圖像反差的大小，標(biāo)準(zhǔn)差越大，則圖像灰度級分布越分散，包含更多的信息。其定義為:

3)圖像信息熵(Entropy):反映其包含的信息量的多少，熵值越大，信息量越多，定義為:

其中，p(i)表示灰度值為i的像素數(shù)目與圖像總像素數(shù)之比。

4)清晰度，又稱為平均梯度，反映圖像細(xì)節(jié)反差程度和紋理變化特征，一般來說，平均梯度值越大，表明圖像越清晰。定義如下:

其中，ΔFx ΔFy分別為融合圖像F在X與Y方向上的差分。

3 結(jié)束語

根據(jù)圖像多分辨率分析理論，在小波域采用低頻系數(shù)取平均、高頻系數(shù)模值取大的融合規(guī)則進(jìn)行融合實驗，相對于空域直接融合，取得了較好的效果，但基于區(qū)域特性的融合規(guī)則大大增加了運算量。近年來，多尺度幾何分析開始應(yīng)用于圖像處理領(lǐng)域[7]，該理論的發(fā)展將為圖像融合技術(shù)提供新的思路和解決方案。

參考文獻(xiàn):

[1] 敬忠良，肖剛，李振華.圖像融合—理論與應(yīng)用[M].北京:高等教育出版社，2007:10.

[2] Shi Wenzhong， Zhu ChangQing， Tian Yan.Wavelet-based Image Fusion and Quality Assessment[J].International Journal of Applied Earth Observation and Geoinformation，2005(6):241-251.

[3] 李偉.像素級圖像融合方法及應(yīng)用研究[D].華南理工大學(xué)博士論文，2006.

[4] 孫延奎著，小波分析及其應(yīng)用[M].北京:機(jī)械工業(yè)出版社，2005.

[5] 胥妍.基于小波變換技術(shù)的圖像融合方法的研究與應(yīng)用[D].山東師范大學(xué)碩士論文，2008.

[6] 晃銳，張科，李言俊.一種基于小波變換的圖象融合算法[J].電子學(xué)報，2004(5):750-753.

[7] 焦李成，侯彪，王爽，等.圖像多尺度幾何分析理論及應(yīng)用[M].西安:西安電子科技大學(xué)出版社，2008.