劉 娟，蔡演軍，王 瑾

(1.中國(guó)科學(xué)院地球環(huán)境研究所黃土與第四紀(jì)地質(zhì)國(guó)家重點(diǎn)實(shí)驗(yàn)室，西安 710075;2.中國(guó)科學(xué)院大學(xué)，北京 100049)

0 引言

土地這一界面是遙感圖像上反映最直接的環(huán)境信息，也是研究其他環(huán)境要素的基礎(chǔ)[1]。土壤遙感是指通過(guò)對(duì)遙感影像的解譯，識(shí)別和劃分出土壤類(lèi)型，制作土壤類(lèi)型圖，分析其分布規(guī)律，為改良和合理利用土壤服務(wù)。在常規(guī)土壤調(diào)查方法受限或需要消耗大量人力、物力、時(shí)間的地區(qū)，基于遙感的土壤調(diào)查方法具有不可替代的優(yōu)勢(shì)。

目前，國(guó)內(nèi)外已開(kāi)展了很多采用衛(wèi)星數(shù)據(jù)和地形數(shù)據(jù)進(jìn)行土壤資源調(diào)查和制圖的研究。Odeh等[2]以AVHRR和DEM為數(shù)據(jù)源，采用普通克里金等方法對(duì)土壤性質(zhì)進(jìn)行了研究，并結(jié)合地形數(shù)據(jù)，評(píng)價(jià)了各種方法的優(yōu)缺點(diǎn);Dobos等[3-4]基于AVHRR 500 m，1 000 m空間分辨率的圖像，輔以DEM數(shù)據(jù)對(duì)匈牙利土壤進(jìn)行分類(lèi)研究，發(fā)現(xiàn)地形數(shù)據(jù)的加入對(duì)2種分辨率圖像的分類(lèi)精度都有顯著提高，且高分辨率圖像的分類(lèi)精度明顯高于低分辨率圖像。近年來(lái)，國(guó)內(nèi)也有不少學(xué)者借助遙感圖像進(jìn)行土壤調(diào)查與分類(lèi)。沙晉明等[5]用TM圖像，以浙江省龍游縣為例，對(duì)我國(guó)植被覆蓋率高的東南山區(qū)進(jìn)行了土壤調(diào)查;亢慶等[6-7]基于MODIS遙感數(shù)據(jù)和地形數(shù)據(jù)，在第二次全國(guó)土壤普查數(shù)據(jù)庫(kù)的支持下，探討了遙感技術(shù)在常規(guī)土壤調(diào)查工作受限制的干旱地區(qū)進(jìn)行土壤調(diào)查的效果和適用性。

本文根據(jù)青海湖流域內(nèi)土壤的分布特點(diǎn)，選取一個(gè)具有代表性的區(qū)域?yàn)樵囼?yàn)區(qū)，利用TM圖像和地形數(shù)據(jù)，結(jié)合青海省土壤圖和高分辨率遙感圖像，綜合分析了實(shí)地調(diào)查數(shù)據(jù)、遙感信息特征及植被分布情況，并在此基礎(chǔ)上建立了適合試驗(yàn)區(qū)的土壤遙感分類(lèi)系統(tǒng)，探討了遙感技術(shù)在青海湖流域土壤分類(lèi)中的可行性，為該區(qū)域土壤調(diào)查和分類(lèi)提供新的方法參考。

1 試驗(yàn)區(qū)概況

