陳玉，王欽軍

(中國科學(xué)院遙感與數(shù)字地球研究所，北京 100094)

0 引言

遙感影像分析的目的是反映、解釋、量化與描述地表模式，如何從遙感影像的空間與光譜特征出發(fā)構(gòu)建與現(xiàn)實(shí)地物相匹配的特征參量是遙感影像分析的關(guān)鍵[1-2]。多年來，地物的光譜特征參量(波形、特征吸收位置、特征吸收深度等[3-5])，空間特征參量(大小、形狀、紋理、位置、布局等[6-8])均已被有效應(yīng)用于遙感影像分析中，并獲得了成功[9-11]。但是，隨著遙感影像空間分辨率與光譜分辨率的逐步提高，數(shù)據(jù)挖掘的廣度與深度也在不斷增強(qiáng)。有些特定領(lǐng)域的應(yīng)用還需要通過進(jìn)一步挖掘目標(biāo)的空間與光譜信息，構(gòu)建在其所在環(huán)境背景干擾下保持一定不變性的特征參量。比如熱液或油氣滲漏的蝕變暈，點(diǎn)源污染導(dǎo)致的植被脅迫長勢分帶等。該類對象具有環(huán)形分帶或漸變的特征，可表達(dá)為沿著某一位置向周邊或某一側(cè)，光譜特征呈現(xiàn)出有規(guī)律的變化。該類目標(biāo)具有外源侵入的性質(zhì)，而被侵入體類型多樣，形成不同背景干擾著該類目標(biāo)的識別。由于該類目標(biāo)的光譜特征及形狀多變，因而難以用傳統(tǒng)的空間或光譜特征參量來表達(dá)。

1 空間光譜差比參量的提出

對于該類目標(biāo)對象，由于受到各種干擾影響，線性光譜混合模型有：

(1)

(2)

0≤fk i≤1

(3)

式中：Ri b為b波段第i像元的光譜反射率；fk i為對應(yīng)于第i像元的第k個組分所占比例；Ck b為第k個基本組分b波段的反射率；n為像元i所包含的基本組分?jǐn)?shù)目；εi b為誤差。為突出感興趣目標(biāo)地物，將公式(1)進(jìn)行變形有

(4)

Ri b-Rj b=f1i(C1i b-C1j b)

(5)

同樣對于其他波段b’有

Ri b′-Rj b′=f1i(C1i b′-C1j b′)

(6)

將公式(5)與公式(6)比有

(7)

公式(7)與地表干擾物的覆蓋組分無關(guān)，只與感興趣目標(biāo)的空間光譜變化有關(guān)，該參量稱為空間光譜差比參量Vijbb’，即不同波段空間相鄰像元反射率差值的比值。

2 空間光譜差比參量的模擬應(yīng)用

2.1 圖像模擬

以某類熱液蝕變?yōu)槔?，模擬存在如圖1(a)所示的熱液蝕變目標(biāo)，其地表蝕變暈呈三層分帶特征，各分帶光譜特征曲線如圖1(b)所示(假設(shè)有4個特征波段)。目標(biāo)物之外為取值0到1之間均勻分布的隨機(jī)像元。

圖1 某類模擬熱液蝕變目標(biāo)圖像及光譜特征

對于這一目標(biāo)對象，若不存在混合像元現(xiàn)象，可以采用多種方法直接識別出目標(biāo)地物。但是自然界地物環(huán)境復(fù)雜，受地質(zhì)環(huán)境影響，一些熱液蝕變目標(biāo)在地表可能局部缺失、形狀各異。此外，由于受到不同地表覆蓋類型及強(qiáng)度的影響，這些目標(biāo)的光譜特征會產(chǎn)生變異。為了簡化分析，我們以同一類型干擾物不同強(qiáng)度下影響為例，假定存在9個該類熱液蝕變目標(biāo)均受到某一光譜曲線如圖2所示地表覆蓋的影響。該覆蓋物在9個目標(biāo)中強(qiáng)度占比分別從0增加到80%。并且這9個目標(biāo)形狀各異還有部分有缺失現(xiàn)象。模擬的4個波段在干擾下各波段圖像如圖3所示。為了考察算法去噪性能在模擬圖像中加入了0～0.1之間的隨機(jī)噪聲。

圖2 模擬的干擾地物光譜曲線

注：各個圖像上目標(biāo)物從左到右，從上到下干擾依次占比以10%為梯度從0增加到80%。圖3 不同強(qiáng)度干擾下目標(biāo)地物各波段圖像

2.2 傳統(tǒng)分類方法

采用傳統(tǒng)的監(jiān)督分類的方法選擇圖1(b)的標(biāo)準(zhǔn)曲線作為參考光譜，并利用ENVI現(xiàn)有軟件分類功能，分別采用支持向量機(jī)法、光譜角填圖法、最大似然法、平行管道法、最小距離法、光譜信息散度法對目標(biāo)按照分類的方式進(jìn)行識別，結(jié)果如圖4所示。其中紅色、黃色、綠色分別代表識別出的目標(biāo)物的內(nèi)環(huán)、中環(huán)、外環(huán)。黑色表示不屬于目標(biāo)地物。

