陳全+周忠發(fā)+魏小島+

摘要：由于喀斯特山區(qū)地表破碎，環(huán)境異質(zhì)性強(qiáng)，影像上混合像元現(xiàn)象嚴(yán)重，增加了草地地上生物量遙感估測的難度。以Ｌａｎｄｓａｔ ＥＴＭ＋為數(shù)據(jù)源，將光譜混合分析技術(shù)應(yīng)用于喀斯特山區(qū)草地地上生物量的遙感估測中，分析了基于線性混合分解模型分解的草地分量與實測草地地上生物量的相關(guān)關(guān)系，進(jìn)而構(gòu)建了喀斯特山區(qū)草地地上生物量估測模型為■＝２ ４５１．１５８ｘ－２８０．４６１（ｒ２＝０．８０１ ２），驗證結(jié)果總體精度可達(dá)８５％以上。基于光譜混合分析建立的草地地上生物量估測模型提高了喀斯特山區(qū)草地地上生物量遙感估測的精度，為喀斯特山區(qū)草地生長狀況調(diào)查與監(jiān)測提供了技術(shù)支持。

關(guān)鍵詞：喀斯特山區(qū)；草地地上生物量；遙感估測；光譜混合分析

中圖分類號：Ｓ１２７；Ｓ８１２．５文獻(xiàn)標(biāo)識碼：Ａ文章編號：0439－８114（２０14）09－2013-04

Ｅｓｔｉｍａｔｉｏｎ ｏｆ Ｋａｒｓｔ Ｍｏｕｎｔａｉｎ Ｇｒａｓｓｌａｎｄ Ａｂｏｖｅｇｒｏｕｎｄ Ｂｉｏｍａｓｓ Ｂａｓｅｄ ｏｎ Ｌｉｎｅａｒ Ｓｐｅｃｔｒａｌ Ｍｉｘｔｕｒｅ Ｍｏｄｅｌ

ＣＨＥＮ Ｑｕａｎ１，２， ＺＨＯＵ Ｚｈｏｎｇ－ｆａ１，２， ＷＥＩ Ｘｉａｏ－ｄａｏ１，２

（１． Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｓｏｕｔｈ Ｃｈｉｎａ Ｋａｒｓｔ， Ｇｕｉｚｈｏｕ Ｎｏｒｍａｌ Ｕｎｉｖｅｒｓｉｔｙ， Ｇｕｉｙａｎｇ ５５０００１， Ｃｈｉｎａ；

２． Ｔｈｅ Ｓｔａｔｅ Ｋｅｙ Ｌａｂｏｒａｔｏｒｙ Ｉｎｃｕｂａｔｉｏｎ Ｂａｓｅ ｆｏｒ Ｋａｒｓｔ Ｍｏｕｎｔａｉｎ Ｅｃｏｌｏｇｙ Ｅｎｖｉｒｏｎｍｅｎｔ ｏｆ Ｇｕｉｚｈｏｕ Ｐｒｏｖｉｎｃｅ， Ｇｕｉｙａｎｇ ５５０００１， Ｃｈｉｎａ）

Ａｂｓｔｒａｃｔ：Ｔｈｅ ｇｒａｓｓｌａｎｄ ｉｓ ａｎ ｉｍｐｏｒｔａｎｔ ｐａｒｔ ｏｆ ｔｈｅ Ｋａｒｓｔ ｍｏｕｎｔａｉｎ ｅｃｏｓｙｓｔｅｍｓ in Ｓｏｕｔｈern Ｃｈｉｎａ and ｔｈｅ ｂａｓｉｓ ｏｆ protecting ｔｈｅ ｋａｒｓｔ ｍｏｕｎｔａｉｎｏｕｓ ｅｃｏｌｏｇｉｃａｌ ｅｎｖｉｒｏｎｍｅｎｔ ａｎｄ ｄｅｖｅｌｏｐing ｏｆ ｇｒａｓｓｌａｎｄ ａｎｉｍａｌ ｈｕｓｂａｎｄｒｙ． Ｉｎ ｔｙｐｉｃａｌ ｋａｒｓｔ ｍｏｕｎｔａｉｎ ａｒｅａ，ｔｈｅ ｃｏｍｐｌｅｘ ａｎｄ ｂｒｏｋｅｎ ｇｒｏｕｎｄ ｓｕｒｆａｃｅ ｌｅａｄｓ ｔｏ ｓｅｒｉｏｕｓ ｍｉｘｅｄ ｐｉｘｅｌ ｐｈｅｎｏｍｅｎｏｎ． Lｉｎｅａｒ ｓｐｅｃｔｒａｌ ｍｉｘｔｕｒｅ model was used ｔｏ ｅｓｔｉｍａｔｅ ｇｒａｓｓｌａｎｄ ａｂｏｖｅｇｒｏｕｎｄ ｂｉｏｍａｓｓ ｂａｓｅｄ ｏｎ Ｌａｎｄｓａｔ ＥＴＭ＋ ｉｍａｇｅ．Ｔｈｅ results showed ｔｈａｔ ｔｈｅ ｇｒａｓｓｌａｎｄ ｆｒａｃｔｉｏｎ ｅｘｔｒａｃｔｅｄ with ｔｈｅ ｌｉｎｅａｒ ｍｉｘｅｄ ｍｏｄｅｌ ｈａd ｏｂｖｉｏｕｓ ｌｉｎｅａｌ ｃｏｒｒｅｌａｔｉｏｎ ｗｉｔｈ ｔｈｅ ｆｉｅｌｄ ｄａｔａ ｍｅａｓｕｒｅｄ． Ｇｒａｓｓｌａｎｄ ａｂｏｖｅｇｒｏｕｎｄ ｂｉｏｍａｓｓ ｅｓｔｉｍａｔｉｏｎ ｍｏｄｅｌ ｃａｎ ｂｅ ｅｘｐｒｅｓｓｅｄ ｂｙ ｔｈｅ ｅｑｕａｔｉｏｎ ■＝２ ４５１．１５８ｘ－２８０．４６１（ｒ２＝０．８０１ ２）． Ｔｈｅ ａｃｃｕｒａｃｙ ｖａｌｉｄａｔｉｏｎ ｓｈｏｗed ｔｈａｔ ｔｈｅ ｏｖｅｒａｌｌ ａｃｃｕｒａｃｙ was ａｂｏｖｅ ８５％． Ｔｈｅ ｅｓｔｉｍａｔｉｏｎ ｍｏｄｅｌ of ｇｒａｓｓｌａｎｄ ａｂｏｖｅｇｒｏｕｎｄ ｂｉｏｍａｓｓ ｗｈｉｃｈ ｂａｓｅｄ ｏｎ ｓｐｅｃｔｒａｌ ｍｉｘｔｕｒｅ model ｃａｎ ｉｍｐｒｏｖｅ ｔｈｅ ｅｓｔｉｍａｔｉｏｎ ａｃｃｕｒａｃｙ ｉｎ ｋａｒｓｔ ｍｏｕｎｔａｉｎｏｕｓ ａｒｅａ． Ｔｈｉｓ ｍｅｔｈｏｄ ｃａｎ ｐｒｏｖｉｄｅ ｔｅｃｈｎｉｃａｌ ｓｕｐｐｏｒｔ ｆｏｒ ｓｕｒｖｅｙing ａｎｄ ｍｏｎｉｔｏｒｉｎｇ the growth of ｋａｒｓｔ ｍｏｕｎｔａｉｎｏｕｓ ｇｒａｓｓ．

