張?bào)汴?楊 桄 楊永波 黃俊華

(空軍航空大學(xué) 長春 130022)

高光譜圖像融合算法研究與進(jìn)展

張?bào)汴?楊 桄 楊永波 黃俊華

(空軍航空大學(xué) 長春 130022)

高光譜圖像融合是改善其圖像質(zhì)量的重要途徑。將高光譜圖像融合分為高光譜自身波段間的融合及其與高空間分辨率的其他遙感圖像如全色圖像等融合兩大類。從分析高光譜圖像融合特點(diǎn)入手,回顧了近年來出現(xiàn)的融合算法,闡述了其基本原理,進(jìn)行了歸納、分類,指出各自的優(yōu)缺點(diǎn),并總結(jié)了現(xiàn)行高光譜圖像融合算法的評價(jià)方法。最后,探討了高光譜圖像融合下一步的研究方向。

高光譜; 圖像融合; 融合算法評價(jià)

Class Number TP751

1 引言

高光譜遙感能夠獲取地物豐富的光譜信息,在目標(biāo)檢測、精細(xì)分類、人工解譯等領(lǐng)域有著廣闊的應(yīng)用前景[1]。然而,高光譜又存在嚴(yán)重的數(shù)據(jù)冗余問題,其信息量與波段數(shù)目并不成正比；此外,由于傳感器能量有限,高光譜分辨率的獲取往往以犧牲空間分辨率為代價(jià),圖像通常存在空間分辨率較低的問題。圖像融合是提升高光譜圖像質(zhì)量的重要方法。高光譜圖像融合與通常意義下的多源遙感圖像融合有著細(xì)微的區(qū)別,除了參與融合的圖像數(shù)量通常多于兩幅,融合可能與數(shù)據(jù)降維聯(lián)系在一起外,融合的結(jié)果也隨需求不同而呈現(xiàn)出多種形式,既可能是匯集了多幅圖像細(xì)節(jié)特征的灰度或合成彩色圖像,又可能是空間分辨率得到提升的高光譜圖像。

根據(jù)融合源圖像類型不同,高光譜圖像融合可以分為高光譜圖像波段間融合與高光譜圖像與他源圖像融合兩大類。高光譜圖像波段間融合,是指對高光譜圖像去除冗余、合并互補(bǔ)信息到單一圖像[2]的過程。而高光譜與他源圖像融合,包括高空間分辨率的全色圖像、多光譜圖像以及SAR圖像等多源遙感圖像,目的是整合幾類圖像的優(yōu)勢,改善高光譜圖像的空間特性。本文對目前這兩類高光譜圖像融合的算法研究進(jìn)行了回顧,總結(jié)了現(xiàn)有成果,概括分析了研究中所存在的問題,并對高光譜圖像融合算法的發(fā)展進(jìn)行了展望。

2 高光譜圖像自身波段融合算法

傳統(tǒng)的多源圖像融合,如紅外圖像與可見光圖像融合,針對的多是兩類差異較大、互補(bǔ)信息明顯的圖像,而高光譜圖像各波段圖像成像條件相同,波段間相關(guān)性強(qiáng),光譜特征差異不那么明顯,一方面融合前基本不需要配準(zhǔn),可以方便實(shí)施像素級融合,另一方面,又要求算法具有善于保留細(xì)節(jié)信息、多波段輸入等特性。這類融合有兩大目的:一是集中展示圖像信息,為目視解譯服務(wù)；二是一種有效的降維手段。根據(jù)融合過程是否需將圖像變換到頻率域進(jìn)行,融合算法可分為基于空間域以及基于變換域的兩大類。前者直接在空間域上進(jìn)行融合,后者則先將圖像從空域轉(zhuǎn)到頻域再進(jìn)行融合。

2.1 基于空間域的融合算法

