陳潔， 杜磊， 李京， 韓亞超， 高子弘

(1.中國國土資源航空物探遙感中心，北京 100083； 2.中國科學(xué)院遙感與數(shù)字地球研究所，北京 100101)

基于噪聲白化的高光譜數(shù)據(jù)子空間維數(shù)算法

陳潔1,2， 杜磊1， 李京1， 韓亞超1， 高子弘1

(1.中國國土資源航空物探遙感中心，北京 100083； 2.中國科學(xué)院遙感與數(shù)字地球研究所，北京 100101)

高光譜影像數(shù)據(jù)的相鄰波段間相關(guān)性較強(qiáng)，信號與噪聲共存，根據(jù)最小二乘原理，使觀測數(shù)據(jù)與噪聲的投影誤差之和最小化的HySime (hyperspectral signal identification by minimum error)算法，通過數(shù)據(jù)觀測值減去噪聲估計值后得到信號的估計值，進(jìn)而可以計算信號相關(guān)矩陣的估計值。該算法在準(zhǔn)確估計噪聲的情況下是可行的，但實際上經(jīng)光譜降維去相關(guān)后得到的各像元噪聲估計值往往并不準(zhǔn)確，因此，原始的HySime算法得到的結(jié)果可能并不理想。提出一種基于噪聲白化的HySime改進(jìn)算法，它不必進(jìn)行逐像元的噪聲去除，而是先對原始數(shù)據(jù)進(jìn)行噪聲白化處理，然后準(zhǔn)確獲取噪聲的協(xié)方差矩陣估計值，再利用HySime算法進(jìn)行信號相關(guān)矩陣計算，實現(xiàn)了提高算法精度的目的。通過模擬和實驗數(shù)據(jù)的驗證，改進(jìn)的算法結(jié)果更準(zhǔn)確穩(wěn)定，與經(jīng)典的NSP (noise subspace projection)算法在不同情況下所得結(jié)果有很好的一致性，通過引入噪聲白化的過程，提高了算法對非白噪聲的適應(yīng)性。

高光譜降維; HySime; 噪聲白化; 信號相關(guān)矩陣

0 引言

以小于20 nm光譜分辨率進(jìn)行數(shù)據(jù)采集的高光譜數(shù)據(jù)的相鄰波段間通常具有較強(qiáng)的相關(guān)性，它們的DN值以及顯示出來的視覺效果往往很相似。若把高光譜圖像表達(dá)成一個矩陣，矩陣的每一列代表其對應(yīng)波段的光譜響應(yīng)數(shù)據(jù)值，可以通過對該矩陣進(jìn)行一系列運算，使得變化后的矩陣前若干列波段數(shù)據(jù)就可以表示原圖像的大部分信息，而后幾列的數(shù)據(jù)則主要以噪聲為主。通過這樣的方法可以用低維數(shù)的數(shù)據(jù)代表高維數(shù)據(jù)，減小了數(shù)據(jù)處理量的同時，也抑制了噪聲的影響。數(shù)據(jù)降維結(jié)果的好壞將直接影響到下一步端元提取算法的最終效果，在數(shù)據(jù)降維的過程中，如何選擇合適的去相關(guān)方法和怎樣準(zhǔn)確確定所需低維數(shù)據(jù)的維數(shù)，是2個重要的關(guān)鍵步驟。

