劉萬(wàn)軍　楊秀紅　曲海成,2　孟　煜　姜慶玲

1(遼寧工程技術(shù)大學(xué)軟件學(xué)院　遼寧 葫蘆島 125105)2(哈爾濱工業(yè)大學(xué)電子與信息工程學(xué)院　黑龍江 哈爾濱 150006)3(鐵嶺師范高等專(zhuān)科學(xué)?！∵|寧 鐵嶺 112000)

?

基于端元子集優(yōu)選的高光譜解混算法研究

劉萬(wàn)軍1楊秀紅1曲海成1,2孟煜1姜慶玲3

1(遼寧工程技術(shù)大學(xué)軟件學(xué)院遼寧 葫蘆島 125105)2(哈爾濱工業(yè)大學(xué)電子與信息工程學(xué)院黑龍江 哈爾濱 150006)3(鐵嶺師范高等專(zhuān)科學(xué)校遼寧 鐵嶺 112000)

摘要針對(duì)采用最大體積單體MVS(Maximization Volume Simplex)端元提取算法進(jìn)行端元初選時(shí)存在相似端元光譜問(wèn)題，提出一種光譜信息散度SID(Spectral Information Divergence)和光譜梯度角SGA(Spectral Gradient Angle)相結(jié)合以區(qū)分兩個(gè)相似端元光譜的方法。該方法對(duì)經(jīng)過(guò)端元初選之后的端元子集進(jìn)行端元的二次選擇，采用以SID_SG作為最相似端元選擇的判據(jù)，除去相似端元，降低相似端元對(duì)解混精度的影響，利用全約束最小二乘法進(jìn)行豐度估計(jì)。實(shí)驗(yàn)結(jié)果表明，提出的優(yōu)化方法與傳統(tǒng)方法相比，提高了端元的選擇精度，重構(gòu)影像與原始影像之間的均方根誤差RMSE(Root Mean Square Error)也有所降低，分布更加均勻。該方法對(duì)高光譜遙感影像進(jìn)行深度解譯具有十分重要的意義。

關(guān)鍵詞光譜解混端元初選二次提取除去端元解混算法

0引言

高光譜遙感技術(shù)已經(jīng)廣泛應(yīng)用于軍事和民用等各個(gè)領(lǐng)域。然而，阻礙高光譜遙感技術(shù)向定量化方向深入發(fā)展的主要因素是混合像元的普遍存在[1]。因?yàn)楦吖庾V傳感器有限的空間分辨率，每個(gè)像元對(duì)應(yīng)著較大的區(qū)域，所以每個(gè)像元一般包含多種物質(zhì)，多種物質(zhì)光譜混合的混合光譜就是傳感器探測(cè)的光譜[2]。解決混合像元問(wèn)題最有效的方法就是混合像元分解[3]。高光譜混合像元分解一般包括兩個(gè)方面，分別是端元提取[4]和豐度估計(jì)[5]。當(dāng)前端元提取方法有N-FINDR[6]、頂點(diǎn)成分分析[7]VCA(Vertex Component Analysis)、迭代誤差分析IEA(Iterative Error Analysis)[8]和純像元指數(shù)PPI(Pure Pixel Index)等。傳統(tǒng)的豐度估計(jì)算法主要有非限制性最小二乘法，非負(fù)限制性最小二乘法，和為一限制性最小二乘法和全限制性最小二乘法[9]。在進(jìn)行高光譜解混時(shí)需要盡可能地選擇出最精確的端元，降低解混時(shí)出現(xiàn)的誤差。

本文結(jié)合兩種常用的光譜區(qū)分方法，提出了一種將光譜信息散度和光譜梯度角相結(jié)合的方法用以區(qū)分兩個(gè)相似端元光譜。該方法不僅能夠從形狀上反映兩條光譜曲線的匹配程度，而且能夠?qū)蓷l有較小差異的光譜曲線做出較大的響應(yīng)。首先利用最大體積單體端元提取算法進(jìn)行端元的初次提取，然后再利用光譜信息散度和光譜梯度角SID_SG(Spectral Information Divergence_ Spectral Gradient)相結(jié)合的算法找出端元子集中存在的相似端元，去除相似端元，最后利用全約束最小二乘法FCLS(Full Constraint Least Square)進(jìn)行端元的豐度估計(jì)，從而實(shí)現(xiàn)高光譜的混合像元分解。