圖4 基于傳統(tǒng)各分類算法識別結(jié)果

2.3 利用空間光譜差比參量的方法

對于該類目標(biāo)對象，在不同強(qiáng)度或類型的背景干擾下其空間光譜差比參量保持不變。由公式(7)，對于這一目標(biāo)對象，若不存在混合像元現(xiàn)象，可以采用多種方法直接識別出目標(biāo)地物。但是自然界地物環(huán)境復(fù)雜，受地質(zhì)環(huán)境影響，一些熱液蝕變目標(biāo)在地表可能局部缺失、形狀各異。此外，由于受到不同地表覆蓋類型及強(qiáng)度的影響，這些目標(biāo)的光譜特征會產(chǎn)生變異。為了簡化分析，我們以同一類型干擾物不同強(qiáng)度下影響為例，假定存在9個該類熱液蝕變目標(biāo)均受到某一光譜曲線如圖2所示地表覆蓋的影響。該覆蓋物在9個目標(biāo)中強(qiáng)度占比分別從0增加到80%，并且這9個目標(biāo)形狀各異還有部分有缺失現(xiàn)象。模擬的4個波段在干擾下各波段圖像如圖3所示。為了考察算法去噪性能在模擬圖像中加入了0～0.1之間的隨機(jī)噪聲。

假設(shè)Vijbb′-(Rib-Rjb)/(Rib′-Rjb′)表示空間上相鄰i、j位置處光譜波段b、b’的空間光譜差比參量。以字母下標(biāo)n代表內(nèi)環(huán)，z代表中環(huán)，w代表外環(huán)，數(shù)字下標(biāo)1～4分別代表相應(yīng)的波段。Vnz12代表空間上內(nèi)環(huán)與中環(huán)，波段1與波段2之間的空間光譜差比參量。可知對于圖1所示的模擬圖像有：Vnz12=(Rn1-Rz1)/(Rn2-Rz2)=-1、Vzw12=2/3、Vnz23=-2、Vzw23=-3/2、Vnz34=-1、Vzw34=-2/3這些空間光譜差比參量在9個目標(biāo)地物中保持不變，與覆蓋物的強(qiáng)度與類型無關(guān)，可用于該類目標(biāo)地物識別。為此設(shè)計的算法流程如圖5所示。首先進(jìn)行圖像分割(分割過程中參數(shù)僅考慮光譜特性)，然后，分別對分割后的圖斑計算其與相鄰圖斑的空間光譜差比參量，在滿足給定誤差限的情況下輸出匹配結(jié)果。圖6為圖像分割結(jié)果(借助eCognition軟件實(shí)現(xiàn))，最終識別輸出結(jié)果如圖7所示(基于ArcGIS二次開發(fā)實(shí)現(xiàn))。

圖5 基于空間光譜差比參量的識別流程

圖6 圖像分割結(jié)果

圖7 基于空間光譜差比參量的識別結(jié)果

2.4 精度評價

對比圖4及圖7各方法識別結(jié)果，采用整體精度及Kappa系數(shù)評價各分類方法(表1)可以看出，傳統(tǒng)的分類方法在背景干擾下均會出現(xiàn)識別出錯的情況，由于模擬圖像的波段設(shè)置問題，這種錯誤在不同的算法中有不同程度的體現(xiàn)。而基于空間光譜差比參量的識別方法有效消除了背景干擾，整體精度達(dá)到99.98%，Kappa系數(shù)為0.998，可以很好地將9個目標(biāo)識別出來。

表1 各類方法分類整體精度及Kappa系數(shù)比較

3 結(jié)束語

本文通過理論分析與公式推導(dǎo)首次提出了空間光譜差比參量的概念，該參量可有效用于表達(dá)光譜特征隨空間呈現(xiàn)規(guī)律性變化的一類地物。基于模擬圖像，利用該參量構(gòu)建的識別方法總體識別精度達(dá)到了99.8%，大大高于傳統(tǒng)方法。本次實(shí)驗(yàn)方法是在圖像分割的基礎(chǔ)上，以圖斑為單位開展算法設(shè)計。實(shí)際上，也可不經(jīng)過圖像分割以像元代替圖斑開展應(yīng)用，但是算法的時間與空間復(fù)雜度均大大增強(qiáng)，且誤差限的設(shè)計需要更具有靈活性，這些有待后續(xù)進(jìn)一步的深入研究。此外，該參量也具有一定的局限性。首先，在公式推導(dǎo)過程中假設(shè)了臨近的像元受到的干擾類型及強(qiáng)度差異遠(yuǎn)小于距離較遠(yuǎn)的像元，該假設(shè)要求背景干擾在單個目標(biāo)物范圍內(nèi)應(yīng)該保持盡可能一致，對于不一致的情況則需要我們進(jìn)一步抽析出其中不變的特征參量。其次，該參量在推導(dǎo)過程中忽略了誤差項(xiàng)，誤差項(xiàng)對結(jié)果的影響隨著臨近圖斑光譜反射率差值的增大而減小，這就提醒我們在利用該參量時選擇光譜隨空間變化較大的特征波段組合效果更佳。