經(jīng)過多年的研究，高光譜數(shù)據(jù)在去相關(guān)方面已有相對成熟和穩(wěn)定的方法，并已得到廣泛應(yīng)用。去相關(guān)變換，就是去除高光譜數(shù)據(jù)波段間相關(guān)性的變換，即變換后數(shù)據(jù)波段處于不同的維度空間，代表數(shù)據(jù)的矩陣列相互正交。主成分變換(principal components analysis, PCA)和最小噪聲分離(minimumnoise fractions, MNF)是當(dāng)今使用最頻繁的兩種去相關(guān)算法。PCA變換就是通過一種基于K-L的線性變化，它將原始圖像矩陣中的變量變成一組不相關(guān)的隨機(jī)變量，并按照協(xié)方大小順序排列，組成新的圖像矩陣。一般情況下該矩陣第一個主成分含有所有波段中80%的方差信息。但PCA變換對噪聲比較敏感，即信息量大的主成分分量，信噪比不一定高，當(dāng)某個信息量大的主成分中包含的噪聲的方差大于信號的方差時，該主成分分量形成的圖像質(zhì)量就差。因此，Andrew等[1]提出了MNF變換，該方法本質(zhì)上是兩次PCA的重疊運算，變換后的數(shù)據(jù)向量相互正交，第一分量集中了主要的信息，與PCA變化不同的是向量按照信噪比大小排列，從而克服了噪聲對影像質(zhì)量的影響。隨后，James等[2]提出了NAPC變換方法(noise-adjusted principal component, NAPC)，一般認(rèn)為MNF和NAPC算法等價，稱呼起來也一般將兩者等同[3]。實際應(yīng)用中，往往按照NAPC方法的思路來實現(xiàn)。

但數(shù)據(jù)維數(shù)的確定方法還處在研究階段。PCA和MNF變換方法，在變換之后都有同一個問題需要解決，即“選取多少個主分量”，換句話說，就是低維數(shù)據(jù)究竟選多少維來代表高維數(shù)據(jù)才合適。實際應(yīng)用中需要研究者根據(jù)個人經(jīng)驗進(jìn)行選擇，這使得維數(shù)的確定有很強(qiáng)的主觀性和盲目性，而維數(shù)選取的正確與否，將直接影響后續(xù)處理的精度和數(shù)據(jù)分析，因此，選用適當(dāng)?shù)臏?zhǔn)則和方法確定最優(yōu)的維數(shù)，對接下來的數(shù)據(jù)處理和分析有著重要意義?；谏鲜隹紤]，本文分析和梳理了常用的噪聲白化處理方式和子空間維數(shù)確定方法，提出了一種以噪聲白化為基礎(chǔ)的高光譜數(shù)據(jù)子空間維數(shù)計算方法，然后利用該算法對模擬和實際數(shù)據(jù)進(jìn)行維數(shù)計算，最后通過與不同的算法的對比來檢驗其穩(wěn)定性和準(zhǔn)確性。

1 噪聲估算與白化

在各種因素的影響下，獲取影像信息的同時不可避免地采集了噪聲，而圖像噪聲的影響是不容忽視的，為此在對影像質(zhì)量進(jìn)行評估、解譯、分析和應(yīng)用之前，都需要進(jìn)行噪聲估算。目前，遙感影像噪聲估計主要有3種方法： 試驗室法、暗電流法和圖像法。試驗室法和暗電流法都需要對圖像有一定的先驗知識，需要一系列復(fù)雜的精確測量，在實際應(yīng)用中實現(xiàn)較難； 圖像法則利用遙感影像直接進(jìn)行分析來估算噪聲，因此該方法被廣泛應(yīng)用于圖像的噪聲估算[4]。針對高光譜影像，基于圖像法進(jìn)行噪聲估算，主要利用了高光譜圖像的相同類別和性質(zhì)像元的空間域相關(guān)性和高光譜數(shù)據(jù)波段間的波段域相關(guān)性。

1.1 噪聲估算