1高光譜解混算法

1.1最大體積單體MVS

一般而言，高光譜圖像中的每個(gè)像元都能看作是圖像中各端元的線性混合的像元:

(1)

(2)

0≤ki≤1

(3)

式中，x是高光譜圖像中任一P維的光譜向量(P為圖像的波段數(shù)目)，E=[e1,e2,…,en]是P×N的矩陣，N是端元的數(shù)目，E中的每一列都是端元向量。c=(c1,c2,…,cn)T是系數(shù)向量，ci是像元x中端元ei所占比例，n表示誤差項(xiàng)。

圖1　兩個(gè)波段三個(gè)端元

在誤差項(xiàng)n非常小時(shí)，只有滿足上面3個(gè)式子的點(diǎn)的集合，才能夠構(gòu)成高維空間的凸集，這些端元就坐落在凸面的單形體的各個(gè)頂點(diǎn)[10]。這里采用兩個(gè)波段、三個(gè)端元作為例子，從而清晰地表明它們之間存在的幾何關(guān)系，如圖1所示。從圖1中看出,三角形體的各個(gè)頂點(diǎn)分別是端元a、b和c，那么三角形內(nèi)部的點(diǎn)就對(duì)應(yīng)著圖像中混合像元。所以提取高光譜圖像的端元問(wèn)題就轉(zhuǎn)化成求單形體的頂點(diǎn)問(wèn)題[11]。

從圖1可以看出，圖1中的端元a，端元b以及端元c所構(gòu)成的三角形面積是圖像中的任意三個(gè)像元所構(gòu)成的面積集合中的最大的一個(gè)，表示如下:

S(a,b,c)=max{S(i,j,k)}

(4)

式中，S(i,j,k)表示的是由像元i、j和k所構(gòu)成的三角形面積。所以尋找單形體的頂點(diǎn)的問(wèn)題就轉(zhuǎn)化為求最大體積單體的問(wèn)題，體積公式如下:

(5)

(6)

式中，ei表示第i個(gè)端元的列向量，V表示的是由e1,e2,…,en這N個(gè)端元所構(gòu)成的單形體體積，|·|表示行列式的運(yùn)算符。因?yàn)檫@里用到了求行列式的運(yùn)算，所以要求E必須為方陣。

1.2光譜信息散度SID

Ii(x)=-lgai

(7)

Ii(y)=-lgbi

(8)

通過(guò)式(7)和式(8)，可得y關(guān)于x的相對(duì)熵：

(9)

同樣可以得到x關(guān)于y的相對(duì)熵：

(10)

因而x和y的光譜信息散度為：

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

(11)

光譜信息散度法的優(yōu)點(diǎn)是可以對(duì)兩條光譜進(jìn)行整體上的比較。

1.3光譜梯度角SGA

光譜梯度角是在對(duì)兩條光譜向量一階求導(dǎo)的基礎(chǔ)上求解它們的梯度向量，然后再得出兩個(gè)梯度向量間的廣義夾角大小[13]。

光譜曲線x的梯度向量為：

sGA(x)=[x2-x1,x3-x2,…,xn-xn-1]

(12)

光譜曲線y的梯度向量為：

sGA(y)=[y2-y1,y3-y2,…,yn-yn-1]

(13)

則梯度向量的廣義夾角大小為：

(14)

該方法的優(yōu)勢(shì)在于能夠反映光譜曲線的斜率變化，能反映出光譜曲線的局部特征的變化。

1.4信息散度與梯度角相結(jié)合

兩條光譜曲線的銳角夾角范圍是(0,π/2)，并且梯度向量的銳角夾角范圍同樣是(0,π/2)。因?yàn)檎泻瘮?shù)的特性是當(dāng)夾角是π/4時(shí)值為1；當(dāng)夾角大于π/4時(shí)，其值隨著夾角的增大迅速變大。根據(jù)正切的函數(shù)特征，將光譜梯度夾角控制在[π/4，π/2]范圍內(nèi)，當(dāng)夾角值改變較小時(shí)，其結(jié)果就能有明顯變化。因此為了不僅可以在整體上比較兩條光譜曲線的匹配程度，而且可以從局部上區(qū)分兩條光譜曲線的差異，提出將SID與SGA相結(jié)合的公式如下：

(15)

1.5均方根誤差RMSE

