高光譜遙感圖像目標(biāo)探測(cè)與分類技術(shù)研究目的及意義

高光譜目標(biāo)探測(cè)與分類技術(shù)是高光譜遙感最為重要的應(yīng)用方向，可以有效地應(yīng)用于環(huán)境檢測(cè)、城市規(guī)劃、地質(zhì)巖礦的識(shí)別、海洋水色定量檢測(cè)、植被的精細(xì)分類以及軍事目標(biāo)的探測(cè)等諸多方面。

一般地，目標(biāo)探測(cè)與分類可以從三個(gè)方面入手：圖象空間、光譜空間和特征空間。圖象空間反映了地物的分布和變化以及不同地物之間的空間關(guān)系，圖象是對(duì)地物最直觀的表示。單純基于圖象空間的目標(biāo)探測(cè)與分類屬于一般圖象處理的范疇，處理效果通常取決于圖象的空間分辨率，而往往這正是高光譜圖象的局限所在；光譜空間由景物中的各個(gè)象元對(duì)應(yīng)的光譜曲線所組成，基于光譜空間的分析方法側(cè)重于地物本身的物理屬性，主要利用不同地物有不同光譜特征這一性質(zhì)來(lái)進(jìn)行探測(cè)與處理；景物中的每個(gè)象元在高維特征空間均對(duì)應(yīng)于一個(gè)點(diǎn)，利用高光譜數(shù)據(jù)在高維空間的特性可以更有效的進(jìn)行目標(biāo)探測(cè)與分類的處理，這也是高光譜遙感的巨大內(nèi)涵及獨(dú)到之處。

高光譜遙感是一門新興的科學(xué)，是當(dāng)前遙感技術(shù)的前沿；相對(duì)于高光譜數(shù)據(jù)提供的巨大的數(shù)據(jù)量，當(dāng)前的數(shù)據(jù)處理技術(shù)還有很大的發(fā)展空間。尤其是國(guó)內(nèi)，高光譜圖象處理算法研究尚處于初級(jí)階段，與國(guó)外相比還有很大的差距。因此，本論文期望根據(jù)高光譜數(shù)據(jù)的內(nèi)在特性，綜合高光譜數(shù)據(jù)在圖象空間、光譜空間尤其是特征空間的信息建立起具有一定理論意義和實(shí)用價(jià)值的目標(biāo)探測(cè)與分類算法，進(jìn)一步促進(jìn)高光譜遙感的廣泛應(yīng)用和發(fā)展。

1.2研究現(xiàn)狀

1.2.1成象光譜技術(shù)的發(fā)展概況