通過上述分析可知，高光譜數(shù)據(jù)相鄰波段間往往具有很強(qiáng)的相關(guān)性，噪聲估算的基本思想就以此為基礎(chǔ)： 多元線性回歸算法的運用一方面可以消除波段間的相關(guān)性，另一方面可以將計算得到的波段殘差作為噪聲的估計。目前常用的方法有Roger提出的全圖光譜維去相關(guān)法[5]，即對整個圖像進(jìn)行逐像元的噪聲估計，計算噪聲協(xié)方差矩陣并予以去除，認(rèn)為其結(jié)果即為信號真實值。和Roger及Arnold提出的空間、光譜維去相關(guān)法(spectral and spatial decorrelation method, SSDC)[6]類似，SSDC方法同樣是在進(jìn)行信號去相關(guān)處理的同時，計算信號殘差值，并以此作為噪聲估計進(jìn)行去除，其優(yōu)勢在于它利用高光譜圖像像元在相鄰空間域和光譜域的相關(guān)性，受地物覆蓋類型影響小，自動化程度高，無需人工干預(yù)，是目前較穩(wěn)定的高光譜圖像噪聲評估方法之一[7]。SSDC方法實現(xiàn)步驟如下[8]： ①首先需要對整個高光譜圖像數(shù)據(jù)進(jìn)行分塊處理。數(shù)據(jù)分塊要求數(shù)據(jù)塊連續(xù)、均勻且相互不重疊。理論上分塊越小則均一性越高、數(shù)據(jù)越獨立，但為了使分塊數(shù)據(jù)的回歸分析結(jié)果具有一定的統(tǒng)計意義，數(shù)據(jù)塊的大小一般推薦為15行×15列； ②求取對各數(shù)據(jù)塊進(jìn)行多元線性回歸運算的回歸參數(shù)，并計算噪聲方差； ③針對高光譜影像的任一波段，取所有均勻子塊方差的均值作為該波段圖像的噪聲方差的最優(yōu)估計。實際應(yīng)用中，由于分塊得到的子塊并不總是均勻的，因此可以只選取其中合理的子集進(jìn)行方差檢驗。子集的選取可以通過對方差進(jìn)行排序，然后前后各去掉15%來實現(xiàn)。

1.2 噪聲白化

噪聲白化是指對噪聲數(shù)據(jù)進(jìn)行去相關(guān)，且噪聲協(xié)方差被單位化的過程。若a為m×1維隨機(jī)向量序列(均值=0)，則其協(xié)方差矩陣可表示為一個m×m維的矩陣Ca，如果該協(xié)方差矩陣Ca非奇異，且不為單位矩陣，則稱隨機(jī)向量a為有色(或非白)隨機(jī)向量。令協(xié)方差矩陣特征值分解為

Ca=VDVT,

(1)

且

(2)

記變換

(3)

則b的協(xié)方差矩陣為單位矩陣I，即隨機(jī)向量b是標(biāo)準(zhǔn)白色隨機(jī)向量，換言之，有色的隨機(jī)向量a經(jīng)過線性變換Wa后，變成了白色隨機(jī)向量。一般將線性變換矩陣W稱為隨機(jī)向量a的白化矩陣。

通過研究發(fā)現(xiàn)，對于某一隨機(jī)變量，能使其進(jìn)行白化的變換矩陣并不是唯一存在的[9]，例如上述分析的W矩陣和通常MNF變換中使用的白化矩陣

(4)

都可以達(dá)到使a白化的效果。進(jìn)一步的研究顯示，當(dāng)采用不同的白化矩陣，所得的白化結(jié)果b并不一樣，以變換矩陣W和FT為例進(jìn)行如下實驗說明：

通過以上實驗和結(jié)果可知，雖然變換矩陣W和FT都可以使隨機(jī)變量a“白化”，但白化矩陣W的變換結(jié)果并不改變變量a原有的分布特征，而白化矩陣FT則不論變量a原有的分布特征如何，變化后均呈高斯分布。在對高光譜數(shù)據(jù)進(jìn)行處理時，通常會假設(shè)高光譜信號數(shù)據(jù)與噪聲之間不相關(guān)，即噪聲協(xié)方差矩陣為對角陣，因此當(dāng)對高光譜觀測數(shù)據(jù)進(jìn)行白化時，建議采用噪聲FT作為其白化矩陣。

2 HySime算法

José M. Bioucas-Dias(2008)[10]根據(jù)最小二乘原理，提出了HySime (hyperspectral signal identification by minimum error)算法。HySime是一種估算高光譜數(shù)據(jù)子空間維數(shù)的方法，該方法首先估計信號與噪聲相關(guān)矩陣，然后選擇以最小均方誤差形式最優(yōu)地表達(dá)信號子空間的特征向量子集。

(5)

特征值分解為

(6)

式中

E=[e1,…,eL]

(7)

所張成的L維空間可以分解為2個互相正交的子空間〈Ek〉和〈Ek〉⊥，對應(yīng)的特征向量分別記為

Ek=?ei1,…,eik」,

(8)

(9)

用于記錄特征向量的排列順序L，即

L={i1,i2,…,iL}。

(10)

若子空間〈Ek〉上的投影矩陣記為