青海湖流域地處青藏高原東北部，整個(gè)流域近似織梭形，自西北向東南傾斜，是一個(gè)封閉的內(nèi)陸盆地，周?chē)椒宥嘣诤２? 000 m以上，最高處為西北部海拔5 291 m的崗格爾肖合力山。從相對(duì)高度2 000 m左右的山嶺到湖面之間，環(huán)帶狀發(fā)育著寬窄不一的侵蝕構(gòu)造地貌、堆積地貌和風(fēng)積地貌。青海湖流域氣候類(lèi)型為半干旱的溫帶大陸性氣候。氣溫偏低，寒冷期長(zhǎng)，沒(méi)有明顯四季之分，干旱少雨，太陽(yáng)輻射強(qiáng)烈，氣溫日較差大。流域內(nèi)年平均氣溫在-1.1～4.0℃，呈現(xiàn)出由東南向西北遞減的趨勢(shì)。流域年平均降水量在291～579 mm之間，受地形和湖區(qū)影響，降水分布極不均勻[8]。

試驗(yàn)區(qū)位于青海湖流域西北部，E98°50′27″～99°46′32″，N37°1′16″～37°48′39″，南北長(zhǎng) 70.7 km，東西長(zhǎng)53.5 km，面積約為3 783 km2(圖1)。試驗(yàn)區(qū)內(nèi)有4條河流流經(jīng)，分別為布哈河、峻河、夏日哈曲和吉爾孟曲。布哈河是青海湖流域內(nèi)的最大河流，多年平均徑流量為7.825×108m3，占青海湖每年水源補(bǔ)給量的近一半。試驗(yàn)區(qū)內(nèi)的植被類(lèi)型隨海拔的升高，依次為溫性草原、高寒草原、高寒草甸、高寒沼澤草甸。在高海拔地區(qū)陰坡、半陰坡，還分布有斑塊狀的高寒灌叢。

圖1 試驗(yàn)區(qū)地理位置Fig.1 Location of test area

2 數(shù)據(jù)源

2.1 遙感數(shù)據(jù)

選用青海湖流域2009年8月11日Landsat5 TM圖像。該圖像幾乎沒(méi)有云層覆蓋，且8月青海湖流域的植被信息比較豐富，有利于分類(lèi)特征的提取。同時(shí)還購(gòu)買(mǎi)了試驗(yàn)區(qū)內(nèi)部分地區(qū)2010年7月、2012年9月的GeoEye-1高分辨率圖像，用于分類(lèi)樣本和驗(yàn)證樣本的選擇。由于試驗(yàn)區(qū)受人類(lèi)活動(dòng)干擾少且土壤類(lèi)型變化緩慢，雖然高分辨率圖像獲取時(shí)間比TM圖像晚了1～2 a，但仍可認(rèn)為二者具有土壤類(lèi)型識(shí)別對(duì)比的同步性。

2.2 實(shí)地調(diào)查數(shù)據(jù)

研究人員于2009年8月29日—9月3日和2012年9月23—27日分別對(duì)研究區(qū)進(jìn)行了實(shí)地調(diào)查。調(diào)查內(nèi)容包括地形、地貌、土地利用方式、土壤要素和植被類(lèi)型及覆蓋度。調(diào)查時(shí)選擇的土壤剖面與植被調(diào)查剖面一致，包括多種景觀類(lèi)型。觀察土壤剖面的顏色、質(zhì)地、結(jié)構(gòu)及松緊度等采樣點(diǎn)的基本情況。根據(jù)土壤剖面顏色的不同，分層均勻取土約200 g，用于土壤水分、粒度等性質(zhì)的測(cè)定。

2.3 土壤數(shù)據(jù)

使用了青海省測(cè)繪局1984年制作的1∶100萬(wàn)的青海省土壤圖。該土壤圖是基于土壤地理發(fā)生學(xué)進(jìn)行土壤類(lèi)型區(qū)分的，其基本分類(lèi)單元為亞類(lèi)，涉及10個(gè)土綱，17個(gè)土類(lèi)，48個(gè)亞類(lèi)及3個(gè)非土壤單元。在ArcGIS9.3軟件平臺(tái)下將其矢量化，裁剪出屬于青海湖流域的部分，然后分析青海湖流域內(nèi)的土壤分布，進(jìn)而選擇試驗(yàn)區(qū)的范圍。本試驗(yàn)區(qū)內(nèi)包括了5個(gè)土類(lèi)，10個(gè)亞類(lèi)和1個(gè)非土壤單元。

