徐爽;沈润平;杨晓月
【摘 要】ASD Field Spec Pro FRTM spectroradiometer was used to measure the spectral response of the vegetable and grass at different vegetation coverage levels. The data were applied to calculate six vegetation indices, i. e. , NDVI (normalized difference vegetation index) , DVI (difference vegetation index) , RVI (ratio vegetation index) , MVI (modified vegetation index) , MSAVI (modified soil adjusted vegetation index) and GEMI (global environment monitoring index). Then the best combination of spectral bands was analyzed. Furthermore, the performance of different vegetation indices was investigated when they were used to estimate the vegetation coverage by using the dimidiate pixel model. The results show that, for the green vegetable, the best combinations of bands in the spectral region from 620 to 740 nm and from 780 to 900 nm have the best correlation with the vegetation index, whereas for the grass, the best combinations of bands are from 620 to 750 nm and from 760 to 900 nm, with the correlation coefficients of the two cases being all l
arger than 0. 8. The bands of Landsat7 and HJ -1A CCD1 simulated according to the spectral response function were employed to calculate the six vegetation indices. The average overall accuracy for estimating the vegetation fraction by DVI and MSAVI is 83.7% and 79.5% respectively, indicating that they are superior to the other four vegetation indices as the input of vegetation index for the dimidiate pixel model.%选用蔬菜地和草地2种植被类型,利用ASD光谱仪实测二者在不同覆盖度下的光谱响应,分析了归一化植被指数(NDVI)、差值植被指数(DVI)、比值植被指数(RVI)、修正植被指数(MVI)、修改型土壤调节植被指数(MSAVI)以及全球环境监测植被指数(GEMI)等6种植被指数所用的最佳波段及其组合,进而研究了利用像元二分模型估算植被覆盖度时的不同植被指数的表现.结果表明,与蔬菜地植被指数相关系数较高的波段组合为620 ～ 740 nm谱段和780 ～ 900 nm谱段内波段的组合,与草地植被指数相关系数较高的波段组合为620 ～750 nm谱段和760 ～900 nm谱段内波段的组合,相关系数均达0.8以上；在高光谱数据构建的植被指数和模拟卫星数据构建的植被指数中,用DVI和MSAVI估算植被覆盖度,平均总体精度分别达到83.7％和79.5％,与其他4种植被指数相比,这2种指数更适合于利用像元二分模型进行植被覆盖度的估算.
【期刊名称】《国土资源遥感》
【年(卷),期】2012(000)004
【总页数】6页(P95-100)
【关键词】植被覆盖度;植被指数;像元二分模型
【作 者】徐爽;沈润平;杨晓月
【作者单位】南京信息工程大学气象灾害省部共建教育部重点实验室,南京210044;南京信息工程大学遥感学院,南京 210044;南京信息工程大学气象灾害省部共建教育部重点实验室,南京210044;南京信息工程大学气象灾害省部共建教育部重点实验室,南京210044;南京信息工程大学遥感学院,南京 210044
【正文语种】中 文
【中图分类】TP79
植被作为陆地生态系统的主要成分之一，对气候变化具有重要的调节作用［1］。植被覆盖度(vegetation coverage)是指植被在地面的垂直投影面积占统计区总面积的百分比［2］，是衡量地表植被状况的一个重要指标，也是影响土壤侵蚀与水土流失的主要因子，对于区域环境变化和监测研究具有重要意义［3］。利用遥感数据估算植被覆盖度是测定区域植被覆盖度的主要手段之一，目前应用的遥感测量方法较多，其中应用最广泛的是回归模型法和像元分解模型法。Wittich等［4］和刁兆岩等［5］都曾建立了归一化差值植被指数(normalized difference vegetation index，NDVI)和植被覆盖度的回归模型，并对所选研究区植被覆盖度进行了估算;但回归模型法一般只适用于特定区域和特定植被类型，并且需要大量的实测数据建立模型，有一定的局限性，不利于普遍推广应用。像元分解法则具有理论基础较好，不依赖于实测数据，不受地域限制，可削弱大气、土壤背景和植被类型等影响，较易于推广等优点。其中，像元二分模型(dimidiate pixel model)法得到广泛的认可，具体应用时，需要首先构建植被指数。
目前，应用NDVI来估算植被覆盖度的研究较多［6－8］，有关其估算精度问题，Jiang 等［9］对 NDVI代入模型所存在的误差问题进行了公式推导;Small［10］研究发现，在中等植被覆盖度时利用TM数据计算的NDVI随着叶面积指数(leaf area index，LAI)的增加逐渐趋
于饱和，并且会高估植被密度较高区域周围零散分布的植被区的覆盖度，这种现象主要是由于阴影对NDVI的影响造成的。由于不同植被指数的构建有着不同的理论依据和应用价值，利用像元二分模型法开展不同植被指数应用效果的比较尚缺乏系统、深入的研究。
本文基于地表高光谱实测数据，构建了归一化植被指数(NDVI)、差值植被指数(difference vegetation index，DVI)、比值植被指数(ratio vegetation index，RVI)、修正植被指数(modified vegetation index，MVI)、修改型土壤调节植被指数(modified soiladjusted vegetation index，MSAVI)及全球环境监测植被指数(global environment monitoring index，GEMI)等6种植被指数，采用像元二分模型对植被覆盖度进行估算，对比分析不同植被指数估算植被覆盖度的效果，为植被覆盖度的遥感估算提供理论依据。
选用位于 N 32.203°，E 118.705°处，覆盖度达100%的青菜地为代表的蔬菜地和草地。光谱测定日期为2011－10－14，天气状况良好，晴朗无云。采用ASD Field Spec Pro FRTM野外光谱辐射仪进行采集，光谱采集范围350～2500 nm，采样间隔1.4 nm(350～1000 nm)和2 nm(1000～2500 nm)，重采样间隔1 nm。测定光谱时，传感器探头垂直向下，视场角25°，距离冠层顶1.6 m;每次采集10条光谱，最后取光谱平均值。所测光谱数据经处理后，选择适当谱带构建植被指数，以便利用像元二分模型进行植被覆盖度估算。
植被覆盖度通过数码相机拍摄获得。为了保证植被覆盖度测量的准确性，并保证与光谱测量范围相对应，在拍摄前放置一铁圈于草地周围，铁圈大小尽量与光谱仪的视场角范围一致;拍摄时相机垂直向下，保证整个铁圈(样圆)内的草地在拍摄范围内(测量光谱时将铁圈移走)。为研究不同植被覆盖度下光谱曲线的变化，在利用光谱仪完成样区内植被覆盖度100%下的光谱测量后，分7次在该样圆内进行人工均匀去除一定量植被，调节植被覆盖度，直至样圆内所有植被全部去除;每去除一次植被后，进行一次光谱测量和植被覆盖度拍照。对蔬菜地也做同样处理，对不同覆盖度下的蔬菜地进行了光谱采集和覆盖度测量。
在Photoshop和ENVI软件下完成植被覆盖度的求算。首先在Photoshop软件下进行样圆裁切，然后在ENVI软件下进行监督分类，求算植被所占百分比，得到样地的植被覆盖度。草地植被覆盖度分别为 100%，78%，62%，44%，37%，25%，18% 和0%;蔬菜地的植被覆盖度分别为100%，85%，74% ，62%，58%，52%，42%，39%，33%，25%，10%和0%。利用ViewSpecPro软件剔除2个样地光谱数据中的异常数据，并取光谱数据的平均值得到不同覆盖度下的植被光谱曲线。
首先要确定有效谱段范围和波段组合，波段选择方法主要有逐步回归选择法、简单波段自
相关选择法和主成分系数权重选择法［11－12］等。比较所有波段组合下求算的植被指数与植被覆盖度的相关系数大小，来确定构建植被指数的较优波段。由于运算量较大，所以通过VBA编程实现运算。另外，根据光谱响应函数模拟卫星传感器波段的反射率，利用高光谱数据模拟美国陆地卫星Landsat7 TM数据和我国环境减灾卫星HJ－1A CCD数据，分别计算NDVI，DVI，RVI，MVI，MSAVI和 GEMI等 6 种植被指数，用以分析代入像元二分模型的适宜性。6种植被指数的计算公式见表1。
像元二分模型［19］是一种简单实用的遥感估算模型，它假设一个像元的地表由有植被覆盖地表与无植被覆盖地表共同组成，而传感器观测到的光谱信息也由这2个组分的因子线性加权合成，各因子的权重是各自的面积在像元中所占比例，如植被的权重可以看作是植被覆盖度。
根据像元二分模型的原理，通过遥感传感器所获信息S可以表达为由绿植被组分所贡献的信息Sv和由土壤组分所贡献的信息Ss组成，即
设像元中有植被覆盖的面积比例为fc(即该像元的植被覆盖度)，则无植被覆盖地表(土壤覆盖)的面积比例为1－fc;设全部由植被所覆盖的纯像元所得的遥感信息为Sveg，则混合像元
中的植被组分所贡献的信息Sv可以表示为Sveg与fc的乘积，即
同理，设全部由土壤覆盖的纯像元所得的遥感信息为Ssoil，则混合像元中的土壤成分所贡献的信息Ss可以表示为Ssoil与1－fc的乘积，即
换算率将式(2)(3)代入式(1)，可得
对式(4)进行变换，可得
式中Ssoil与Sveg是模型的2个参数。本文中Ssoil即为植被覆盖度为0%时的植被指数信息;Sveg即为植被覆盖度为100%时的植被指数信息;S则为所求覆盖度下的植被指数信息。将 NDVI代入式(5)［20］，可得
将表示其他5种植被指数的变量(DVI，RVI，MVI，MSAVI和GEMI)也分别代入式(5)进行植被覆盖度估算。
图1为不同植被覆盖度下蔬菜地和草地的光谱反射率曲线图。我眼中的网络世界
从图1可以看出，随着植被覆盖度的增大，可见光谱段内蓝光波段反射率有减小趋势，绿光
胜利电视台
波段反射率有增大趋势(但不是很明显)，红光波段反射率减小较为明显;近红外谱段反射率则逐渐增大，并且增加幅度较大。裸土光谱反射率在红光波段内明显高于有植被覆盖时的反射率，在近红外谱段则明显低于有植被覆盖时的反射率。
寻圣诞老人3.2.1 植被指数最佳波段选择
因植被光谱反射率随植被覆盖度变化的规律在1000 nm后表现不明显，因此本文选取1000 nm前的数据进行了植被指数构建。由于在可见光和近红外谱段植被光谱随植被覆盖度变化的规律性明显，因此首先对这2个谱段内所有波段两两组合，求算其植被指数，然后计算植被指数与植被覆盖度的相关系数，依照相关系数的大小来确定植被覆盖度估算较优的植被指数波段组合。确定的结果是:蔬菜地植被指数相关系数较高的波段组合为620～740 nm谱段和780～900 nm谱段内波段的组合，草地植被指数相关系数较高的波段组合为620～750 nm谱段和760～900 nm谱段内波段的组合，植被指数的相关系数均达0.8以上。本文选取相关系数最高的最佳波段组合构建 NDVI，DVI，RVI，MVI，MSAVI和GEMI(表2)，再分别代入像元二分模型估算对应的植被覆盖度。
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议