基于MODIS的干旱区植被覆盖度反演及植被指数优选
第2期,总第89期2011年6月15日
国土资源遥感
No.2,2011Jun.,2011
REMOTE SENSING FOR LAND &RESOURCES
基于MODIS 的干旱区植被覆盖度反演
及植被指数优选
112
郭玉川,何英,李霞
(1.新疆农业大学水利与土木工程学院,乌鲁木齐830052;2.新疆农业大学草业与环境科学学院,乌鲁木齐830052)
摘要:为进一步掌握塔里木河下游输水后的生态恢复程度,以输水河畔的植被覆盖度为研究对象,采用MODIS 数据构建研究区多种植被指数;结合现场实测植被覆盖度,给出离散坐标下实测植被覆盖度与各种植被指数间的二维散MSAVI 、SAVI 和EVI 点图,据此得出二者相关方程,反演区域植被覆盖度。结果表明,基于MODIS 数据构建的NDVI 、等植被指数均与植被覆盖度有较好的相关关系,采用这些植被指数反演植被覆盖度的精度由高到低依次为NDVI 、EVI 、MSAVI 及SAVI 。建议在干旱荒漠区,使用MODIS 数据进行植被变化监测,最优选择的植被指数为NDVI 。关键词:塔里木河下游;植被覆盖度;植被指数;MODIS 中图法分类号:TP 751.1
文献标识码:A
文章编号:1001-070X (2011)02-0115-04
0引言
自2000年5月实施塔里木河下游生态输水工
1
1.1
数据获取
地面植被覆盖度获取
程以来,输水河段的生态恢复进展情况成为关注热
点,植被覆盖度作为刻画生态恢复的重要因子之一也备受关注。目前对植被覆盖度的研究主要有地表实测和遥感测量两种方法
[1,2]
地面植被覆盖度调查样地的布设主要考虑不同
离河距离、不同地下水埋深以及有胡杨、怪柳等植被生长的地段。在2006年5 6月,在英苏(40ʎ25' 52.3ᵡN ,87ʎ56'18. 7ᵡE )、88ʎ 喀尔达依(40ʎ22'00. 9ᵡN ,09'52.2ᵡE )和阿拉干(40ʎ08'29. 6ᵡN ,88ʎ20'51. 6ᵡE )3
B 和C ,个地点分别布设3个剖面(A 、每个剖面有3条断面),剖面为南北或东西走向,剖面间距2000m ,
每个剖面在河道左、右两岸各设8个样地,样地面积9条断面共设144个样地(图1)
。均为250m ˑ 250m ,
。地表实测准确、客
[1]
观,但耗时较大,只能在小尺度范围内提供植被结构不宜大范围推广和分布状况的变化信息,
;遥感
测量植被覆盖度则具有宏观、时效性强、数据获取容
易等诸多优点,是当前植被覆盖度研究的主流方向。对于干旱区植被覆盖度,目前的研究多集中在通过分析某一植被指数的变化趋势来反映植被生长情况
[3-5]
,或采用混合像元解混法反演区域地表植被
,这些方法较之传统的地面调查可以更
覆盖度
[6-9]
更快地反映区域植被覆盖度的变化趋势,但大多好、
而针对同一区域针对某一植被指数进行分析研究,
不同植被指数与植被覆盖度相关程度的分析研究的
报道尚不多见。
本文结合塔里木河下游输水河段实测样地的植
对该区域基于MODIS 的几种常见植被覆盖度数据,
被指数与实测植被覆盖度的相关性进行了探讨,采用经验统计模型建立了塔里木河下游输水河畔4种植
并对遥感反演植被被指数与植被覆盖度的相关方程,
覆盖度常用的植被指数进行了优选,以期为研究塔里
木河下游输水河段的生态恢复提供科学依据。
图1
Fig.1
地面植被指数调查样地布设示意图
Distribution of the vegetation coverage survey plots
收稿日期:2010-08-23;修订日期:2010-09-27
基金项目:国家自然科学基金项目(编号:40961027)和水利部公益项目(编号:200801050)共同资助。
·116·
国土资源遥感2011年
植被覆盖度通过实测与目估法来确定。对样地柽柳等大型乔灌木的植被覆盖度通过逐棵内胡杨、
实测获得,对中小型灌草类通过实测得到多个1m ˑ 1m 样地内的植被投影面积,再由目估法推求出整个样地的植被覆盖度。1.2取
植被指数获取
结合研究区域的地理环境特征,植被指数的选
[10,11]
MOD09GHK 数据)进行植被指数的计算。对于采用首先将空间分辨率为蓝光波段修正NDVI 的EVI ,
500m 的蓝光波段重采样为250m 的分辨率,以获取250m 分辨率的植被指数;在获取各植被指数对植被指数结果进行投影变换和几何精纠正,使后,
其具有相同的投影系统,并使各像元的植被指数与植被覆盖度调查样地的坐标一一对应;最后根据各样地坐标在配准后的植被指数空间分布图上查出与样地一一对应的植被指数,再由SPSS 统计分析软件分别对各植被指数/植被覆盖度散点图进行拟合曲线优选,从而求出各植被指数与植被覆盖度的相关方程。2.2
植被覆盖度的遥感反演
根据已获得的植被指数与植被覆盖度相关方程,采用ERDAS 的Model Maker 模块对已求出的各植被指数进行代入计算,反演不同植被指数所对应的地表植被覆盖度。
主要考虑研究区土壤背景的影响;对于大气
气溶胶的影响,由于所需气象参数难以获取而不做专门考虑,只是尽可能选择天气晴朗、无风的时段进行植被指数的反演。选取用于对比的植被指数包括考虑应用最为广泛的归一化差值植被指数(NDVI )、土壤背景干扰的土壤调节植被指数(SAVI )和修正以及综合考虑了土的土壤调节植被指数(MSAVI ),
壤和大气共同影响的增强型植被指数(EVI )。各植被指数的计算公式为
NDVI =SAVI =
MSAVI =[2ρN +1-
ρN -ρR
ρN +ρR
(1)(2)
ρN -ρR
(1+L )
ρN +ρR +L
3
3.1
结果分析
植被覆盖度与植被指数相关性分析
通过查找144个样地对应坐标下的植被指数,
2+1)-8(ρN -ρR )]/2ρN (3)
EVI =
ρN -ρR
1+l )
ρN +C 1ρR -C 2ρB +l
(4)
采用SPSS 统计分析软件对样地实测植被覆盖度(Fv )与对应坐标下的植被指数(VI )的相关性进行分析(表1)。结果表明,实测样地的植被覆盖度与可以采用植同坐标下的植被指数相关性非常显著,被指数反演对应的植被覆盖度。
表1Tab.1
实测植被覆盖度与植被指数相关分析①Analysis of correlation for F vs measured
in site and different VI s
覆盖度-指数F v -NDVI F v -SAVI F v -MSAVI F v -EVI
样本数[1**********]4
P 值<0.001<0.001<0.001<0.001
相关系数0.9180.8380.8560.866
相关性非常显著非常显著非常显著非常显著
式中,ρN 、ρR 和ρB 分别为近红外波段、红光波其值段和蓝光波段的反射率;L 为土壤亮度指数,Huete 取决于植被的密度,
[12]
建议L 的最佳值为
0.5;l 为背景调节参数;C 1、C 2为大气修正参数。C 1、C 2和l 分别取值为6.0、7.5对于MODIS 数据,和1.0。
考虑到植被覆盖度遥感反演对时间及空间分辨率的要求以及实际应用的需要,选取与地面实测数根据上述公式即可得据同时段获取的MODIS 数据,
到各植被指数在研究区的空间分布情况。
①F v 为植被覆盖度。
2
2.1
研究方法
植被指数与植被覆盖度相关分析
地表植被覆盖度调查时间为2006年6月22
以植被指数为横坐标、植被覆盖度为纵坐标,绘出二维散点图;采用SPSS 软件分析各植被指数与植被覆盖度的最佳关系拟合曲线(图2)。
结果表明,各种植被指数与植被覆盖度均呈正
2相关关系,且以三次曲线拟合的R 最大,二次曲线2
但三次曲线在植被覆盖度大于拟合的R 次之,
日 8月18日,样地大小为250m ˑ 250m 。为与调查样地的植被覆盖度进行对比分析,选取2006年7月23日获取的MODIS 数据(MOD09GQK 及
15%时拟合效果明显不如二次曲线的拟合效果。
(a )NDVI -F v (b )SAVI -F
v
(c )MSAVI -F v (d )EVI -F v
图2Fig.2
各植被指数与植被覆盖度相关关系
Correlation between different VI s and F vs
2
覆盖度的样本判断系数R 最大,为0.85;采用SAVI 反演植被覆盖度的R 2最小,为0.703。说明在
无论采用何种植被指数反演塔里木河下游研究区,
综合考虑,可选取二次曲线作为该地区植被指
数与植被覆盖度的拟合曲线。
NDVI -F v 拟合曲线为
y =-207.871x 2+153.826x -4.04,R 2=0.850(5)
SAVI -F v 拟合曲线为
y =-43.44x 2+94.410x -5.229,R 2=0.703
MSAVI -F v 拟合曲线为
y =-90.835x 2+108.669x -7.644,R 2=0.738(7)
EVI -F v 拟合曲线为
y =-355.754x 2+215.261x -3.519,R 2=0.750
表2Tab.2
植被覆盖度
F v NDVI
F v SAVI F v MSAVI F v EVI
均值5.0574.7974.8854.921
最小值1.5900.0930.0671.182
最大值19.36019.51019.32020.067
植被覆盖度,均可在一定程度上反映地表植被的覆盖程度,但采用NDVI 反演的植被覆盖度最接近实际地表植被覆盖情况。
3.2区域植被覆盖度的遥感反演
采用ERDAS IMAGINE 的Model make 工具,根据上述相关方程分别反演该区的地表植被覆盖度,结果如表2及图3所示。
(6)
(8)
由式(5) (8)可以看出,采用NDVI 反演植被
不同植被指数反演的植被覆盖度The F vs inversed by different VI s
0%≤F v <5%
52.559.753.753.5
5%≤F v <10%
31.430.231.734.7
10%≤F v <15%
15.09.313.411.0
F v ≥15%1.00.81.20.7
合计
[1**********]0
(a )F v NDVI (b )F v SAVI (c )F v MSAVI (d )F v EVI
图3Fig.3
SAVI 、MSAVI 和EVI 反演的植被覆盖度由NDVI 、
The Fvs inversed by NDVI ,SAVI ,MSAVI and EVI
由表2及图3可以看出,由不同植被指数反演的研究区植被覆盖度的均值及最大值比较接近,最小值略有差异,且空间分布特征相似;植被覆盖度在10% 20%之间的灌草地面积约占研究区面积的10% 16%,主要分布在输水河道两侧600m 范围内及河道附近的低洼地区;植被覆盖度小于5%的面积占研究区总面积的84%以上,分布在远离输水河道的荒漠过渡带和沙漠地区。
[2]邢著荣,冯幼贵,杨贵军,等.基于遥感的植被覆盖度估算方法
.遥感技术与应用,2009,24(6):849-854.述评[J ]
[3]邓铭江.塔里木河下游生态输水及植被恢复遥感监测评价
[J ].冰川冻土,2007,29(3):380-386.
[4]闫正龙,黄强,牛宝茹,等.应急输水工程对塔里木河下游地
.应用生态学报,2008,19(3):65-区植被覆盖度的影响[J ]67.
[5]徐水师,谭克龙,闫正龙.塔里木河下游植被覆盖动态变化遥
J ].西北大学学报,2010,40(2):323-326.感分析[
[6]刘静,银山,张国盛.毛乌素沙地17年间植被覆盖度变化的
2009,23(7):162-167.遥感监测[J ].干旱区资源与环境,
[7]刘静,银山,张国盛.基于高光谱混合像元分解的干旱地区
.应用生态学报,2010,21(7):152-稀疏植被覆盖度估测[J ]158.
[8]戴俣俣,丁贤荣,王文种.基于MODIS 影像的植被覆盖度提取
.遥感应用,2009(2):67-70.研究[J ]
[9]吴云,曾源,吴炳方.基于MODIS 数据的三北防护林工程区
2009,28(9):1712-植被覆盖度提取与分析[J ].生态学杂志,1718.
[10]罗亚,徐建华,岳文泽.基于遥感影像的植被指数研究方法述
J ].生态科学,2005,24(1):75-79.评[
[11]杨嘉,郭铌,贾建华.西北地区MODIS /NDVI与MODIS /EVI
2007,25(1):38-43.对比分析[J ].生态学杂志,
[12]Huete A R.A Soil Adjusted Vegetation Index (SAVI )[J ].Remote
Sens Environ ,1988,25(3):295-309.
4结论
(1)基于MODIS 的NDVI 、MSAVI 、SAVI 和EVI
等4种植被指数均可以较好地反映地表植被的覆盖
情况,并且与植被覆盖度呈显著正相关关系。(2)对于MODIS 数据,4种植被指数与地表植
EVI 、被覆盖度的相关程度由高到低依次为NDVI 、MSAVI 及SAVI 。建议在塔里木河下游,采用基于MODIS 数据的NDVI 进行植被覆盖变化监测。参考文献:
[1]程红芳,章文波,陈锋.植被覆盖度遥感估算方法研究进展
[J ].国土资源遥感,2008(1):13-17.
Remote Sensing Inversion of Vegetation Coverage and Optimization of
Vegetation Index Based on MODIS Data in Arid Area
GUO Yu -chuan 1,HE Ying 1,LI Xia 2
(1.College of Water Conservancy and Civil Engineering ,Xinjiang Agricultural University ,Urumqi 830052,China ;2.College of Prataculture and Environment Science ,Xinjiang Agricultural University ,Urumqi 830052,China )
Abstract :To further understand the ecologic restoration level after water supplying ,the authors studied vegetation coverage (Fv )in the water supply area of the lower reaches of the Tarim River.Correlation equation and two -dimensional scatter plot between Fv and vegetation index (VI )could be obtained based on MODIS data by constructing a variety of VIs in combination with measured Fvs ,and the value of Fv could be figured out based on the relevant equation.The results show that there exist good correlations between the VIs such as NDVI ,MSAVI ,SAVI and EVI and Fvs on the basis of MODIS data.With these VIs ,the inversion accuracy of the Fvs in the order from higher to lower is NDVI ,EVI ,MSAVI and SAVI.The authors suggest using MODIS data to monitor vegetation change in the arid desert area ,with the optimal choice of VI being NDVI.
Key words :Lower reaches of the Tarim River ;Vegetation coverage (Fv );Vegetation Index (VI );MODIS
第一作者简介:郭玉川(1982-),男,博士研究生,主要从事地下水资源利用与保护研究。
通讯作者:李霞(1957-),女,教授,硕士生导师,主要从事生态环境保护与利用研究。E -mail :xjlx782@126.com
(责任编辑:刘心季)