2.4 地形數(shù)據(jù)

本文用到的地形數(shù)據(jù)包括:試驗(yàn)區(qū)1∶5萬(wàn)DEM數(shù)據(jù)，其空間分辨率為25 m，經(jīng)預(yù)處理，重采樣后，用于坡度、坡向等分類(lèi)特征的提取;試驗(yàn)區(qū)1∶10萬(wàn)地形圖，主要用于TM圖像的幾何精糾正。

3 技術(shù)路線

本文采用最大似然方法進(jìn)行分類(lèi)，技術(shù)流程如圖2所示。

圖2 技術(shù)路線Fig.2 Technique flow chart

4 數(shù)據(jù)預(yù)處理與分類(lèi)數(shù)據(jù)集建立

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

在ENVI4.7軟件平臺(tái)的支持下，以青海湖流域1∶10萬(wàn)比例尺地形圖為基準(zhǔn)，對(duì)青海湖流域2009年的TM圖像進(jìn)行幾何精糾正，并保證誤差在1個(gè)像元內(nèi)。選擇Transverse Mercator投影，WGS84坐標(biāo)系，采用雙線性?xún)?nèi)插法進(jìn)行重采樣，分辨率為28.5 m。用試驗(yàn)區(qū)土壤圖的矢量邊界裁剪出試驗(yàn)區(qū)對(duì)應(yīng)的TM圖像和DEM，用于分類(lèi)特征的提取。

4.2 分類(lèi)特征提取

土壤是在母質(zhì)、氣候、地形、生物、時(shí)間等自然因子及人為因素綜合影響下發(fā)生、發(fā)展和演化的，土壤特征反映了各種因素共同作用的結(jié)果。在遙感圖像上，不同土壤類(lèi)型的特征差別不明顯，同時(shí)，由于土壤性狀主要表現(xiàn)在剖面上，而不是表現(xiàn)在土壤表面，因此僅依靠土壤表面電磁波譜的輻射特性來(lái)判別土壤類(lèi)型并不直接。土壤與上述成土因子關(guān)系密切，特別是受主導(dǎo)因素的影響較大[9]，可借助間接解譯標(biāo)志對(duì)其進(jìn)行綜合分析。

4.2.1 主成分變換分量特征

主成分變換能夠?qū)崿F(xiàn)數(shù)據(jù)壓縮和圖像增強(qiáng)，通過(guò)該變換，可提取植被和土壤等主要信息，減少信息冗余和噪聲對(duì)分類(lèi)的影響。一般認(rèn)為T(mén)M的7個(gè)波段可分為 4 組:1，2，3;5，7;4;6。所以，本文中將1，2，3波段和5，7波段分別進(jìn)行主成分變換，取前者的第1和第2主分量，后者的第1主分量以及TM的第4，6波段作為分類(lèi)特征，參與分類(lèi)(分別表示為 PC123，1-2;PC57，1;TM4;TM6)。

4.2.2 纓帽變換分量特征

TM數(shù)據(jù)除熱紅外波段的6個(gè)波段經(jīng)纓帽變換后得到3個(gè)分量。第1分量表征“土壤亮度”，反映了像元總輻射水平，集中了大部分土壤信息，其變化主要反映了不同類(lèi)別土壤的變化;第2分量表征“綠度”，與綠色植被長(zhǎng)勢(shì)、覆蓋度密切相關(guān);第3分量為濕度指數(shù)，對(duì)土壤水分和植物水分都比較敏感[10]，尤其反映的土壤水分狀況有利于土壤分類(lèi)(表示為T(mén)C1-3)[11]。

4.2.3 歸一化植被指數(shù)