Ｋｅｙ ｗｏｒｄｓ： kａｒｓｔ mountain； ｇｒａｓｓｌａｎｄ ａｂｏｖｅｇｒｏｕｎｄ ｂｉｏｍａｓｓ； ｒｅｍｏｔｅ ｓｅｎｓｉｎｇ ｅｓｔｉｍａｔｉｏｎ； ｓｐｅｃｔｒａｌ ｍｉｘｔｕｒｅ ａｎａｌｙｓｉｓ

草地是中國南方喀斯特山區(qū)脆弱生態(tài)系統(tǒng)的重要組成部分，是維持喀斯特山區(qū)生態(tài)平衡的最后一道屏障，同時也是喀斯特山區(qū)生態(tài)環(huán)境保護(hù)與發(fā)展草地畜牧業(yè)的基礎(chǔ)。草地植被地上生物量能客觀反映草地固碳能力與草地載畜能力的大小，是衡量草地資源的重要指標(biāo)，也是進(jìn)行草業(yè)規(guī)劃、生態(tài)評價、草原退化監(jiān)測的基礎(chǔ)［１］。開展喀斯特山區(qū)草地地上生物量的遙感估測，對保護(hù)喀斯特山區(qū)脆弱生態(tài)環(huán)境與實現(xiàn)草地畜牧業(yè)的可持續(xù)發(fā)展，以及推動定量遙感在喀斯特復(fù)雜背景下的應(yīng)用均具有重要的科學(xué)及現(xiàn)實意義。綜合國內(nèi)外已有研究［２－５］可以發(fā)現(xiàn)，草地資源遙感估測的研究多數(shù)都是通過對遙感數(shù)據(jù)的波段組合或利用遙感數(shù)據(jù)計算出的各種植被指數(shù)與草地地上生物量進(jìn)行相關(guān)分析，建立經(jīng)驗?zāi)Ｐ?，并推廣模型以求取大范圍區(qū)域的草地地上生物量。然而，在自然條件復(fù)雜的喀斯特山區(qū)，受遙感影像空間分辨率影響，圖像中的一個像元往往由多個地物組分構(gòu)成，內(nèi)部空間異質(zhì)性很強(qiáng)，草地分布破碎，其植被指數(shù)受巖石、土壤影響強(qiáng)烈，基于植被指數(shù)對喀斯特山區(qū)草地地上生物量進(jìn)行遙感估測誤差明顯。光譜混合分析法能夠較好地解決混合像元的問題，將其引入草地地上生物量的估測，可以有效提高定量遙感的精度。

研究選擇貴州省黔南布依族苗族自治州的獨(dú)山縣兔場鎮(zhèn)作為研究區(qū)域。獨(dú)山縣兔場鎮(zhèn)石漠化退耕還林（草）與草地畜牧業(yè)綜合治理小區(qū)位于東經(jīng)１０７°５１′－１０７°６５′，北緯２５°９１′－２６°０１′之間，地處貴州高原向廣西丘陵盆地過渡的斜坡地帶，境內(nèi)巖溶地貌發(fā)育，地貌類型復(fù)雜多樣，以山地為主，其間交錯分布有少量的丘陵、壩地，面積１０６．２ ｋｍ２，平均海拔１ ０７５ ｍ，屬于中亞熱帶濕潤季風(fēng)氣候區(qū)，年平均氣溫１３．６ ℃，年均降水量１ ４６８ ｍｍ，適宜天然草地生長［６］。選擇Ｌａｎｄｓａｔ ＥＴＭ＋影像為數(shù)據(jù)源，將光譜混合分析技術(shù)應(yīng)用到喀斯特山區(qū)草地地上生物量的估測中，分析了基于線性混合分解模型分解的草地分量與實測草地地上生物量的相關(guān)關(guān)系，進(jìn)而建立了喀斯特山區(qū)草地地上生物量估測模型。

１數(shù)據(jù)獲取與預(yù)處理

１．１影像數(shù)據(jù)獲取與預(yù)處理

選用的遙感影像為２０１１年９月２５日獲取的Ｌａｎｄｓａｔ ＥＴＭ＋影像，軌道號為１２６／４２，太陽高度角５５．５６°，云量較少，影像質(zhì)量較好。使用ＥＮＶＩ ４．８軟件對影像進(jìn)行去條帶處理之后，以１∶５００００地形圖為基準(zhǔn)，均勻選取控制點(diǎn)，利用多項式模型對影像進(jìn)行幾何精校正，校正后圖像均方根誤差在０．４５個像元以內(nèi)，滿足精度要求。精校正后的影像經(jīng)輻射定標(biāo)后，采用Ｌａｎｄｓａｔ ＦＬＡＡＳＨ進(jìn)行大氣校正，得到除熱紅外波段后的其他５個波段的地表反射率［７］。

１．２地面數(shù)據(jù)采集

