李現(xiàn)虎，呂京國，江珊

(北京建筑大學(xué) 測繪與城市空間信息學(xué)院，北京 102616)

0 引言

影像融合是都空間信息和光譜信息的優(yōu)化處理，其可以獲得同時包含有高空間分辨率和高光譜分辨率的影像，為各個領(lǐng)域提供具有更高針對性和準(zhǔn)確性信息的研究影像[1]。根據(jù)影像融合特點可以將影像融合算法分為3個層次：像素級融合、決策級融合和特征級融合。其中像素級影像融合是目前最常用的影像融合方法，其獲取的影像信息豐富、性能效果最佳、適用影像種類最廣。本文研究也是基于像素級影像融合進(jìn)行研究[2]。但由于其在影像融合過程中需要對影像的所有像素信息進(jìn)行操作，所以隨著遙感影像空間分辨率的增加，其數(shù)據(jù)量也越來越大，融合過程對于計算機(jī)性能的要求也在增加，所以影像融合時間和效果問題也越來越突出。

近年來發(fā)展起來的壓縮感知(CS)理論提出，其可以在影像采樣的同時對影像進(jìn)行壓縮，并且不會丟失影像的有用信息，達(dá)到完全重構(gòu)原始影像的目的。Alin Achim[3]、 Atul Diveker[4]等一些學(xué)者結(jié)合壓縮感知理論能夠?qū)τ跋駢嚎s采樣的特點，提出了基于壓縮感知理論的影像融合算法，通過對影像壓縮采樣后的測量值進(jìn)行融合減少了計算量，提高了影像融合效率[5]，在一定程度上解決了大影像數(shù)據(jù)融合的問題。但由于影像中不同影像塊具有不同的信息復(fù)雜度，如果采用統(tǒng)一的融合系數(shù)進(jìn)行計算，無法在融合過程中對于每一塊影像都達(dá)到較為理想的效果。針對該問題，Lu Gan[6]、楊森林[7]等人提出了分塊壓縮感知理論(BCS)，通過分析影像的局部信息特征將影像分為不同的影像塊，對具有不同信息復(fù)雜度的影像塊采用不同的權(quán)值進(jìn)行融合，提高了影像整體的融合效果。但該過程中，影像融合系數(shù)的選擇不能隨不同的影像自動調(diào)整，所以該算法不具有自適應(yīng)性。

本文以分塊壓縮感知影像融合理論為基礎(chǔ)，利用粒子群優(yōu)化算法全局尋優(yōu)的特點，提出基于粒子群優(yōu)化算法的分塊壓縮感知影像融合算法。首先利用分塊壓縮感知理論對待融合影像進(jìn)行分塊壓縮采樣；然后通過粒子群優(yōu)化算法對不同影像塊的壓縮測量值求取最優(yōu)融合系數(shù)，并以該系數(shù)為權(quán)值進(jìn)行加權(quán)融合，得到融合后的測量值；最后對融合后的測量值進(jìn)行影像重構(gòu)，得到融合影像。該算法在一定程度上解決了分塊壓縮感知影像融合過程融合系數(shù)不能自適應(yīng)的問題。

1 相關(guān)理論

1.1 傳統(tǒng)壓縮感知理論

壓縮感知理論最早是在2006年由Candes和Donoho提出的[8]。之后在現(xiàn)有的信號調(diào)和分析、時頻分析、統(tǒng)計概率理論、矩陣分析、泛函分析、統(tǒng)計概率理論等理論的基礎(chǔ)上得到了迅速發(fā)展[9]。壓縮感知是一種新的在采樣的同時實現(xiàn)壓縮目的的理論框架，其壓縮采樣過程如圖1所示。

圖1 壓縮感知理論框架