植被可以作為土壤類(lèi)型劃分的重要標(biāo)志，植被指數(shù)是對(duì)地表植被活動(dòng)最簡(jiǎn)單，有效的度量。目前已經(jīng)定義了40多種植被指數(shù)，其中歸一化植被指數(shù)(normalized differential vegetation index，NDVI)與植被分布密度呈線性關(guān)系，被認(rèn)為是反映植被生長(zhǎng)狀態(tài)及植被覆蓋度的最佳指標(biāo)[12]，即

式中Red和NIR分別為T(mén)M3和TM4圖像像元亮度值。

4.2.4 水體指數(shù)和濕度指數(shù)

Mcfeeters[13]于 1996 年提出了水體指數(shù)(normalized differential water index，NDWI)，它是基于綠波段和近紅外波段的的歸一化差異水體指數(shù)，通常被用來(lái)提取圖像中的水體信息，如河流、湖泊等，即

式中Green為T(mén)M2圖像像元亮度值。

TM5波段對(duì)植物和土壤的水分含量敏感，而TM2對(duì)水體的反射較強(qiáng)，基于這2個(gè)波段的特點(diǎn)，將它們經(jīng)標(biāo)準(zhǔn)化處理得到濕度指數(shù)(normalized differential moisture index，NDMI)[14]，即

式中SWIR為T(mén)M5圖像像元亮度值。

4.2.5 地形特征

地形是5大成土因素之一，對(duì)土壤的形成有重要意義。試驗(yàn)區(qū)內(nèi)多河流和高山，相對(duì)高差較大，所以海拔(Ele)可以參與分類(lèi)。此外，坡度(Slo)、坡向(Asp)和表面曲率(Cur)可改變水、熱及植被狀況，從而影響土壤的發(fā)育，因此它們也可以參與分類(lèi)。

4.3 分類(lèi)特征數(shù)據(jù)集的建立

將以上得到的15個(gè)分類(lèi)特征進(jìn)行組合，得到分類(lèi)特征數(shù)據(jù)集。在ArcGIS9.3中，將GPS中的實(shí)地采樣點(diǎn)生成矢量樣點(diǎn)圖，配準(zhǔn)之后與遙感圖像進(jìn)行疊加，并分析它們?cè)谶b感圖像上的特點(diǎn)。

5 土壤遙感分類(lèi)及分類(lèi)特征分析

5.1 訓(xùn)練樣本采集

將裁剪出的研究區(qū)域的1∶100萬(wàn)土壤矢量圖以相同的地理坐標(biāo)系疊加到TM圖像上，結(jié)合實(shí)地調(diào)查數(shù)據(jù)、高分辨率遙感圖像，在各個(gè)土壤類(lèi)型中按照約3∶1的比例采集訓(xùn)練數(shù)據(jù)和驗(yàn)證數(shù)據(jù)。采集樣本過(guò)程中應(yīng)注意以下3個(gè)問(wèn)題:

1)分布在不同區(qū)域的相同類(lèi)型的土壤，應(yīng)分別對(duì)各個(gè)區(qū)域采集樣本。原因在于它們?cè)谶b感圖像上表現(xiàn)出的景觀特征可能存在差異，即使在同一區(qū)域內(nèi)，景觀特征可能也不同。例如:在研究區(qū)中有2個(gè)區(qū)域是草甸沼澤土，但是它們?cè)诰坝^特征上卻有很大差異，且在其中一個(gè)較大的區(qū)域中，也存在2種不同的景觀類(lèi)型，所以將草甸沼澤土分為3個(gè)區(qū)域分別采集樣本。此外，普通高山草甸土也有同樣問(wèn)題，根據(jù)實(shí)際情況，分為3個(gè)區(qū)域分別采集樣本。

2)避免在土壤類(lèi)型的交界處選擇樣本。

3)河流附近的土壤類(lèi)型比較復(fù)雜，故不宜在近河流處采集樣本。

5.2 可分離性分析

