基于光譜重構(gòu)的三步高光譜圖像融合方法
【專利摘要】基于光譜重構(gòu)的三步高光譜圖像融合方法，即通過波段選擇、光譜重構(gòu)的處理方法對主成分分析圖像融合進行完善，得到更好的融合效果。該方法如下：1）通過人機交互接口模塊獲得待融合的高光譜圖像數(shù)據(jù)和高空間分辨率圖像數(shù)據(jù)，并初始化相關(guān)參數(shù)；2）通過高光譜圖像波段選擇模塊，設計波段選擇矩陣，對初始高光譜圖像進行波段選擇；3）通過高光譜圖像-高空間分辨率圖像融合模塊，對輸入的圖像進行主成分分析融合；4）通過光譜重構(gòu)模塊，對初始融合圖像進行高精度光譜重構(gòu)；5）通過融合結(jié)果輸出模塊，輸出最終的融合結(jié)果。本發(fā)明的優(yōu)點在于減少融合數(shù)據(jù)量，降低光譜畸變，提升光譜保持性，增強圖像質(zhì)量，該方法融合質(zhì)量好，使用范圍廣。
【專利說明】基于光譜重構(gòu)的三步高光譜圖像融合方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及一種基于光譜重構(gòu)的三步高光譜圖像融合方法，適用于高光譜圖像處理系統(tǒng)，屬于高光譜遙感圖像處理領(lǐng)域。
【背景技術(shù)】
[0002]隨著20世紀80年代初期成像光譜概念的出現(xiàn)，光學遙感進入了 一個嶄新的階
段-高光譜遙感(Hyperspectral Remote Sensing)階段。高光譜遙感指集探測器技術(shù)、
精密光學機械、微弱信號檢測、計算機技術(shù)、信息處理技術(shù)于一體的，能夠獲取具有高的光譜分辨率圖像的遙感科學和技術(shù)，是當前模式識別及遙感信息處理研究領(lǐng)域中的一個研究執(zhí)占。
[0003]與常規(guī)遙感圖像相比，高光譜遙感圖像具有“圖譜合一”的突出特點，既能表達空間信息，又能表達光譜信息，以三維數(shù)據(jù)立方體的形式呈現(xiàn)。在遙感應用中，有時會要求圖像同時具有高的空間分辨率和高的光譜分辨率。然而，由于當前技術(shù)條件的限制，現(xiàn)有高光譜成像儀器很難提供這樣的數(shù)據(jù)。高光譜圖像擁有較高的光譜分辨率，可達10nm，甚至更高，但其空間分辨率較低，即空間的細節(jié)表現(xiàn)能力較差；而普通遙感圖像(如全色圖像等)具有相對較高的空間分辨率，但不能表達光譜信息。因此，將高光譜遙感圖像和對應的具有高的空間分辨率的遙感圖像進行融合，可優(yōu)化圖像質(zhì)量，得到新的高光譜圖像，既具有較高的空間細節(jié)表現(xiàn)能力，同時又保留圖像的光譜特性。
[0004]目前很多傳統(tǒng)多源遙感數(shù)據(jù)融合方法只適用于普通遙感圖像間的融合以及多光譜圖像與全色數(shù)據(jù)的融合，由于高光譜圖像的巨大數(shù)據(jù)量，使得傳統(tǒng)融合方法效率低、難適用。在圖像融合領(lǐng)域，使用多源遙感數(shù)據(jù)進行融合的方法主要包括直接加權(quán)法、主成分分析(Principal Component Analysis, PCA)法、高通濾波(High-Pass Filter, HPF)算法和小波變換(Wavelet Translation, WT)法等。盡管前三種方法提高了高光譜數(shù)據(jù)各波段圖像的空間分辨率，但融合后的圖像光譜失真嚴重，容易增加噪聲，甚至影響視覺效果；而小波變換法保持了原始光譜特性，但空間分辨率效果較差。因此有必要采用新的方式解決高光譜圖像融合問題。
[0005]PCA是一種統(tǒng)計相關(guān)分析技術(shù)，在多數(shù)實際問題中，由于指標較多及指標間有一定的相關(guān)性，勢必增加分析問題的復雜性，PCA就是設法將原有指標重新組合成一組新的相互無關(guān)的幾個綜合指標來代替原來指標，是研究如何通過少數(shù)幾個主成分來解釋多變量的協(xié)方差結(jié)構(gòu)的一種方法。PCA變換后的新子空間中，圖像間的相關(guān)性基本消除，其中的第一主分量包含了絕大部分的信息。
[0006]壓縮感知是當前信號重構(gòu)研究的熱門理論，該理論指出，在一定條件下，可用少量壓縮信號高精度重構(gòu)原始信號，由此衍生了一系列信號重構(gòu)方法，如正交匹配追蹤(Orthogonal Matching Pursuit, OMP)算法等,利用最小二乘求解的貪婪方法快速求解原始信號。
[0007]本發(fā)明將壓縮感知波段重構(gòu)技術(shù)引入融合過程進行圖像優(yōu)化，首先對高光譜圖像進行波段選擇，再進行高光譜圖像與高空間分辨率圖像的PCA融合，最后再對其進行波段維的重構(gòu)過程，得到同時包含高空間分辨率空間信息和光譜信息的圖像，形成一套新的高光譜圖像融合方法體系。