研究于２０１１年９月初至9月中旬在研究區(qū)開展草地資源地面調(diào)查工作，考慮到喀斯特山區(qū)草地分布的生境條件與植被類型，選取具有代表性的面積為２０ ｍ×２０ ｍ的不同草地類型樣區(qū)２０個，在其中共布設(shè)了８０塊樣地，樣地大?。?ｍ×１ ｍ，樣地均選取在樣區(qū)內(nèi)草地植被長勢較均一的位置。齊地面刈割樣地內(nèi)草本植物，稱取鮮重，統(tǒng)計作為該樣地草地地上生物量。樣地采用ＧＰＳ定位，經(jīng)緯度定位數(shù)據(jù)精度優(yōu)于２ ｍ，高程數(shù)據(jù)精度優(yōu)于５ ｍ。樣地基礎(chǔ)信息調(diào)查主要錄入樣地經(jīng)緯度、海拔、地形、坡度、土壤巖性、草地類型等信息，同時還采集了研究區(qū)典型地物的光譜數(shù)據(jù)，為光譜混合分析提供數(shù)據(jù)。使用的光譜儀型號為ＧＥＲ１５００，波長范圍３５０～１ ０５０ ｎｍ。

２光譜混合分析的草地信息提取

２．１光譜混合模型的選擇

國內(nèi)外對遙感光譜混合分析（Ｓｐｅｃｔｒａｌ Ｍｉｘｔｕｒｅ Ａｎａｌｙｓｉｓ，ＳＭＡ）的研究已較成熟，研究和發(fā)展了多種混合光譜分解模型，其中線性光譜混合模型（Ｌｉｎｅａｒ Ｓｐｅｃｔｒａｌ Ｍｉｘｔｕｒｅ Ｍｏｄｅｌ，ＬＳＭＭ）因具有一定的理論依據(jù)、試驗驗證基礎(chǔ)、一定的精度保證和可操作性強(qiáng)而被廣泛應(yīng)用，是最常用的一種模型［８］。

線性光譜混合模型（ＬＳＳＭ）定義像元在某一光譜波段的反射率（亮度值）是由構(gòu)成像元的端元的反射率（亮度值）以其所占像元面積比例為權(quán)重系數(shù)的線性組合?？捎靡韵鹿接嬎悖?/p>

Ｒｉλ＝■ｆkｉＣｋλ＋εｉλ（１）

■ｆkｉ＝１（ｋ＝１，２，…，ｎ）（２）

０≤ｆkｉ≤１（３）

式中，Ｒｉλ為第λ波段第ｉ像元的光譜反射率； ｆkｉ為對應(yīng)于ｉ像元的第ｋ個基本組分（端元）所占的分量值；Ｃｋλ為第ｋ個基本組分在第λ波段的光譜反射率； εｉλ為殘余誤差值（即光譜的非模型化部分）；ｎ為基本組分的數(shù)目。公式（２）、（３）為ＬＳＳＭ的限制條件。模型計算的結(jié)果表現(xiàn)為各端元豐度圖像和以均方根誤差表示的殘余誤差圖像。

評價模型用殘差εｉλ或均方根誤差ＲＭＳＥ表示：

ＲＭＳＥ＝［■（εｉλ）２／ｎ］１／２ （４）

２．２端元選取

在ＬＳＭＭ模型中，端元（Ｅｎｄ－ｍｅｍｂｅｒ）的選取是成功進(jìn)行光譜混合分析的關(guān)鍵。不當(dāng)?shù)亩嗽x取將導(dǎo)致噪聲、不穩(wěn)定解和不正確的豐度圖像［９］。

端元組分光譜的確定有多種途徑，可以從地物光譜數(shù)據(jù)庫中選取，從圖像自身的像元光譜中獲得或者從外部數(shù)據(jù)源中提取端元等［１０］?？λ固厣絽^(qū)海拔較高，常年多云雨難以及時獲取高質(zhì)量的遙感影像數(shù)據(jù)，影像上相同物質(zhì)的光譜曲線與實驗室光譜曲線匹配難度較大。而從影像自身選出的組分更易于獲得且具有與數(shù)據(jù)相同的度量尺度，同時參考對地面調(diào)查所獲得的草地光譜信息的分析。本研究選擇通過分析影像的光譜特征并結(jié)合地面調(diào)查光譜數(shù)據(jù)的方法來進(jìn)行純凈像元的提取。首先，對多光譜影像進(jìn)行最小噪聲分離變換（ＭＮＦ）分離噪聲，對基于ＭＮＦ變換結(jié)果的前3個分量進(jìn)行純凈像元指數(shù)（ＰＰＩ）計算，將可能是純凈像元的像元集導(dǎo)入Ｎ維可視化分析器進(jìn)行分析，參考地面光譜數(shù)據(jù)，確定純凈像元，最終將各類端元對應(yīng)像元的平均反射率作為該類端元的反射率。通過光譜分析，并參考喀斯特山區(qū)土地覆蓋特點(diǎn)，最終選擇的端元分別為：①草地端元②林地端元③裸土端元④水體等低反射率端元⑤居民點(diǎn)、裸巖等高反射率端元。各類端元的光譜反射率特征曲線如圖１所示。

２．３線性光譜混合分解結(jié)果

基本組分端元及其光譜特征確定以后，基于獲取的端元光譜特征，采用全受限的ＬＳＭＭ對Landsat ＥＴＭ＋影像進(jìn)行分解，得到草地分量圖。圖２為光譜混合分析的圖像結(jié)果，圖２ａ為草地端元分量圖。在光譜分解結(jié)果中，白色或較亮顏色表示像元中該端元組分比例較高，相反，黑色或較暗顏色則表示像元中該端元組分的比例較低。均方根誤差ＲＭＳＥ是評價ＬＳＭＭ 模型的一個重要指標(biāo)，均方根誤差越小分類計算的準(zhǔn)確性就越高。圖２ｆ為應(yīng)用線性模型分解求解過程產(chǎn)生的各像元ＲＭＳＥ圖像，其中最大的均方根誤差僅有０．００５，小于０．００３的均方根誤差占到了９６．２％，表明此次光譜混合分解較成功［１１］。

３草地地上生物量遙感估測結(jié)果

３．１建立估測模型