基于空間域的融合算法主要有加權(quán)平均法、偏袒法、主成分分析法(Principal Component Analysis, PCA)等,其中PCA變換是處理相關(guān)性較大的高維數(shù)據(jù)的經(jīng)典方法之一[3～5]。其基本思想是利用幾個不相關(guān)、信息量大的綜合變量來概括原始變量的信息。

(1)

易知,通過線性變換得到的矩陣Z可以有多種。若要用z1來代替原來的p個波段的信息,以方差來表征信息量,則增加限制條件,使

1)z1方差Var(z1)最大；

2)z1與其余波段zi不相關(guān),即:

(2)

這樣可以唯一確定變換矩陣A,從而求出Z。滿足這些條件的z1稱為第一主成分；同理進(jìn)行迭代可以繼續(xù)求出第二、第三主成分。通常前幾個主成分可以包含原圖像85%以上的信息。

PCA變換可以認(rèn)為是原始高光譜圖像方差意義下的波段最優(yōu)融合方案。但是,這種最優(yōu)是針對全局的,并沒有考慮圖像局部感興趣信息,較為籠統(tǒng)。因此學(xué)者們作出各種改進(jìn)。例如針對PCA線性分析未必適用于高光譜圖像一些非線性特征的問題,Jia等提出分段主成分變換算法[6],將整個光譜空間分解成若干子空間,在各子空間中分別進(jìn)行PCA特征提取,使高光譜數(shù)據(jù)特性得到充分利用；G. A. Licciardi等則提出非線性PCA算法(Non-linear PCA, NLPCA)[7],快速降低數(shù)據(jù)維數(shù)的同時盡量避免原始PCA算法可能造成的光譜失真；為改進(jìn)原始PCA將整個波段圖像視為一維向量造成空間上下文信息的浪費(fèi),Ch. Theoharatos等將二維主成分分析(two-dimensional PCA,2DPCA)工具引入高光譜圖像融合[8],在不破壞空間信息的前提下使特征提取更為直觀、迅速；王秀朋等則聯(lián)合投影尋蹤(Projection Pursuit,PP)算法[9],將PCA在全局最優(yōu)分析特性與PP分析良好的局部特性相結(jié)合,改善了高光譜圖像的特征提取融合效果；朱院院等針對復(fù)雜背景數(shù)據(jù)容易導(dǎo)致PCA融合丟失弱小目標(biāo)的問題,選用歐式距離(Euclidean Distance, ED)、光譜角度(Spectral Angle, SA)、光譜信息散度(Spectral Information Divergence, SID)、正交投影散度(Orthogonal Projection Divergence, OPD)作為非監(jiān)督分類相似性測度對PCA進(jìn)行改進(jìn)[10],提高了融合圖像的對比度。

總體來說,PCA變換法原理簡單、運(yùn)算效率高,且為方差意義下的全局最優(yōu)方案,適用于大目標(biāo)、背景較均勻的高光譜圖像融合。此外,在融合應(yīng)用中PCA還常與其他方法結(jié)合使用。

2.2 基于多分辨率分析的融合算法

基于變換域的多分辨率分析的方法能夠更有針對性的改善圖像空間特性。Piella等提出的多分辨率分析框架[11]基本思想是采用某種多分辨率分析工具將融合源圖像分解為低頻與高頻子帶,采用合適的規(guī)則將兩類子帶分別融合,從而得到融合圖像的高頻與低頻子帶系數(shù),最后通過逆變換將其還原為空間域上的圖像。1983年Burt P.J和Adelson E.H提出的拉普拉斯金字塔變換[12]、90年代以來發(fā)展迅速的小波變換工具以及后來出現(xiàn)的多小波變換、脊波變換、曲波變換和輪廓波變換等各種多分辨率分析工具都被應(yīng)用到圖像融合領(lǐng)域,取得令人滿意的效果。