(11)

則y在〈Ek〉上的投影為

(12)

(13)

式中： 第一項表示投影誤差能量，隨著k的增大，該項減?。?第二項表示噪聲能量，隨著k的增大，該項增大。

由于公式(11)和

(14)

將式(13)進(jìn)一步表示為

(15)

(16)

由此，式(15)等號右側(cè)最小值的求取，只需找出δi的所有負(fù)值即可，這樣對應(yīng)的(k,L)即為所求的值。

José M. Bioucas-Dias在介紹HySime算法時，將端元數(shù)目估算過程與信號子空間的確定聯(lián)系起來，通過該算法確定的子空間能最近似地表達(dá)原始數(shù)據(jù)信息，所以該子空間的維數(shù)等價于需要估算的端元數(shù)目。

3 算法改進(jìn)

分析式(15)和式(16)可知，Hysime算法中，噪聲相關(guān)矩陣的估算和信號相關(guān)矩陣的估算是非常重要的2個步驟，它們所對應(yīng)的特征向量值計算也直接影響著最終的結(jié)果。在原始的Hysime算法中，噪聲估算采用的是逐像元的進(jìn)行噪聲估計的全圖光譜維去相關(guān)法，得到每個像元的噪聲估計值后，用原始數(shù)據(jù)觀測值減去噪聲估計值后得到近似真實信號的估計值，進(jìn)而可以計算信號相關(guān)矩陣的估計值。也就是說，原始HySime算法中處理噪聲采取的策略是“去除”[11]，該方法對于準(zhǔn)確估計圖像噪聲的情況下是可行的，但實際應(yīng)用中，要對整幅圖像的每個像元噪聲值都進(jìn)行精準(zhǔn)估計的計算量和難度很大，而且往往不太可能真正實現(xiàn)，因此，利用原始HySime算法進(jìn)行噪聲估計和去除的結(jié)果精度可能并不理想。

(17)

通過上式得解算可知，當(dāng)噪聲為白噪聲，不論噪聲大小如何變化，HySime算法所求得的特征向量保持不變，即特征向量所確定的子空間不隨噪聲的變化而變化。由此可見，若某一觀測數(shù)據(jù)的噪聲為有色噪聲時，可以首先對其噪聲進(jìn)行白化變換，然后再利用HySime算法對數(shù)據(jù)進(jìn)行噪聲估計，從而消除該算法對噪聲變化的影響，提高對噪聲變化的適應(yīng)能力。

(18)

對數(shù)據(jù)進(jìn)行噪聲白化，將白化后的觀測數(shù)據(jù)記為yw，并將yw作為原始Hysime算法的輸入數(shù)據(jù)，此時噪聲為白噪聲，所以其協(xié)方差矩陣為單位矩陣，即

(19)

觀測數(shù)據(jù)相關(guān)矩陣記為Ryw，則信號相關(guān)矩陣為觀測數(shù)據(jù)減去噪聲數(shù)據(jù)，即

(20)

(21)

故式(15)可表示為

(22)

即噪聲白化之后，(2-pij)的所有負(fù)值對應(yīng)的(k,L)即為所求的子空間維數(shù)。

4 驗證與分析

4.1 模擬數(shù)據(jù)實驗

模擬數(shù)據(jù)是在已知的信號數(shù)據(jù)中人工添加噪聲構(gòu)成的。信號數(shù)據(jù)由若干端元波譜乘以各自對應(yīng)的豐度得到，其中端元波譜從美國地質(zhì)調(diào)查局(United States Geological Survey,USGS)光譜庫中隨機(jī)選擇，并將光譜曲線參照AVIRIS的響應(yīng)函數(shù)和中心波長重采樣到224個波段； 豐度數(shù)據(jù)按照Dirichlet分布隨機(jī)生成。選擇Dirichlet分布的原因在于，它可以保證各端元的豐度值為非負(fù)且所有的豐度值相加和為1，而且通過改變Dirichlet分布的參數(shù)，可以方便地改變各端元的分布形狀，使得模擬數(shù)據(jù)的生成有很強(qiáng)的靈活性[13]。在獲得信號數(shù)據(jù)之后，可以根據(jù)實驗分析的需要加入一定量的噪聲值以得到特定的信噪比數(shù)據(jù)。另外，為模擬真實的高光譜數(shù)據(jù)中不存在只含有某一種端元的純像元情況，可通過設(shè)置最大豐度閾值ρ來控制像元“純度”，即在得到服從Dirichlet分布的豐度樣本后，根據(jù)需要對所得候選樣本進(jìn)行篩選。