本研究于２０１１年９月初至９月中旬在研究區(qū)開展草地資源地面調(diào)查工作，地面調(diào)查時間與影像獲取時間基本同步，共布設(shè)８０塊樣地，考慮到之后估測模型的精度驗證，１６塊樣地數(shù)據(jù)不參與建模。利用ＳＰＳＳ１９．０軟件對研究區(qū)６４個樣地調(diào)查所得草地植被地上生物量數(shù)據(jù)與線性分解得到的草地分量進(jìn)行相關(guān)分析發(fā)現(xiàn)，二者相關(guān)系數(shù)為０．８９６，呈顯著相關(guān)關(guān)系，因而可通過對得到的草地分量進(jìn)行變換得到草地地上生物量。

為了進(jìn)一步分析草地地上生物量與草地分量之間的定量關(guān)系，采用一元線性回歸模型、指數(shù)函數(shù)、對數(shù)函數(shù)、冪函數(shù)等對二者進(jìn)行模擬預(yù)測分析。對不同模型測算結(jié)果的比較與分析結(jié)果表明，基于草地分量的一元線性回歸模型擬合精度最高。喀斯特山區(qū)草地地上生物量最優(yōu)擬合模型為：

■＝２ ４８１．１５８ｘ－２８０．４６１（r２＝０．８０１ ２）（５）

式中，ｙ為草地秋季當(dāng)年地上生物量（ｇ／ｍ２），ｘ為草地分量。二者的線性擬合圖見圖3。

３．２模型精度驗證

為評價草地地上生物量估測模型的應(yīng)用精度，通過野外同期采樣的另外１６個樣地數(shù)據(jù)與模型估算的數(shù)據(jù)進(jìn)行比較，并利用誤差統(tǒng)計公式［公式（６）］對回歸模型的精度進(jìn)行綜合評價［１２］。

Ｆ＝（Ｌ■－Ｌ■／Ｌ■）×１００%（６）

式中Ｆ為相對誤差，Ｌ■和Ｌ■分別代表地面實測的草地地上生物量和經(jīng)模型反演得到的數(shù)值。結(jié)果（表１）表明，模擬值同實測值之間的相對誤差為５．３６％～２４．６２％，平均相對誤差為１３．２９％，模型的總體精度達(dá)到８５％以上，說明利用線性光譜分解提取的草地分量建立的一元線性回歸模型能較好地反映草地地上植被生物量狀況，能夠滿足喀斯特山區(qū)中尺度的草地地上生物量估算。

４結(jié)論

受喀斯特山區(qū)地形破碎、環(huán)境差異大、遙感影像空間分辨率和混合像元的限制等因素影響，利用傳統(tǒng)的植被指數(shù)方法對喀斯特山區(qū)草地地上生物量進(jìn)行遙感估測誤差明顯。將光譜混合分析技術(shù)應(yīng)用于喀斯特山區(qū)草地地上生物量的估測研究，通過線性光譜分解，得到符合研究區(qū)環(huán)境特征的端元光譜和各分量的豐度，利用得到的草地分量與地面樣地調(diào)查所得草地地上生物量數(shù)據(jù)進(jìn)行回歸，建立了草地地上生物量估測的線性回歸模型■＝２４８１．１５８ｘ－２８５０．４６１（ｒ２＝０．８０１ ２），基于地面同期調(diào)查數(shù)據(jù)的精度驗證也顯示利用光譜混合模型進(jìn)行喀斯特山區(qū)草地地上生物量估測的精度可達(dá)８５％以上，能夠滿足喀斯特山區(qū)的草地地上生物量的估算。因此，基于光譜混合分析對喀斯特山區(qū)草地地上生物量進(jìn)行遙感估測是可行的。

參考文獻(xiàn)：

［１］ 董永平，吳新宏，戎郁萍，等．草原遙感監(jiān)測技術(shù)［Ｍ］．北京：化學(xué)工業(yè)出版社，２００５．

［２］ 查勇，ＧＡＯ Ｊ，倪紹祥．國際草地資源遙感研究新進(jìn)展［Ｊ］．地理科學(xué)進(jìn)展，２００３，２２（６）：６０７－６１７．

［３］ 徐斌，楊秀春，陶偉國，等．中國草原產(chǎn)草量遙感監(jiān)測［Ｊ］．生態(tài)學(xué)報，２００７，２７（２）：４０５－４１３．

［４］ 楊秀春，徐斌，朱曉華，等．北方農(nóng)牧交錯帶草原產(chǎn)草量遙感監(jiān)測模型［Ｊ］．地理研究，２００７，２６（２）：２１３－２２１．

［５］ 王靜，郭鈮，王振國，等．甘南草地地上部生物量遙感監(jiān)測模型［Ｊ］．干旱氣象，２０１０，２８（２）：１２８－１３３．

［６］ 熊康寧，陳永畢，陳滸，等．點(diǎn)石成金——貴州石漠化治理技術(shù)與模式［Ｍ］．貴陽：貴州科學(xué)技術(shù)出版社，２０１１．

［７］ 趙英時．遙感應(yīng)用分析原理與方法［Ｍ］．北京：科學(xué)出版社，２００３．

［８］ ＡＳＮＥＲ Ｇ， ＬＯＢＥＬＬ Ｄ． Ａ ｂｉｏｇｅｏｐｈｙｓｉｃａｌ ａｐｐｒｏａｃｈ ｆｏｒ ａｕｔｏｍａｔｅｄ ＳＷＩＲ ｕｎｍｉｘｉｎｇ ｏｆ ｓｏｉｌｓ ａｎｄ ｖｅｇｅｔａｔｉｏｎ［Ｊ］． Ｒｅｍｏｔｅ Ｓｅｎｓｉｎｇ ｏｆ Ｅｎｖｉｒｏｎｍｅｎｔ，２０００，７４（１）：９９－１１２．

［９］ ＧＩＬＡＢＥＲＴ Ｍ Ａ． Ａｎ ａｔｍｏｓｐｈｅｒｉｃ ｃｏｒｒｅｃｔｉｏｎ ｍｅｔｈｏｄ ｆｏｒ ｔｈｅ ａｕｔｏｍａｔｉｃ ｒｅｔｒｉｅｖａｌ ｏｆ ｓｕｒｆａｃｅ ｒｅｆｌｅｃｔａｎｃｅ ｆｒｏｍ ＴＭ ｉｍａｇｅｓ［Ｊ］． Ｉｎｔ Ｊ Ｒｅｍｏｔｅ Ｓｅｎｓｉｎｇ，１９９４，１５（１０）：２０６５－２０８６．