小波變換在處理圖像的過程中表現(xiàn)出良好的方向性、時頻分析特性和各尺度上的獨(dú)立性,使之成為圖像融合應(yīng)用最廣泛的多分辨率分析工具[13]。直接應(yīng)用小波變換加權(quán)平均的方法融合會導(dǎo)致融合圖像出現(xiàn)分塊現(xiàn)象,學(xué)者們紛紛對融合規(guī)則和小波工具進(jìn)行了改良。針對小波融合中依據(jù)單一特征較為局限的問題,張鈞萍等提出一種基于局部方差、局部信息熵、局部能量、局部梯度等多特征聯(lián)合指標(biāo)的確定低頻與高頻系數(shù)權(quán)值的融合規(guī)則[14],提高了融合圖像的分類精度；陳瀚孜[2]將第二代提升小波應(yīng)用到高光譜圖像融合,使小波函數(shù)的構(gòu)造更具靈活性,提高了運(yùn)算速度并使恢復(fù)圖像品質(zhì)更好；董廣軍等則使用多個小波函數(shù)構(gòu)成的小波包作為多分辨率分析工具[15],克服了傳統(tǒng)小波在處理高頻細(xì)節(jié)信息的不足。但是,小波變換反映的是信號的零維奇異性,在分析“點(diǎn)”的奇異性時是最優(yōu)的[16],而實(shí)際圖像地物輪廓更多的表現(xiàn)為“線”的奇異性,因此又有學(xué)者將更接近自然二維圖像的最佳稀疏表達(dá)的脊波、曲波、輪廓波等變換引入高光譜圖像融合。田養(yǎng)軍等利用曲波變換(Curvelet)作為多分辨率分析工具[17]研究了高光譜圖像融合在土地調(diào)查中的應(yīng)用；常戚戚等將非下采樣輪廓波變換(Nonsub Sampled Contourlet Transform,NSCT)應(yīng)用到高光譜特征圖像融合中[18]。與小波融合相比,這些多分辨率分析工具能夠更好地保持高光譜圖像的空間和光譜特性。

此外,將多分辨率分析與主成分分析相結(jié)合,綜合這兩類算法的特長,也是一種行之有效的融合方法。常戚戚、郭雷等提出了基于混合Contourlet變換和PCA變換的融合方法[19],首先對將選出的多波段高光譜圖像進(jìn)行Contourlet變換,利用PCA 變換對得到的系列多尺度多方向的子帶系數(shù)分別進(jìn)行自適應(yīng)融合處理,再通過Contourlet逆變換得到結(jié)果,該方法能夠有效抑制原始圖像中的噪聲；朱衛(wèi)東等對多波段圖像二代脊波(Bandelet)變換得到的系數(shù)和幾何流進(jìn)行PCA變換[20],得到其主成分后再經(jīng)過Bandelet逆變換重構(gòu)圖像,融合圖像性能也優(yōu)于直接使用這兩種方法得到的圖像。

2.3 其他方法

除了上述提到的兩大類方法,耦合脈沖神經(jīng)網(wǎng)絡(luò)(Pulse Coupled Neural Networks,PCNN)等也被成功引入到高光譜圖像融合中。文獻(xiàn)[21]采用多通道PCNN模型來對輸入的多波段圖像進(jìn)行非線性融合處理,通過記錄點(diǎn)火時刻的賦時矩陣等方法增強(qiáng)了融合結(jié)果；付朝陽則結(jié)合小波變換[22]、趙春暉等結(jié)合二代曲波變換[23],對低頻、高頻系數(shù)分別采取多通道PCNN模型按照不同融合規(guī)則得到融合結(jié)果,豐富了融合圖像的紋理等細(xì)節(jié)信息。但是,PCNN參數(shù)要經(jīng)過多次實(shí)驗(yàn)驗(yàn)證得到,不能自適應(yīng)確定,限制了算法的通用性。

3 高光譜與高空間分辨率圖像融合算法

除了自身融合,光譜信息豐富的高光譜圖像與高空間分辨率的其他類遙感圖像融合,各取所長,各補(bǔ)所需,也是研究的一大熱點(diǎn)。同時根據(jù)融合目的不同,融合結(jié)果可以是空間分辨率提高的高光譜圖像,也可以是提高了目視效果的灰度或假彩色圖像。