壓縮感知從全局出發(fā)探究信號的本質(zhì)結(jié)構(gòu)，擺脫了信號頻率等物理量度之間的聯(lián)系[10]。只要信號具有某個可壓縮的稀疏域，那么就可以利用與變換矩陣非相干的測量矩陣將變換系數(shù)線性投影為低維觀測向量，從而把信號投影到低維空間。投影空間采樣值中包含足夠的信息來重構(gòu)信號，使用稀疏最優(yōu)化理論就能夠精確地、高概率地從少量的采樣值中重構(gòu)原始信號。壓縮感知理論主要包括信號的稀疏表示、測量矩陣的設(shè)計與重構(gòu)算法3個部分，其優(yōu)點就是針對可稀疏的信號，將傳統(tǒng)的數(shù)據(jù)釆集與數(shù)據(jù)壓縮合二為一，在獲取信號同時對數(shù)據(jù)進(jìn)行壓縮，且大大降低了傳統(tǒng)信號獲取和處理流程中的潛在消耗。

1.2 分塊壓縮感知理論

然而，在實際影像壓縮處理過程中，圖像中不同位置的信息量即紋理信息往往不同，如果采用統(tǒng)一的采樣率或統(tǒng)一的融合規(guī)則，往往達(dá)不到最為理想的效果。根據(jù)BCS 理論[6]，一幅Ir×Ic大小的如下，劃分為一個個B×B大小的小塊，對其中每一塊采用相同的測量矩陣φB進(jìn)行CS 測量。

yi=φBxi

(1)

式中：yi為CS 測量值矩陣；xi為第i塊圖像經(jīng)矢量化后的信號；φB為測量矩陣(φB∈RnB×B2)；nB=Ir×Ic，為總像素數(shù)，n為要獲得的壓縮感知測量值的個數(shù)。 其是一個局部壓縮采樣，相當(dāng)于對整幅圖像采用具有如下的分塊對角結(jié)構(gòu)的測量矩陣進(jìn)行全局壓縮采樣：

由于BCS 中只需要存儲一個nB×B2的測量矩陣，而不是全局測量的n×N測量矩陣，因此，BCS 具有更高的存儲效率。

另外，由于影像在不同分塊中具有不同的信息復(fù)雜度特征，在加權(quán)影像融合中具有不同的最優(yōu)融合系數(shù)，如果采用同一個融合系數(shù)對影像進(jìn)行加權(quán)融合，則不能夠使影像的不同信息復(fù)雜度位置得到最優(yōu)的融合效果。本文采用分塊壓縮感知具有2個方面的優(yōu)勢：采用局部壓縮采樣的方法提高了存儲效率；針對不同分塊影像在融合中選用不同的最優(yōu)融合系數(shù)，從而使影像不同部位都能得到最優(yōu)的融合效果。

1.3 粒子群優(yōu)化算法

粒子群優(yōu)化算法PSO (particle swarm optimization)是一種進(jìn)化計算方法(evolutionary computation)。 其基本思想是通過群體中個體之間的協(xié)作和信息共享來尋找最優(yōu)解[11]。 在該算法中，鳥被抽象為沒有體積和質(zhì)量的點(粒子)，并將其延伸到N維空間，粒子i在N維空間的位置矢量表示為Xi=(x1，x2，L，xN)，飛行速度的矢量表示為Vi=(v1，v2，L，vN)。每個粒子都有一個由目標(biāo)函數(shù)決定的適應(yīng)值(fitness value)，并知道其到已經(jīng)發(fā)現(xiàn)的最好位置(pbest)和目前所在位置Xi。此過程可以看作粒子自身的飛行經(jīng)驗。另外，每個粒子還知道整個群體當(dāng)前所發(fā)現(xiàn)的最好位置gbest(gbest是pbest中的最好值)，此過程可以看作群體粒子的飛行經(jīng)驗。粒子通過自身與群體中同伴的最好的經(jīng)驗來決定下一步的運(yùn)動方向與速度。在該過程中，粒子群優(yōu)化算法(PSO)首先被初始化為一群隨機(jī)粒子(隨機(jī)解)，然后通過迭代尋找最優(yōu)解。在每一次的迭代中，粒子通過兩個極值(pbest，gbest)來更新自己的位置與速度。

Vi=Vi+c1×rand()×(pbest-Xi)+
c2×rand()×(gbest-Xi)