［１０］ 李曉松，李增元，吳波，等．基于光譜混合分析的毛烏素沙地油蒿群落覆蓋度提取［Ｊ］．遙感學(xué)報，２００７，１１（６）：９２３－９３０．

［１１］ 李素，李文正，周建軍，等．遙感影像混合像元分解中的端元選擇方法綜述［Ｊ］．地理與地理信息科學(xué)，２００７，２３（５）：３５－３８．

［１２］ 梁天剛，崔霞，馮琦勝，等．２００１－２００８年甘南牧區(qū)草地地上生物量與載畜量遙感動態(tài)監(jiān)測［Ｊ］．草業(yè)學(xué)報，２００９，１８（６）：１２－２２．

２光譜混合分析的草地信息提取

２．１光譜混合模型的選擇

國內(nèi)外對遙感光譜混合分析（Ｓｐｅｃｔｒａｌ Ｍｉｘｔｕｒｅ Ａｎａｌｙｓｉｓ，ＳＭＡ）的研究已較成熟，研究和發(fā)展了多種混合光譜分解模型，其中線性光譜混合模型（Ｌｉｎｅａｒ Ｓｐｅｃｔｒａｌ Ｍｉｘｔｕｒｅ Ｍｏｄｅｌ，ＬＳＭＭ）因具有一定的理論依據(jù)、試驗驗證基礎(chǔ)、一定的精度保證和可操作性強(qiáng)而被廣泛應(yīng)用，是最常用的一種模型［８］。

線性光譜混合模型（ＬＳＳＭ）定義像元在某一光譜波段的反射率（亮度值）是由構(gòu)成像元的端元的反射率（亮度值）以其所占像元面積比例為權(quán)重系數(shù)的線性組合。可用以下公式計算：

Ｒｉλ＝■ｆkｉＣｋλ＋εｉλ（１）

■ｆkｉ＝１（ｋ＝１，２，…，ｎ）（２）

０≤ｆkｉ≤１（３）

式中，Ｒｉλ為第λ波段第ｉ像元的光譜反射率； ｆkｉ為對應(yīng)于ｉ像元的第ｋ個基本組分（端元）所占的分量值；Ｃｋλ為第ｋ個基本組分在第λ波段的光譜反射率； εｉλ為殘余誤差值（即光譜的非模型化部分）；ｎ為基本組分的數(shù)目。公式（２）、（３）為ＬＳＳＭ的限制條件。模型計算的結(jié)果表現(xiàn)為各端元豐度圖像和以均方根誤差表示的殘余誤差圖像。

評價模型用殘差εｉλ或均方根誤差ＲＭＳＥ表示：

ＲＭＳＥ＝［■（εｉλ）２／ｎ］１／２ （４）

２．２端元選取

在ＬＳＭＭ模型中，端元（Ｅｎｄ－ｍｅｍｂｅｒ）的選取是成功進(jìn)行光譜混合分析的關(guān)鍵。不當(dāng)?shù)亩嗽x取將導(dǎo)致噪聲、不穩(wěn)定解和不正確的豐度圖像［９］。

端元組分光譜的確定有多種途徑，可以從地物光譜數(shù)據(jù)庫中選取，從圖像自身的像元光譜中獲得或者從外部數(shù)據(jù)源中提取端元等［１０］?？λ固厣絽^(qū)海拔較高，常年多云雨難以及時獲取高質(zhì)量的遙感影像數(shù)據(jù)，影像上相同物質(zhì)的光譜曲線與實驗室光譜曲線匹配難度較大。而從影像自身選出的組分更易于獲得且具有與數(shù)據(jù)相同的度量尺度，同時參考對地面調(diào)查所獲得的草地光譜信息的分析。本研究選擇通過分析影像的光譜特征并結(jié)合地面調(diào)查光譜數(shù)據(jù)的方法來進(jìn)行純凈像元的提取。首先，對多光譜影像進(jìn)行最小噪聲分離變換（ＭＮＦ）分離噪聲，對基于ＭＮＦ變換結(jié)果的前3個分量進(jìn)行純凈像元指數(shù)（ＰＰＩ）計算，將可能是純凈像元的像元集導(dǎo)入Ｎ維可視化分析器進(jìn)行分析，參考地面光譜數(shù)據(jù)，確定純凈像元，最終將各類端元對應(yīng)像元的平均反射率作為該類端元的反射率。通過光譜分析，并參考喀斯特山區(qū)土地覆蓋特點(diǎn)，最終選擇的端元分別為：①草地端元②林地端元③裸土端元④水體等低反射率端元⑤居民點(diǎn)、裸巖等高反射率端元。各類端元的光譜反射率特征曲線如圖１所示。

２．３線性光譜混合分解結(jié)果

基本組分端元及其光譜特征確定以后，基于獲取的端元光譜特征，采用全受限的ＬＳＭＭ對Landsat ＥＴＭ＋影像進(jìn)行分解，得到草地分量圖。圖２為光譜混合分析的圖像結(jié)果，圖２ａ為草地端元分量圖。在光譜分解結(jié)果中，白色或較亮顏色表示像元中該端元組分比例較高，相反，黑色或較暗顏色則表示像元中該端元組分的比例較低。均方根誤差ＲＭＳＥ是評價ＬＳＭＭ 模型的一個重要指標(biāo)，均方根誤差越小分類計算的準(zhǔn)確性就越高。圖２ｆ為應(yīng)用線性模型分解求解過程產(chǎn)生的各像元ＲＭＳＥ圖像，其中最大的均方根誤差僅有０．００５，小于０．００３的均方根誤差占到了９６．２％，表明此次光譜混合分解較成功［１１］。

３草地地上生物量遙感估測結(jié)果

３．１建立估測模型

本研究于２０１１年９月初至９月中旬在研究區(qū)開展草地資源地面調(diào)查工作，地面調(diào)查時間與影像獲取時間基本同步，共布設(shè)８０塊樣地，考慮到之后估測模型的精度驗證，１６塊樣地數(shù)據(jù)不參與建模。利用ＳＰＳＳ１９．０軟件對研究區(qū)６４個樣地調(diào)查所得草地植被地上生物量數(shù)據(jù)與線性分解得到的草地分量進(jìn)行相關(guān)分析發(fā)現(xiàn)，二者相關(guān)系數(shù)為０．８９６，呈顯著相關(guān)關(guān)系，因而可通過對得到的草地分量進(jìn)行變換得到草地地上生物量。