3.1 基于分量替換的融合算法

假設(shè)融合前的高光譜圖像為I1,高空間分辨率圖像為I2。經(jīng)過筆者歸納,發(fā)現(xiàn)許多算法的本質(zhì)思想是“替換”,即對兩類圖像進(jìn)行某種變換,然后以I2某分量去替換高光譜圖像經(jīng)變換后對應(yīng)的分量,從而將I2的優(yōu)良特性傳遞給I1,經(jīng)過逆變換便可得到融合結(jié)果。例如經(jīng)典的基于顏色空間變換的IHS方法[24],選取高光譜圖像的三個波段賦值給R、G、B三通道,然后通過以下公式進(jìn)行亮度與光譜信息分離:

(3)

(4)

(5)

其中I代表了高光譜圖像的亮度信息,以I2取代I再進(jìn)行逆變換,即可將I2的空間亮度等信息注入到融合圖像中。類似的基于顏色空間變換還有Brovey變換、HSV變換[25]等方法,但是,這些方法一次只能對三個波段進(jìn)行處理,且對融合圖像的光譜空間造成一定扭曲,因此在實(shí)際應(yīng)用中通常結(jié)合其他方法使用。如陸東華等利用高通濾波法來修正Brovey變換中使用的全色圖像[26]、楊可明等將HSV與小波包變換結(jié)合實(shí)施融合[27]。此外,還有基于PCA變換的方法[28],以I2代替高光譜PCA變換的第一主成分；小波變換方法,以I2小波分解后的高頻細(xì)節(jié)信息代替高光譜各波段圖像小波分解的高頻信息等；正交(Gram Schmidt)變換法[29],將高光譜圖像變換到正交空間,以I2取代其變換后的第一分量等。此外,楊可明等還引入電力系統(tǒng)中常用的諧波分析思路[30～31],提取高光譜影像諧波分解后的一系列諧波能量譜特征分量中對光譜波形無貢獻(xiàn)的諧波余項(xiàng),將其替換為I2,再經(jīng)過逆變換得到融合結(jié)果。

此外,英國學(xué)者Liu等提出了基于平滑濾波的亮度調(diào)節(jié)(SFIM)的融合方法,將圖像DN值表示為太陽輻射值與地表反射率的乘積,采用平滑濾波的方法模擬低分辨率圖像的DN值并進(jìn)行分量替換,這種方法能夠較好地保持光譜信息,但是會存在邊緣模糊的問題；寒冰等對此進(jìn)行了改進(jìn),結(jié)合瞬時視場角理論進(jìn)行DN值重新分配,提高了融合圖像的空間特性。

這些方法運(yùn)算簡單,能夠提高高光譜圖像的空間特性,但是不能很好地解決融合后光譜保持的問題,因此沒有得到廣泛認(rèn)可。

3.2 基于光譜解混模型的融合算法

除了替換的思想,基于混合光譜分析(Spectral Mixture Analysis,SMA)的方法也被應(yīng)用于高光譜圖像融合。這種方法的思路是在對高光譜圖像進(jìn)行光譜解混的基礎(chǔ)上,結(jié)合高空間分辨率圖像的空間信息,不在破壞光譜特征的前提下將分解端元“定位”以提高高光譜圖像空間分辨率。

非負(fù)矩陣分解(Nonnegative Matrix Factorization,NMF)是一種有效的光譜解混方法。它將高光譜數(shù)據(jù)視為矩陣,通過非負(fù)矩陣分解將其分解為端元光譜矩陣與豐度矩陣的乘積,然后經(jīng)過一定處理、重構(gòu)達(dá)到融合目的。崔艷榮等交替利用MNF方法對高光譜數(shù)據(jù)與高空間分辨率的多光譜數(shù)據(jù)進(jìn)行分解[32],得到高光譜分辨率的端元矩陣以及高空間分辨率的豐度矩陣,相乘還原為融合圖像；而安振宇等則引入光譜約束項(xiàng)對高光譜數(shù)據(jù)非負(fù)矩陣分解后的基底進(jìn)行增強(qiáng)[33],保證了光譜信息的保持。

