摘要：由于喀斯特山区地表破碎，环境异质性强，影像上混合像元现象严重，增加了草地地上生物量遥感估测的难度。以ＬａｎｄｓａｔＥＴＭ＋为数据源，将光谱混合分析技术应用于喀斯特山区草地地上生物量的遥感估测中，分析了基于线性混合分解模型分解的草地分量与实测草地地上生物量的相关关系，进而构建了喀斯特山区草地地上生物量估测模型为■＝２４５１．１５８ｘ－２８０．４６１（ｒ２＝０．８０１２），验证结果总体精度可达８５％以上。基于光谱混合分析建立的草地地上生物量估测模型提高了喀斯特山区草地地上生物量遥感估测的精度，为喀斯特山区草地生长状况调查与监测提供了技术支持。

 

　　关键词：喀斯特山区；草地地上生物量；遥感估测；光谱混合分析

 

　　草地是中国南方喀斯特山区脆弱生态系统的重要组成部分，是维持喀斯特山区生态平衡的最后一道屏障，同时也是喀斯特山区生态环境保护与发展草地畜牧业的基础。草地植被地上生物量能客观反映草地固碳能力与草地载畜能力的大小，是衡量草地资源的重要指标，也是进行草业规划、生态评价、草原退化监测的基础［１］。开展喀斯特山区草地地上生物量的遥感估测，对保护喀斯特山区脆弱生态环境与实现草地畜牧业的可持续发展，以及推动定量遥感在喀斯特复杂背景下的应用均具有重要的科学及现实意义。综合国内外已有研究［２－５］可以发现，草地资源遥感估测的研究多数都是通过对遥感数据的波段组合或利用遥感数据计算出的各种植被指数与草地地上生物量进行相关分析，建立经验模型，并推广模型以求取大范围区域的草地地上生物量。然而，在自然条件复杂的喀斯特山区，受遥感影像空间分辨率影响，图像中的一个像元往往由多个地物组分构成，内部空间异质性很强，草地分布破碎，其植被指数受岩石、土壤影响强烈，基于植被指数对喀斯特山区草地地上生物量进行遥感估测误差明显。光谱混合分析法能够较好地解决混合像元的问题，将其引入草地地上生物量的估测，可以有效提高定量遥感的精度。

 

　　研究选择贵州省黔南布依族苗族自治州的独山县兔场镇作为研究区域。独山县兔场镇石漠化退耕还林（草）与草地畜牧业综合治理小区位于东经１０７°５１′－１０７°６５′，北纬２５°９１′－２６°０１′之间，地处贵州高原向广西丘陵盆地过渡的斜坡地带，境内岩溶地貌发育，地貌类型复杂多样，以山地为主，其间交错分布有少量的丘陵、坝地，面积１０６．２ｋｍ２，平均海拔１０７５ｍ，属于中亚热带湿润季风气候区，年平均气温１３．６℃，年均降水量１４６８ｍｍ，适宜天然草地生长［６］。选择ＬａｎｄｓａｔＥＴＭ＋影像为数据源，将光谱混合分析技术应用到喀斯特山区草地地上生物量的估测中，分析了基于线性混合分解模型分解的草地分量与实测草地地上生物量的相关关系，进而建立了喀斯特山区草地地上生物量估测模型。

 

　　１数据获取与预处理

 

　　１．１影像数据获取与预处理

 

　　选用的遥感影像为２０１１年９月２５日获取的ＬａｎｄｓａｔＥＴＭ＋影像，轨道号为１２６／４２，太阳高度角５５．５６°，云量较少，影像质量较好。使用ＥＮＶＩ４．８软件对影像进行去条带处理之后，以１∶５００００地形图为基准，均匀选取控制点，利用多项式模型对影像进行几何精校正，校正后图像均方根误差在０．４５个像元以内，满足精度要求。精校正后的影像经辐射定标后，采用ＬａｎｄｓａｔＦＬＡＡＳＨ进行大气校正，得到除热红外波段后的其他５个波段的地表反射率［７］。

 

　　１．２地面数据采集

 