(2)

Xi=Xi+Vi

(3)

式中：i=1，2，3L，N；N是此群中粒子的總數(shù)；Vi是粒子的速度；rand()是介于(0，1)之間的隨機(jī)數(shù)；Xi是粒子的當(dāng)前位置；c1和c2是學(xué)習(xí)因子，通常取值為c1=c2=2。

PSO 的算法流程如下：

(1)初始化一群粒子(粒子群規(guī)模為N)，包括隨機(jī)位置和速度；

(2)評價每個粒子的適應(yīng)值；

(3)將每個粒子的當(dāng)前適應(yīng)值與經(jīng)過的最好位置pbest作比較，如果其當(dāng)前適應(yīng)值較好，則將其當(dāng)前適應(yīng)值作為最好位置pbest；

(4)將每個粒子的當(dāng)前適應(yīng)值與群體粒子經(jīng)過的最好位置gbest作比較，如果其當(dāng)前適應(yīng)值較好，則將其作為當(dāng)前群體最好值gbest；

(5)根據(jù)公式(2)、公式(3) 調(diào)整粒子的位置和速度；

(6)判斷是否達(dá)到結(jié)束條件，若達(dá)到則結(jié)束，否則轉(zhuǎn)向步驟(2)。

迭代終止條件根據(jù)具體問題進(jìn)行選擇。一般選為最大迭代次數(shù)Gk，或粒子群目前為止尋找到的滿足預(yù)定最小適應(yīng)閾值的最優(yōu)位置[12]，粒子群優(yōu)化算法流程如圖2所示。

圖2 粒子群優(yōu)化算法流程圖

2 基于粒子群優(yōu)化的分塊壓縮感知影像融合算法

本文綜合利用了分塊壓縮感知融合算法存儲需求小、融合決策過程簡單、有利于大數(shù)據(jù)量遙感影像融合[9]的優(yōu)點，同時也應(yīng)用了粒子群優(yōu)化算法全局尋優(yōu)的特點計算分塊壓縮感知融合過程中每塊影像融合時的最優(yōu)融合系數(shù)，提出了基于粒子群優(yōu)化的分塊壓縮感知影像融合算法。此算法利用粒子群優(yōu)化算法確定不同分塊影像融合權(quán)重系數(shù)ω(其中ω為m×n矩陣，待融合影像分塊維度為m×n)，避免了傳統(tǒng)經(jīng)驗值不能隨不同影像的不同特點自適應(yīng)調(diào)整的問題。該算法中采用粒子群優(yōu)化算法選取融合參數(shù)ω(ω1，ω2，Lωm×n)，在粒子群優(yōu)化算法迭代過程中，將影像融合系數(shù)ωi(i=1，2，3，L，m×n)作為粒子，將信息熵和平均梯度函數(shù)作為優(yōu)化目標(biāo)函數(shù)，目標(biāo)是使融合的測量值能夠使信息熵和平均梯度值最大，從而得到限定優(yōu)化目標(biāo)函數(shù)之下的最優(yōu)值。將該最優(yōu)值作為目標(biāo)融合權(quán)重系數(shù)，然后在分塊壓縮感知影像融合算法中采用加權(quán)融合規(guī)則進(jìn)行影像融合。該算法具備一定的自適應(yīng)性，從而保證了融合效果的最優(yōu)化選擇?；诹Ｗ尤簝?yōu)化的分塊壓縮感知影像融合算法流程如圖3所示。

圖3 基于粒子群優(yōu)化的分塊壓縮感知影像融合算法流程圖

具體影像融合步驟描述如下：

(1)原始影像預(yù)處理。將全色波段影像與多光譜影像進(jìn)行空間配準(zhǔn)，得到空間信息高度匹配的兩幅影像，并根據(jù)具體研究區(qū)域進(jìn)行拼接、裁剪等工作；

(2)分別對全色波段影像與多光譜影像進(jìn)行分塊；

(4)將圖像測量值進(jìn)行融合：y′=ωy1+(1-ω)y，ω∈(0，1)；