Estiman等提出運(yùn)用極大后驗(yàn)估計(jì)(MAP)的方法,通過挖掘高光譜影像以及高空間分辨率影像之間的關(guān)系,基于隨機(jī)光譜解混模型(SMM)實(shí)現(xiàn)高空間分辨率高光譜數(shù)據(jù)的模型估計(jì),實(shí)現(xiàn)了高光譜影像整體空間信息的整體增強(qiáng)。豐明博等提出的基于投影和小波分解的多、高光譜遙感圖像融合算法[34]也反映了光譜解混的思想。首先確定純地物光譜,然后根據(jù)源圖像利用投影的方法得到混合地物的比例,基于多光譜圖像得到模擬的高光譜圖像,再通過小波融合方法將模擬與真實(shí)的高光譜圖像融合,一定程度上修正了融合結(jié)果的光譜特征；此外,Yeji Kim、Jaewan Choi等利用光譜解混技術(shù)來模擬擴(kuò)充多光譜圖像波段、進(jìn)而提取相關(guān)波段的高頻信息使之與高光譜圖像融合[35]。此外,還有學(xué)者提出了基于邊緣信息混合像元分解的光譜信息保持融合算法[36]。

基于光譜解混的高光譜圖像融合能夠很好的減少光譜失真現(xiàn)象,同時保持源圖像的光譜和空間分辨率,提高了高光譜圖像的質(zhì)量,最充分的體現(xiàn)了高光譜圖像融合的意義。但是這類算法對端元提取算法結(jié)果有較大依賴性,尤其是在地物光譜庫不完善、缺乏先驗(yàn)知識的情況下不能正確提取端元會容易造成端元定位錯誤。

3.3 其他算法

除了基于替換和光譜解混兩大思路,還有基于遺傳算法、基于粗糙集理論、基于光譜空間復(fù)原等方法[37～39]。表1簡單總結(jié)了這幾種算法的優(yōu)缺點(diǎn)。

4 融合算法評價(jià)

算法評價(jià)是算法研究的重要環(huán)節(jié)。除了評價(jià)融合算法自身的運(yùn)算復(fù)雜度、穩(wěn)定性等性能,對融合圖像質(zhì)量的評價(jià)也是判斷融合算法好壞的一項(xiàng)重要內(nèi)容。圖像質(zhì)量評價(jià)可分為主觀評價(jià)和客觀評價(jià)兩大類。前者根據(jù)人眼對融合圖像的視覺感受對色彩(灰度)、紋理信息、清晰度、對比度等作出判斷[40],具有簡單直觀的特點(diǎn),但是易受個人主觀影響,且不能判斷光譜信息的保持；客觀評價(jià)則通過一系列指標(biāo)對融合圖像進(jìn)行量化分析。對高光譜融合圖像質(zhì)量的客觀評價(jià)要從其幾何空間特性以及光譜空間特性入手[41]。

通常評價(jià)單幅圖像空間特性的指標(biāo)[21,27,35]包括均值、標(biāo)準(zhǔn)差、平均梯度、信息熵等,適用于融合圖像單一波段的評價(jià)。這些指標(biāo)越大說明融合圖像對比度越高,細(xì)節(jié)信息越豐富。但是,由于高光譜圖像會存在噪聲污染,細(xì)節(jié)紋理信息與噪聲都會提高圖像的標(biāo)準(zhǔn)差、平均梯度等指標(biāo),在應(yīng)用時也要注意結(jié)合圖像實(shí)際情況。在評價(jià)高光譜與高空間分辨率圖像融合算法時還涉及光譜特征保持度指標(biāo)。這類指標(biāo)包括圖像相關(guān)系數(shù)[32](融合圖像各波段與原波段像元光譜相關(guān)性)、波形保持度[31](融合圖像像元波譜曲線與原光譜波形形態(tài)的相似程度)、波譜角誤差[8](融合圖像像元光譜與真實(shí)光譜之間的角度誤差)等。光譜保持指標(biāo)多散見于各文獻(xiàn),缺乏普適性的系統(tǒng)研究。此外,融合圖像分類精度等應(yīng)用性指標(biāo)也常用于對融合算法質(zhì)量評價(jià)。