RMSE[14]是均方根誤差，是一個(gè)度量準(zhǔn)則，用來(lái)評(píng)價(jià)兩幅圖像的相似性，均方根誤差越小，說(shuō)明二者越相近。

(16)

式中，假設(shè)端元豐度矩陣為S=(s11,s12,…,s1n;…;sp1,sp2,…,spn),真實(shí)豐度矩陣為Z=(z11,z12,…,z1n;…;zp1,zp2,…,zpn),sij和zij分別為矩陣S和Z的第i行第j列元素。

1.6豐度估計(jì)誤差

(17)

式中，m是圖像中混合光譜條數(shù)；n是端元組中端元個(gè)數(shù)；qij是第i條光譜第j個(gè)端元的計(jì)算豐度值；pij是第i條光譜第j個(gè)端元的假定豐度值；fa是豐度估計(jì)誤差。

2基于端元子集優(yōu)選的高光譜解混

利用最大體積單體端元提取算法進(jìn)行端元的初次提取，先去除部分端元，再進(jìn)行端元的二次提取迭代循環(huán)，去掉相似的端元，剩下的端元為最佳端元，最后再利用全約束最小二乘法進(jìn)行豐度的估計(jì)，最終實(shí)現(xiàn)高光譜的混合像元分解。端元初次提取的目的是為了去除部分多余的端元，端元二次選擇的目的是進(jìn)一步從原始高光譜影像中找到純像元，從而準(zhǔn)確求解各純物質(zhì)在混合像元中所占的比例。以SID_SG作為最相似的端元選擇準(zhǔn)則，能夠選擇出最相似的端元組，去除均方根誤差較大的端元，保留均方根誤差較小的端元，從而降低相似端元對(duì)解混精度的影響。

算法流程如圖2所示。

圖2　基于端元子集優(yōu)選的高光譜解混算法流程圖

端元二次選擇優(yōu)選算法步驟：

步驟1找出一組相似端元。將端元集合利用光譜匹配算法SID_SG進(jìn)行第一次循環(huán)，找出最相似的一組端元，分別是端元A和端元B。

步驟2計(jì)算出去除端元A后的均方根誤差。去除端元集合中的端元A，利用剩下的端元計(jì)算出均方根誤差。

步驟3計(jì)算出去除端元B后的均方根誤差。去除端元集合中的端元B，利用剩下的端元計(jì)算出均方根誤差。

步驟4比較步驟2和步驟3中的均方根誤差。將步驟2和步驟3中計(jì)算出的均方根誤差進(jìn)行對(duì)比，留下均方根誤差比較小的端元。

步驟5重復(fù)步驟1～4，直至選出均方根誤差最小并且求解豐度所運(yùn)行的時(shí)間是最小的端元組，最終形成新的端元集合。

選出優(yōu)化后的端元組之后，再采用全約束最小二乘法求解高光譜最終的豐度。

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

3.1端元數(shù)目不同時(shí)的對(duì)比實(shí)驗(yàn)

該實(shí)驗(yàn)采用的模擬圖像大小是64×64，端元數(shù)目是3至8個(gè)，并且添加的信噪比SNR(Signal to Noise Ratio)為30 dB。在端元數(shù)目不同的情況下對(duì)比N-FINDR、單體增長(zhǎng)算法SGA(Simplex Growing Algorithm)和最大體積單體(MVS)這三種端元提取算法的優(yōu)劣，見(jiàn)表1所示。

表1　端元數(shù)不同時(shí)的實(shí)驗(yàn)對(duì)比

從表1中可以看出，N-FINDR、單體增長(zhǎng)算法(SGA)、最大體積單體(MVS)這三種算法在端元數(shù)目不同的情況下對(duì)端元提取的結(jié)果進(jìn)行相似性度量的對(duì)比。隨著端元數(shù)目的增多， MVS算法的效果是最好的，SGA算法要優(yōu)于N-FINDR算法，所以采用MVS算法進(jìn)行端元初次選擇的效果要優(yōu)于SGA算法和N-FINDR算法。

3.2信噪比不同時(shí)的對(duì)比實(shí)驗(yàn)

該實(shí)驗(yàn)采用的模擬圖像是64×64，端元數(shù)目是4條，并且添加的信噪比為10至50 dB。在不同的信噪比下對(duì)比N-FINDR，單體增長(zhǎng)算法(SGA)，最大體積單體(MVS)這三種端元提取算法的魯棒性，見(jiàn)表2所示。

