郭 佳，李 輝，王 成，李思眾

(西安應(yīng)用光學(xué)研究所，西安710065)

隨著傳感器技術(shù)的迅猛發(fā)展，圖像融合技術(shù)作為多傳感器信息融合的一個非常重要的分支——可視信息的融合，近二十年來，引起了世界范圍內(nèi)的廣泛關(guān)注和研究熱潮［1－2］。圖像配準(zhǔn)是圖像融合技術(shù)的基本環(huán)節(jié)和首要問題，它是對取自不同時間、不同傳感器或不同視角的同一景物的兩幅圖像或多幅圖像進(jìn)行匹配、疊加的過程，只有經(jīng)過配準(zhǔn)后的圖像才能進(jìn)行有效的融合［3－6］。

在多傳感器數(shù)據(jù)融合中，兩種常見的傳感器是紅外和可見光傳感器。研究高效準(zhǔn)確的紅外與可見光圖像配準(zhǔn)方法具有廣闊的應(yīng)用前景，特別是在軍事領(lǐng)域中，可以獲得更加全面的精確可靠的目標(biāo)信息，提升對目標(biāo)的精確打擊以及全天候作戰(zhàn)的能力。

傳統(tǒng)的圖像配準(zhǔn)方法一般是采用選取控制點的方法，該方法在相同傳感器下是可行的，但在多傳感器圖像配準(zhǔn)時出現(xiàn)了一定的困難。這是由于多傳感器圖像反映的是物體不同的物理特性:其中一幅傳感器圖像中出現(xiàn)的特征在另一種傳感器圖像中未必會出現(xiàn)，如在紅外圖像出現(xiàn)的熱油箱未必在可見光圖像中出現(xiàn)。另一方面，選取控制點需要人工的干預(yù)，難以實現(xiàn)自動配準(zhǔn)，費時費力。此外，由于紅外與可見光處于不同波段，圖像間的相關(guān)性較小，所以紅外與可見光圖像的配準(zhǔn)相對來說，難度更大。目前，國內(nèi)外有關(guān)紅外與可見光圖像的配準(zhǔn)尚無成熟的算法。本文針對紅外與可見光的圖像特性，采用互信息作為兩幅圖像的統(tǒng)計相關(guān)性測度，并使用仿射模型實現(xiàn)兩者的配準(zhǔn)，得到了較好的效果。

1 互信息理論

互信息(Mutual Information)［7－8］是信息論中的一個基本概念，是兩個隨機(jī)變量統(tǒng)計相關(guān)性的測度。原理上，當(dāng)兩幅基于共同物理結(jié)構(gòu)的圖像達(dá)到最佳配準(zhǔn)時，他們對應(yīng)像素的灰度互信息應(yīng)為最大，由于該測度不需要對不同成像模式下圖像灰度間的關(guān)系作任何假設(shè)，也不需要對圖像進(jìn)行分割或任何預(yù)處理，因此，該測度可以被廣泛的應(yīng)用于多模圖像配準(zhǔn)。

熵［9］是互信息理論中的一個重要概念，它是用來測量一個信源所含信息量的測度。假設(shè)一個信源A輸出N個消息，其中有n個不同的消息，第i個消息(i=1，2，…，n)重復(fù) hi次，則 hi/N 為每個輸出消息的重復(fù)頻率，故可以用概率替換，即pi=hi/N。則該信源的平均信息量(熵)為:

熵表達(dá)的是一個系統(tǒng)的復(fù)雜性或不確定性。對灰度圖像來說，可以將圖像的灰度看作是一個隨機(jī)變量，圖像中每個點的灰度取值為該隨機(jī)變量的一個事件，則可以根據(jù)圖像的灰度信息計算出每級灰度發(fā)生的概率pi=hi/N，其中hi為圖像中灰度值為i的像素點的總數(shù)，N為圖像中的像素總和。如果圖像中的灰度級越多，像素灰度值分布越分散，則每級灰度的概率值很接近，或者說圖像中任一點灰度值具有很大的不確定性，我們所獲得信息量也就越大，該圖像的熵值也較大。反之，如果圖像中的灰度值比較集中，則一些灰度級的概念值也較大，不確定性減少，熵值較小。