5 結(jié)語

本文總結(jié)了目前國內(nèi)外高光譜圖像融合算法的研究進(jìn)展,對各類算法進(jìn)行了初步的分類和歸納。總體而言,高光譜波段間融合注重保留各波段的細(xì)節(jié)信息,因此能夠進(jìn)行針對性融合的多分辨率分析方法更受青睞；而對于高光譜與高分辨率圖像融合,目前許多算法都由通用遙感圖像融合方法基礎(chǔ)上發(fā)展而來,而能夠保持光譜特征的算法正在成為研究重點(diǎn)。筆者認(rèn)為今后應(yīng)在以下方面作進(jìn)一步研究:

1) 對于高光譜自身波段間的融合,新提出的曲波、脊波等多分辨率分析工具雖然提升了融合圖像質(zhì)量,但是運(yùn)算復(fù)雜度較高,對數(shù)據(jù)量大的高光譜圖像來說更是增加了融合時間。因此,改善多尺度分析工具、合理快速選擇融合波段都是值得研究的問題；

2) 高空間分辨率的遙感圖像可以輔助高光譜圖像進(jìn)行校正、去噪、混合像元分解等處理,全面提高高光譜圖像質(zhì)量。因此,在研究融合算法時不能僅僅局限于提高其空間分辨率、避免光譜扭曲,要充分利用兩類圖像的信息,挖掘融合的價(jià)值；

3) 目前對高光譜融合算法的評價(jià)多是沿用傳統(tǒng)多源圖像融合中的方法,缺乏廣泛認(rèn)可的能夠全面反映高光譜融合數(shù)據(jù)特性的評價(jià)指標(biāo)。建立合理有效的高光譜圖像融合評價(jià)體系也是下一步值得關(guān)注的問題；

4) 目前對高光譜圖像融合更多的停留在理論層面,通過融合改善了空間特性的高光譜圖像沒有得到很好的應(yīng)用檢驗(yàn)。下一步研究中應(yīng)加強(qiáng)理論實(shí)踐相結(jié)合,推動高光譜遙感的應(yīng)用。

[1] 張兵,高連如.高光譜遙感圖像分類與目標(biāo)探測[M].北京:科學(xué)出版社,2011.

[2] 陳瀚孜.高光譜圖像融合算法研究[D].哈爾濱:哈爾濱工程大學(xué)碩士學(xué)位論文,2008:4-37.

[3] Wen Dou, Yun Hao,Chen, et al. A general framework for component substitution image fusion: An implementation using the image fusionmethod[J]. Comput ers and G eosciences,2007,33(2):219-228.

[4] Zheng Du, Myong K J. Band selection of hyperspectral images for automatic detection of poult ry skin tumors[J]. IEEE Trans. on Automation Science and Engineering,2007,4(1):332-339.

[5] Ch. Theoharatos, V. Tsagaris, N. Fragoulis. Hyperspectral image fusion using 2-D principal component analysis[C]//2ndInternational Conference on Space Technology. IEEE,2011:1-4.

[6] Jia Xiu-Ping, Richards John A. Segmented principal components transformation for efficient hyperspectral remote- sensing image display and classification[J]. IEEE Trans. On G. R. S.,1999,37(1):538-542.

[7] G. A. Licciardi, M. M. Khan, J. Chanussot. Fusion of hyperspectral and panchromatic images: a hybrid use of indusion and nonlinear PCA[C]//ICIP IEEE,2012:2133-2137.