表2　信噪比不同時(shí)的實(shí)驗(yàn)對(duì)比

從表2中可以看出，N-FINDR、單體增長(zhǎng)算法(SGA)、最大體積單體(MVS)這三種算法在信噪比不同的情況下對(duì)端元提取的結(jié)果進(jìn)行相似性度量的對(duì)比。對(duì)比表2中數(shù)據(jù)發(fā)現(xiàn)N-FINDR、SGA和MVS這三種算法的魯棒性從高到低依次是MVS算法、SGA算法、N-FINDR算法，所以采用MVS算法進(jìn)行端元的初次選擇的穩(wěn)定性要優(yōu)于SGA算法和N-FINDR算法。

3.3分類(lèi)精度對(duì)比實(shí)驗(yàn)

實(shí)驗(yàn)數(shù)據(jù)集為高光譜研究中經(jīng)常采用的Pavia數(shù)據(jù)[15]，該研究區(qū)域位于意大利北部的帕維亞，該地區(qū)地形復(fù)雜，主要有灌木、水體、森林、農(nóng)田、草地和城鎮(zhèn)等地物類(lèi)型。

為了更好地證明SID_SG方法的優(yōu)越性，利用實(shí)驗(yàn)數(shù)據(jù)將SID_SG方法與光譜角匹配SAM(Spectral Angle Mapping)法，光譜信息散度(SID)法，光譜梯度角(SGA)法這三種方法做對(duì)比實(shí)驗(yàn)，將它們?cè)诜诸?lèi)精度上進(jìn)行對(duì)比驗(yàn)證，見(jiàn)表3所示。

表3　不同分類(lèi)方法的精度分析

從表3中可以看出，分類(lèi)精度最高的是SID_SG，達(dá)到78.96%，它對(duì)所調(diào)查的樣本點(diǎn)正確分類(lèi)的個(gè)數(shù)是30個(gè)。分類(lèi)精度最低的是SAM，它的精度是60.51%，它對(duì)所調(diào)查的樣本點(diǎn)正確分類(lèi)的個(gè)數(shù)是24個(gè)。SID的分類(lèi)精度低于SID_SG的分類(lèi)精度，高于SAM的分類(lèi)精度，它對(duì)所調(diào)查的樣本點(diǎn)正確分類(lèi)的個(gè)數(shù)是27個(gè)，其精度是69.75%。SGA的分類(lèi)精度也低于SID_SG的分類(lèi)精度，高于SAM的分類(lèi)精度，它對(duì)所調(diào)查的樣本點(diǎn)正確分類(lèi)的個(gè)數(shù)是25個(gè)，其精度是62.81%。從這些數(shù)據(jù)可以看出，SID_SG的分類(lèi)精度要高于SAM、SGA和SID，所以采用SID_SG算法提取相似端元的效果要優(yōu)于SAM、SGA和SID。

3.4均方根誤差及豐度估計(jì)時(shí)間對(duì)比

實(shí)驗(yàn)數(shù)據(jù)集為高光譜解混研究中經(jīng)常采用的美國(guó) Nevada 州南部沙漠地區(qū)的Cuprite數(shù)據(jù)[16]。該數(shù)據(jù)中共有224個(gè)波段，波長(zhǎng)范圍為0.37到2.48 μm，光譜分辨率是10 nm。

RMSE是均方根誤差，是一個(gè)度量準(zhǔn)則，用來(lái)評(píng)價(jià)兩幅圖像的相似性，均方根誤差越小，說(shuō)明二者越相近。

利用實(shí)驗(yàn)數(shù)據(jù)，根據(jù)本文優(yōu)選算法分別提取出不同的端元數(shù)，得到的RMSE隨著端元個(gè)數(shù)增加的變化情況，如圖3所示。利用全約束最小二乘法解混求豐度運(yùn)行的時(shí)間，見(jiàn)表4所示。根據(jù)圖3中不同端元數(shù)目得出的均方根誤差值以及表4中給出的在不同端元數(shù)目下利用全約束最小二乘法解混求解豐度所運(yùn)行的時(shí)間值可以看出去除十個(gè)相似的端元是最適宜的，所以在端元二次提取時(shí)提取出12個(gè)端元是最合適的。

