基于PCA變換和光譜補(bǔ)償?shù)倪b感影像融合方法

包 磊，徐其志

(1.北京航空航天大學(xué)計(jì)算機(jī)學(xué)院，北京100191;2.北京航空航天大學(xué)數(shù)字媒體北京市重點(diǎn)實(shí)驗(yàn)室，北京100191)

隨著遙感技術(shù)的不斷發(fā)展，人們獲取到大量同時(shí)相全色和多光譜遙感圖像。通常，在實(shí)際應(yīng)用中單一傳感器獲得的遙感圖像存在較大局限，例如多光譜圖像的分辨率低，而全色圖像中目標(biāo)和背景的色彩辨識(shí)度過低，二者分開處理時(shí)解譯人員和機(jī)器難以準(zhǔn)確地分析和判讀。為此，研究人員將全色和多光譜圖像的信息綜合在一起形成光譜信息豐富、空間細(xì)節(jié)清晰的融合圖像，降低圖像分析和解譯的不確定性。迄今為止，人們提出了基于矩陣變換［1－2］、基于比率變換［3］和基于多分辨分析［4－5］等多種圖像融合方法。

其中，PCA融合方法簡(jiǎn)單，可極大地提高圖像的空間分辨率，改善視覺效果。但是，該方法在提高圖像空間分辨率的同時(shí)，光譜往往嚴(yán)重失真。為解決這個(gè)問題，本文首先分析導(dǎo)致融合圖像失真的根本原因，并在主成分分析(Principal Component Analysis，PCA)變換的基礎(chǔ)上，結(jié)合小波變換的分頻和方向?yàn)V波特性，提出了一種基于PCA變換的遙感圖像高保真融合方法。同時(shí)，本文給出了一種新型融合流程，將低分辨率圖像逐層上采樣并與相應(yīng)的高分辨率圖像融合，獲得較好的視覺效果。最后，對(duì)IKONOS和QuickBird圖像進(jìn)行實(shí)驗(yàn)驗(yàn)證，結(jié)果表明本方法生成融合圖像的空間細(xì)節(jié)清晰，光譜保真度好。

1 PCA融合圖像的失真原因分析

PCA方法是一種數(shù)據(jù)壓縮和信息抽取的多元統(tǒng)計(jì)分析技術(shù)。在用PCA融合方法對(duì)低分辨率多光譜圖像與高分辨率圖像進(jìn)行融合時(shí)，首先對(duì)低分辨率多光譜圖像進(jìn)行PCA變換，獲得其成分變量;然后，對(duì)高分辨率圖像進(jìn)行線性拉伸，使之與低分辨率多光譜圖像第一主成分分量具有相同的均值與方差;最后，用拉伸后圖像替換第一主成分，再通過PCA反變換到RGB空間，即得到最終的融合圖像。

在融合過程中，光譜畸變的主要原因是傳感器的光譜響應(yīng)范圍不一致。理想情況下，各個(gè)多光譜波段應(yīng)近似是全色波段覆蓋范圍的一個(gè)劃分。但實(shí)際情況往往不是這樣，例如QuickBird和IKONOS衛(wèi)星的光譜響應(yīng)(圖來源于digitalglobe網(wǎng)站)，如圖1所示。在圖1中，全色與多光譜波段的覆蓋范圍存在一定的差異:在可見光高頻部分，藍(lán)波段大部分沒有落在全色波段內(nèi);在可見光低頻部分，全色波段的光譜響應(yīng)覆蓋范圍延伸至近紅外波段外側(cè)，甚至在波長(zhǎng)大于950μm的情況下，只有全色波段覆蓋。此外，多光譜圖像光譜曲線彼此間存在部分交叉，并不連續(xù)［6］。因此，在部分區(qū)域，多光譜圖像第一主成分與全色圖像二者的亮度存在較大差異，該差異是產(chǎn)生PCA融合失真的根源。