模擬數(shù)據(jù)實驗的目的主要有2個方面，一是對比不同情況下改進(jìn)算法與原算法的計算結(jié)果，二是以NSP(noise subspace projection)算法[14]作為參考，評價改進(jìn)算法的優(yōu)劣。NSP方法只利用樣本協(xié)方差矩陣或者相關(guān)矩陣的特征值估計端元數(shù)目，并且在樣本數(shù)據(jù)量較少和缺少先驗知識的情況下，仍能較好地進(jìn)行端元估算，具有較好的穩(wěn)定性和實用性。

實驗策略如下： 分別取5,10,15,20為不同的端元數(shù)p，端元豐度滿足Dirichlet分布，構(gòu)成信號數(shù)據(jù)，對以上4種信號數(shù)據(jù)分別加入一定量的噪音，使它們的信噪比SNR分別為15 dB，25 dB和35 dB，加入的噪聲又分為白噪聲和有色噪聲，豐度閾值ρ取1，上述每種組合都生成100組模擬數(shù)據(jù)，利用原始Hysime算法、改進(jìn)的Hysime算法和NSP算法對上述數(shù)據(jù)進(jìn)行測試。由于人工加入的噪聲信號相對簡單，所以噪聲估算方法采用逐像元估算的全圖光譜維去相關(guān)法。這里將虛警概率取值為10e-3和10e-4的NSP算法記為NSP_10e-3和NSP_10e-4。針對每種參數(shù)組合下的100組模擬數(shù)據(jù)進(jìn)行信號子空間估計，將計算結(jié)果的均值和標(biāo)準(zhǔn)差作為評價指標(biāo)。將各種算法最接近端元真實值的均值，和最小的標(biāo)準(zhǔn)差用粗體字體標(biāo)出，若這兩個評價標(biāo)準(zhǔn)同時出現(xiàn)在同一個算法中，則表明該算法結(jié)果較其他算法更準(zhǔn)確和穩(wěn)定。所得結(jié)果如表1所示。

表1 不同參數(shù)條件下不同算法所得結(jié)果的均值和方差比較

通過表1的分析，可得出： ①在大多數(shù)情況下，改進(jìn)的HySime和NSP算法估算的子空間維數(shù)更接近真實端元數(shù)目，不同情況下該算法的穩(wěn)定性也較高； ②當(dāng)噪聲為白噪聲時，改進(jìn)的HySime與原始的HySime算法在不同參數(shù)條件下還具有一定程度的一致性，但噪聲為有色噪聲時，改進(jìn)的算法在端元估計的準(zhǔn)確性和穩(wěn)定性上都要明顯優(yōu)于原算法; ③改進(jìn)的Hysime與NSP算法在不同參數(shù)條件下對端元的估算結(jié)果基本一致，這是因為兩算法在確定最優(yōu)子空間之前，都對原始數(shù)據(jù)的噪聲進(jìn)行了白化處理; ④白噪聲情況下，隨著信噪比的增加，各算法對信號子空間維數(shù)估計的精度也隨著提高。因為信噪比越高，信號對噪聲的抑制能力越好，噪聲對算法的影響也越小，這也從一個側(cè)面說明了噪聲估算和去除的重要性; ⑤在低信噪比條件下，噪聲為有色噪聲時Hysime_White與NSP算法所得結(jié)果與真實值吻合更好。因為當(dāng)噪聲趨向于單波段噪聲時，噪聲主要存在于圖像的少數(shù)幾個波段，此時噪聲更容易與信號分離，因而估計的信號子空間也更準(zhǔn)確。

4.2 真實數(shù)據(jù)實驗