聯(lián)合熵H(A，B)是檢測隨機(jī)變量A和B相關(guān)性的統(tǒng)計量。對于兩個隨機(jī)變量A和B，它們的概率分布分別為pA(i)和pB(j)，聯(lián)合概率分布為pAB(i，j)，則它們的聯(lián)合熵H(A，B)為:

如果H(A，B)=H(A)+H(B)，則A和B是相互獨立的。反之，可用兩個變量的聯(lián)合分布與完全獨立時的概率分布的廣義距離來作為變量之間的相似性測度，即互信息:

互信息是用來度量A、B兩個信源之間的依賴程度。對于待配準(zhǔn)的兩幅多傳感器圖像，可以認(rèn)為它們是關(guān)于圖像灰度的兩個隨機(jī)變量集A和B，并且它們存在較強的依賴關(guān)系。i和j分別為兩幅圖像中的對應(yīng)部分的像素灰度值，i、j通過坐標(biāo)變換T相聯(lián)系，則兩幅圖像的互信息I(A，B)為:

其中pA(i)、pB(j)和pAB(i，j)分別為兩幅圖像灰度的邊緣概率分布和它們的聯(lián)合概率分布。

由式(3)、式(4)可以看出，當(dāng)圖像對準(zhǔn)時，A、B相互包含的信息總量應(yīng)為最大。式(4)實際上是測量兩幅圖像灰度的聯(lián)合概率分布 pAB(i，j)和當(dāng)pA(i)與pB(j)獨立時的聯(lián)合概率分布pA(i)pB(j)之間的距離，它是兩幅圖像灰度相關(guān)性的一個測度。當(dāng)采用(4)定義互信息時，互信息法配準(zhǔn)圖像基于如下假設(shè)，當(dāng)兩幅圖像對準(zhǔn)時，兩幅圖像之間的灰度相關(guān)性應(yīng)為最大;當(dāng)誤配準(zhǔn)發(fā)生時，該測度值會降低。

2 仿射變換

仿射變換［10－16］定義為:若變換 C:Rn→ Rn，C(X)=T(X)+a，T 是非奇異線性變換，a∈Rn，則變換C稱為仿射變換。

仿射變換是一種線性變換，即具有平行線變換成平行線，有限點映射到有限點的一般特性，它可以描述3D空間中的目標(biāo)到2D平面的成像過程。具體表現(xiàn)可以是各個方向尺度變換系數(shù)一致的均勻尺度變換或變換系數(shù)不一致的非均勻尺度變換及剪切變換等。它具有平移、旋轉(zhuǎn)、和縮放不變性。仿射變換一般應(yīng)用于場景圖的配準(zhǔn)，或者是簡單/復(fù)雜背景中擴(kuò)展目標(biāo)的配準(zhǔn)。

本文使用的仿射變換模型可以表示為:

其中，(x，y)和(x'，y')分別是兩幅圖像中對應(yīng)點的坐標(biāo)，s為尺度縮放因子，θ為旋轉(zhuǎn)角度，tx，ty為平移量，這三者為圖像配準(zhǔn)的配準(zhǔn)參數(shù)。

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

本文選取分辨率均為180×240的紅外和可見光圖像，對基于圖像互信息的配準(zhǔn)方法進(jìn)行了仿真，實驗結(jié)果如圖1和圖2所示，并給出了實驗相關(guān)數(shù)據(jù)，如表1所示。從圖1和表1中的實驗數(shù)據(jù)中可以看出，圖1中兩幅圖像(b)和(c)間有明顯的尺度、旋轉(zhuǎn)和平移變化，并且，由圖2可以看出，圖1中兩幅圖像(a)和(b)之間的尺度、旋轉(zhuǎn)和平移等差別得到了很好的修正。實驗結(jié)果表明，采用互信息作為紅外和可見光圖像的統(tǒng)計相關(guān)性測度，運用仿射變換模型進(jìn)行成像差異較大的紅外與可見光圖像間的配準(zhǔn)，效果較為理想，具有較高的實用價值。