[8] 吳學(xué)明,楊武年.基于2DPCA的遙感圖像融合算法研究[J].地質(zhì)科技情報(bào),2008,27(5):97-103.

[9] 王秀朋,張洪才,趙永強(qiáng),等.基于投影尋蹤的高光譜圖像特征融合算法[J].計(jì)算機(jī)測量與控制,2006,14(11):1539-1541.

[10] 朱院院,高教波,高澤東,等.基于相似性分類的高光譜主成分融合方法比較[J].光學(xué)技術(shù),2016,42(2):97-103.

[11] G. Piella. A general framework for multiresolution image fusion: From pixels to regions[J]. Information Fusion,2003,4(4):259-280.

[12] Burt P J, Adelson E H. The Laplacian pyramid as a compact image code[J]. IEEE Transaction on Communication,1983,31(4):532-540.

[13] Tania Stathaki. Image Fusion: Algorithms and Applications[M]. Elsevier SingaporePte Ltd.,2014.

[14] 張鈞萍,張曄.基于多特征多分辨率融合的高光譜圖像分類[J].紅外與毫米波學(xué)報(bào),2004,23(5):345-348.

[15] 董廣軍,張永生,范永弘.PHI高光譜數(shù)據(jù)和高空間分辨率遙感圖像融合技術(shù)研究[J].紅外與毫米波學(xué)報(bào),2006,25(2):123-127.

[16] 亓辰.高光譜與高空間分辨率遙感圖像融合算法研究[D].哈爾濱:哈爾濱工業(yè)大學(xué)碩士學(xué)位論文,2008.

[17] 田養(yǎng)軍,薛春紀(jì),馬智民,等.曲波變換的高光譜遙感圖像融合方法在土地利用調(diào)查中的應(yīng)用[J].遙感學(xué)報(bào),2009,13(2):313-320.

[18] 常威威,郭雷,付朝陽.基于非采樣Contourlet變換的高光譜圖像融合算法[J].火力與指揮控制,2009,34(10):121-125.

[19] 常威威,郭雷,付朝陽,等.基于混合Contourlet變換和主成分變換的高光譜圖像融合[J].宇航學(xué)報(bào),2009,30(6):2360-2366.

[20] 朱衛(wèi)東,李全海,徐克科,等.基于二代Bandelet和主成分變換的高光譜遙感圖像融合[J].同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,39(7):1068-1074.

[21] 常威威,郭雷,付朝陽,等.利用脈沖耦合神經(jīng)網(wǎng)絡(luò)的高光譜多波段圖像融合方法[J].紅外與毫米波學(xué)報(bào),2010,29(3):205-210.

[22] 付朝陽,郭雷,常威威.基于小波變換和多通道脈沖耦合神經(jīng)網(wǎng)絡(luò)的高光譜圖像融合[J].吉林大學(xué)學(xué)報(bào)(工學(xué)版),2011,41(3):838-844.

[23] 趙春暉,季亞新.基于二代曲波變換和PCNN的高光譜圖像融合[J].哈爾濱工程大學(xué)學(xué)報(bào),2008,29(7):729-735.

[24] 王克成.基于小波理論的超光譜圖像融合研究[D].哈爾濱:哈爾濱工程大學(xué)碩士學(xué)位論文,2005:34-58.

[25] Kim Y, Eo Y, Kim Y, et al. Generalized IHS-Based Satellite Imagery Fusion Using Spectral Response Functions[J]. Etri Journal,2011,33(4):497-505.

[26] 陸冬華,趙英俊.基于改進(jìn)型Brovey算法的高光譜數(shù)據(jù)融合技術(shù)[J].世界核地質(zhì)科學(xué),2006,23(3):177-182.

[27] 楊可明,錢小麗,王林偉,等.結(jié)合HSV和小波包變換的高光譜影像融合及信息識別研究[J].測繪通報(bào),2014(5):11-15.