圖3　端元數(shù)目和RMSE值的關(guān)系

端元數(shù)目求解豐度所運(yùn)行時(shí)間(s)221156211097208681969318668176341660315578145421352112503

3.5整體豐度估計(jì)誤差對(duì)比

為了研究噪聲對(duì)本文算法性能的影響, 本文在混合光譜中加入了隨機(jī)噪聲。加入噪聲的公式為：

(18)

式中，r′是加噪聲之后的光譜；r是無(wú)噪聲的光譜；u(0,1)是均值為0，方差為1的隨機(jī)噪聲；SNR是信噪比；X是噪聲的假想發(fā)射率，假定為0.5。

本文實(shí)驗(yàn)中將本文算法、MVS算法以及FCLS算法進(jìn)行對(duì)比，這三種算法的豐度估計(jì)誤差隨信噪比的變化情況，如圖4所示。

圖4　三種算法的豐度估計(jì)誤差隨信噪比變化情況

從圖4中可以看出, 當(dāng)信噪比SNR的范圍是[12，100]時(shí), 本文算法豐度估計(jì)誤差fa是最小的，并且當(dāng)信噪比達(dá)到100時(shí)，本文算法的豐度估計(jì)誤差小于0.1，綜上所述，本文算法要優(yōu)于最大體積單體算法和全約束最小二乘算法。

3.6優(yōu)選前后端元及誤差分布對(duì)比

利用Cuprite數(shù)據(jù)進(jìn)行優(yōu)選前后的誤差分布圖對(duì)比。在端元初選時(shí)，利用最大體積單體端元提取算法進(jìn)行端元的一次提取處理，提取出22個(gè)端元，得到其誤差分布圖，如圖5所示。對(duì)22個(gè)端元進(jìn)行端元的二次提取處理后，得到12個(gè)端元的誤差分布圖，如圖6所示。

圖5　初次提取誤差分布圖　　圖6　二次提取誤差分布圖

對(duì)比圖5和圖6，可以看出優(yōu)選后誤差分布更加均勻，降低了誤差，提高了端元選擇的精度，減少了相似端元對(duì)解混精度的影響。

4結(jié)語(yǔ)

本文提出了一種將光譜信息散度和光譜梯度角相結(jié)合的方法用以區(qū)分兩個(gè)相似端元光譜。在利用最大體積單體端元提取算法進(jìn)行端元初選之后，再采用以SID_SG作為最相似端元選擇的判據(jù)進(jìn)行端元的二次選擇，除去相似的端元，降低相似端元對(duì)解混精度的影響。從實(shí)驗(yàn)的結(jié)果可以看出, 提出的基于端元子集優(yōu)選的高光譜解混算法比傳統(tǒng)算法提高了端元的選擇精度，分布也更加均勻，同時(shí)也減少了重構(gòu)影像與原始影像之間的均方根誤差(RMSE)。該方法對(duì)高光譜遙感影像進(jìn)行深度解譯具有十分重要的意義。

參考文獻(xiàn)

[1] Chen G,Song F,Wang J,et al.FRET spectral unmixing:a ratiometric fluorescent nanoprobe for hypochlorite[J].Chemical Communications,2011,24(4):2949-2951.

[2] Roth K L,Dennison P E,Roberts D A.Comparing endmember selection techniques for accurate mapping of plant species and land cover using imaging spectrometer data[J].Remote Sensing of Environment,2012,32(8):139-152.

[3] Gu Y,Wang S,Jia X.Spectral unmixing in multiple-kernel hilbert space for hyperspectral imagery[J].Gon and Rmo Nng Ranaon on,2013,51(10):3968-3981.

[4] Martin G,Plaza A.Spatial-spectral preprocessing prior to endmember identification and unmixing of remotely sensed hyperspectral data[J].Ld O N Ald Arh Obrvaon and Rmo Nng Jornal of,2012,5(12):380-395.

[5] Gonzalez C,Resano J,Plaza A,et al.FPGA implementation of abundance estimation for spectral unmixing of hyperspectral data using the image space reconstruction algorithm[J].Ld O N Ald Arh Obrvaon and Rmo Nng Jornal of,2012,5(3):248-261.

[6] Xiong W,Chang C,Wu C,et al.Fast algorithms to implement N-FINDR for hyperspectral endmember extraction[J].Ld O N Ald Arh Obrvaon and Rmo Nng Jornal of,2011,4(3):545-564.