(5)采用粒子群優(yōu)化算法選取融合參數(shù)ω(ω1，ω2，Lωm×n)。在粒子群優(yōu)化算法迭代過程中，將影像融合系數(shù)ωi(i=1，2，3，L，m×n)作為粒子，將信息熵和平均梯度函數(shù)作為優(yōu)化目標(biāo)函數(shù)，目標(biāo)是使融合的測量值的信息熵和平均梯度值最大；

(6)將上述得到的新的壓縮信息y′通過OMP(正交匹配追蹤)算法進(jìn)行壓縮感知的影像重構(gòu)，得到影像分別與對應(yīng)的每塊多光譜影像進(jìn)行直方圖匹配，得到最終的融合影像。

本文在影像分塊的過程中采取了均勻分塊的方法，但由于不同分塊數(shù)對影像融合結(jié)果也會產(chǎn)生影響。所以實驗先對分塊壓縮感知影像融合在不同分塊數(shù)時進(jìn)行了融合結(jié)果分析(圖4)，可以發(fā)現(xiàn)融合結(jié)果隨分塊數(shù)的增加，影像融合評價指標(biāo)中AG、EN、PSNR、STD 都在增加，指標(biāo)DD 在減小，表示融合效果越來越好。但隨著分塊數(shù)的增加，其指標(biāo)變化越來越慢，當(dāng)分塊數(shù)為6×6時，AG、EN、PSNR 和STD 均達(dá)到最大值，且當(dāng)分塊數(shù)大于6×6時，上述各指標(biāo)數(shù)開始有變小的趨勢，DD 減小的幅度越來越小，并且影像融合時間急速增加。綜合各種指標(biāo)分析，可以發(fā)現(xiàn)在該實驗中，當(dāng)分塊數(shù)為6×6時，影像融合效果最好。

注：各指標(biāo)數(shù)值為多組實驗結(jié)果的平均值，并且各組實驗結(jié)果單獨分析，與上述結(jié)果也較為相似。圖4 BCS影像融合算法在不同分塊時的融合結(jié)果評價分析

3 影像融合實驗結(jié)果評價與分析

本文所用數(shù)據(jù)為512像素×512像素的WorldView影像(全色波段分辨率為0.5 m，多光譜波段分辨率為1.8 m)、512像素×512像素的高分2號影像(全色波段分辨率為0.81 m，多光譜波段分辨率為3.24 m)、512×512像素的Landsat-8影像(全色波段分辨率為15 m，多光譜波段分辨率為30 m)，其中多光譜波段采用標(biāo)準(zhǔn)假彩色組合對上述算法進(jìn)行試驗驗證，并與傳統(tǒng)的IHS、PCA、CS、IHS-CS影像融合算法結(jié)果進(jìn)行比較。不同算法融合結(jié)果如圖5所示。

通過對圖5實驗數(shù)據(jù)中不同算法結(jié)果影像分析，可以發(fā)現(xiàn)各算法融合影像的空間分辨率均得到了明顯改善，但不同算法的融合效果之間又有差異；另外，不同算法融合的影像在光譜保持方面也具有不同的效果。其中，CS 融合算法與本文影像融合算法在光譜保持方面效果較好，但在空間分辨率提升方面，CS 融合結(jié)果較本文算法融合結(jié)果較差，如圖6(a4)、圖6(b4) 黃色框圖標(biāo)出的部分，其紋理顯示方面較圖6(a3)、圖6(b3) 明顯更為清晰。另外，IHS變換在影像空間分辨率的提高方面也具有較好的優(yōu)勢，但其采用了替換融合的思想，所以導(dǎo)致其融合結(jié)果的光譜保持方面的效果較差，如圖6(a2)、圖6(b2) 與原始影像圖6(a1)、圖6(b1) 光譜的對比發(fā)現(xiàn)，IHS變換融合影像光譜較暗，而CS影像融合圖6(a3)、圖6(b3) 與本文影像融合算法的結(jié)果圖6(a4)、圖6(b4)與原始多光譜影像對比，其光譜效果視覺效果上變化不大。綜合光譜特征的保持及空間分辨率的提升效果發(fā)現(xiàn)，本文算法具有較其他方法具有更大的優(yōu)勢，其融合效果改善較為明顯。

