楊鵬飛，廖秀英，程 輝，帥慕容，謝貽文

(1.湖南科技大學(xué) 資源環(huán)境與安全工程學(xué)院，湖南 湘潭 411201；2.湖南科技大學(xué) 先進(jìn)礦山裝備教育部工程研究中心，湖南 湘潭 411201)

如何使不同類型的遙感信息在該領(lǐng)域發(fā)揮更大的作用成為遙感應(yīng)用重點(diǎn)探討的問題。通過遙感數(shù)據(jù)融合可以提高遙感數(shù)據(jù)的使用率、減少目標(biāo)信息識(shí)別的不確定性，為目標(biāo)信息的準(zhǔn)確提取奠定基礎(chǔ)[1-3]。近年來隨著對(duì)高分辨率遙感影像應(yīng)用研究的深入[4-6]，尤其是在地物識(shí)別和特征提取方面。針對(duì)高分辨率遙感影像分類[7-8]的各種智能算法是目前研究的熱點(diǎn)，而大多數(shù)智能算法都是基于影像地物光譜信息，因此保持光譜信息在遙感影像融合中顯得格外重要。高分辨率遙感數(shù)據(jù)的全色影像具有較好的空間特征信息，而多光譜影像則擁有較好的光譜信息,融合處理能夠提高后續(xù)影像的使用。目前對(duì)遙感影像數(shù)據(jù)融合評(píng)價(jià)的研究頗多,但大多數(shù)對(duì)遙感影像融合評(píng)價(jià)指標(biāo)都是基于融合前后影像的空間細(xì)節(jié)和信息量等方面進(jìn)行的，雖然這些指標(biāo)在一定程度上能夠反映不同融合方法對(duì)影像空間分辨率的保持能力，但是無法體現(xiàn)融合方法對(duì)地物光譜特征和光譜信息的保持能力。針對(duì)以上問題，以光譜信息散度[9-12]和光譜梯度角[13-15]作為評(píng)價(jià)光譜信息保持度的指標(biāo)，分析在高分辨率遙感影像融合中應(yīng)用最為廣泛的幾種融合方法，評(píng)價(jià)和比較不同融合方法對(duì)SPOT6影像地物光譜信息的保持能力，為其他同類型的高分辨遙感影像融合方法的選擇提供經(jīng)驗(yàn)借鑒。

1 融合方法

1.1 高通濾波融合

高通濾波融合法(HPF)采用空間高通濾波器去除全色圖像空間信號(hào)中低頻信號(hào)成分，提取全色影像空間信息特征，相應(yīng)的采用低通濾波器去除多光譜高頻成分，提取影像的光譜信息。最后經(jīng)過處理后的具有空間信息特征和光譜信息特征的圖像相加以提高新的融合影像。

1.2 小波變換融合

絕大多數(shù)傳統(tǒng)的融合方法都是基于圖像的時(shí)間域進(jìn)行運(yùn)算處理的，并沒有從圖像的頻率域進(jìn)行分析。小波變換融合的思想是將小波視作一種介于時(shí)間域和頻率域的函數(shù)表示。融合過程中首先是將全色影像和多光譜影像分別采用小波基進(jìn)行二維小波分解，從分解后的結(jié)果中提取小波系數(shù)，最后對(duì)提取小波系數(shù)采用逆變換達(dá)到影像融合的目的。

1.3 Gram-Schmidt變換融合

Gram-Schmidt融合是當(dāng)下一種比較新穎基于像元級(jí)別的融合，是對(duì)遙感影像采用正交化處理。其基本流程是首先利用原多光譜影像模擬全色影像，然后將模擬出的高分辨率全色影像作為Gram-Schmidt變換的初始分量，接著用高空間分辨率的影像代替Gram-Schmidt 變換后得到的第一分量，最后將替換的數(shù)據(jù)進(jìn)行反變換得到融合后圖像。

1.4 Pan-Sharpening變換融合

Pan-Sharpening是一種基于最小二乘逼近法來分別擬合全色灰度值、多光譜灰度值以及融合后的灰度值并計(jì)算多光譜影像和全色影像間灰度值的關(guān)系。這種方法使得大多數(shù)影像都能擬合全色波段數(shù)據(jù)，因此融合效果具有獨(dú)立性，同時(shí)在信息保留上更有優(yōu)勢(shì)。