本文選擇 ENVI中 Jeffries-Matusita，即 J-M距離來(lái)衡量訓(xùn)練樣本的可分離性。其取值在0.0～2.0之間，值越大說(shuō)明類(lèi)對(duì)之間的可分性越好，一般認(rèn)為大于1.9的類(lèi)對(duì)之間具有較好的可分性。本研究以1.9作為類(lèi)對(duì)合并的閾值，對(duì)J-M距離小于1.9的類(lèi)對(duì)進(jìn)行合并(表1)。

表1 J－M距離小于1.9的部分土壤類(lèi)對(duì)Tab.1 Soil pairs which J－M distances is less than 1.9

5.3 類(lèi)別合并

根據(jù)J-M距離，對(duì)可分離性不高的土壤類(lèi)別進(jìn)行合并。將草甸沼澤土2，3與泥炭沼澤土進(jìn)行合并，合并后統(tǒng)稱(chēng)為草甸泥炭沼澤土。草甸沼澤土1與其他亞類(lèi)均可分，還稱(chēng)為草甸沼澤土。把普通高山草甸土和碳酸鹽高山草甸土也進(jìn)行合并，統(tǒng)稱(chēng)為高山草甸土，且合并后的高山草甸土與草甸沼澤土、草甸土可分，J-M距離為1.90。進(jìn)行類(lèi)別合并后試驗(yàn)區(qū)有9個(gè)土壤亞類(lèi)和一個(gè)非土壤單元。

5.4 分類(lèi)特征分析

將不同的分類(lèi)特征進(jìn)行組合，并分析較難區(qū)分的類(lèi)對(duì)之間的可分離性，具體情況見(jiàn)表2，3。

表2 分類(lèi)特征組合Tab.2 Groups of classification feature

表3 不同分類(lèi)特征組合的J－M距離分析Tab.3 Analysis of J－M distances with different groups of classification feature

從表3可以看出:隨著參與分類(lèi)的特征數(shù)量的增加，土壤類(lèi)對(duì)之間的J-M距離也隨之增大，說(shuō)明所選擇的分類(lèi)特征是有效的。組合3的效果比組合2略好，即對(duì)于分類(lèi)NDMI比NDWI更有效。但是為了更好地區(qū)分土壤類(lèi)對(duì)之間的可分性，這2個(gè)分類(lèi)特征都將參與分類(lèi)。組合5中加入了海拔，大大提高了土壤類(lèi)對(duì)之間的可分性，使其中6個(gè)土壤類(lèi)對(duì)由不可分到可分。尤其是草甸沼澤土、草甸土-草甸泥炭沼澤土和草甸沼澤土-沼澤鹽土，它們的J-M 距離分別由 1.68和 1.52提高到 1.98和2.00。組合6—8是在組合5的基礎(chǔ)上分別加入了坡度、坡向、曲率，它們的加入都不同程度地提高了土壤類(lèi)對(duì)之間的可分性。但是對(duì)于栗鈣土-暗栗鈣土來(lái)說(shuō)，坡度和曲率的效果比坡向要好。

5.5 土壤遙感分類(lèi)

雖然有些分類(lèi)特征對(duì)土壤類(lèi)對(duì)可分性的貢獻(xiàn)效率不高，但是為了更有效地區(qū)分土壤類(lèi)對(duì)，將以上15個(gè)分類(lèi)特征都參與到試驗(yàn)區(qū)的遙感土壤分類(lèi)中。圖3為采用最大似然方法得到的分類(lèi)圖。

圖3 試驗(yàn)區(qū)土壤遙感分類(lèi)圖Fig.3 Soil map of the test area

6 分類(lèi)精度檢驗(yàn)及分析

采用混淆矩陣進(jìn)行精度評(píng)價(jià)，結(jié)果見(jiàn)表4。

表4 混淆矩陣Tab.4 Confusion matric