【發(fā)明內(nèi)容】

[0008]本發(fā)明的目的在于提供一種基于光譜重構(gòu)的三步高光譜圖像融合方法，發(fā)明了一種運用壓縮感知光譜重構(gòu)技術(shù)進行波段選擇，從而以少量數(shù)據(jù)進行快速、精確融合，最終再高質(zhì)量重構(gòu)光譜維信息的多源圖像融合方法。本方法既提高融合過程的效率，又增強了高光譜圖像融合的光譜保持性，融合圖像性能好，使用范圍廣。
[0009]本發(fā)明所涉及的方法流程包括以下步驟:(1)獲得待融合的高光譜圖像和高空間分辨率圖像數(shù)據(jù)，初始化相關(guān)參數(shù)；(2)利用波段選擇矩陣對高光譜圖像進行波段選擇，以減少融合數(shù)據(jù)量；(3)計算高光譜圖像的主成分，并用高空間分辨率圖像與其進行PCA融合；(4)通過壓縮感知光譜重構(gòu)方法進行光譜的高精度重構(gòu)；(5)輸出最終的高光譜圖像融合結(jié)果。
[0010]下面對該方法流程各步驟進行詳細說明:
[0011](1)獲得待融合的高光譜圖像和高空間分辨率圖像數(shù)據(jù)，初始化相關(guān)參數(shù):待融合數(shù)據(jù)為波段數(shù)為n的高光譜圖像X= (X1, X2,, Xn)T和高空間分辨率圖像X，Xi (i=l, 2，…，n)表示第i波段的圖像，設置合適的波段采樣比r。
[0012](2)利用波段選擇矩陣對高光譜圖像進行波段選擇:按照合理的方式和采樣比對高光譜圖像進行波段選擇，通常選擇等間隔波段選擇，采樣比為大于I的正整數(shù)，即r=n/m=2, 3，…，其中m為選擇后的高光譜圖像波段數(shù)，通過波段選擇矩陣O對高光譜圖像進行波段選擇，即得到s=Ox。
[0013](3)計算高光譜圖像的主成分，并用高空間分辨率圖像與其進行PCA融合:對采樣后的高光譜圖像s進行PCA變換，得到第一主成分PC;，將X與PC;作直方圖匹配得到PCnew,用PCnew代替/5C1 ):與其余主成分一起做PCA逆變換，得到初始融合圖像S。
[0014](4)通過壓縮感知光譜重構(gòu)方法進行光譜的高精度重構(gòu):
[0015]
【權(quán)利要求】
1.基于光譜重構(gòu)的三步高光譜圖像融合方法，其特征在于:通過波段選擇和重構(gòu)提高融合過程的效率，并增強光譜保持性、穩(wěn)定性，具有人機交互接口模塊、高光譜圖像波段選擇模塊、高光譜圖像-高空間分辨率圖像融合模塊、光譜重構(gòu)模塊、融合結(jié)果輸出模塊這五個功能模塊。包括如下步驟: (1)獲得待融合的高光譜圖像數(shù)據(jù)和高空間分辨率圖像數(shù)據(jù)，初始化相關(guān)參數(shù):假設X為待融合的高空間分辨率圖像，X= (X1, X2, , xn)T為待融合的低空間分辨率的高光譜圖像，其中n是高光譜數(shù)據(jù)的波段總數(shù)，Xi, i=l,…，n代表第i波段的圖像響應數(shù)值，完成初始參數(shù)設置步驟； (2)設置合適的采樣比r，根據(jù)設置的采樣比和起始波段設計得到波段選擇矩陣O的表達，通過計算式S=Ox獲得降維后的高光譜圖像S，完成波段選擇步驟； (3)利用降維的高光譜圖像s與空間分辨率較高的遙感圖像X進行PCA融合:對采樣后的高光譜圖像s進行PCA變換，得到第一主成分，將高空間分辨率圖像X與PC: ?作直方圖匹配得到PCnew,用PCnew代替尸C1后與其余主成分一起做PCA逆變換，得到初始融合圖像S ； (4)基于壓縮感知光譜重構(gòu)方法對圖像S進行重構(gòu)，根據(jù)初始參數(shù)設計正交的稀疏變換基屯，通過S=OX’=O Wy求解圖像的稀疏表示y，并經(jīng)過正交反變換得到重構(gòu)圖像
【文檔編號】G06T5/50GK103714572SQ201310656056
【公開日】2014年4月9日 申請日期:2013年12月6日 優(yōu)先權(quán)日:2013年12月6日
【發(fā)明者】尹繼豪, 孫建穎, 姜志國, 余萬科, 朱紅梅 申請人:北京航空航天大學