4 結(jié)論

本文針對紅外與可見光圖像的特性及其配準(zhǔn)的特殊性，提出了一種基于互信息的圖像配準(zhǔn)算法，并進(jìn)行了算法仿真。實驗結(jié)果表明，此算法是有效的，能夠較好的進(jìn)行圖像間尺度、旋轉(zhuǎn)和平移差別的修正，可以很好地解決成像差異較大的紅外與可見光圖像間的配準(zhǔn)問題。

圖1 紅外與可見光圖像基于圖像互信息的配準(zhǔn)結(jié)果

表1 紅外與可見光圖像基于圖像互信息配準(zhǔn)的相關(guān)實驗數(shù)據(jù)

圖2 配準(zhǔn)后的圖像

［1］ 王利平，劉濱，陳錢，等.多傳感器目標(biāo)識別的圖像數(shù)據(jù)融合技術(shù)研究［J］.光電子技術(shù)，1998，18(9):239－242.

［2］ 佘二永.多源圖像融合方法研究［D］.湖南:國防科學(xué)技術(shù)大學(xué)，2005.

［3］ 廖惜春，丘敏，麥漢榮.基于參數(shù)估計的多傳感器數(shù)據(jù)融合算法研究［J］.傳感技術(shù)學(xué)報，2007，20(1):193－197.

［4］ 趙杰，江晶，蓋旭剛.多無源傳感器的在線配準(zhǔn)［J］.傳感技術(shù)學(xué)報，2007，20(8):1894－1898.

［5］ 倪國強，劉瓊.多源圖像配準(zhǔn)技術(shù)分析與展望［J］.光電工程，2004，31(9):1－6.

［6］ 陶冰潔，王敬儒，張啟衡.采用仿射變換的紅外與可見光圖像配準(zhǔn)方法［J］.光電工程，2004，31(11):39－41.

［7］ 江靜，張雪松，蘭西柱.紅外與可見光圖像互信息法自動配準(zhǔn)算法研究［J］.華北科技學(xué)院學(xué)報，2006，3(4):73－78.

［8］ 高鵬.基于邊緣和互信息法的紅外圖像配準(zhǔn)［J］.紅外，2013，34(1):30－36.

［9］ 黃杰賢，楊冬濤，龔昌來.互信息熵與區(qū)域特征結(jié)合的圖像匹配研究［J］.激光與紅外，2013，43(1):98－103.

［10］曾文鋒，李樹山，王江安.基于仿射變換模型的圖像配準(zhǔn)中的平移、旋轉(zhuǎn)和縮放［J］.紅外與激光工程，2001，30(1):18－20.

［11］王金輝，陳冰，王建莊.實時圖像仿射變換系統(tǒng)的研究與實現(xiàn)［J］.自動控制與檢測，2012(2):59－62.

［12］ Mehfuza Holia，Thakar V K.Image Registration for Recovering Affine Transformation Using Nelder Mead Simplex Method for Optimization［J］.International Journal of Image Processing(IJIP)，2009，3(5):218－228.

［13］胡津翔，陳祺，王宏強.基于仿射變換的雙波段紅外圖像配準(zhǔn)算法研究［J］.科技廣場，2011，(1):103－106.

［14］張斌，周楊，陳書文，等.基于遺傳算法與仿射變換的遙感圖像配準(zhǔn)［J］.科學(xué)技術(shù)與工程，2010，10(21):5168－5170.

［15］ Noppadol Chumchob，Ke Chen.A Robust Affine Image Registration Methood［J］.International Journal of Numerical Analysis and Modeling，2008，6(2):311－334.

［16］ Tanner C，Carter T，Hawkes D J，et al.Cylindrical Affine Transformation Model for Image Registration［J］.Proc SPIE Medical Imaging:Image Processing，2010，7623(2):1－6.