[7] Krafft C,Diderhoshan M A,Recknagel P,et al.Crisp and soft multivariate methods visualize individual cell nuclei in Raman images of liver tissue sections[J].Vibrational Spectroscopy,2011,55(1):90-100.

[8] Fang Y,Yan H.Error analysis for remote nonlinear iterative learning control system with wireless channel noise[J].Journal of Shanghai University (English Edition),2011,15(3):7-11.

[9] Chang C I.An information-theoretic approach to spectral variability,similarity,and discrimination for hyperspectral image analysis[J].IEEE Transactionson InformationTheory,2000,45(5):1927-1932.

[10] 耿修瑞,張兵,張霞,等.一種基于高維空間凸面單行體體積的高光譜圖像解混算法[J].自然科學(xué)進(jìn)展,2004,7(14):810-814.

[11] 耿修瑞,趙永超,周冠華.一種利用單行體體積自動(dòng)提取高光譜圖像端元的算法[J].自然科學(xué)進(jìn)展,2006,9(16):1196-1200.

[12] Jarosik N,Bennett C L,Dunkley J,et al.Seven-year wilkinson microwave anisotropy probe (WMAP) observations: sky maps,systematic errors,and basic results[J].The Astrophysical Journal Supplement Series,2011,192(2):14.

[13] Liu J,Zhang J.A new maximum simplex volume method based on householder transformation for endmember extraction[J].Ranaon on Gon and Rmo Nng,2012,50(6):104-118.

[14] Chai T,Draxler R R.Root mean square error (RMSE) or mean absolute error (MAE)?[J].Geoscientific Model Development Discussions,2014,7(1):1525-1534.

[15] Brasco L,Pratelli A.Sharp stability of some spectral inequalities[J].Geometric and Functional Analysis,2012,22(1):107-135.

[16] Chen B,Vodacek A,Cahill N D.A novel adaptive scheme for evaluating spectral similarity in high-resolution urban scenes[J].Selected Topics in Applied Earth Observations and Remote Sensing,IEEE Journal of,2013,9(3):1376-1385.

收稿日期：2014-09-28。國(guó)家高技術(shù)研究發(fā)展計(jì)劃項(xiàng)目(2012AA 12A405)；國(guó)家自然科學(xué)基金項(xiàng)目(61172144)。劉萬(wàn)軍，教授，主研領(lǐng)域：數(shù)字圖像處理，運(yùn)動(dòng)目標(biāo)檢測(cè)與跟蹤。楊秀紅，碩士。曲海成，博士。孟煜，碩士。姜慶玲，博士。

中圖分類(lèi)號(hào)TP75

文獻(xiàn)標(biāo)識(shí)碼A

DOI:10.3969/j.issn.1000-386x.2016.07.058

HYPERSPECTRAL UNMIXING ALGORITHM BASED ON ENDMEMBER SUBSET OPTIMISATION

Liu Wanjun1Yang Xiuhong1Qu Haicheng1,2Meng Yu1Jiang Qingling3

1(SchoolofSoftware,LiaoningTechnicalUniversity,Huludao125105,Liaoning,China)2(SchoolofElectronicsandInformationEngineering,HarbinInstituteofTechnology,Harbin150006,Heilongjiang,China)3(TielingNormalCollege,Tieling112000,Liaoning,China)

AbstractFor the problem of maximisation volume simplex (MVS)endmember extraction algorithm that when applying in primary endmembers selection it would have similar endmembers spectra, we proposed an algorithm which distinguishes two similar endmember spectra by combining the spectral information divergence (SID) and the spectral gradient angle (SGA). This algorithm carries out secondary selection on the endmember subset derived from primary endmembesr selection, and adopts SID_SG rule as the criteria for selecting the most similar endmembers to remove the similar endmembers and to reduce the influence of similar endmembers on unmixing accuracy, and uses the full constraint least square for abundance estimation. Experimental results showed that proposed optimisation algorithm improved the accuracy of endmember selection than the traditional method, the root mean square error (RMSE) between reconstruction images and original image was rather reduced, and the distribution was more evenly as well. The algorithm is of very great significance to the deep interpretation on hyperspectral remote sensing image.

KeywordsSpectral unmixingPrimary endmembers selectionSecondary extractionRemoval of endmemberUnmixing algorithm