由表4可知，由本文方法計(jì)算得到總體分類(lèi)精度為91.76%，Kappa系數(shù)為0.895 7。各土壤類(lèi)型的錯(cuò)分誤差和漏分誤差見(jiàn)表5。

表5 分類(lèi)誤差分析Tab.5 Error analysis of classification (%)

結(jié)合表4和表5發(fā)現(xiàn)，暗栗鈣土的錯(cuò)分誤差是將屬于自己的像元錯(cuò)分為栗鈣土，而栗鈣土的漏分誤差則是暗栗鈣土的錯(cuò)分造成的。栗鈣土和暗栗鈣土是栗鈣土下面的2個(gè)亞類(lèi)，它們?cè)谒帧⒅脖环矫嬗兄嗨菩?圖3)，造成了界線的模糊，所以在邊緣部分遙感區(qū)分難度較大。前面的分類(lèi)特征分析時(shí)也曾提到，坡度和曲率對(duì)栗鈣土-暗栗鈣土的可分性貢獻(xiàn)度很大，而坡度和曲率的變化在2個(gè)土壤亞類(lèi)的邊界也是一個(gè)漸變的過(guò)程。泥炭土的錯(cuò)分誤差為27.59%，原因是把本該屬于自己的像元錯(cuò)分為草甸泥炭沼澤土。試驗(yàn)區(qū)內(nèi)泥炭土的分布面積很小，而草甸泥炭沼澤土的分布面積很大，泥炭土與草甸泥炭沼澤土中的某些區(qū)域具有及其相似的景觀特征，經(jīng)調(diào)查它們的植被覆蓋均為高寒沼澤草甸。

此外，草甸沼澤土、草甸土的部分像元被誤分為高山草甸土，導(dǎo)致錯(cuò)分誤差也很大。這是因?yàn)楸狙芯克鶇⒖嫉耐寥缊D上，草甸沼澤土、草甸土本來(lái)就是1個(gè)合并類(lèi)，所以比較容易混淆。但歸根結(jié)底還是因?yàn)樗鼈兎植挤秶鷥?nèi)的生物和地理因素比較相似。

7 結(jié)論

本研究以TM圖像和地形數(shù)據(jù)為主，高分辨率遙感圖像和土壤圖為輔，參照土壤發(fā)生分類(lèi)系統(tǒng)，對(duì)試驗(yàn)區(qū)土壤進(jìn)行了遙感分類(lèi)研究，并對(duì)分類(lèi)結(jié)果進(jìn)行了精度驗(yàn)證和分析。

1)研究中所采用的從遙感圖像和DEM中提取的土壤分類(lèi)特征，對(duì)于該尺度的土壤類(lèi)型區(qū)分是有效的，能夠較好地區(qū)分出大部分的土壤亞類(lèi);但對(duì)于某些光譜特征和地形條件相似的亞類(lèi)，區(qū)分度不高，需要進(jìn)行類(lèi)別的合并等調(diào)整。這也說(shuō)明了土壤發(fā)生分類(lèi)系統(tǒng)并不直接適用于遙感分類(lèi)。

2)在青海湖流域這種地形復(fù)雜的區(qū)域，高分辨率遙感圖像的參與對(duì)樣本選擇具有一定的作用。

3)試驗(yàn)區(qū)范圍有限，訓(xùn)練樣本的采集相對(duì)較簡(jiǎn)單，如果對(duì)整個(gè)青海湖流域進(jìn)行土壤遙感分類(lèi)，訓(xùn)練樣本的采集將會(huì)比較困難，而且分類(lèi)精度將會(huì)有所降低，所以本文所研究的土壤遙感分類(lèi)系統(tǒng)在整個(gè)青海湖流域的推廣還有待進(jìn)一步研究和驗(yàn)證。

[1]趙英時(shí).遙感應(yīng)用分析原理與方法[M].北京:科學(xué)出版社，2003.Zhao Y S.The principles and methods for analysis and application of remote sensing[M].Beijing:Science Press，2003.