為了進(jìn)一步分析草地地上生物量與草地分量之間的定量關(guān)系，采用一元線性回歸模型、指數(shù)函數(shù)、對數(shù)函數(shù)、冪函數(shù)等對二者進(jìn)行模擬預(yù)測分析。對不同模型測算結(jié)果的比較與分析結(jié)果表明，基于草地分量的一元線性回歸模型擬合精度最高?？λ固厣絽^(qū)草地地上生物量最優(yōu)擬合模型為：

■＝２ ４８１．１５８ｘ－２８０．４６１（r２＝０．８０１ ２）（５）

式中，ｙ為草地秋季當(dāng)年地上生物量（ｇ／ｍ２），ｘ為草地分量。二者的線性擬合圖見圖3。

３．２模型精度驗證

為評價草地地上生物量估測模型的應(yīng)用精度，通過野外同期采樣的另外１６個樣地數(shù)據(jù)與模型估算的數(shù)據(jù)進(jìn)行比較，并利用誤差統(tǒng)計公式［公式（６）］對回歸模型的精度進(jìn)行綜合評價［１２］。

Ｆ＝（Ｌ■－Ｌ■／Ｌ■）×１００%（６）

式中Ｆ為相對誤差，Ｌ■和Ｌ■分別代表地面實測的草地地上生物量和經(jīng)模型反演得到的數(shù)值。結(jié)果（表１）表明，模擬值同實測值之間的相對誤差為５．３６％～２４．６２％，平均相對誤差為１３．２９％，模型的總體精度達(dá)到８５％以上，說明利用線性光譜分解提取的草地分量建立的一元線性回歸模型能較好地反映草地地上植被生物量狀況，能夠滿足喀斯特山區(qū)中尺度的草地地上生物量估算。

４結(jié)論

受喀斯特山區(qū)地形破碎、環(huán)境差異大、遙感影像空間分辨率和混合像元的限制等因素影響，利用傳統(tǒng)的植被指數(shù)方法對喀斯特山區(qū)草地地上生物量進(jìn)行遙感估測誤差明顯。將光譜混合分析技術(shù)應(yīng)用于喀斯特山區(qū)草地地上生物量的估測研究，通過線性光譜分解，得到符合研究區(qū)環(huán)境特征的端元光譜和各分量的豐度，利用得到的草地分量與地面樣地調(diào)查所得草地地上生物量數(shù)據(jù)進(jìn)行回歸，建立了草地地上生物量估測的線性回歸模型■＝２４８１．１５８ｘ－２８５０．４６１（ｒ２＝０．８０１ ２），基于地面同期調(diào)查數(shù)據(jù)的精度驗證也顯示利用光譜混合模型進(jìn)行喀斯特山區(qū)草地地上生物量估測的精度可達(dá)８５％以上，能夠滿足喀斯特山區(qū)的草地地上生物量的估算。因此，基于光譜混合分析對喀斯特山區(qū)草地地上生物量進(jìn)行遙感估測是可行的。

參考文獻(xiàn)：

［１］ 董永平，吳新宏，戎郁萍，等．草原遙感監(jiān)測技術(shù)［Ｍ］．北京：化學(xué)工業(yè)出版社，２００５．

［２］ 查勇，ＧＡＯ Ｊ，倪紹祥．國際草地資源遙感研究新進(jìn)展［Ｊ］．地理科學(xué)進(jìn)展，２００３，２２（６）：６０７－６１７．

［３］ 徐斌，楊秀春，陶偉國，等．中國草原產(chǎn)草量遙感監(jiān)測［Ｊ］．生態(tài)學(xué)報，２００７，２７（２）：４０５－４１３．

［４］ 楊秀春，徐斌，朱曉華，等．北方農(nóng)牧交錯帶草原產(chǎn)草量遙感監(jiān)測模型［Ｊ］．地理研究，２００７，２６（２）：２１３－２２１．

［５］ 王靜，郭鈮，王振國，等．甘南草地地上部生物量遙感監(jiān)測模型［Ｊ］．干旱氣象，２０１０，２８（２）：１２８－１３３．

［６］ 熊康寧，陳永畢，陳滸，等．點(diǎn)石成金——貴州石漠化治理技術(shù)與模式［Ｍ］．貴陽：貴州科學(xué)技術(shù)出版社，２０１１．

［７］ 趙英時．遙感應(yīng)用分析原理與方法［Ｍ］．北京：科學(xué)出版社，２００３．

［８］ ＡＳＮＥＲ Ｇ， ＬＯＢＥＬＬ Ｄ． Ａ ｂｉｏｇｅｏｐｈｙｓｉｃａｌ ａｐｐｒｏａｃｈ ｆｏｒ ａｕｔｏｍａｔｅｄ ＳＷＩＲ ｕｎｍｉｘｉｎｇ ｏｆ ｓｏｉｌｓ ａｎｄ ｖｅｇｅｔａｔｉｏｎ［Ｊ］． Ｒｅｍｏｔｅ Ｓｅｎｓｉｎｇ ｏｆ Ｅｎｖｉｒｏｎｍｅｎｔ，２０００，７４（１）：９９－１１２．

［９］ ＧＩＬＡＢＥＲＴ Ｍ Ａ． Ａｎ ａｔｍｏｓｐｈｅｒｉｃ ｃｏｒｒｅｃｔｉｏｎ ｍｅｔｈｏｄ ｆｏｒ ｔｈｅ ａｕｔｏｍａｔｉｃ ｒｅｔｒｉｅｖａｌ ｏｆ ｓｕｒｆａｃｅ ｒｅｆｌｅｃｔａｎｃｅ ｆｒｏｍ ＴＭ ｉｍａｇｅｓ［Ｊ］． Ｉｎｔ Ｊ Ｒｅｍｏｔｅ Ｓｅｎｓｉｎｇ，１９９４，１５（１０）：２０６５－２０８６．

［１０］ 李曉松，李增元，吳波，等．基于光譜混合分析的毛烏素沙地油蒿群落覆蓋度提?。郏剩荩b感學(xué)報，２００７，１１（６）：９２３－９３０．

［１１］ 李素，李文正，周建軍，等．遙感影像混合像元分解中的端元選擇方法綜述［Ｊ］．地理與地理信息科學(xué)，２００７，２３（５）：３５－３８．