實驗數(shù)據(jù)1997年6月采集的Nevada州Cuprite采礦區(qū)AVIRIS數(shù)據(jù)。該數(shù)據(jù)的DN值代表反射率，它含有400 nm至2 500 nm 的224個波段，輻射分辨率10 nm，空間分辨率20 m。通過初步分析，去掉了由于水汽吸收和信噪比較低的波段，僅使用余下的192個波段。

圖1為USGS發(fā)布的利用AVIRIS數(shù)據(jù)對Cuprite礦區(qū)的礦物填圖結(jié)果。

圖1 Cuprite采礦區(qū)典型礦物分布圖(來源USGS網(wǎng)站)

Cuprite采礦區(qū)有詳細(xì)的地面調(diào)查數(shù)據(jù)及背景資料。Gregg Swayze 等[15-16]利用1990年AVIRIS獲取的數(shù)據(jù)，從中確定了18種礦物類型，并通過對實地采樣樣本的X射線衍射礦物類別分析，驗證了分析結(jié)果; 利用1995年AVIRIS獲取的數(shù)據(jù)，分析了25種礦物的分布情況，并對23種礦物進(jìn)行了礦物填圖。

從上述已有研究的結(jié)果中可以總結(jié)出，對Cuprite采礦區(qū)進(jìn)行“維數(shù)估算”所得結(jié)果大致在18～30之間[14-15]。本文使用與模擬數(shù)據(jù)相同的四種算法，即HySime、改進(jìn)HySime、虛警概率值分別為10e-3和10e-4的NSP算法對該實驗區(qū)域進(jìn)行最優(yōu)子空間估算，其中噪聲估算分別采用全圖光譜維去相關(guān)法和SSDC方法進(jìn)行，所得結(jié)果如表2所示。

表2 Cuprite采礦區(qū)域維數(shù)估算結(jié)果

表2中的維數(shù)估算結(jié)果為26和27，與前人所估計的18-30維[14-15]有較好的一致性。同時也略大于18種和25種端元的結(jié)果，這是由于地物混合中非線性因素的影響、未知地物類型的存在以及光譜變異，得出的光譜維數(shù)估計值往往會大于實際的端元數(shù)[17]。

為了對本次維數(shù)估算結(jié)果的有效性進(jìn)行評估，對Cuprite采礦區(qū)高光譜圖像采用MNF變換，選擇其結(jié)果的不同波段組合進(jìn)行純像元指數(shù)(pure pixel index,PPI)計算，投影的迭代次數(shù)設(shè)置為20 000，閾值系數(shù)設(shè)置為2，在結(jié)果中統(tǒng)計PPI>1的像元數(shù)，結(jié)果如下表3所示：

表3 不同參數(shù)下的PPI計算結(jié)果

PPI越大說明擁有的純像元可能性越大，光譜空間維的可分性越強(qiáng)。表中的數(shù)據(jù)說明，當(dāng)波段數(shù)為26和30時，像元數(shù)目比較接近。通過查看MNF變化后的第40分量影像，其所含有的原影像信息已十分有限，主要體現(xiàn)為噪聲。因此可以推斷，經(jīng)MNF變換后的影像主要信息大多集中在前40個波段中，隨著波段數(shù)的增加，其所含有的噪聲也隨之增多，而此時選取前26和30個波段進(jìn)行PPI的計算所得像元數(shù)目比較接近，說明前26和前30個波段子集已經(jīng)可以對原始高光譜影像信息進(jìn)行較完整地表達(dá)，兩個波段數(shù)下所得的像元數(shù)接近，兩者具有較好的一致性。而波段數(shù)為10和20時，像元數(shù)較少，表明此時選擇的波段子集過小，導(dǎo)致原影像信息表達(dá)得不完全。40波段之后主要為噪聲信息，40和50波段數(shù)時像元數(shù)的增多可能是受噪聲數(shù)據(jù)的影響。