商光譜分辨率遙感1(#vpRTBPPFPL 將很妙糧略的尾指利用很多很窄的電磁波波段從感興趣的物體獲取有關(guān)數(shù)據(jù)。高光譜遙感具有102λ的光譜分辨率，它是在電磁波譜的紫外、可見光、近紅外和中紅外區(qū)域，獲取許多非常窄且光譜連續(xù)的圖象數(shù)據(jù)的技術(shù)[Lillesand &Kiefer,1994],所以高光譜遙感技術(shù)通常又被稱為成像光譜技術(shù)。由于成像光譜儀能夠得到上百通道、連續(xù)波段的圖像，從而可以從每個(gè)圖像像元中提取一條完整的光譜曲線[Goetzet al,1981,1985](如圖1.1)

高光譜遙感的出現(xiàn)是遙感界的一場(chǎng)革命。它使本來(lái)在寬波段遙感中不可探測(cè)的物質(zhì)，在高光譜遙感中能被探測(cè)，研究表明許多地表物質(zhì)的吸收特征在吸收峰深度一般處的寬度為20~40nm[Hunt,1980]。而成象光譜系統(tǒng)獲得的連續(xù)波段寬度一般都在10nm以內(nèi)，因此這種數(shù)據(jù)能以足夠的光譜分辨率區(qū)分出那些具有診斷性光譜特征的地表物質(zhì)，這也是高光譜技術(shù)用于探測(cè)甚至識(shí)別地物的物質(zhì)基礎(chǔ)。由于高光譜遙感將確定物質(zhì)或地物性質(zhì)的光譜與把握其空間和幾何關(guān)系的圖象革命性地結(jié)合在一起，它所獲取的地球表面圖象包含了豐富的空間、輻射和光譜三重信息，因而在地質(zhì)、農(nóng)業(yè)、環(huán)境、軍事、水文、大氣等方面都有著巨大的應(yīng)用前景[童慶禧，1999;童慶禧，1995,1996;鄭蘭芬，1995;Cloutis,1996;Green,1998]。自從上世紀(jì)70年代末美國(guó)噴氣推進(jìn)實(shí)驗(yàn)室(JPL)在NASA支持下首先對(duì)成象光譜儀進(jìn)行概念設(shè)計(jì)[Goetz,1992]以來(lái)，高光譜遙感技術(shù)就受到了極大的關(guān)注，世界上一些有條件的國(guó)家，如美、加、澳、法、德等都競(jìng)相投入了大量資本進(jìn)行成象光譜儀器的研制和應(yīng)用[童慶禧，1998]。比如，在NASA支持下，美國(guó)對(duì)內(nèi)華達(dá)州Comstock Lode地區(qū)進(jìn)行了長(zhǎng)達(dá)15年的系統(tǒng)高光譜遙感研究，除獲得了大量的地面光譜資料外，還先后用AIS、AVIRIS、Probe-1和HYMAP等航空高光譜成象儀獲取了多光譜范圍、多時(shí)期、多分辨率的高光譜圖象、取得了一系列成果，從而直接推動(dòng)了高光譜在地質(zhì)遙感應(yīng)用中的發(fā)展[Smailbegovic,2000]。特別是近20年來(lái)，在美國(guó)、加拿大、澳大利亞、日本和中國(guó)等國(guó)的重視下，高光譜遙感技術(shù)不斷發(fā)展。最早的航空成象儀AIS是上世界80年代才出現(xiàn)的[Kruse,1988],隨后從80年代后期至90年代，各國(guó)競(jìng)相研制，出現(xiàn)了一系列航空成象光譜儀。如美國(guó)的HYDICE、SEBASS[Calvin,et al,2000]和AVIRIS[Green,etal,1998],加拿大的FLI、CASI和SFSI,德國(guó)的ROSIS,荷蘭的CAESAR[童慶禧，1995],澳大利亞的HYMAP[童慶禧，1998]。其中HYMAP主要針對(duì)可見光-短波紅外，共有128波段，信噪比達(dá)1000:1,SEBASS主要針對(duì)中紅外(3~5.4μm)和熱紅外(7.8~13.5μm),各分128波段，信噪比達(dá)2000:1,溫度敏感度0.05K[Calvin,2000]。我國(guó)在九五和十五期間，中科院上海技術(shù)物理研究所成功研制了PHI推帚式高光譜成像儀和OMIS實(shí)用性模塊化高光譜成像儀，也達(dá)到了很高的水平。在航空遙感技術(shù)進(jìn)步的同時(shí)，航天高光譜遙感也蓬勃發(fā)展起來(lái)。近年來(lái)，國(guó)外相繼發(fā)射了一系列高光譜衛(wèi)星，如美國(guó)EO-1衛(wèi)星中的HYPERION,MightySat-Ⅱ中的FTHSI,EOS系統(tǒng)中的MODIS、ASTER,美國(guó)以及歐空局ENVISAT系統(tǒng)中的MERIS,歐空局PROBA衛(wèi)星中的CHRIS等。除美國(guó)和歐洲以外，澳大利亞和日本也各有完整的航天高光譜遙感系統(tǒng)發(fā)展計(jì)劃，澳大利亞的ARIES系統(tǒng)，日本也將在今后5年內(nèi)發(fā)射高光譜衛(wèi)星。在我國(guó)的神舟三號(hào)無(wú)人飛船中搭載了一個(gè)中分辨率的成象光譜儀(CMODIS),另外，在我國(guó)的環(huán)境與減災(zāi)小衛(wèi)星星座中，也包括有高光譜遙感器。

1.2.2目標(biāo)探測(cè)與分類技術(shù)研究現(xiàn)狀

隨著成象光譜儀的迅速發(fā)展，高光譜圖象分析與處理算法也日益蓬勃發(fā)展。

目標(biāo)探測(cè)與分類技術(shù)是其中一個(gè)最為重要的分支。一般的，目標(biāo)探測(cè)與分類可以從圖象空間、光譜空間和特征空間來(lái)進(jìn)行。理想情況下，如果成象光譜儀能獲取沒有噪聲和外界干擾的具有無(wú)限空間分辨率的圖象，那么可以從圖象空間上探測(cè)目標(biāo)，然后再根據(jù)所研究區(qū)域的先驗(yàn)知識(shí)從光譜庫(kù)中找出與目標(biāo)光譜具有最佳匹配效果的參考光譜，從而達(dá)到目標(biāo)識(shí)別與分類的目的。然而空間分辨率與光譜分辨率往往是一對(duì)矛盾，目前高光譜遙感數(shù)據(jù)的空間分辨率尚不夠高，使得某些感興趣的小目標(biāo)在圖象中只能以子象元的形式存在，即圖象中只存在感興趣目標(biāo)與別的地物的混合象元，而沒有只含有此目標(biāo)的純粹象元，這無(wú)疑是對(duì)傳統(tǒng)的目標(biāo)探測(cè)與分類技術(shù)的一個(gè)巨大的挑戰(zhàn)。但隨著光譜分辨率的提高、波段數(shù)的增多，高光譜圖象有足夠的維數(shù)容納地物在其特征空間的單形體結(jié)構(gòu)，而傳統(tǒng)的一般光學(xué)圖象處理手段以及多光譜圖象處理方法顯然已經(jīng)不能完全適應(yīng)或者不能充分利用高光譜數(shù)據(jù)在特征空間中的這種內(nèi)在屬性，各種針對(duì)高光譜數(shù)據(jù)這種內(nèi)在特

點(diǎn)的處理算法應(yīng)運(yùn)而生。

目標(biāo)探測(cè)的四個(gè)方面

1.低維空間情況：

高光譜圖象通常都具有幾十甚至上百個(gè)波段，對(duì)如此大的數(shù)據(jù)量進(jìn)行處理往往需要以時(shí)間為代價(jià)，尤其是當(dāng)這些數(shù)據(jù)需要傳輸?shù)臅r(shí)候還要求傳輸介質(zhì)有足夠的帶寬(比如星上數(shù)據(jù)下傳),所以很多情況下，在保證目標(biāo)與背景有足夠分離度的基礎(chǔ)上先對(duì)原始數(shù)據(jù)進(jìn)行降維操作，這其實(shí)也是一個(gè)特征提取的過(guò)程。在遙感圖象分析中，特征提取可以從兩個(gè)意義上來(lái)實(shí)施：一種是按照一定的準(zhǔn)則直接從原始特征空間中選出一個(gè)子集(即子空間),光學(xué)遙感中的波段選擇即屬于這一類。另一種是找到一個(gè)映射關(guān)系P:P:x→y。將原始特征空間X={x?x?,…,x。}映射到維數(shù)降低了的特征空間Y上去，其中Y={y?,y,,…,y,},n<m。高光譜遙感圖象特征提取包含的內(nèi)容非常廣泛，提取方法也很多，光譜維特征提取和空間維特征提取是表現(xiàn)圖象特征提取的兩個(gè)方面。主成分分析(PCA)是一種把原來(lái)多個(gè)指標(biāo)化為少數(shù)幾個(gè)互不相關(guān)的綜合指標(biāo)的一種分析技術(shù)。應(yīng)用PCA可以在基本保持?jǐn)?shù)據(jù)信息量的基礎(chǔ)上把海量的高光譜數(shù)據(jù)投影到僅有少數(shù)幾個(gè)主分量所構(gòu)成的正交子空間中。Jia &Richards[1999]發(fā)展的分塊主成分分析法用于特征提取，取得了一定的分類和顯示效果，他們將全部波段的相關(guān)陣按照相鄰波段的相關(guān)性分成若干塊。一般高度相關(guān)的塊沿對(duì)角線分布，而相關(guān)性低的遠(yuǎn)離對(duì)角線。塊矩陣本身集合了相鄰波段間相關(guān)性高的波段。因此沿相關(guān)矩陣對(duì)角線可分成若干塊(波段組),并對(duì)每組進(jìn)行主成分變換，最后將每組的重要特征(主成分)再重新組合在一起作為進(jìn)一步主成分分析與特征選擇之用。典范分析也是將較多的變量化為少數(shù)幾個(gè)典范變量，通過(guò)這較少的典范變量之間的相關(guān)性來(lái)綜合地描述兩個(gè)多元隨機(jī)變量之間關(guān)系的一種數(shù)學(xué)方法。PCA主要想最大限度地將不同類別分開，而典范分析則是在低維變量間尋找能代表高維變量的相關(guān)性，以達(dá)到分類、識(shí)別目標(biāo)的目的。劉建貴[1999]在分析K-L變換性質(zhì)的基礎(chǔ)上，根據(jù)高光譜數(shù)據(jù)用于城市目標(biāo)識(shí)別提取的特點(diǎn)，提出了面向分類的特征提取的規(guī)范分析(CA)改進(jìn)方法，設(shè)法使原特征空間的各類的樣本點(diǎn)在光譜維上的投影能使類間距離與類內(nèi)距離的比值達(dá)到最大。根據(jù)這一原則來(lái)決定變換矩陣的選擇。劉建貴用這種改進(jìn)的方法實(shí)施對(duì)北京沙河鎮(zhèn)城市地物特征的提取，取得了較好的結(jié)果。