以上4種融合方法都是基于像素級(jí)的融合，融合效果相較于傳統(tǒng)的融合方法效果較好并且一次性能對(duì)所有波段進(jìn)行融合，但以上融合方法原理不同，優(yōu)缺點(diǎn)也不盡相同，如表1所示。由于以上4種融合方法的優(yōu)缺點(diǎn)及效果不同，因此采用這4種融合方法對(duì)SPOT衛(wèi)星的多光譜影像和全色影像進(jìn)行融合對(duì)比實(shí)驗(yàn)。

表1 4種融合方法優(yōu)缺點(diǎn)對(duì)比

2 評(píng)價(jià)指標(biāo)

光譜梯度角是一種類似于光譜角的光譜區(qū)分方法，也是一種常用的評(píng)價(jià)光譜曲線形狀相似性的指標(biāo)。它是將兩條光譜曲線視為二維空間的向量，首先對(duì)兩條光譜曲線進(jìn)行一階求導(dǎo)，獲取其對(duì)應(yīng)的梯度向量，然后計(jì)算兩梯度向量之間廣義夾角。對(duì)于以上兩光譜曲線x,y,其梯度向量分別是

SG(x) = (x2-x1,x3-x2,

x4-x3,x5-x4,…xn-xn-1),

(1)

SG(y) = (y2-y1,y3-y2,y4-y3,

y5-y4…yn-yn-1).

(2)

兩梯度向量構(gòu)成的廣義夾角表示為

(3)

光譜梯度角作為評(píng)價(jià)指標(biāo)在反映光譜的局部特征變化上有較大優(yōu)勢(shì)，尤其是針對(duì)光譜曲線斜率的變化。

光譜信息散度主要是從信息學(xué)理論角度對(duì)光譜曲線加以區(qū)分。按照信息論的概念，自信息量I(x)是指某消息x出現(xiàn)后其自身所攜帶的信息量，一般用概率的負(fù)對(duì)數(shù)來計(jì)算，因此首先定義曲線x和曲線y的自信息(self-information)定義為

I(xi/yi)=-logp(xi/yi).

(4)

兩條曲線的相關(guān)熵按照定義為

(5)

則光譜信息散度為

SID(x，y)=D(x||y)+D(y‖x)=

(6)

光譜梯度角和光譜信息散度分別從兩個(gè)不同方面反映光譜特性。光譜梯度角類似于光譜角，從局部形狀特征的相似度上區(qū)分光譜曲線；光譜信息散度從整體上對(duì)光譜曲線比較，從信息論的角度予以區(qū)別。雖然光譜梯度角和光譜信息散度的側(cè)重點(diǎn)不同，但二者優(yōu)缺點(diǎn)互補(bǔ)。將光譜梯度角和光譜信息散度結(jié)合能較好地體現(xiàn)融合方法對(duì)原始影像光譜信息的保持能力。

3 實(shí)驗(yàn)結(jié)果與分析

3.1 研究區(qū)域

本次實(shí)驗(yàn)數(shù)據(jù)選擇的株洲市株洲縣部分地區(qū)SPOT6影像，數(shù)據(jù)軌道號(hào)是202～238，數(shù)據(jù)的接收日期為2014-01-23。影像參數(shù)如表2所示。

株洲縣屬于典型的亞熱帶季風(fēng)氣候，境內(nèi)多屬丘陵地帶，影像上的地物主要以耕地、林地、道路、水域、建設(shè)用地等，原始多光譜影像空間分辨率為6 m(如圖1所示)，原始全色影像空間分辨率為1.5 m(如圖2所示)。

圖1 原始多光譜影像

圖2 原始全色影像

SPOT6發(fā)射時(shí)間2012-09-09軌道高度/km695 軌道形式太陽同步軌道重訪周期/d2～3幅寬/km60 動(dòng)態(tài)范圍/bits12 空間分辨率全色/m1.5 多光譜/m6 波譜范圍/μm1全色0.45～0.75 2藍(lán)色0.45～0.52 3綠色0.53～0.59 4紅色0.62～0.69 5近紅外0.76～0.89

3.2 數(shù)據(jù)預(yù)處理

在對(duì)影像進(jìn)行融合處理前，首先要對(duì)影像預(yù)處理，包括影像的輻射校正、大氣校正、空間配準(zhǔn)等預(yù)處理。數(shù)據(jù)預(yù)處理流程如圖3所示。

圖3 數(shù)據(jù)處理流程圖