為了進(jìn)一步驗證上述分析，選取上述6種情況下所得像元集合的交集，共得到475個像元。表明這475個像元在6種情況下都被標(biāo)記為純像元，它們受噪聲影響較小，代表了最有可能成為端元的像元，記這一集合為Φ。光譜角表征的是未知光譜的響應(yīng)曲線構(gòu)成的N維空間矢量，與已知類別的光譜響應(yīng)曲線構(gòu)成的N維空間矢量之間的廣義夾角，夾角越小說明越相似。因此，可以通過計算光譜角來分析上述6種情況中除去交集的475個像元，剩余的像元相對于集合Φ的“獨立性”，光譜夾角越大，表明像元越“獨立”于集合Φ，像元屬于純凈像元的可能性越高。上述6種情況下計算所得的光譜夾角的最大值中，波段數(shù)為10時夾角最大值最小，為2.756°； 波段數(shù)為50時夾角最大值最大，為4.188°，而波段數(shù)為26和30時，夾角最大值基本一致，為4.114°，且與其他波段數(shù)的夾角最大值有明顯區(qū)別。這表明維數(shù)估算方法得到的信號子空間能很好地表達(dá)信息，也表明估算得到的維數(shù)應(yīng)用于數(shù)據(jù)分析是有效的。

5 結(jié)論

1)改進(jìn)的HySime維數(shù)估算算法，不需要逐像元地進(jìn)行噪聲估算，一方面減少了計算量、提高了解算效率； 另一方面，消除了因噪聲估算不準(zhǔn)確帶來的誤差效果，提高了結(jié)果精度。

2)改進(jìn)的HySime維數(shù)估算算法不依賴于噪聲均值為零的假設(shè)條件，使原算法有了進(jìn)一步的普適性，針對白化和有色噪聲數(shù)據(jù)，都能獲得接近真實維數(shù)值得結(jié)果。

3) 噪聲估計在實際的高光譜數(shù)據(jù)處理和分析中占有重要作用，不論噪聲為白色或有色的，對它的準(zhǔn)確估計和消除，都會顯著提高高光譜數(shù)據(jù)的解譯精度。

[1] Green A A,Berman M,Switzer P,et al.A transformation for ordering multispectral data in terms of image quality with implications for noise removal[J].IEEE Transactions on Geoscience and Remote Sensing,1988,26(1):65-74.

[2] Lee J B,Woodyatt A S,Berman M.Enhancement of high spectral resolution remote-sensing data by a noise-adjusted principal components transform[J].IEEE Transactions on Geoscience and Remote Sensing,1990,28(3):295-304.

[3] Chang C I,Du Q.Interference and noise-adjusted principal components analysis[J].IEEE Transactions on Geoscience and Remote Sensing,1999,37(5):2387-2396.

[4] 張兵,高連如.高光譜圖像分類與目標(biāo)探測[M].北京:科學(xué)出版社,2011. Zhang B,Gao L R.Hyperspectral Image Classification and Target Detection[M].Beijing:Science Press,2011.

[5] Roger R E.Principal components transform with simple, automatic noise adjustment[J].International Journal of Remote Sensing,1996,17(14):2719-2727.

[6] Roger R E,Arnold J F.Reliably estimating the noise in AVIRIS hyperspectral images[J].International Journal of Remote Sensing,1996,17(10):1951-1962.

[7] 洪波.高光譜遙感圖像信噪比估算方法研究[D].北京:中國科學(xué)院大學(xué),2013. Hong B.Study on Methods for SNR Estimation of Hyperspectral Remote Sensing Images[D].Beijing:University of Chinese Academy of Sciences,2013.

[8] 高連如.高光譜遙感目標(biāo)探測中的信息增強(qiáng)與特征提取研究[D].北京:中國科學(xué)院遙感應(yīng)用研究所,2007. Gao L R.Research on Information Enhancement and Feature Extraction in Hyperspectral Remote Sensing Object Detection[D].Beijing:Institute of Remote Sensing and Digital Earth,Chinese Academy of Sciences,2007.

[9] Chang C I,Du Q.Estimation of number of spectrally distinct signal sources in hyperspectral imagery[J].IEEE Transactions on Geoscience and Remote Sensing,2004,42(3):608-619.

[10]Bioucas-Dias J M,Nascimento J M P.Hyperspectral subspace identification[J].IEEE Transactions on Geoscience and Remote Sensing,2008,46(8):2435-2445.

[11]Eldar Y C,Oppenheim A V.MMSE whitening and subspace whitening[J].IEEE Transactions on Information Theory,2003,49(7):1846-1851.