［１２］ 梁天剛，崔霞，馮琦勝，等．２００１－２００８年甘南牧區(qū)草地地上生物量與載畜量遙感動態(tài)監(jiān)測［Ｊ］．草業(yè)學(xué)報，２００９，１８（６）：１２－２２．

２光譜混合分析的草地信息提取

２．１光譜混合模型的選擇

國內(nèi)外對遙感光譜混合分析（Ｓｐｅｃｔｒａｌ Ｍｉｘｔｕｒｅ Ａｎａｌｙｓｉｓ，ＳＭＡ）的研究已較成熟，研究和發(fā)展了多種混合光譜分解模型，其中線性光譜混合模型（Ｌｉｎｅａｒ Ｓｐｅｃｔｒａｌ Ｍｉｘｔｕｒｅ Ｍｏｄｅｌ，ＬＳＭＭ）因具有一定的理論依據(jù)、試驗驗證基礎(chǔ)、一定的精度保證和可操作性強(qiáng)而被廣泛應(yīng)用，是最常用的一種模型［８］。

線性光譜混合模型（ＬＳＳＭ）定義像元在某一光譜波段的反射率（亮度值）是由構(gòu)成像元的端元的反射率（亮度值）以其所占像元面積比例為權(quán)重系數(shù)的線性組合?？捎靡韵鹿接嬎悖?/p>

Ｒｉλ＝■ｆkｉＣｋλ＋εｉλ（１）

■ｆkｉ＝１（ｋ＝１，２，…，ｎ）（２）

０≤ｆkｉ≤１（３）

式中，Ｒｉλ為第λ波段第ｉ像元的光譜反射率； ｆkｉ為對應(yīng)于ｉ像元的第ｋ個基本組分（端元）所占的分量值；Ｃｋλ為第ｋ個基本組分在第λ波段的光譜反射率； εｉλ為殘余誤差值（即光譜的非模型化部分）；ｎ為基本組分的數(shù)目。公式（２）、（３）為ＬＳＳＭ的限制條件。模型計算的結(jié)果表現(xiàn)為各端元豐度圖像和以均方根誤差表示的殘余誤差圖像。

評價模型用殘差εｉλ或均方根誤差ＲＭＳＥ表示：

ＲＭＳＥ＝［■（εｉλ）２／ｎ］１／２ （４）

２．２端元選取

在ＬＳＭＭ模型中，端元（Ｅｎｄ－ｍｅｍｂｅｒ）的選取是成功進(jìn)行光譜混合分析的關(guān)鍵。不當(dāng)?shù)亩嗽x取將導(dǎo)致噪聲、不穩(wěn)定解和不正確的豐度圖像［９］。

端元組分光譜的確定有多種途徑，可以從地物光譜數(shù)據(jù)庫中選取，從圖像自身的像元光譜中獲得或者從外部數(shù)據(jù)源中提取端元等［１０］?？λ固厣絽^(qū)海拔較高，常年多云雨難以及時獲取高質(zhì)量的遙感影像數(shù)據(jù)，影像上相同物質(zhì)的光譜曲線與實驗室光譜曲線匹配難度較大。而從影像自身選出的組分更易于獲得且具有與數(shù)據(jù)相同的度量尺度，同時參考對地面調(diào)查所獲得的草地光譜信息的分析。本研究選擇通過分析影像的光譜特征并結(jié)合地面調(diào)查光譜數(shù)據(jù)的方法來進(jìn)行純凈像元的提取。首先，對多光譜影像進(jìn)行最小噪聲分離變換（ＭＮＦ）分離噪聲，對基于ＭＮＦ變換結(jié)果的前3個分量進(jìn)行純凈像元指數(shù)（ＰＰＩ）計算，將可能是純凈像元的像元集導(dǎo)入Ｎ維可視化分析器進(jìn)行分析，參考地面光譜數(shù)據(jù)，確定純凈像元，最終將各類端元對應(yīng)像元的平均反射率作為該類端元的反射率。通過光譜分析，并參考喀斯特山區(qū)土地覆蓋特點(diǎn)，最終選擇的端元分別為：①草地端元②林地端元③裸土端元④水體等低反射率端元⑤居民點(diǎn)、裸巖等高反射率端元。各類端元的光譜反射率特征曲線如圖１所示。

２．３線性光譜混合分解結(jié)果

基本組分端元及其光譜特征確定以后，基于獲取的端元光譜特征，采用全受限的ＬＳＭＭ對Landsat ＥＴＭ＋影像進(jìn)行分解，得到草地分量圖。圖２為光譜混合分析的圖像結(jié)果，圖２ａ為草地端元分量圖。在光譜分解結(jié)果中，白色或較亮顏色表示像元中該端元組分比例較高，相反，黑色或較暗顏色則表示像元中該端元組分的比例較低。均方根誤差ＲＭＳＥ是評價ＬＳＭＭ 模型的一個重要指標(biāo)，均方根誤差越小分類計算的準(zhǔn)確性就越高。圖２ｆ為應(yīng)用線性模型分解求解過程產(chǎn)生的各像元ＲＭＳＥ圖像，其中最大的均方根誤差僅有０．００５，小于０．００３的均方根誤差占到了９６．２％，表明此次光譜混合分解較成功［１１］。

３草地地上生物量遙感估測結(jié)果

３．１建立估測模型

本研究于２０１１年９月初至９月中旬在研究區(qū)開展草地資源地面調(diào)查工作，地面調(diào)查時間與影像獲取時間基本同步，共布設(shè)８０塊樣地，考慮到之后估測模型的精度驗證，１６塊樣地數(shù)據(jù)不參與建模。利用ＳＰＳＳ１９．０軟件對研究區(qū)６４個樣地調(diào)查所得草地植被地上生物量數(shù)據(jù)與線性分解得到的草地分量進(jìn)行相關(guān)分析發(fā)現(xiàn)，二者相關(guān)系數(shù)為０．８９６，呈顯著相關(guān)關(guān)系，因而可通過對得到的草地分量進(jìn)行變換得到草地地上生物量。

為了進(jìn)一步分析草地地上生物量與草地分量之間的定量關(guān)系，采用一元線性回歸模型、指數(shù)函數(shù)、對數(shù)函數(shù)、冪函數(shù)等對二者進(jìn)行模擬預(yù)測分析。對不同模型測算結(jié)果的比較與分析結(jié)果表明，基于草地分量的一元線性回歸模型擬合精度最高?？λ固厣絽^(qū)草地地上生物量最優(yōu)擬合模型為：