輻射校正目的是將影像上地物的像元值(DN值)轉(zhuǎn)為輻射亮度值，以便于分析地物的光譜曲線。大氣校正目的是進(jìn)一步消除外部因素如大氣、光照等對(duì)地物反射的影響，獲取地物真實(shí)的反射數(shù)據(jù)。影像經(jīng)過配準(zhǔn)之后進(jìn)行融合處理可以有效的降低多光譜影像和全色影像之間的像元誤差。

3.3 融合效果評(píng)價(jià)

3.3.1 主觀評(píng)價(jià)

圖4—圖7為SPOT6衛(wèi)星采用4種不同融合方法得到效果圖。從主觀的視覺上看4種融合方法明顯的提高影像的空間分辨率，色彩細(xì)節(jié)與原始多光譜影像較為接近。從對(duì)地物細(xì)節(jié)的保持能力看4種融合在影像紋理特性和細(xì)節(jié)特征上均得到不同程度的提高。

圖4 小波變換融

圖5 Gram-schmidit融合

圖6 高通變換融合

圖7 Pan-Sharping融合

3.3.2 客觀評(píng)價(jià)

為了獲取融合影像各類地物的光譜曲線特征，提高融合后影像在地物分類上的應(yīng)用。實(shí)驗(yàn)針對(duì)融合后影像與原始多光譜影像分別選取了林地、耕地、水體、建設(shè)用地、道路5種地物并獲得5類地物在不同波段上的光譜特征曲線，如圖8—圖12所示。

圖8 道路光譜曲線

圖9 耕地光譜曲線

圖10 建設(shè)用地光譜曲線

圖11 林地光譜曲線

圖12 水體光譜曲線

通過分析以上5種地物光譜曲線可以得出4種融合方法都較好地保持了地物的光譜曲線，融合影像的光譜曲線變化趨勢(shì)與原多光譜曲線變化趨勢(shì)基本一致。在林地及建設(shè)用地上4種融合方法光譜曲線的變化趨勢(shì)相當(dāng)一致，而在道路、人工建筑和水體的光譜曲線的變化趨勢(shì)基本保持一致。對(duì)地物光譜信息的保持度采用光譜信息梯度角和光譜信息散度定量評(píng)價(jià)。

3.3.3 光譜梯度角和光譜信息散度分析

根據(jù)5類地物的光譜曲線，結(jié)合光譜梯度角與光譜信息散度優(yōu)缺點(diǎn)互補(bǔ)的特點(diǎn)，對(duì)于4種融合方法對(duì)地物光譜信息的保持能力做進(jìn)一步的定量評(píng)價(jià)。選取5種地物類型，通過4種不同的融合方法得到的光譜曲線與原始多光譜影像光譜曲線進(jìn)行光譜梯度角與光譜信息散度的計(jì)算，綜合計(jì)算結(jié)果如表3所示。

通過表3統(tǒng)計(jì)的數(shù)據(jù)可以看出，對(duì)于道路，Gram-Schmidt變換融合后的光譜梯度角最小， 高通濾波光譜信息散度最?。粚?duì)于耕地、建設(shè)用地、林地，高通濾波變換融合后的光譜梯度角最小，且光譜信息散度最小；對(duì)于水體，小波變換融合后的光譜梯度角最小，高通濾波的光譜信息散度最小。

綜上所述，高通變換相對(duì)于其它三種融合方法而言能夠較好的保持原始多光譜影像的光譜信息，同時(shí)也能夠最大程度保持多光譜影像的光譜特性。針對(duì)于其他方法，小波變換的光譜信息保持能力優(yōu)于Gram-Schmidt變換，光譜信息保持效果最差的是Pan-Sharpening變換。

表3 4種融合方法的光譜梯度角和光譜信息散度

4 結(jié) 論

從光譜梯度角與光譜信息散度兩個(gè)角度,利用小波變換、高通濾波、Pan-Sharpening、Gram-Schmidt 4 種不同的融合方法，選取5種典型地物作為研究對(duì)象，對(duì)比融合前后影像同一地物的光譜曲線之間的光譜梯度角以及光譜信息散度數(shù)值，通過這兩個(gè)指標(biāo)對(duì)不同融合方法的光譜信息保持能力做了詳細(xì)評(píng)價(jià)。主觀和客觀的評(píng)價(jià)表明：小波變換、高通濾波、Pan-Sharpening融合和Gram-Schmidt4種融合方法都能有效提高影像空間分辨率，但在光譜信息的保持能力上，高通濾波變換最大程度地保持了原始多光譜影像的光譜信息，在對(duì)比實(shí)驗(yàn)研究過程中融合效果最好的方法，能較好地提高SPOT影像的應(yīng)用精度。