在傳統(tǒng)的PCA融合方法中，直接用全色波段替代占有多光譜圖像90%以上能量的主成分，但通過以上分析可知，全色波段與各多光譜亮度分量是一個(gè)復(fù)雜的非線性關(guān)系。二者直接替換后，在亮度存在較大差異的區(qū)域會(huì)產(chǎn)生嚴(yán)重的光譜失真。現(xiàn)有方法往往通過亮度匹配減小第一主成分與全色圖像的亮度差異，但效果較差，部分區(qū)域的亮度差異甚至變大，從而導(dǎo)致更大的光譜失真。

圖1 QuickBird與IKONOS光譜響應(yīng)示意圖Fig.1 QuickBird and IKONOS relative spectral response

小波變換方法提供了多尺度和多分辨率的分析手段，在圖像融合處理中發(fā)揮著重要作用。小波變換融合是將圖像分解成一系列不同分辨率、頻率、方向的子帶。通常，小波融合圖像存在較嚴(yán)重的“振鈴”模糊，因此，該方法往往與其他方法結(jié)合提高融合圖像質(zhì)量。

2 基于PCA變換的高保真融合方法

在充分分析PCA方法不足的基礎(chǔ)上，本文將從以下兩個(gè)方面對(duì)其進(jìn)行優(yōu)化:(1)選取PCA融合光譜失真的區(qū)域;(2)在光譜失真區(qū)域，引入小波融合避免圖像的光譜失真。

由分析可知，當(dāng)全色圖像與第一主成分二者局部均值之差大于閾值時(shí)，融合將產(chǎn)生嚴(yán)重的光譜失真。在光譜失真區(qū)域，本文采用對(duì)全色圖像和第一主成分進(jìn)行小波融合的方法，而非直接用全色圖像替換第一主成分。這樣既不損失圖像的細(xì)節(jié)信息，又盡可能的保留多光譜圖像原有的光譜信息。算法流程如圖2所示，具體步驟如下:

(1)對(duì)多光譜圖像進(jìn)行PCA變換，得到第一主成分P;

(2)分別對(duì)全色圖像和P進(jìn)行二級(jí)小波變換，其低頻子帶分別記為和;

(5)對(duì)P*替換第一主成分并進(jìn)行PCA反變換，得到融合圖像。

圖2 新型算法流程示意圖Fig.2 Flow chart of the new algorithm

3 彌補(bǔ)走樣缺陷的新型融合流程

在傳統(tǒng)的圖像融合過程中，通常是把低分辨率圖像重采樣到高分辨率圖像的尺寸，而由于不同傳感器圖像源的分辨率有較大的差異，經(jīng)常會(huì)出現(xiàn)走樣現(xiàn)象，即由于拉伸或采樣不充分重建后造成的信息失真，或多或少呈現(xiàn)鋸齒狀或臺(tái)階狀外觀，影響了融合影像的視覺效果。

針對(duì)這種情況，本文采用新型融合流程，若低分辨率圖像M尺寸為m，則高分辨率圖像P尺寸為4m，得到分辨率為4m的融合結(jié)果，如圖3所示。

圖3 多級(jí)融合流程圖Fig.3 Flow chart ofmulti-level fusion

圖4為傳統(tǒng)的和本文的融合流程的對(duì)比實(shí)驗(yàn)結(jié)果圖(放大32倍)，容易看出圖4(a)具有更光滑的邊緣。此外，這種流程具有廣泛的適用性，隨著圖像分辨率差異的增大，視覺效果的提升也更為顯著。

圖4 圖像融合效果對(duì)比Fig.4 Com parison of the different fusion results

4 實(shí)驗(yàn)結(jié)果及分析

由于篇幅的限制，本文僅展示兩組具有代表性的遙感圖像的實(shí)驗(yàn)結(jié)果，以論證本文方法的有效性和適應(yīng)性。第一組實(shí)驗(yàn)數(shù)據(jù)為QuickBird衛(wèi)星(簡(jiǎn)稱SQ)拍攝的寧波地區(qū)遙感圖像;第二組實(shí)驗(yàn)數(shù)據(jù)為IKONOS衛(wèi)星(簡(jiǎn)稱SI)拍攝的上海地區(qū)遙感圖像，如圖5所示。本文采用主觀和客觀相結(jié)合的方法對(duì)融合圖像的質(zhì)量進(jìn)行評(píng)價(jià)，其中，客觀評(píng)價(jià)指標(biāo)主要包括光譜扭曲度、圖像清晰度和圖像信息熵三個(gè)指標(biāo)。