2.已知目標(biāo)、已知背景：

在目標(biāo)端元與背景端元都知道的情況下，最簡(jiǎn)單的目標(biāo)提取方法是直接利用目標(biāo)光譜對(duì)圖象中的象元逐一匹配，從而得到目標(biāo)在圖象中的分布情況。通常的整波形匹配算法主要有以下四種：根據(jù)光譜之間的距離，如最小距離匹配；根據(jù)光譜之間的夾角，代表性的算法是光譜角度填圖(SAM);根據(jù)光譜之間的相似系數(shù)，如交叉相關(guān)光譜匹配技術(shù)(CCSM);編碼匹配算法，如二值編碼匹配算法。整波形匹配的缺陷在于它對(duì)噪聲非常敏感，因而要求圖象光譜有很高的信噪比。實(shí)際上，對(duì)圖象光譜定標(biāo)和反射率轉(zhuǎn)換的精度往往很難達(dá)到光譜匹配的要求。注意到，上述算法并沒有用到已知的背景信息，從而只能在沒有抑制背景的情況下定性得到目標(biāo)地物的大致覆蓋區(qū)域，如果要定量的得到感興趣目標(biāo)在圖象中各個(gè)象元中所占的比例，需要引入混合象元模型。線性混合模型是其中應(yīng)用最為廣泛的模型，它假定：在一定的條件下，高光譜圖象中每個(gè)象元都可以近似認(rèn)為是圖象中各個(gè)端元的線性混合象元。因此，在已經(jīng)知道圖象想所有端元(包括目標(biāo)與背景)的前提下，可以利用線性解混定量地得到目標(biāo)地物在圖象中每個(gè)象元中