[12]Cawse K,Robin A,Sears M.The effect of noise whitening on methods for determining the intrinsic dimension of a hyperspectral image[C]//Proceedings of the 3rd Workshop on Hyperspectral Image and Signal Processing:Evolution in Remote Sensing.Lisbon,Portugal:IEEE Computer Society,2011.

[13]Nascimento J M P,Dias J M B.Vertex component analysis:A fast algorithm to unmix hyperspectral data[J].IEEE Transactions on Geoscience and Remote Sensing,2005,43(4):898-910.

[14]Chang C I,Du Q.Noise subspace projection approaches to determination of intrinsic dimensionality of hyperspectral imagery[C]//Proceedings of the SPIE 3871,Image and Signal Processing for Remote Sensing V.Florence,Italy:SPIE,1999:34-44.

[15]Swayze G,Clark R N,Kruse F,et al.Ground-truthing AVIRIS mineral mapping at Cuprite,Nevada[C]//Summaries of the Third Annual JPL Airborne Geoscience Workshop.Denver,CO:JPL Publication,1992:47-49.

[16]Swayze G A.The Hydrothermal and Structural History of the Cuprite Mining District,Southwestern Nevada:An Integrated Geological and Geophysical Approach[D].Boulder,Colorado:University of Colorado,1997.

[17]Chang C I,Xiong W,Liu W M,et al.Linear spectral mixture analysis based approaches to estimation of virtual dimensionality in hyperspectral imagery[J].IEEE Transactions on Geoscience and Remote Sensing,2010,48(11):3960-3979.

(責(zé)任編輯： 邢宇)

Hyperspectral data subspace dimension algorithm based on noise whitening

CHEN Jie1,2, DU Lei1, LI Jing1， HAN Yachao1, GAO Zihong1

(1.ChinaAeroGeophysicalSurveyandRemoteSensingCenterforLandandResources,Beijing100083,China; 2.InstituteofRemoteSensingandDigitalEarth,ChineseAcademyofSciences,Beijing100101,China)

The correlation between adjacent bands of hyperspectral image data is relatively strong. However, signal coexists with noise. The HySime (hyperspectral signal identification by minimum error) algorithm which is based on the principle of least squares is designed to calculate the estimated noise value and the estimated signal correlation matrix value. The algorithm is effective with accurate noise value but ineffective with estimated noise value obtained from spectral dimension reduction and decorrelation process. This paper proposes an improved HySime algorithm based on noise whitening process. Instead of removing noise pixel by pixel, the algorithm carries out the noise whitening process on the original data first, obtains the noise covariance matrix estimated value accurately, and uses the HySime algorithm to calculate the signal correlation matrix value so as to improve the precision of the resultant value. Simulation and experiment have reached some conclusions: Firstly, the improved HySime algorithm is more accurate and stable than the original HySime algorithm; Secondly, the improved HySime algorithm results have better consistency under different conditions compared with the classic NSP (noise subspace the projection) algorithm; Finally, the improved HySime algorithm improves the adaptability of non-white data noise with the noise whitening process.

Hyperspectral dimension reduction； HySime； noise whitening； signal correlation matrix

10.6046/gtzyyg.2017.02.09

陳潔，杜磊，李京,等.基于噪聲白化的高光譜數(shù)據(jù)子空間維數(shù)算法[J].國土資源遙感,2017,29(2):60-66.(Chen J,Du L,Li J,et al.Hyperspectral data subspace dimension algorithm based on noise whitening[J].Remote Sensing for Land and Resources,2017,29(2):60-66.)

2015-12-22；

2016-01-15

中國地質(zhì)調(diào)查局地質(zhì)調(diào)查項目“天山-北山重要成礦區(qū)帶遙感調(diào)查”(編號： 121201003000150008)和“高光譜地質(zhì)調(diào)查技術(shù)方法研究”(編號： 1212031513012)共同資助。

陳潔(1980-),男,高級工程師,主要從事航空攝影技術(shù)方法、數(shù)據(jù)處理及遙感地質(zhì)應(yīng)用的研究。Email： 6592296@qq.com。

TP 751.1

A

1001-070X(2017)02-0060-07