全 文 ：书新疆典型荒漠草地的高光谱特征提取和分析研究
钱育蓉１，２，于炯１，贾振红３，杨峰４，帕力旦·吐尔逊１
（１．新疆大学软件学院，新疆 乌鲁木齐８３０００８；２．南京大学生命科学学院，江苏 南京２１００９３；３．新疆大学信息科学与工程学院，
新疆 乌鲁木齐８３００４６；４．四川农业大学农学院，四川 成都６１１１３０）
摘要：本研究以新疆阜康地区的典型荒漠植被为对象，通过微分和ＳａｖｉｔｚｋｙＧｏｌａｙ平滑滤波处理，从野外高光谱仪
采集的高光谱谱线中，提取了典型荒漠植被的高光谱特征（红边、绿峰、红谷），结果显示，荒漠植被虽然也具有红边
和绿峰等吸收特征，但是枝叶稀疏、覆盖度较低的荒漠植被高光谱受下垫面影响，不具有明显绿峰，红边特征也随
着植被盖度的降低而减弱；由于荒漠植被１０月份进入凋落枯萎时期，其光谱反射率受土壤光谱干扰而数值偏大，
红边位置较５月份呈现出“红移”的规律，ＮＤＶＩ和ＲＶＩ指数也呈现出１０月份低于５月份的现象；不同草地类型、不
同季相的荒漠植被，其红边、绿峰、红谷特征和植被指数均具有显著差异，但是ＮＤＶＩ和ＲＶＩ具有较高的数据相关
性。新疆典型荒漠草地的光谱特征提取和分析研究不但丰富了高光谱信息挖掘手段，还为基于高光谱信息的草地
类型分类研究提供了数据基础和科学依据。
关键词：一阶导数；ＳａｖｉｔｚｋｙＧｏｌａｙ平滑滤波；红移；绿峰特征；植被指数
中图分类号：Ｓ８１２　　文献标识码：Ａ　　文章编号：１００４５７５９（２０１３）０１０１５７１０
　　高光谱遥感的发展是２０世纪末的最后２０年中人类在对地观测方面所取得的重大技术突破之一，与传统的
常规遥感相比，高光谱仪可以把光谱分离成几十甚至数百个很窄的波段来接收信息，光谱范围从可见光到热红外
的电磁辐射波谱，所有波段排列在一起能形成一条连续的、完整的光谱曲线［１］。由于其高分辨率等特性在对植被
特征参数进行计算方面大大优于常规多光谱遥感，甚至能完成运用常规遥感方法所不能完成的某些植被特征的
计算，如红边特征、绿峰特征和导数光谱等［２］，因此利用高光谱技术对地表植被无损伤探测的理论与技术得到了
长足的发展，例如：提高遥感信息的信噪比（ＳＮＲ）、改进遥感数据的分析方法、植被生物量的估测等，但目前目标
植被多集中于玉米（犣犲犪犿犪狔狊）、水稻（犗狉狔狕犪狊犪狋犻狏犪）、茶叶等经济作物［３］，针对荒漠植被生化参数的遥感估测研
究还不多见。
草地作为主要陆地生态系统之一，曾几度遭受农垦过牧和不良管理地摧残，在整个生物圈中处于压力的最底
层。任继周［４］利用综合顺序分类法分析了１９５０－２０００年和２００１－２０５０年间的草原类型演替及碳汇动态后，指
出我国潜在草地面积最大的是冻原和高山草地，主要分布在西藏、青海、新疆和四川，其次是冷荒漠草地类型，主
要分布在新疆、内蒙古和甘肃，而且热荒漠草地大类几乎全部分布在新疆。由此可见，新疆草地具有面积广阔、类
型丰富等特点，其荒漠草地类型极具典型性。另一方面，草地资源是重要的国土资源，也是有生命的可更新自然
资源，具有极其重要的生态、经济和社会价值［５］。然而，由于自然社会等因素的长期作用，新疆草地生态严重失
衡，干旱、风沙、盐碱等环境问题日趋严峻，环境质量逐渐恶化，草地生产力持续下降，畜牧业的可持续发展受到严
重影响［６］。
新疆地处我国西北干旱半干旱区域，荒漠植被覆盖度较低，大面积土壤裸露，植被下垫面对野外植被的光谱
特征均有不同程度的影响，消除大气、土壤背景、凋落物等低频光谱成分对植被光谱的影响，可以突出目标植被的
光谱特征，有利于揭示荒漠植被的组织、结构，叶面积指数、蛋白质含量等［７］。本研究以新疆阜康地区的典型荒漠
植被为对象，依据不同海拔划分了３个样区、２３个样带，利用野外高光谱仪采集的典型荒漠植被野外高光谱数
第２２卷　第１期
Ｖｏｌ．２２，Ｎｏ．１
草　业　学　报
ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ 　　
１５７－１６６
２０１３年２月
收稿日期：２０１１０８１７；改回日期：２０１１１０１２
基金项目：新疆高校科研计划项目（ＸＪＥＤＵ２０１２Ｉ１０），新疆大学博士启动基金（ＢＳ１００１２８）和国家自然科学基金项目（６１０６３０４２和６１２６２０８８）资
助。
作者简介：钱育蓉（１９８０），山东武城人，副教授，博士。Ｅｍａｉｌ：ｑｙｒ＠ｘｊｕ．ｅｄｕ．ｃｎ
通讯作者。Ｅｍａｉｌ：ｙｕｊｉｏｎｇ＠ｘｊｕ．ｅｄｕ．ｃｎ；ｊｚｈｈ＠ｘｊｕ．ｅｄｕ．ｃｎ
据，对荒漠植被的高光谱数据进行了特征参量法和植
图１　研究区及采样样区示意图
犉犻犵．１　犜犺犲狅狏犲狉狏犻犲狑狅犳狊犪犿狆犾犻狀犵狕狅狀犲狑犻狋犺犻狀狊狋狌犱狔犪狉犲犪
被指数法的分析，提取出典型植被的光谱内在特征，为
基于高光谱的植被分类提供了理论依据，也为荒漠草
地的监测及“三化”防治奠定了科学根据。
１　材料与方法
１．１　研究区概况
阜康市位于新疆维吾尔自治区中北部，天山东段
北麓，准噶尔盆地南缘，地理坐标为北纬４３°４５′～
４５°３０′、东经８７°４６′～８８°４４′。阜康东界吉木萨尔县，
西与米泉市接壤，南至博格达峰与乌鲁木齐市相连，北
部伸入准噶尔盆地与富蕴县毗邻，市域东西相距７６
ｋｍ，南北绵长１９８ｋｍ，行政区总面积８８４４．８５ｋｍ２，
总人口１６．９４万人，是汉、回、维、哈等２０多个民族的
聚居地（图１）。
研究区地势自东南向西北缓缓倾斜，海拔高程自
５４４５ｍ至４５０ｍ，地势起伏大，其地貌南部为天山支
脉博格达山、中部为山前冲积平原、北部大部分为古尔
班通古特沙漠，从山区过渡为平原再至沙漠（图１），构
成典型的干旱半干旱的自然景观。冬季寒冷漫长，夏
季炎热，春秋季节不明显，且春秋季气温变化剧烈，降
水量分布不均，山区多，平原次之，沙漠最少。春夏多
风，光照充足，热量丰富。年均气温在平原区为６．７℃，
在山区为２．５４℃。年均降水量２０５ｍｍ，南部山区年
平均降水量５３０ｍｍ，中部平原区年平均降水量１８７．５
ｍｍ，北部沙漠边缘年平均降水量１３２ｍｍ，降水的年季变化和季节变化大，年均无霜期１７４ｄ。
１．２　高光谱数据
本研究使用由美国分析光谱仪器公司（ＡＳＤ）生产的便携式ＦｉｅｌｄＳｐｅｃＲＦＲ野外光谱辐射仪采集野外植被高
光谱数据，其光谱范围覆盖３８４．７～１０７５．５ｎｍ，共５１２个波段，分辨率为１．３～１．５ｎｍ。野外光谱测量的时间分
别为２０１０年５月和１０月，光谱采集尽量在晴朗、干燥、微风、无云的天气下进行，同时为了减少不同太阳高度角
对反射率影响，测定时间为北京时间１２：００～１３：００（太阳高度角为５０°～５５°）。光谱仪视场角（ＦＯＶ）１５°，测量中
始终保持探头铅垂向下，距冠层顶垂直高度约１ｍ，地面视场范围直径约为０．５ｍ，测量条件：温度１０～２４℃（５
月份）和３～１３℃（１０月份）。每个观测点记录２７个采样光谱，即在某类植被的１个位置正上方测３次，旋转探头
约１２０°后测３次，然后再旋转探头约１２０°测３次，取平均值作为１个重复的光谱反射率值，最后以３个重复的光
谱反射率的平均值作为处理的光谱反射率值。在室外太阳光下测量植被光谱反射率时，反射率大小与太阳高度
角有关，并且由于植被冠层不是绝对的漫反射体，而校正白板是标准漫反射体，测量时需要及时进行标准白板校
正，即每测完１个处理，用标准白板校正１次（标准白板反射率为１，因此所得目标物光谱为无量纲的相对反射率）。
荒漠植被的研究主要集中于７种当地典型植被：生长在沙漠边缘和农牧交错带的琵琶柴（犚犲犪狌犿狌狉犻犪狊狅狅狀
犵狅狉犻犮犪）、柽柳（犜犪犿犪狉犻狓犮犺犻狀犲狀狊犻狊）、矮梭梭（犎犪犾狅狓狔犾狅狀犪犿犿狅犱犲狀犱狉狅狀）、水柏枝（犕狔狉犻犮犪狉犻犪犫狉犪犮狋犲犪狋犪）和木地肤
（犓狅犮犺犻犪狆狉狅狊狋狉犪狋犪）（图２），以及生长在低山带和泉水溢出带上的芨芨草（犃犮犺狀犪狋犺犲狉狌犿狊狆犾犲狀犱犲狀狊）和锦鸡儿（犆犪
狉犪犵犪狀犪犫狅犵犱犪狊犺犪狀犻犮犪）（图３）。研究仅测高度一般在０．２～１．４ｍ的禾草和灌木层，并在生物量较大的月份（５－６
月）测量了不同空间分布的典型荒漠植被的光谱值。２０１０年５月底获得研究区内不同空间分布下７种荒漠植被
的高光谱有效数据２１６条，２０１０年１０月获得其中荒漠植被的高光谱有效数据６６条。
８５１ ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ（２０１３） Ｖｏｌ．２２，Ｎｏ．１
图２　５种荒漠植被的高光谱谱线
犉犻犵．２　犜犺犲犺狔狆犲狉狊狆犲犮狋狉犪犾犮狌狉狏犲狅犳犳犻狏犲犱犲狊犲狉狋狏犲犵犲狋犪狋犻狅狀
　“植被名狀”表明在图１中样点狀上所测植被 “Ｖｅｇｅｔａｔｉｏｎｎａｍｅ狀”ｍｅａｎｓｔｈｅｍｅａｓｕｒｅｄ
ｖｅｇｅｔａｔｉｏｎｉｎｓａｍｐｌｅａｒｅａ狀ｆｒｏｍｆｉｇ．１．下同Ｔｈｅｓａｍｅｂｅｌｏｗ．
２　基于高光谱的植被特征提取
２．１　获取光谱反射率
高光谱原始数据的校正由ＡＳＤ光谱仪自带软件包自动完成，主要包括：原始数据辐射强度中暗电流去除，各
光谱数据通道的增益校准，波长订正等。获取光谱反射率的基本原理是认为在极短的时间内测定的植被和校正
白板的辐射强度之比即为两者的反射率之比，则植被的反射率为：
犚犼＝犚犼０×犐犼／犐犼０，犼＝１，２，……，狀
式中，狀为通道数或波段数，此处为５１２，犼为波段号，犚犼 为植被在犼波段的反射率，犚犼０ 是校正白板在犼波段的反
射率，不考虑二向反射时应是校正白板的半球反射率，犐犼和犐犼０ 则分别是植被和校正白板在犼波段响应的辐射强
度。
２．２　一阶导数光谱
导数光谱也被称为光谱的微分技术，它可以迅速地确定光谱弯曲点及最大最小反射率的波长位置，能部分消
９５１第２２卷第１期 草业学报２０１３年
除背景、光照等因素的影响，这对于植被稀疏、混合光谱成分多的干旱半干旱地区来说，无疑是一种消除背景噪
声、提取光谱特征的好方法［８］。相对于植被的光谱曲线，土壤的光谱曲线是接近线性的，其导数曲线值整体上很
小，而植被由于光合色素产生的效应，其导数值在７００～７５０ｎｍ附近为一尖峰，这一位置导数值较大，同时土壤
光谱在这个范围内的一阶导数值接近于０。利用这一特点就可用来对混合光谱中的土壤信号进行压缩，因此从
理论上讲，导数可以非常好地消除土壤背景信号。
另外，采用导数光谱技术可以消除光谱数据之间的系统误差、减弱大气辐射、散射和吸收对目标光谱的影响，
以便提取可识别地物的光谱吸收峰参数（波长位置、深度、宽度和吸收光谱指数等）；光谱的一阶导数定义和计算
可以参照钱育蓉等［９］归纳的高光谱特征变量的定义及描述，图４是对图２中梭梭和琵琶柴在不同盖度和土壤背
景下所采集的土壤光谱一阶导数的平滑光谱。
结果（图４）表明，１）光谱曲线两端的噪声比较大，有效光谱波段大致范围为５００～９００ｎｍ；２）植被盖度间差
异在红边幅值上有较明显显示，在琵琶柴和梭梭的红边位置（７１０～７２０ｎｍ）上，低盖度植被红边幅值较低，高盖
度植被红边幅值较高。
图３　２种低山荒漠植被的高光谱谱线
犉犻犵．３　犜犺犲犺狔狆犲狉狊狆犲犮狋狉犪犾犮狌狉狏犲狅犳狋狑狅犾狅狑犿狅狌狀狋犪犻狀犱犲狊犲狉狋狏犲犵犲狋犪狋犻狅狀狊
２．３　光谱平滑处理
实测的光谱数据包含３种误差：１）由于光谱仪是由３个不同波段范围的探头组成，它们之间会出现台阶；２）
随机误差，导致反射率光谱曲线出现一些小抖动的现象；３）水气吸收带的杂乱跳动。这些误差不仅影响曲线的美
观，更主要的是影响即将要提取的光谱特征值。因此，必须减少和限制这种随机噪声水平，对反射率曲线进行光
滑处理消除光谱曲线的抖动、杂乱的水气吸收带以及“台阶”跳跃现象，使光谱曲线变得光滑、连续［１０］。图５ａ为７
种典型荒漠植被在１０月份的地面实测光谱的一阶导数，可以看出一阶导数对噪声较敏感，红边幅值处出现的抖
动造成了部分导数光谱呈现红边“双峰”的现象（如芨芨草）。鉴于光谱导数技术对光谱信噪比非常敏感，本研究
采用基于光谱导数数据常用的ＳａｖｉｔｚｋｙＧｏｌａｙ平滑处理方法，处理结果如图５ｂ、５ｃ和５ｄ。
　　ＳａｖｉｔｚｋｙＧｏｌａｙ平滑滤波因为可保留光谱的一些细微特征（如光谱吸收峰），因此在光谱学中应用非常广泛，
它是对数据窗口（包含奇数个点，２犿＋１）中的数据点进行狀次多形式的拟合，应用最小二乘法法则，解方程组求
得数据点的拟合方程，把经过拟合处理之后的窗口中间点的拟合值代替原始数据，窗口每次向前移动一个点直至
信号点结束，从而实现平滑处理。本研究采用 Ｍａｔｌａｂ工具包进行了ＳａｖｉｔｚｋｙＧｏｌａｙ平滑滤波处理，图５ｂ、５ｃ、５ｄ
是对实验信号图５ａ进行窗口参数分别为“５”、“７”、“９”的ＳａｖｉｔｚｋｙＧｏｌａｙ平滑以后的结果。对比处理结果可以看
出，该方法能有效的去除光谱中存在的噪声，并能保留光谱原有光谱特征，是地物实测光谱和高光谱图像预处理
的一种有效手段。
０６１ ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ（２０１３） Ｖｏｌ．２２，Ｎｏ．１
图４　琵琶柴和梭梭的高光谱一阶导数平滑谱线
犉犻犵．４　犜犺犲狊犿狅狅狋犺犲犱犳犻狉狊狋犱犲狉犻狏犪狋犻狏犲狊狆犲犮狋狉犪狅犳犚犲犪狌犿狌狉犻犪犪狀犱犎犪犾狅狓狔犾狅狀
图５　典型荒漠植被高光谱谱线经犛犪狏犻狋狕犽狔犌狅犾犪狔光滑处理后的曲线
犉犻犵．５　犜犺犲犛犪狏犻狋狕犽狔犌狅犾犪狔狊犿狅狅狋犺犲犱犺狔狆犲狉狊狆犲犮狋狉犪犾犮狌狉狏犲狅犳狋狔狆犻犮犪犾犱犲狊犲狉狋狏犲犵犲狋犪狋犻狅狀
　ａ：原始一阶导数曲线Ｆｉｒｓｔｄｅｒｉｖａｔｉｏｎｏｆｏｒｉｇｉｎａｌｒｅｆｌｅｃｔａｎｃｅ；ｂ：经参数为“５”的ＳａｖｉｔｚｋｙＧｏｌａｙ光滑处理后的曲线ＦｉｒｓｔｄｅｒｉｖａｔｉｏｎｏｆＳａｖｉｔｚｋｙＧｏ
ｌａｙｓｍｏｏｔｈｅｄｒｅｆｌｅｃｔａｎｃｅｗｉｔｈｐａｒａｍｅｔｅｒ“５”；ｃ：经参数为“７”的ＳａｖｉｔｚｋｙＧｏｌａｙ光滑处理后的曲线ＦｉｒｓｔｄｅｒｉｖａｔｉｏｎｏｆＳａｖｉｔｚｋｙＧｏｌａｙｓｍｏｏｔｈｅｄｒｅ
ｆｌｅｃｔａｎｃｅｗｉｔｈｐａｒａｍｅｔｅｒ“７”；ｄ：经参数为“９”的ＳａｖｉｔｚｋｙＧｏｌａｙ光滑处理后的曲线ＦｉｒｓｔｄｅｒｉｖａｔｉｏｎｏｆＳａｖｉｔｚｋｙＧｏｌａｙｓｍｏｏｔｈｅｄｒｅｆｌｅｃｔａｎｃｅｗｉｔｈｐａ
ｒａｍｅｔｅｒ“９”．
１６１第２２卷第１期 草业学报２０１３年
３　典型荒漠植被光谱特征分析
３．１　光谱吸收特征提取
植被盖度、组成和形态等因素影响了植物的光谱特征，因此每种草地植被在绿光反射带、红光吸收带和“红
边”高差等方面均存在差异。根据高光谱特征变量对各吸收峰及吸收谷的定义及描述［９］，编写程序求得荒漠草地
７种典型植被在５月和１０月的高光谱吸收特征：红边幅值、红边位置、红谷反射率、红谷位置、绿峰反射率和绿峰
位置，结果如图６所示，图中７种典型植被类型１、２、３、４、５、６、７依次为木地肤、琵琶柴、柽柳、矮梭梭、水柏枝、芨
芨草和锦鸡儿。
图６　７种典型荒漠植被不同季相下的高光谱特征变量
犉犻犵．６　犜犺犲犺狔狆犲狉狊狆犲犮狋狉犪犾犮犺犪狉犪犮狋犲狉犻狊狋犻犮狊狅犳狊犲狏犲狀狋狔狆犻犮犪犾犱犲狊犲狉狋狏犲犵犲狋犪狋犻狅狀狌狀犱犲狉犱犻犳犳犲狉犲狀狋狊犲犪狊狅狀狊
　１：木地肤犓狅犮犺犻犪狆狉狅狊狋狉犪狋犪；２：琵琶柴犚犲犪狌犿狌狉犻犪狊狅狅狀犵狅狉犻犮犪；３：柽柳犜犪犿犪狉犻狓犮犺犻狀犲狀狊犻狊；４：矮梭梭犎犪犾狅狓狔犾狅狀犪犿犿狅犱犲狀犱狉狅狀；５：水柏枝犕狔狉犻犮犪狉犻犪
犫狉犪犮狋犲犪狋犪；６：芨芨草犃犮犺狀犪狋犺犲狉狌犿狊狆犾犲狀犱犲狀狊；７：锦鸡儿犆犪狉犪犵犪狀犪犫狅犵犱犪狊犺犪狀犻犮犪．下同Ｔｈｅｓａｍｅｂｅｌｏｗ．
２６１ ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ（２０１３） Ｖｏｌ．２２，Ｎｏ．１
由于低山草甸草地（锦鸡儿和芨芨草）和荒漠草地类型植被（木地肤、琵琶柴、柽柳、梭梭、水柏枝）在５月和
１０月份具有明显的季相差异，导致其不同季相的红边、红谷、绿峰参数差异较大（图６）。１）不同植被类型的光谱
特征参数差异很大，但是在红边位置上都出现了１０月份较５月份呈现出“红移”的规律（图６ｂ），主要原因是：１０
月份已经进入植被的凋落枯萎时期，而５月份正是植被的返青生长期；２）在１０月份，除琵琶柴和柽柳还可以找
出绿峰位置，其余５种植被均无明显绿锋（图２、图３和图６ｆ）。结合实地调查现场照片，１０月份众多荒漠植被均
为过冬准备，芨芨草呈现枯黄状态，锦鸡儿、木地肤和水柏枝等荒漠植被的地上部分或者凋落或者枯黄，绿色叶片
大量减少，因此无明显绿峰；其对应绿峰反射率偏高，是由于植被稀疏、土壤裸露，植被光谱反射率受到下垫面反
射光谱影响所致，因此该部分数值不具有可参考性；３）１０月份的绿峰反射率和红谷反射率均高于５月份的，其原
因主要是由于１０月份的植被盖度较５月份低，而盐生木质灌木植被下垫面的光谱反射率较高（图２），使得植被
盖度较低的野外植被光谱数据混合了土壤的光谱，导致光谱反射率数值偏高；４）以琵琶柴为例，虽然它的高光谱
形态特征在５月和１０月比较接近（图２），从特征参数上看，反射率特征也没有较大季相差异，但是不同季相的红
边、绿峰和红谷特征具有显著差异（图６ａ、图６ｃ、图６ｅ），位置特征上的差异也较明显（图６ｂ、图６ｄ、图６ｆ），因此利
用荒漠植被的光谱特征参数可以对连续的高光谱谱线降维，并数量化其吸收峰、反射峰，使得利用高光谱数据进
行草地资源的分类研究成为可能［２］；５）野外水柏枝的外观在不同的生长生活条件下呈现出不同颜色，使得光谱
的蓝绿光谱段和红光谱段的光谱反射受到影响，致使其红谷位置和绿峰位置的参数误差较大，这将会导致其吸收
特征参数的不稳定，也会增加基于高光谱的植被分类的难度。
总之，荒漠植被不同于普通的绿色植被，具有许多荒漠植被所特有的属性，如植被多枝少叶、叶片多呈针状或
刺状、植被覆盖度低等的特点，不同的植被类型高光谱特征复杂，就是同一种植被类型在不同的季节也呈现出迥
异的光谱特征。因此，利用７种植被的高光谱数据提取的个别特征参数误差较大，例如水柏枝的红边幅值和红谷
位置误差较大，芨芨草的红谷反射率和绿峰反射率误差也较大。
３．２　光谱指数特征提取
在草地遥感领域，植被指数（ｖｅｇｅｔａｔｉｏｎｉｎｄｅｘ）作为一种遥感手段已广泛应用于植被覆盖度评价、产量估测以
及自然灾害预测预报等方面。李建龙等［１１］、黄敬峰等［１２］于１９９２－１９９４年在新疆乌鲁木齐、阜康和阿勒泰试验
区，利用试验区牧草产量资料和ＮＯＡＡ／ＡＶＨＲＲ的光谱资料，进行了天然草地牧草产量与草地类型的光谱监测
研究，结果显示天然草地牧草产量与光谱植被指数、遥感植被指数关系密切，可用光谱植被指数和遥感植被指数
对天然草地牧草产量进行监测。同时还指出处于干旱区的新疆，由于降水量少植被盖度和牧草产量都低，用
ＮＤＶＩ（ｔｈｅｎｏｒｍａｌｉｚｅｄｄｉｆｆｅｒｅｎｃｅｖｅｇｅｔａｔｉｏｎｉｎｄｅｘ，归一化植被指数）和ＲＶＩ（ｒａｔｉｏｖｅｇｅｔａｔｉｏｎｉｎｄｅｘ，比值植被指
数）建立的监测模型没有显著差异，一般情况下线性模型足以反映牧草产量的动态变化，自变量选用 ＮＤＶＩ和
ＲＶＩ均可［１１１４］。
植被的光谱指数是根据植被反射波段的特性计算出来的，对地表植被生长、覆盖状况、地表生物参数具有一
定指示意义的各种数值，广泛应用于植被的时空动态监测以及一些生物参数的估算，容易获取且计算方便。可见
光波段与近红外波段数值的不同形式组合，是植被指数的核心。遥感估算植被覆盖研究中最常用的植被指数为
ＮＤＶＩ和ＲＶＩ，它们对绿色植被变化灵敏，对土壤或者枯草不灵敏，能够较好地反映出草地的覆盖度、生物量和叶
面积指数变化，与地上现存净初级生物量有较好的相关性。
ＮＤＶＩ被定义为近红外波段（ＮＩＲ）与可见光红波段（Ｒ）数值之差和这２个波段数值之和的比值，即ＮＤＶＩ＝
（ＴＭ４－ＴＭ３）／（ＴＭ４＋ＴＭ３），它是目前表征植被覆盖特征应用最广泛的指标，同时也被广泛地应用到初级生
产力（ＮＰＰ）、植被覆盖、生物量以及物候等研究领域；ＲＶＩ被定义为近红外波段（ＮＩＲ）与可见光红波段（Ｒ）数值
的比值，即ＲＶＩ＝ＮＩＲ／Ｒ，它是绿色植物的灵敏指示参数，与叶面积指数、叶干生物量、叶绿素含量相关性高，可
用于检测和估算植物生物量。其中ＮＩＲ为７６０～９００ｎｍ近红外波段的反射，相当于ＴＭ／ＥＴＭ＋遥感影像的第
４通道接受的地物光谱辐射能量，是植被高反射率波段，便于进行植被种类的识别；Ｒ为６３０～６９０ｎｍ红光波段
的反射率，相当于ＴＭ／ＥＴＭ＋遥感影像的第３通道接受的地物光谱辐射能量，该波段对植物覆盖度、植物生长
状况敏感。
３６１第２２卷第１期 草业学报２０１３年
图７　７种典型荒漠植被不同季相下的植被指数
犉犻犵．７　犜犺犲狏犲犵犲狋犪狋犻狅狀犻狀犱犲狓狅犳狊犲狏犲狀狋狔狆犻犮犪犾犱犲狊犲狉狋狏犲犵犲狋犪狋犻狅狀狊狌狀犱犲狉犱犻犳犳犲狉犲狀狋狊犲犪狊狅狀狊
对比分析７种典型荒漠植被（木地肤、琵琶柴、柽柳、矮梭梭、水柏枝、芨芨草和锦鸡儿）在５月和１０月份的植
被指数，结果（图７）表明，１）不同月份下的植被指数数值不同。除了琵琶柴的植被指数季相间差异不明显，其余
６种典型荒漠植被５月份的ＮＤＶＩ数值和ＲＶＩ数值普遍高于１０月份的。这种季相间的植被指数差异主要是由
于５月份雪山融水和雨水充沛，是荒漠植被的生长季，而１０月份气候干旱、雨量减少，大部分荒漠植被的地上部
分枯黄凋落；２）不同草地类型上的ＮＤＶＩ和ＲＶＩ数值规律变异较大。以５月份为例，ＮＤＶＩ指数呈现出柽柳＞
锦鸡儿＞芨芨草＞水柏枝＞琵琶柴＞木地肤＞梭梭的规律，但是个别类型差异不显著，如水柏枝和芨芨草、木地
肤和琵琶柴。草原化荒漠植被由于植被盖度低，土壤背景对植被的光谱反射影响较大，造成植被覆盖度相近的草
地类型呈现出植被指数相似的现象；３）对比图７中ＮＤＶＩ和ＲＶＩ可以发现：２组数据的规律非常相似，生物量
大、覆盖度高的植被ＮＤＶＩ和ＲＶＩ均较高，反之则较低，且５月份ＮＤＶＩ和ＲＶＩ数值的相关系数高达０．９２５３
（狀＝７２），１０月份的相关系数更高（犚＝０．９７２０，狀＝２２）。
４　结论与讨论
１）新疆典型荒漠植被的高光谱线经微分和平滑处理后，具有红边和绿峰等光谱吸收特征，但受下垫面土壤
光谱噪声影响，秋季中个别植被不具有明显绿峰，红边特征也随着覆盖度的降低而减弱，其中１０月份典型荒漠植
被的红边位置较５月份呈现出“红移”的规律，ＮＤＶＩ和ＲＶＩ指数也呈现出１０月份低于５月份的现象。还有学者
利用高光谱谱线中水分子特征吸收波段进行草地可燃物含水率的估算，并以此快速、便捷、动态的预测草原火险
发生率［１５］，这些研究都充分利用了高光谱谱线中特征波段上所包含的丰富的吸收特征。
２）７种典型荒漠植被在１０月份已经进入凋落枯萎时期，地上部分大多凋落或枯黄，绿色叶片大量减少，反射
率受到土壤光谱干扰而呈现数值偏高的现象，大多无明显绿峰；同时，虽然各种荒漠植被的高光谱谱线在不同季
相中显示出相同的曲线形态特征，但是其红边、绿峰和红谷特征具有显著差异（图６）。因此利用荒漠植被的光谱
特征参数可以对连续的高光谱谱线降维，并数量化其吸收峰、反射峰，使得利用高光谱数据进行草地资源的分类
研究成为可能［２］，还可以识别植被混合群落中的草地退化指示植被［１６］，为基于高光谱的天然草地退化遥感监测
研究提供了科学依据。
３）虽然不同草地类型上的植被指数数值变异较大，但是ＮＤＶＩ和ＲＶＩ具有较高的数据相关性。７种典型荒
漠植被在不同季相的植被指数数值差异较大：ＮＤＶＩ为０．２１～０．７７，ＲＶＩ为１．８～７．９（图７），但是荒漠植被由于
具有植被盖度低，土壤背景对植被的光谱反射影响较大等共性，使得它们均呈现出植被指数相关的现象：春秋两
季荒漠植被ＮＤＶＩ和ＲＶＩ数值的相关系数分别为０．９２５３（狀＝７２）和０．９７２０（狀＝２２）。考虑到基于影像提取的
诸多光谱指数均与草地生物量有较好的定量关系［１７］，那么对荒漠草地的高光谱遥感图像进行特征提取不但可以
及时、准确的掌握草地资源的变化状况，指导人们正确开发利用草地资源、保护草地生态环境，还可以成为定性、
４６１ ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ（２０１３） Ｖｏｌ．２２，Ｎｏ．１
定量探测草地资源和发展精细草业的重要保证。
草地高光谱研究以其快速、便捷、准确等优势受到世界各国学者的普遍关注，在区域和全球尺度上已经广泛
应用于草地植被生化参数的无损伤估算、牧草品质测定、含水率估测、退化草地指示草种识别等领域［１３，１５］。目
前，在我国天然草地植被监管中，高光谱遥感技术的应用还不甚广泛。随着高光谱遥感技术的不断发展和成熟，
各种高光谱遥感卫星影像［１８，１９］（例如：Ｈｙｐｅｒｉｏｎ、ＣＨＲＩＳ／ＰＲＯＢＡ和 ＭＯＤＩＳ）和满足不同需要的应用软件的出
现，高光谱遥感技术会更多地应用于草地植被测定和生态环境监测中［２０］。
伴随着荒漠植被高光谱遥感应用的逐步深化，高光谱数据信息挖掘技术的逐步提高，高光谱遥感可以由实验
研究阶段逐步转向实际应用阶段，并在防治草地退化、加快草地植被恢复与重建、监测草地演替动态、合理保护和
利用草地资源中提供宝贵的信息，进一步实现对干旱半干旱地区植被的类别划分、大面积实时监测，并对荒漠草
地的退化进行及时预警，也必将推动草地管理进入一个全新的阶段。
参考文献：
［１］　杨可明，郭达志．植被高光谱特征分析及其病害信息提取研究［Ｊ］．地理与地理信息科学，２００６，２２（４）：３１３４．
［２］　钱育蓉，贾振红，于炯，等．ＢＰＡＮＮ在荒漠草地高光谱分类研究中的应用［Ｊ］．计算机工程与应用，２０１１，４７（１２）：２２５
２２８．
［３］　赵杰文，王开亮，欧阳琴，等．高光谱技术分析茶树叶片中叶绿素含量及分布［Ｊ］．光谱学与光谱分析，２０１１，３１（２）：５１２
５１５．
［４］　任继周，梁天刚，林慧龙，等．草地对全球气候变化的响应及其碳汇潜势研究［Ｊ］．草业学报，２０１１，２０（２）：１２２．
［５］　刘兴元，龙瑞军，尚占环．草地生态系统服务功能及其价值评估方法研究［Ｊ］．草业学报，２０１１，２０（１）：１６７１７４．
［６］　杨峰，钱育蓉，李建龙，等．天山北坡典型荒漠草地退化特征及其成因［Ｊ］．自然资源学报，２０１１，２６（８）：１３０６１３１４．
［７］　范燕敏，武红旗，靳瑰丽．新疆草地类型高光谱特征分析［Ｊ］．草业科学，２００６，２３（６）：１５１８．
［８］　丁建丽，张飞，塔西甫拉提·特依拜．塔里木盆地南缘典型植被光谱特征分析［Ｊ］．干旱区资源与环境，２００８，２２（１１）：１６０
１６６．
［９］　钱育蓉，杨峰，李建龙，等．利用高光谱数据快速估算高羊茅牧草光合色素的研究［Ｊ］．草业学报，２００９，１８（４）：９４１０２．
［１０］　陈文霞，陈安升，蔡之华．基于高光谱吸收特征参数的分类研究［Ｊ］．计算机工程与应用，２００８，４４（２８）：２３０２３２．
［１１］　李建龙，黄敬峰，王秀珍．草地遥感［Ｍ］．北京：气象出版社，１９９７：１４６１５１．
［１２］　黄敬峰，王秀珍，胡新博．新疆北部不同类型天然草地产草量遥感监测模型［Ｊ］．中国草地，１９９９，１：７１１．
［１３］　黄敬峰，桑长青，金杰．天山北坡天然草场光谱植被指数的基本特征［Ｊ］．遥感技术与应用，１９９４，９（１）：２９３３．
［１４］　李建龙，蒋平．遥感技术在大面积天然草地估产和预报中的应用探讨［Ｊ］．武汉测绘科技大学学报，１９９８，２３（２）：１５３１５７．
［１５］　包玉龙，张继权，赵云升，等．基于高光谱的草地可燃物含水率估测方法［Ｊ］．红外，２０１１，３２（７）：３８４２．
［１６］　王焕炯，范闻捷，崔要奎，等．草地退化的高光谱遥感监测方法［Ｊ］．光谱学与光谱分析，２０１０，３０（１０）：２７３４２７３８．
［１７］　陈鹏飞，王卷乐，廖秀英，等．基于环境减灾卫星遥感数据的呼伦贝尔草地地上生物量反演研究［Ｊ］．自然资源学报，２０１０，
２５（７）：１１２２１１３１．
［１８］　李新辉，宋小宁，冷佩．利用ＣＨＲＩＳ／ＰＲＯＢＡ数据定量反演草地ＬＡＩ方法研究［Ｊ］．国土资源遥感，２０１１，３：６１６６．
［１９］　陈健飞，林征，陈颖彪．基于高光谱线形混合光谱分解识别人工地物［Ｊ］．应用基础与工程科学学报，２００９，１７（２）：２０６
２１８．
［２０］　许冬梅，王．近红外光谱技术在天然草地植被管理中的应用［Ｊ］．光谱学与光谱分析，２００７，２７（１０）：２０１３２０１６．
５６１第２２卷第１期 草业学报２０１３年
犈狓狋狉犪犮狋犻狅狀犪狀犱犪狀犪犾狔狊犻狊狅犳犺狔狆犲狉狊狆犲犮狋狉犪犾犱犪狋犪犳狉狅犿狋狔狆犻犮犪犾犱犲狊犲狉狋犵狉犪狊狊犾犪狀犱犻狀犡犻狀犼犻犪狀犵
ＱＩＡＮＹｕｒｏｎｇ１，２，ＹＵＪｉｏｎｇ１，ＪＩＡＺｈｅｎｈｏｎｇ３，ＹＡＮＧＦｅｎｇ４，ＰａｌｉｄａｎＴｕｅｒｘｕｎ１
（１．ＳｃｈｏｏｌｏｆＳｏｆｔｗａｒｅ，ＸｉｎｊｉａｎｇＵｎｉｖｅｒｓｉｔｙ，Ｕｒｕｍｑｉ８３０００８，Ｃｈｉｎａ；２．ＳｃｈｏｏｌｏｆＬｉｆｅＳｃｉｅｎｃｅ，Ｎａｎｊｉｎｇ
Ｕｎｉｖｅｒｓｉｔｙ，Ｎａｎｊｉｎｇ２１００９３，Ｃｈｉｎａ；３．ＣｏｌｅｇｅｏｆＩｎｆｏｒｍａｔｉｏｎＳｃｉｅｎｃｅａｎｄＥｎｇｉｎｅｅｒｉｎｇ，Ｘｉｎｊｉａｎｇ
Ｕｎｉｖｅｒｓｉｔｙ，Ｕｒｕｍｑｉ８３００４６，Ｃｈｉｎａ；４．ＣｏｌｅｇｅｏｆＡｇｒｉｃｕｌｔｕｒａｌ，ＳｉｃｈｕａｎＡｇｒｉｃｕｌｔｕｒａｌ
Ｕｎｉｖｅｒｓｉｔｙ，Ｃｈｅｎｇｄｕ６１１１３０，Ｃｈｉｎａ）
犃犫狊狋狉犪犮狋：Ｔｈｅｈｙｐｅｒｓｐｅｃｔｒａｌｆｅａｔｕｒｅｓ（ｒｅｄｅｄｇｅ，ｇｒｅｅｎｐｅａｋ，ｒｅｄｖａｌｅｙ）ｏｆｔｙｐｉｃａｌｄｅｓｅｒｔｇｒａｓｓｌａｎｄｓｉｎＸｉｎｊｉａｎｇ
ＦｕｋａｎｇＣｉｔｙｗｅｒｅｅｘｔｒａｃｔｅｄｆｒｏｍｈｙｐｅｒｓｐｅｃｔｒａｌｌｉｎｅｓｔｈｒｏｕｇｈｄｉｆｆｅｒｅｎｔｉａｌａｎｄＳａｖｉｔｚｋｙＧｏｌａｙｓｍｏｏｔｈｉｎｇｆｉｌｔｅｒ
ｐｒｏｃｅｓｓｉｎｇ．Ｔｈｅｄｅｓｅｒｔｖｅｇｅｔａｔｉｏｎｓｈａｖｅａｂｓｏｒｐｔｉｏｎｆｅａｔｕｒｅｓ，ｓｕｃｈａｓｒｅｄｅｄｇｅｓａｎｄｇｒｅｅｎｐｅａｋｓｂｕｔｔｈｅｈｙｐｅｒ
ｓｐｅｃｔｒａｌｃｕｒｖｅｓａｒｅｎｏｔｃｌｅａｒｂｅｃａｕｓｅｏｆｔｈｅｅｆｆｅｃｔｏｆｔｈｅｕｎｄｅｒｌｙｉｎｇｓｕｒｆａｃｅｓｈｏｗｉｎｇｔｈｒｏｕｇｈｔｈｅｓｐａｒｓｅｆｏｌｉａｇｅ
ｏｆｔｈｅｌｏｗｃｏｖｅｒａｇｅｄｅｓｅｒｔｖｅｇｅｔａｔｉｏｎ；ＤｕｅｔｏｔｈｅｗｉｔｈｅｒｅｄｌｉｔｔｅｒｄｅｓｅｒｔｖｅｇｅｔａｔｉｏｎｉｎＯｃｔｏｂｅｒ，ｔｈｅｓｐｅｃｔｒａｌｒｅ
ｆｌｅｃｔａｎｃｅｅｆｆｅｃｔｏｆｓｏｉｌｉｎｔｅｒｆｅｒｅｎｃｅｗａｓｔｏｏｌａｒｇｅ，ａｎｄｃｏｍｐａｒｅｄｗｉｔｈｔｈａｔｉｎＭａｙ，ｔｈｅｒｅｄｅｄｇｅｐｏｓｉｔｉｏｎｓｈｏｗｅｄ
ｔｈｅｐｈｅｎｏｍｅｎｏｎｏｆ“ｒｅｄｅｄｇｅｓｈｉｆｔ”．ＴｈｅＮＤＶＩ（ｔｈｅｎｏｒｍａｌｉｚｅｄｄｉｆｆｅｒｅｎｃｅｖｅｇｅｔａｔｉｏｎｉｎｄｅｘ）ａｎｄＲＶＩ（ｒａｔｉｏ
ｖｅｇｅｔａｔｉｏｎｉｎｄｅｘ）ｉｎｄｅｘｖａｌｕｅｓｗｅｒｅｌｏｗｅｒｉｎＯｃｔｏｂｅｒ．Ｔｈｅｒｅｗｅｒｅｓｉｇｎｉｆｉｃａｎｔｄｉｆｆｅｒｅｎｃｅｓｉｎｒｅｄｅｄｇｅ，ｇｒｅｅｎ
ｐｅａｋ，ｒｅｄｐｅａｋａｎｄｖｅｇｅｔａｔｉｏｎｉｎｄｅｘｃｈａｒａｃｔｅｒｉｓｔｉｃｓｂｅｔｗｅｅｎｇｒａｓｓｌａｎｄｔｙｐｅｓ，ｄｉｆｆｅｒｅｎｔｓｅａｓｏｎａｌｄｅｓｅｒｔｖｅｇｅｔａ
ｔｉｏｎｓ，ｂｕｔＮＤＶＩａｎｄＲＶＩｗｅｒｅａｓｓｏｃｉａｔｅｄｗｉｔｈｓｔｒｏｎｇｃｏｒｒｅｌａｔｉｏｎｓ．Ｆｉｎａｌｙ，ｔｈｅｆｅａｔｕｒｅｅｘｔｒａｃｔｉｏｎａｎｄａｎａｌｙｓｉｓ
ｏｆｈｙｐｅｒｓｐｅｃｔｒａｌｄａｔａｆｒｏｍｔｙｐｉｃａｌｄｅｓｅｒｔｇｒａｓｓｌａｎｄｓｉｎＸｉｎｊｉａｎｇｎｏｔｏｎｌｙｅｎｒｉｃｈｅｓｔｈｅｍｅａｎｓｏｆｈｙｐｅｒｓｐｅｃｔｒａｌ
ｄａｔａｍｉｎｉｎｇ，ｂｕｔａｌｓｏｐｒｏｖｉｄｅｓａｓｃｉｅｎｔｉｆｉｃｂａｓｉｓｆｏｒｇｒａｓｓｌａｎｄｃｌａｓｓｉｆｉｃａｔｉｏｎｂａｓｅｄｏｎｈｙｐｅｒｓｐｅｃｔｒａｌｄａｔａ．
犓犲狔狑狅狉犱狊：ｆｉｒｓｔｄｅｒｉｖａｔｉｖｅ；ＳａｖｉｔｚｋｙＧｏｌａｙｓｍｏｏｔｈｉｎｇ；ｒｅｄｅｄｇｅｓｈｉｆｔ；ｇｒｅｅｎｐｅａｋｆｅａｔｕｒｅ；ｖｅｇｅｔａｔｉｏｎｉｎｄｅｘ
６６１ ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ（２０１３） Ｖｏｌ．２２，Ｎｏ．１