　　研究于２０１１年９月初至9月中旬在研究区开展草地资源地面调查工作，考虑到喀斯特山区草地分布的生境条件与植被类型，选取具有代表性的面积为２０ｍ×２０ｍ的不同草地类型样区２０个，在其中共布设了８０块样地，样地大小１ｍ×１ｍ，样地均选取在样区内草地植被长势较均一的位置。齐地面刈割样地内草本植物，称取鲜重，统计作为该样地草地地上生物量。样地采用ＧＰＳ定位，经纬度定位数据精度优于２ｍ，高程数据精度优于５ｍ。样地基础信息调查主要录入样地经纬度、海拔、地形、坡度、土壤岩性、草地类型等信息，同时还采集了研究区典型地物的光谱数据，为光谱混合分析提供数据。使用的光谱仪型号为ＧＥＲ１５００，波长范围３５０～１０５０ｎｍ。

 

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

 

　　２．１光谱混合模型的选择

 

　　国内外对遥感光谱混合分析（ＳｐｅｃｔｒａｌＭｉｘｔｕｒｅＡｎａｌｙｓｉｓ，ＳＭＡ）的研究已较成熟，研究和发展了多种混合光谱分解模型，其中线性光谱混合模型（ＬｉｎｅａｒＳｐｅｃｔｒａｌＭｉｘｔｕｒｅＭｏｄｅｌ，ＬＳＭＭ）因具有一定的理论依据、试验验证基础、一定的精度保证和可操作性强而被广泛应用，是最常用的一种模型［８］。

 

　　线性光谱混合模型（ＬＳＳＭ）定义像元在某一光谱波段的反射率（亮度值）是由构成像元的端元的反射率（亮度值）以其所占像元面积比例为权重系数的线性组合。可用以下公式计算：

 

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

 

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

 

　　０≤ｆkｉ≤１（３）

 

　　式中，Ｒｉλ为第λ波段第ｉ像元的光谱反射率；ｆkｉ为对应于ｉ像元的第ｋ个基本组分（端元）所占的分量值；Ｃｋλ为第ｋ个基本组分在第λ波段的光谱反射率；εｉλ为残余误差值（即光谱的非模型化部分）；ｎ为基本组分的数目。公式（２）、（３）为ＬＳＳＭ的限制条件。模型计算的结果表现为各端元丰度图像和以均方根误差表示的残余误差图像。

 

　　评价模型用残差εｉλ或均方根误差ＲＭＳＥ表示：

 

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

 

　　２．２端元选取

 

　　在ＬＳＭＭ模型中，端元（Ｅｎｄ－ｍｅｍｂｅｒ）的选取是成功进行光谱混合分析的关键。不当的端元选取将导致噪声、不稳定解和不正确的丰度图像［９］。

 

　　端元组分光谱的确定有多种途径，可以从地物光谱数据库中选取，从图像自身的像元光谱中获得或者从外部数据源中提取端元等［１０］。论文答辩喀斯特山区海拔较高，常年多云雨难以及时获取高质量的遥感影像数据，影像上相同物质的光谱曲线与实验室光谱曲线匹配难度较大。而从影像自身选出的组分更易于获得且具有与数据相同的度量尺度，同时参考对地面调查所获得的草地光谱信息的分析。本研究选择通过分析影像的光谱特征并结合地面调查光谱数据的方法来进行纯净像元的提取。首先，对多光谱影像进行最小噪声分离变换（ＭＮＦ）分离噪声，对基于ＭＮＦ变换结果的前3个分量进行纯净像元指数（ＰＰＩ）计算，将可能是纯净像元的像元集导入Ｎ维可视化分析器进行分析，参考地面光谱数据，确定纯净像元，最终将各类端元对应像元的平均反射率作为该类端元的反射率。通过光谱分析，并参考喀斯特山区土地覆盖特点，最终选择的端元分别为：①草地端元②林地端元③裸土端元④水体等低反射率端元⑤居民点、裸岩等高反射率端元。各类端元的光谱反射率特征曲线如图１所示。

 

　　２．３线性光谱混合分解结果

 

　　基本组分端元及其光谱特征确定以后，基于获取的端元光谱特征，采用全受限的ＬＳＭＭ对LandsatＥＴＭ＋影像进行分解，得到草地分量图。图２为光谱混合分析的图像结果，图２ａ为草地端元分量图。在光谱分解结果中，白色或较亮颜色表示像元中该端元组分比例较高，相反，黑色或较暗颜色则表示像元中该端元组分的比例较低。均方根误差ＲＭＳＥ是评价ＬＳＭＭ模型的一个重要指标，均方根误差越小分类计算的准确性就越高。图２ｆ为应用线性模型分解求解过程产生的各像元ＲＭＳＥ图像，其中最大的均方根误差仅有０．００５，小于０．００３的均方根误差占到了９６．２％，表明此次光谱混合分解较成功［１１］。