圖5 不同融合方法的結(jié)果圖像注：第1、2組為WorldView影像；第3組為高分2號影像；第4組為Landsat-8影像。

圖6 不同融合方法的結(jié)果圖像細(xì)節(jié)效果圖注：(a1)、(b1) 原始多光譜影像；(a2)、(b2) IHS融合影像；(a3)、(b3) CS融合影像；(a4)、(b4) 本文方法融合影像。

算法AGENDDPSNRSTDTime/sPCA6.801 811.130 155.976 821.546 350.500 03.083 046IHS8.862 313.148 853.358 421.931 657.065 90.088 041CS6.239 612.580 233.601 020.751 957.410 80.006 999IHS-CS8.826 910.589 053.210 320.795 056.728 60.011 347本文算法9.371 113.827 744.344 325.035 057.770 40.009 711

表1給出了兩組數(shù)據(jù)各種算法影像融合結(jié)果的平均梯度(average gradient，AG)、圖像信息熵(entropy，EN)、光譜扭曲度(degree of distortion，DD)、峰值信噪比(peak signal to noise ratio，PSNR)、標(biāo)準(zhǔn)差(standard deviation，STD)的平均值。從表1實驗數(shù)據(jù)中可以看出，本文算法的結(jié)果影像AG、EN、PSNR、STD值均為各算法結(jié)果影像的最大值，其指標(biāo)參數(shù)為最優(yōu)；DD值較CS算法結(jié)果較大，光譜扭曲度比CS算法稍差，但較其他幾種算法結(jié)果值較小，其光譜扭曲度也比其他幾種算法小，融合效果較好。而在計算的數(shù)據(jù)量方面，對比不同算法在影像融合算法運(yùn)行時間發(fā)現(xiàn)，本文算法通過對采樣后的少數(shù)測量值進(jìn)行融合減少了影像融合過程的計算量，所以其運(yùn)行所需時間比傳統(tǒng)影像融合算法明顯減少，但由于融合過程需要粒子群優(yōu)化獲取最優(yōu)融合系數(shù)，所以其需要時間較CS影像融合算法稍微較多，但其影像融合結(jié)果比CS 算法要好。

4 結(jié)論

本文用粒子群優(yōu)化算法對分塊壓縮感知影像融合算法進(jìn)行改進(jìn)。結(jié)果發(fā)現(xiàn)通過對采樣后的少數(shù)測量值進(jìn)行融合減少了影像融合過程的計算量，在基于粒子群優(yōu)化算法的遙感影像融合中，通過將融合過程中的一系列融合權(quán)重系數(shù)ωi(i=1，2，3，…，m×n)作為粒子群，并在其中尋找最優(yōu)值來確定分塊壓縮感知影像融合中最優(yōu)的融合權(quán)重系數(shù)ωi(ω1，ω2，ω3…，ωm×n)，彌補(bǔ)了傳統(tǒng)經(jīng)驗值非自適應(yīng)的缺點。另外也可以通過控制步長，使計算結(jié)果達(dá)到最優(yōu)，從而在一定程度上平衡了融合影像空間分辨率的提高與光譜特征的保持能力，改善了影像融合效果。

5 結(jié)束語

本文基于分塊粒子群優(yōu)化的壓縮感知影像融合算法在主觀評價和客觀指標(biāo)上都表現(xiàn)出了較大優(yōu)勢，但在光譜扭曲度方面較CS影像融合效果略差。另外由于目前沒有確定比較好的影像分塊數(shù)方法，本文通過一系列實驗分別分析了原始影像在分為不同行列數(shù)的實驗結(jié)果，發(fā)現(xiàn)在分塊數(shù)為6行6列時效果較好，而對于怎么更好地確定影像分塊數(shù)將是下一步研究的地方。