■＝２ ４８１．１５８ｘ－２８０．４６１（r２＝０．８０１ ２）（５）

式中，ｙ為草地秋季當(dāng)年地上生物量（ｇ／ｍ２），ｘ為草地分量。二者的線性擬合圖見圖3。

３．２模型精度驗證

為評價草地地上生物量估測模型的應(yīng)用精度，通過野外同期采樣的另外１６個樣地數(shù)據(jù)與模型估算的數(shù)據(jù)進(jìn)行比較，并利用誤差統(tǒng)計公式［公式（６）］對回歸模型的精度進(jìn)行綜合評價［１２］。

Ｆ＝（Ｌ■－Ｌ■／Ｌ■）×１００%（６）

式中Ｆ為相對誤差，Ｌ■和Ｌ■分別代表地面實測的草地地上生物量和經(jīng)模型反演得到的數(shù)值。結(jié)果（表１）表明，模擬值同實測值之間的相對誤差為５．３６％～２４．６２％，平均相對誤差為１３．２９％，模型的總體精度達(dá)到８５％以上，說明利用線性光譜分解提取的草地分量建立的一元線性回歸模型能較好地反映草地地上植被生物量狀況，能夠滿足喀斯特山區(qū)中尺度的草地地上生物量估算。

４結(jié)論

受喀斯特山區(qū)地形破碎、環(huán)境差異大、遙感影像空間分辨率和混合像元的限制等因素影響，利用傳統(tǒng)的植被指數(shù)方法對喀斯特山區(qū)草地地上生物量進(jìn)行遙感估測誤差明顯。將光譜混合分析技術(shù)應(yīng)用于喀斯特山區(qū)草地地上生物量的估測研究，通過線性光譜分解，得到符合研究區(qū)環(huán)境特征的端元光譜和各分量的豐度，利用得到的草地分量與地面樣地調(diào)查所得草地地上生物量數(shù)據(jù)進(jìn)行回歸，建立了草地地上生物量估測的線性回歸模型■＝２４８１．１５８ｘ－２８５０．４６１（ｒ２＝０．８０１ ２），基于地面同期調(diào)查數(shù)據(jù)的精度驗證也顯示利用光譜混合模型進(jìn)行喀斯特山區(qū)草地地上生物量估測的精度可達(dá)８５％以上，能夠滿足喀斯特山區(qū)的草地地上生物量的估算。因此，基于光譜混合分析對喀斯特山區(qū)草地地上生物量進(jìn)行遙感估測是可行的。

參考文獻(xiàn)：

［１］ 董永平，吳新宏，戎郁萍，等．草原遙感監(jiān)測技術(shù)［Ｍ］．北京：化學(xué)工業(yè)出版社，２００５．

［２］ 查勇，ＧＡＯ Ｊ，倪紹祥．國際草地資源遙感研究新進(jìn)展［Ｊ］．地理科學(xué)進(jìn)展，２００３，２２（６）：６０７－６１７．

［３］ 徐斌，楊秀春，陶偉國，等．中國草原產(chǎn)草量遙感監(jiān)測［Ｊ］．生態(tài)學(xué)報，２００７，２７（２）：４０５－４１３．

［４］ 楊秀春，徐斌，朱曉華，等．北方農(nóng)牧交錯帶草原產(chǎn)草量遙感監(jiān)測模型［Ｊ］．地理研究，２００７，２６（２）：２１３－２２１．

［５］ 王靜，郭鈮，王振國，等．甘南草地地上部生物量遙感監(jiān)測模型［Ｊ］．干旱氣象，２０１０，２８（２）：１２８－１３３．

［６］ 熊康寧，陳永畢，陳滸，等．點(diǎn)石成金——貴州石漠化治理技術(shù)與模式［Ｍ］．貴陽：貴州科學(xué)技術(shù)出版社，２０１１．

［７］ 趙英時．遙感應(yīng)用分析原理與方法［Ｍ］．北京：科學(xué)出版社，２００３．

［８］ ＡＳＮＥＲ Ｇ， ＬＯＢＥＬＬ Ｄ． Ａ ｂｉｏｇｅｏｐｈｙｓｉｃａｌ ａｐｐｒｏａｃｈ ｆｏｒ ａｕｔｏｍａｔｅｄ ＳＷＩＲ ｕｎｍｉｘｉｎｇ ｏｆ ｓｏｉｌｓ ａｎｄ ｖｅｇｅｔａｔｉｏｎ［Ｊ］． Ｒｅｍｏｔｅ Ｓｅｎｓｉｎｇ ｏｆ Ｅｎｖｉｒｏｎｍｅｎｔ，２０００，７４（１）：９９－１１２．

［９］ ＧＩＬＡＢＥＲＴ Ｍ Ａ． Ａｎ ａｔｍｏｓｐｈｅｒｉｃ ｃｏｒｒｅｃｔｉｏｎ ｍｅｔｈｏｄ ｆｏｒ ｔｈｅ ａｕｔｏｍａｔｉｃ ｒｅｔｒｉｅｖａｌ ｏｆ ｓｕｒｆａｃｅ ｒｅｆｌｅｃｔａｎｃｅ ｆｒｏｍ ＴＭ ｉｍａｇｅｓ［Ｊ］． Ｉｎｔ Ｊ Ｒｅｍｏｔｅ Ｓｅｎｓｉｎｇ，１９９４，１５（１０）：２０６５－２０８６．

［１０］ 李曉松，李增元，吳波，等．基于光譜混合分析的毛烏素沙地油蒿群落覆蓋度提取［Ｊ］．遙感學(xué)報，２００７，１１（６）：９２３－９３０．

［１１］ 李素，李文正，周建軍，等．遙感影像混合像元分解中的端元選擇方法綜述［Ｊ］．地理與地理信息科學(xué)，２００７，２３（５）：３５－３８．

［１２］ 梁天剛，崔霞，馮琦勝，等．２００１－２００８年甘南牧區(qū)草地地上生物量與載畜量遙感動態(tài)監(jiān)測［Ｊ］．草業(yè)學(xué)報，２００９，１８（６）：１２－２２．