1、植被光谱分析与植被指数计算在遥感中，常常结合不同波长范围的反射率来增强植被特征，如植被指数（vegetationindicesVI）的计算，植被指数（VI）是两个或多个波长范围内的地物反射率组合运算，以增强植被某一特性或者细节。目前，在科学文献中发布了超过150种植被指数模型，这些植被指数中只有极少数是经过系统的实践检验。本文总结现有植被指数，根据对植被波谱特征产生重要影响的主要化学成份：色素（Pigments）、水分（Water）、碳（Carbon）、氮（Nitrogen）,总结了7大类实用性较强的植被指数，即：宽带绿度、窄带绿度、光利用率、冠层氮、干旱或碳衰减、叶色素、冠层水分含量。这些植
3、四个部分：可见光（Visible）:400nmto700nm近红外（Near-infraredNR:700nmto1300nm短波红外1（Shortwaveinfrared1SWIR-1:1300nmto1900nm短波红外2（Shortwaveinfrared2SWIR-2:1900nmto2500nm其中NIR和SWIR-1的过渡区（1400nm附近）是大气水的强吸收范围，卫星或者航空传感器一般不获取这范围的反射值。SWIR-1和SWIR-2的过渡区（1900nm附近）也是大气水的强吸收范围。植被可分为三个部分组成：植物叶片（PlantFoliage）植被冠层（PlantCanopies非
4、光合作用植被（Non-PhotosyntheticVegetation）这三个部分是植被分析的基础，下面对他们详细介绍。植物叶片（PlantFoliage植物叶片包括叶、叶柄以及其他绿色物质，不同种类的叶片具有不同的形状和化学成份。对波谱特征产生重要影响的主要化学成份包括：色素（Pigments）、水分（Water）、碳（Carbon）、氮（Nitrogen）,这也是遥感反演的基础，如用植被指数来估算叶子的化学成份。色素（Pigments叶色素主要包括叶绿素、叶黄素和花青素。这些都是植被的健康的指标，比如含高浓度叶绿素的植被一般很健康，相反，叶黄素和花青素常常出现在健康较差的植被，濒临死亡的植
5、被出现红色、黄色或棕色。叶色素只影响可见光部分（400nm700nm）,图1为几种叶色素在可见光范围的相对光谱吸收特征。-020.1O0.504004006O55005506Wavelength（nm）图i部分叶色素的相对光谱吸收特征水分（Water）叶子的几何特性、冠层结构和对水的需求影响植被的水分含量。水分对植被反射率的影响波段范围在NIR和SWIR（图2）。在1400nm和1900nm附近有吸收波谷，但是传感器一般会避开这两个波段范围。在970nm和1190nm附近也有强吸收特征，可利用这两个波段范围监测植被水分。碳（Carbon）植物中的碳是以很多形式存在，包括糖，淀粉，纤维素和木质素
6、等。纤维素和木质素的吸收特征表现在短波光谱范围内容（图3）1000100100.140060080010001200140016001800200022002400Wavelength（nm）图2叶片水和碳（纤维素和木质素）相对光谱吸收特征氮（Nitrogen）叶子中的氮元素一般包含在叶绿素、蛋白质以及其他分子中。植被指数（VI）对包含在叶绿素中的氮元素很敏感（大约含6%氮）。包含在蛋白质中的氮元素在1500nm1720nm范围内对叶片波谱特征影响比较大。从上可以看出，植被与辐射的相互作用主要体现在叶片的波谱特征，因此，在可见光谱段内，主要太阳辐射的吸收来自叶绿素、叶黄素和花青素，形成450n
7、m和670nm附近的吸收谷；在近红外谱段内，主要太阳辐射的吸收来自水分，形成970nm和1190nm两个水吸收带；在短波红外谱段内，除了水分，各种形式存在的碳和氮也对太阳辐射的吸收有一定的贡献，形成1400nm和1900nm吸收谷。图3是叶片反射率与透射光谱（TransmittanceSpectra）对比例子，木本植被和草本植被在色素、水分、氮等含量不一样，反射率与透射光谱关系也不一样。VISNIRSWIR1SWIR206.4008001200160020002400Waveiength（nmt图3木本植物（A）和草本植物（B）的叶片反射率与透射光谱植被冠层（PlantCanopieS单片叶子
8、的反射特性对植被冠层光谱特征是重要的，止匕外，叶子数量和冠层结构对植被冠层的散射、吸收也有重要的影响。比如不同的生态系统中，森林、草原、或农业用地等都具有不同的反射特性，虽然它们单个叶子很类似。有很多植被模型用于描述冠层光谱特征。两个最重要的是叶面积指数（LAI）和叶倾叶角分布（LAD）oLAI指每单位面积地上绿叶面积，这表现了冠层中绿色植被的总数；LAD描述了树叶所有类型的定向，常用平均叶倾角（MLA）近似。MLA表示冠层中的每个叶片角度与水平方向的差值的平均值。图4表示LAI和LAD对植被冠层的影响效果，MLA近似LAD。在近红外谱段内，植被强反射太阳辐射，植被冠层在可见光和SWIR-2是
9、强吸收。使用可见光和SWIR-2的植被指数对上层林冠非常敏感图4LAI(A)和MLA(B)的增减对植被冠层的影响非光合作用植被(Non-PhotosyntheticVegetation)在自然界里，还包括占了全球植被覆盖一半的衰老或死亡植被，它们被称为非光合作用植被(简称NPV)NPV的冠层也具有木本森林结构，如树干，茎，和树枝等。NPV主要包含碳元素，以淀粉，纤维素和木质素形式存在，NPV的光谱特征主要受这些物质支配。在短波红外内的波动比较大，与绿色植被相反，SWIR-1和SWIR-2范围内散射占主导。图5显示了绿色植被和NPV冠层光谱特征。V1SNIRSWIR1SWIR2Wavelengt
10、h(nm)图5透射绿色植被和干植被的冠层反射特性的变化(400nm2500nm).植被指数植被指数(VI)是两个或多个波长范围内的地物反射率组合运算，以增强植被某一特性或者细节。所有的植被指数要求从高精度的多光谱或者高光谱反射率数据中计算。未经过大气校正的辐射亮度或者无量纲的DN值数据不适合计算植被指数。下面是7大类27种植被指数的说明，这些植被指数都是经过严格生物条件下测试的。宽带绿度BroadbandGreenness(5中)宽带绿度指数可以简单度量绿色植被的数量和生长状况，它对植物的叶绿素含量、叶子表面冠层、冠层结构比较敏感，这些都是植被光合作用的主要物质，与光合有效辐射(fAPAR也有
11、关系。宽带绿度指数常用于植被物候发育的研究，土地利用和气候影响评估，植被生产力建模等。宽带绿度指数选择的波段范围在可见光和近红外，一般的多光谱都包含这些波段。下面的公式中规定波段的中心波长：pNIR=800nm,pRED=680nm,pBLUE=450nm。表1宽带绿度指数植被指数基本描述归一化植被指数(NormalizedDifferenceVegetationIndex)增加在近红外波段范围绿叶的散射与红色波段范围叶绿素吸收的差异。比值植被指数(SimpleRatioIndex)在近红外波段范围绿叶的散射与红色波段范围叶绿素吸收的比值。增强植被指数(EnhancedVegetationIn
12、dex)增弓5NDVI,解决土壤背景和大气气溶胶对茂密植被的影响大气阻抗植被指数(AtmosphericallyResistantVegetationIndex)增弓HNDVI,更好地解决大气散射的影响。绿波段总和指数(SumGreenIndex绿色波段范围的整体光散射对植被冠层间隙的敏感度。1)归一化植被指数(NormalizedDifferenceVegetationIndexNDVI)NDVI众所周知的一种植被指数，在LAI值很高，即植被茂密时其灵敏度会降低。其计算公式为:PnirPredPIRPREDNDVI=(式1)值的范围是-11,一般绿色植被区的范围是0.20.8。2)比值植被指
13、数(SimpleRatioIndexSRSR指数也是众所周知的一种植被指数，在LAI值很高，即植被茂密时其灵敏度会降低。其计算公式为:Rm食PrfdSR=(式2)值的范围是030+,一般绿色植被区的范围是2&3)增强植被指数(EnhancedVegetationIndexEVI)EVI通过加入蓝色波段以增强植被信号，矫正土壤背景和气溶胶散射的影响。EVI常用于LAI值高，即植被茂密区。其计算公式为：EVI=2.5)PnirPred7Pslue+上/PnirPred(式3)值的范围是-11,一般绿色植被区的范围是0.20.8。4）大气阻抗植被指数（AtmosphericallyResistant
14、VegetationIndeARVDARVI是NDVI的改进，它使用蓝色波段矫正大气散射的影响（如气溶胶），ARVI常用于大气气溶胶浓度很高的区域,如烟尘污染的热带地区或原始刀耕火种地区。其计算公式为：EVI=PnirPblue)(式4)值的范围是-1-1,一般绿色植被区的范围是0.20.8。5）绿度总和指数（SumGreenIndexSG）SG指数对稀疏SG指数是用于探测绿色植被变化最简单的植被指数。由于在可见光范围内，绿色植被对光强吸收,植被的小变化非常敏感。SG指数是500nm600nm范围内平均波谱反射率。2.2窄带绿度总和最后会被转化回反射率。值的范围是050+,一般植被区域是102
15、5。NarrowbandGreenness(7中)窄带绿度指数对叶绿素含量、叶子表面冠层、叶聚丛、冠层结构非常敏感。它使用了红色与近红外区域部分一一红边，红边是介于690nm740nm之间区域，包括吸收与散射。它比宽带绿度指数更加灵敏，特别是对于茂密植被。表2窄带绿度指数植被指数基本描述红边归一化植被指数(RedEdgeNormalizedDifferenceVegetationIndex)使用红边波段的改进型NDVIo改进红边比值植被指数(ModifiedRedEdgeSimpleRatioIndex)使用红边和蓝色波段比值。改进红边归一化植被指数(ModifiedRedEdgeNormal
16、izedDifferenceVegetationIndex)使用蓝色波段，补偿了光散射。Vogelmann红边指数1(VogelmannRedEdgeIndex1标示红色至近红外过渡的交接处，指示树冠压力。Vogelmann红边指数2(VogelmannRedEdgeIndex2标示近红外过渡形状，指示树冠胁迫性和衰老迹象。Vogelmann红边指数3(VogelmannRedEdgeIndex3标示近红外过渡形状，指示树冠胁迫性和衰老迹象。在红外过渡定位最大衍生物，对叶绿素浓度非常敏红边位置指数(RedEdgePositionIndeR感。1)红边归一化植被指数(RedEdgeNormali
17、zedDifferenceVegetationInde，NDVI705)NDVI705是NDVI的改进型，它对叶冠层的微小变化、林窗片断和衰老非常灵敏。它可用于精细农业、森林监测、植被胁迫性探测等。其计算公式为78:NDVI705=（式5）值的范围是-11,一般绿色植被区的范围是0.20.9。2)改进红边比值植被指数(ModifiedRedEdgeSimpleRatioIndemSR705)mSR705改正了叶片的镜面反射效应，可它可用于精细农业、森林监测、植被胁迫性探测等。其计算公式为6:mSR705=(式6)值的范围是030,一般绿色植被区的范围是28。3)改进红边归一化植被指数(Modi
18、fiedRedEdgeNormalizedDifferenceVegetationIndemNDVI705)mNDVI705是NDVI705的改进型，它考虑了叶片的镜面反射效应。它对叶冠层的微小变化、林窗片断和衰老非常灵敏。它可用于精细农业、森林监测、植被胁迫性探测等。其计算公式为：P7OSP7O5-2PmNDVI705=一（式7）值的范围是-11,一般绿色植被区的范围是0.20.7。4)Vogelmann红边指数1(VogelmannRedEdgeIndexVOG1)VOG1指数对叶绿素浓度、叶冠层和水分含量的综合非常敏感。它可应用于植物物候变化研究、精细农业和植被生产力建模。其计算公式为：
19、VOG1=,（式8）值的范围是020,一般绿色植被区的范围是48。5）Vogelmann红边指数2（VogelmannRedEdgeInder2VOG2）VOG2指数对叶绿素浓度、叶冠层和水分含量的综合非常敏感。它可应用于植物物候变化研究、精细农业和植被生产力建模。其计算公式为:P715+P7Z6VOG2=(式9)值的范围是020,一般绿色植被区的范围是48。6）Vogelmann红边指数3（VogelmannRedEdgeIndeZ-VOG3）VOG3指数对叶绿素浓度、叶冠层和水分含量的综合非常敏感。它可应用于植物物候变化研究、精细农业和植被生产力建模。其计算公式为：VOG3=(式10)值的
20、范围是020,一般绿色植被区的范围是487）红边位置指数（RedEdgePositionIndexREPRE可旨数对植被叶绿素浓度变化、叶绿素浓度增加使得吸收特征变宽及红边向长波段方向移动非常敏感。红边位置在690nm740nm范围内急剧倾斜波长范围，一般植被在700nm730nm。RE对旨数的结果输出是在0.69微米0.74微米光谱范围内，植被红边区域内的反射率的最大导数的波长。常用于农作物监测和估产，生态系统干扰探测，光合作用模型，和由气候或其他因素产生的冠层胁迫性。利用率LightUseEfficiency（3种）光利用率指数是用来度量植被在光合作用中对入射光的利用效率。光的利用效率直接
21、与碳吸收效率、植被生长速度和光合有效辐射（fAPAR）有很大的关系。表3光利用率指数植被指数基本描述光化学植被指数(PhotochemicalReflectanceIndex)对估算叶片类胡萝卜素（尤其黄色色素）、叶片胁迫性和碳吸收效率非常有用结构不敏感色素指数(StructureInsensitivePigmentIndex)作为总体冠层结构和植被容量变化的叶片色素含量指标红绿比值指数(RedGreenRatioIndex)对花青素与叶绿素的比率非常敏感1）光化学植被指数（PhotochemicalReflectanceIndexPRI）PRI对活植物的类胡萝卜素（尤其黄色色素）变化非常敏感
22、，类胡萝卜素可标识光合作用光的利用率，或者碳吸收效率。可用于研究植被生产力和胁迫性，常绿灌木植被的健康，森林以及农作物的衰老。其计算公式为:Pst二P第70一目53+2370PRI=MI（式11）值的范围是-11,一般绿色植被区的范围是-0.20.2。2）结构不敏感色素指数（StructureInsensitivePigmentIndexSIPI）SIPI用来最大限度地提高类胡萝卜素（例如a-胡萝卜素和0-胡萝卜素）与叶绿素比率在冠层结构（如叶面积指数）减少时的敏感度，SIPI的增加标识冠层胁迫性的增加。可用于植被健康监测、植物生理胁迫性检测和作物生产和产量分析。其计算公式为：peoo+pes
24、，冠层氮指数使用短波红外测量植被冠层中氮的相对含量。归一化氮指数（NormalizedDifferenceNitrogenIndexNDNI）NDNI是用于估算植被冠层中氮的相对含量。在1510nm的反射率主要取决于叶片氮的含量，以及冠层总体叶生物量。结合叶片氮含量和冠层叶生物量在1520nm范围内预测叶片氮的含量，在1680nm波长范围作为参考反射率，冠层叶生物量这个波长范围具有与1520nm波长范围类似的反射特性，而且1680nm波长范围内没有氮吸收影响。NDNI在植被还是绿色以及覆盖浓密时候，对氮含量的变化非常敏感，它用于精细农业、生态系统分析和森林管理。其计算公式为：NDNI=（式13
25、）值的范围是01,一般绿色植被区的范围是0.020.1。干旱或碳衰减DryorSenescentCarbon（3中）干旱或碳衰减指数是用来估算纤维素和木质素干燥状态的碳含量。干碳分子大量存在于木质材料和衰老、死亡、或休眠的植被，可以使用这些指数可以做植被着火性分析和检测森林的枯枝落叶层。干旱或碳衰减指数是基于纤维素和木质素在短波红外波段吸收特性而计算。表4干旱或碳衰减指数植被指数基本描述归一化木质素指数(NormalizedDifferenceLigninIndex)检测叶片木质素在1754nm特征与1680nm冠层结构区域相互关系产生的木质素增加纤维素吸U指数(CelluloseAbsorp
26、tionIndex)检测纤t素在2000nm波长范围引起的吸收特征植被衰减指数(PlantSenescenceReflectanceIndex使用类胡萝卜素与叶绿素的比值来检测植物开始和衰老程度1）归一化木质素指数（NormalizedDifferenceLigninIndexNDL|NDLI是用来估算植被冠层木质素的相对含量，应用生态系统分析和检测森林的枯枝落叶层。其计算公式为：I胆L/p1754NDLI=（式14）值的范围是01,一般绿色植被区的范围是0.0050.052）纤维素吸U指数（CelluloseAbsorptionIndexCAI）CAI可以指示地表含有干燥植被，纤维素在200
27、0nm2200nm范围内吸收特征非常敏感。应用于农作物残留监测，植物冠层衰老，生态系统中的着火条件和放牧管理。其计算公式为：MnilS502000+02200）02100NDLI=（式15）值的范围是-34+,一般绿色植被区的范围是-24。3）植被衰减指数（PlantSenescenceReflectanceIndPSRIPSR用来最大限度地提高类胡萝卜素（例如a-胡萝卜素和B-胡萝卜素）与叶绿素比率的灵敏度，PSRI的增加预示冠层胁迫性的增加、植被衰老的开始和植物果实的成熟。可用于植被健康监测、植物生理胁迫性检测和作物生产和产量分析。其计算公式为：P75GPSRI=(式16)值的范围是-1-
29、波长花青素反射指数1(AnthocyaninReflectanceIndex1)叶片花青素在绿色波段相对红色波段的吸收特征的变化花青素反射指数2(AnthocyaninReflectanceIndex2)与ARI1类似1)类胡萝卜素反射指数1(CarotenoidReflectanceIndex-CRI1)CRI1对叶片中的类胡萝卜素非常敏感，高的CRI1值意味类胡萝卜素含量相比叶绿素含量多。其计算公式为:CRI1=Psio(式17)值的范围是015+,一般绿色植被区的范围是112。2)类胡萝卜素反射指数2(CarotenoidReflectanceIndex-CRI2)CRI2是CRI1的改
30、进型，在类胡萝卜素浓度高时更加有效，高的CRI2值意味类胡萝卜素含量相比叶绿素含量多。其计算公式为：(式18)值的范围是015+,一般绿色植被区的范围是1113）花青素反射指数1（AnthocyaninReflectanceIndexARI1）ARI1对叶片中的花青素非常敏感，ARI1值越大表明植被冠层增长或者死亡。其计算公式为:ARI1=（式19）值的范围是00.2+,一般绿色植被区的范围是0.0010.1。4）花青素反射指数2（AnthocyaninReflectanceIndexARI2）ARI2对叶片中的花青素非常敏感，ARI2值越大表明植被冠层增长或者死亡。ARI2是ARI1的改进，
31、当花青素浓度高时更加有效。其计算公式为：P800）P7（J0ZARI2=（式20）值的范围是00.2+,一般绿色植被区的范围是0.0010.1。-层水分含量CanopyWaterContent（4种）冠层水分含量指数用于度量植被冠层中水分含量。水分含量是一个重要的植物指标，较高的水含量表明健康植被、生长快及不易着火。冠层水分含量指数基于水在近红外和短波红外范围内的吸收特征，以及光在近红外范围的穿透性，综合起来度量总的水柱含量。表6冠层水分含量指数植被指数基本描述水波段指数(WaterBandIndex)度量在900nm处吸收强度随着冠层水量增加而增加归一化水指数(NormalizedDiffe
32、renceWaterIndex)度量植被水分含量在857nm吸收相对1241nm吸收的增率水分胁迫指数(MoistureStressIndex)检测在1599nm植被水分变化归一化红夕卜指数(NormalizedDifferenceInfraredIndex)度量在1649nm处吸收强度随着冠层水量增加而增加1）水波段指数（WaterBandIndex-WBI）WBI对冠层水分状态的变化非常敏感，随着植被冠层水分的增加，970nm附近吸收强度相比900nm处有所增强。应用包括冠层胁迫性分析，生产力预测与建模，着火威胁条件分析，农作物管理，以及生态系统生理机能研究。其计算公式为：WBI=（式21
33、）一般绿色植被区的范围是0.81.2.2)归一化水指数(NormalizedDifferenceWaterInde-NDWI)NDWI对冠层水分含量的变化非常敏感，因为在857nm和1241nm具有相似的反射率，但是又不同于液态水的吸收特性。应用于冠层胁迫性分析，在浓密叶型植被的叶面积指数的研究，植被生产力模型，着火性研究。其计算公式为：NDWI=（式22）值的范围是-11,一般绿色植被区的范围是-0.10.43)水分胁迫指数(MoistureStressInde-MSI)MSI对叶片水分含量的增加非常敏感。当叶片水分含量的增加，在1599nm处的吸收强度也增加，而在819nm处的吸收强度没有
34、影响。应用于冠层胁迫性分析，生产力预测与建模，着火威胁条件分析，以及生态系统生理机能研究。与其他水指数相反，MSI值越大，水分胁迫性越严重和水分含量越少。其计算公式为：MS（式23）值的范围是03+,一般绿色植被区的范围是0.42。4)归一化红外指数(NormalizedDifferenceInfraredIndeNDII)NDII对农作物冠层的水分含量变化非常敏感，NDII的值越大表示水分含量越多。应用于农作物管理，森林冠层监测，植被胁迫性探测。其计算公式为：Pg工9一P1存P9191nl号NDII=（式24）值的范围是-11,一般绿色植被区的范围是0.0206.环境小卫星高光谱成像仪的植被
35、指数计算环境与灾害监测预报小卫星星座A、B星(HJ-1A/1B星)于2008年9月6日上午11点25分成功发射，HJ-1-A星搭载了CCD相机和超光谱成像仪(HSI),HJ-1-B星搭载了CCD相机和红外相机(IRS。在HJ-1-A卫星装载有一台超光谱成像仪，完成又日!刈宽为50公里、地面像元分辨率为100米、110128个光谱谱段的推扫成像，回访周期约为96ho该数据可免费下载：oENV提供植被指数计算器，包括了上述27种植被指数。它可以根据输入图像波段情况，自动从27种中列出能计算的植被指数。并提供了还提供了生物物理学交叉检验功能，能够提高植被指数的计算精度。我们使用的HJ-1-A-HI激据包含115个波段，波长覆盖459.00956.00nm,图彳t以HDF5格式储存，HDF5图像格式中除了图像文件外，还包括了中心波长、定标文件、成像参数等信息。下面我们实验在ENVI中能计算几种植被指数。直接利用ENVI_HJ1A1B_Tool环境卫星数据读