所占的百分比。常用的光譜解混算法包括最小二乘法、正交子空間投影算法(0SP)[Harsanyi&Chang,1994]、前景背景分析算法(FBA,Foreground/backgound analysis)[Smith etal.,1994;Smith etal.1996]、濾波廊最篡法成 R)本AARedesso,etal.,1995]等。

3.已知目標(biāo)、未知背景：

在僅知道目標(biāo)端元的光譜而未知背景的情況下，可以通過(guò)三條途徑來(lái)進(jìn)行目

標(biāo)提取，第一種就是前面所講的簡(jiǎn)單的匹配算法，由于這種方法既沒有用到目標(biāo)與背景的信息量分布的差別，也沒有利用高光譜圖象在其特征空間目標(biāo)與背景的相對(duì)位置的差別，因而一般得不到好的目標(biāo)探測(cè)效果；第二種途徑就是利用一組正交基底代替背景，然后再通過(guò)光譜解混、OSP、FBA、濾波向量等方法得到目標(biāo)在圖象中的定量分布結(jié)果，這就歸結(jié)為第二種情形。比較常用的正交變換方法包括PCA及奇異值分解(SVD,Singular value decomposition )[Boardman,189];

第三種途徑則是利用樣本相關(guān)矩陣(或者協(xié)方差矩陣)的性質(zhì)對(duì)目標(biāo)進(jìn)行半解混(Parial Unmixing)提取，匹配濾波(MF,Matched Filter)、約束能量最小化算法(CEM,Constrained energy niiRration(Harsanyi,1993,

1994;Farrand,1997]就是應(yīng)用最為廣泛的小目標(biāo)提取算法，與基于混合象元模型側(cè)重于從目標(biāo)與背景在特征空間相對(duì)位置的差異入手的目標(biāo)提取算法相比，CEM(與MF形式上是相同的)主要從能量(信息量)的角度來(lái)壓制背景繼而提取目標(biāo)的?；旌险{(diào)制匹配濾波器(MTMF,Mixture-Tuned Matched Filtering)[Boardman,1998]則結(jié)合了兩種算法的優(yōu)點(diǎn)。

4.未知目標(biāo)、未知背景：

在圖象中的象元滿足線性混合的條件下，其中的每個(gè)象元都可以由圖象中的所有端元線性混合而成，因此高光譜數(shù)據(jù)端元的提取是理解高光譜數(shù)據(jù)、繼而對(duì)數(shù)據(jù)進(jìn)行進(jìn)一步分析(比如解混、填圖等)的前提條件。如何從高光譜圖象中提取端元一直是一個(gè)熱點(diǎn)問題，在這方面，已經(jīng)有很多成熟的方法被提出和應(yīng)用，Boardman[1993]創(chuàng)造性的提出了利用凸面幾何學(xué)分析的方法提取圖象端元的雛形，他認(rèn)為高光譜圖象的所有數(shù)據(jù)在其特征空間中均由圖象中所有地物所對(duì)應(yīng)的純粹象元(端元)為頂點(diǎn)的單形體所包圍，并與Kruse、Green[1995]在一起發(fā)展了純象元指數(shù)(PPI)提取端元的算法。Craig[1994]提出了最小體積變換(MVT):通過(guò)求取包圍整個(gè)高光譜數(shù)據(jù)“云團(tuán)”的具有最小體積的單形體來(lái)獲得端元。