[2]Odeh I O A，McBratney A B，Chittleborough D J.Further results on prediction of soil properties from terrain attributes:Heterotopicco kriging and regression-kriging[J].Geoderma，1995，67(3/4):215-226.

[3]Dobos E，Micheli E，Baumgardner M F，et al.Use of combined digital elevation model and satellite radiometric data for regional soil mapping[J].Geoderma，2000，97(3/4):367-391.

[4]Dobos E，Montanarella L，Negre T，et al.A regional scale soil mapping approach using integrated AVHRR and DEM Data[J].International Journal of Applied Earth Observation and Geoinformation，2001，3(1):30-42.

[5]沙晉明，史 舟，王人潮，等.東南山區(qū)土壤遙感監(jiān)測(cè)的圖像處理及分類(lèi)[J].水土保持學(xué)報(bào)，2000，14(1):38-43，47.Sha J M，Shi Z，Wang R C，et al.Image processing and classification of soil surveying by remote sensing in southeast mountain area[J].Journal of Soil and Water Conservation，2000，14(1):38-43，47.

[6]亢 慶，張?jiān)鱿椋w曉麗.基于MODIS產(chǎn)品的區(qū)域土壤遙感分類(lèi)研究——以新疆為例[J].遙感技術(shù)與應(yīng)用，2007，22(6):690-695.Kang Q，Zhang Z X，Zhao X L.A regional scale soil classification approach based on MODIS data[J].Remote Sensing Technology and Application，2007，22(6):690-695.

[7]亢 慶，張?jiān)鱿?，趙曉麗.基于遙感技術(shù)的干旱區(qū)土壤分類(lèi)研究[J].遙感學(xué)報(bào)，2008，12(1):159-167.Kang Q，Zhang Z X，Zhao X L.A study of soil classification based on remote sensing in arid area[J].Journal of Remote Sensing，2008，12(1):159-167.

[8]陳桂琛，陳孝全，茍新京.青海湖流域生態(tài)環(huán)境保護(hù)與修復(fù)[M].西寧:青海人民出版社，2008:18-39.Chen G S，Chen X Q，Gou X J.Protection and restoration of the ecological environment at Qinghai Lake basin[M].Xining:The Peoples Press of Qinghai，2008:18-39.

[9]梅安新，彭望琭，秦其明，等.遙感導(dǎo)論[M].北京:高等教育出版社，2001.Mei A X，Peng W L，Qin Q M，et al.Introduction to remote sensing[M].Beijing:Higher Education Press，2001.

[10]羅紅霞，龔健雅，蹇代君.亞熱帶低山丘陵區(qū)土壤遙感監(jiān)測(cè)圖像處理[J].長(zhǎng)江流域資源與環(huán)境，2006，15(1):41-47.Luo H X，Gong J Y，Jian D J.Image processing of soil survey by remote sensing in subtropical hilly lands[J].Resources and Environment in the Yangtze basin，2006，15(1):41-47.

[11]Crist E P，Cicone R C.A physically-based transformation of thematic mapper data-the TM tasseled cap[J].IEEE Transactions on Geoscience and Remote Sensing，1984，22(3):256-263.

[12]鄧書(shū)斌.ENVI遙感圖像處理方法[M].北京:科學(xué)出版社，2010.Deng S B.Image processing methods of remote sensing by ENVI[M].Beijing:Science Press，2010.

[13]Mcfeeters S K.The use of normalized difference water index(NDWI)in the delineation of open water feature[J].International Journal of Remote Sensing，1996，17(7):1425-1423.

[14]趙英時(shí).美國(guó)中西部沙山地區(qū)環(huán)境變化的遙感研究[J].地理研究，2001，20(2):213-219.Zhao Y S，A study on environmental change analysis in sand hill of Nebraska using remote sensing[J].Geographical Research，2001，20(2):213-219.