[28] Hassan Kousha, Hassan Gessemian. High Spatial Resolution and Hyperspectral Image Fusion[C]//19thIranian Conference on Electrical Engineering. IEEE,2011.

[29] 于海洋,閆柏琨,甘甫平.基于Gram Schmidt變換的高光譜遙感圖像改進(jìn)融合方法[J].地理與地理信息科學(xué),2007,23(5):39-43.

[30] 楊可明,張濤,王立博.高光譜影像的諧波分析融合算法研究[J].中國礦業(yè)大學(xué)學(xué)報(bào),2014,43(3):547-554.

[31] 楊可明,張濤,王立博,等.諧波分析法高光譜影像融合及其光譜信息保真度評價(jià)[J].光譜學(xué)與光譜分析,2013,23(9):2496-2502.

[32] 崔艷榮,何彬彬,張瑛,等.非負(fù)矩陣分解融合高光譜和多光譜數(shù)據(jù)[J].遙感技術(shù)與應(yīng)用,2015,30(1):82-91.

[33] 安振宇,史振威.基于sc-NMF的高光譜圖像融合[J].紅外與激光工程,2013,42(10):2718-2724.

[34] 豐明博,劉學(xué),趙冬.多/高光譜遙感圖像的投影和小波融合算法[J].測繪學(xué)報(bào),2014,42(2):158-163.

[35] Yeji Kim, Jaewan Choi, Dongyeob Han. Block-Based Fusion Algorithm With Simulated Band Generation for Hyperspectral and Multispectral Images of Partially Different Wavelength Ranges[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing,2015,8(6):2997-3008.

[36] 余旭初,等.高光譜遙感應(yīng)用分析[M].北京:科學(xué)出版社,2013.

[37] 馬一薇.高光譜遙感圖像融合技術(shù)與質(zhì)量評價(jià)方法研究[D].鄭州:信息工程大學(xué)碩士學(xué)位論文,2010.

[38] 董廣軍.高光譜與高空間分辨率遙感信息融合技術(shù)研究[D].鄭州:信息工程大學(xué)碩士學(xué)位論文,2004.

[39] 邵軍勇,潘泉.基于光譜復(fù)原的高光譜高空間遙感數(shù)據(jù)融合研究[J].微電子學(xué)與計(jì)算機(jī),2005,22(8):108-111.

[40] 張良,邵琳.圖像融合在高光譜遙感數(shù)據(jù)處理中的應(yīng)用[J].計(jì)算機(jī)與數(shù)字工程,2010,244(2):118-122.

[41] 滿旺,陳紹杰.高空間分辨率與高光譜分辨率遙感數(shù)據(jù)融合研究[J].測繪,2010,33(6):243-247.

Research Progress of Hyperspectral Image Fusion Algorithms

ZHANG Xiaohan YANG Guang YANG Yongbo HUANG Junhua

(Aviation University of Air Force, Changchun 130022)

Hyperspectral images share rich spectral information while it also has the problem of data redundancy and low spatial resolution. Image fusion is an effective way to improve its quality. According to the resource images of fusion, hyperspectral fusion can be classified into fusion among bands of hyperspectral data and fusion with images having higher spatial resolution such as panchromatic image. In this paper, main fusion methods of hyperspectral image are reviewed. The basic theory of those algorithms is analyzed and their advantages and disadvantages are also discussed. Besides, current fusion algorithm evaluation including indexes and methods is also summarized. At last,the direction of future research of this field is given.

hyperspectral imagery, image fusion, fusion algorithm evaluation

2016年7月11日,

2016年8月28日

吉林省科技發(fā)展計(jì)劃資助項(xiàng)目(編號:20140101213JC);吉林省教育廳“十二五”科研項(xiàng)目(編號:2015448)資助。

張?bào)汴?女,碩士,研究方向:高光譜遙感、圖像解譯。楊桄,男,博士后,教授,研究方向:遙感圖像解譯、地理信息系統(tǒng)。

TP751

10.3969/j.issn.1672-9730.2017.01.007