觀察圖6可知:(1)PCA融合圖像的光譜失真較嚴(yán)重，出現(xiàn)了全局性的色偏(IKONOS衛(wèi)星的融合圖像色彩偏黃，而QuickBird衛(wèi)星的融合圖像色彩偏藍(lán));(2)小波融合圖像的細(xì)節(jié)較模糊;(3)本文方法的融合圖像無明顯的光譜失真，圖像的紋理細(xì)節(jié)清晰。總之，本文方法無論是光譜保持還是空間細(xì)節(jié)保持均較好。

圖5 QuickBird和IKONOS的全色和多光譜圖像Fig.5 Original image of QuickBird and IKONOS

圖6 圖像融合效果對(duì)比Fig.6 Comparison of the different fusion results

表1和表2分別列出了QuickBird和IKONOS融合圖像的客觀評(píng)價(jià)指標(biāo)。由表1，表2中的統(tǒng)計(jì)數(shù)據(jù)可知:本文融合方法的光譜扭曲度具有明顯的優(yōu)勢(shì)(光譜扭曲度越小，光譜失真越少);在空間細(xì)節(jié)保持方面，本文方法略勝于傳統(tǒng)PCA方法，且融合圖像的清晰度僅略低于全色圖像。由此可見，本文方法基本達(dá)到了高保真融合的目標(biāo)。由于沒有利用傳感器的光譜響應(yīng)特性，因此，本文方法的適應(yīng)性較好。

表1 QuickBird融合圖像客觀評(píng)價(jià)Table 1 The assessment of the QuickBird fusion images

表2 IKONOS融合圖像客觀評(píng)價(jià)Table 2 The assessment of the IKONOS fusion images

6 結(jié)束語

本文提出了一種基于PCA變換的遙感圖像高保真融合方法。它通過預(yù)先找到PCA融合方法的光譜失真區(qū)域，對(duì)失真區(qū)域采用PCA與小波變換相結(jié)合的方式進(jìn)行融合，對(duì)其他區(qū)域直接采用PCA方法進(jìn)行融合，從而確保融合圖像不失真。本文分別采用IKONOS和QuickBird遙感圖像進(jìn)行驗(yàn)證，實(shí)驗(yàn)表明:本文方法可以實(shí)現(xiàn)全色與多光譜圖像的高保真融合，且方法的適應(yīng)性好。

［1］Chibani Y，Houseline A.The joint use of HIS transform and redundantwave decomposition for fusingmultispectral and Panchromatic images［J］.International Journal Remote Sensing，2002，23(18):3821－3533.

［2］Kumar S S，Muttan S.PCA－based image fusion［C］// Algorithms and Technologies for Multispectral，Hyperspectral，and Ultraspectral Imagery XII.Proceedings of the SPIE，2006:V6233－62331T.

［3］Liu JG.Smoothing filter－based intensity modulation:a spectral preserve image fusion technique for improving spatial details［J］.International Journal of Remote Sensing，2003，V21(18):3461－3472.

［4］Yijian Pei，Jiang Yu.The improved wavelet transform based image fusion algorithm and the quality assessment［C］//2010 3rd International Congress on Image and Signal Processing，2010.

［5］Bhutada G G，Anand R S，Saxena SC.Edge preserved image enhancement using adaptive fusion of images denoised by wavelet and curvelet transform［J］.Digital Signal Processing，2011，V21(1):118－130.

［6］王明遠(yuǎn).空間對(duì)地觀測(cè)技術(shù)導(dǎo)論［M］.北京:軍事誼文出版社，2002:37－54.

［7］Yang S，Wang M，Jiao L.Fusion of multispectral and panchromatic images based on support value transform and adaptive principal componentanalysis［J］.Information Fusion，2012，13(3):177－184.