计算ndvi值需要的数据_大气校正前后山地流域NDVI变化特征对比

摘要

大气校正是遥感影像处理的重要环节,探究大气校正前后归一化植被指数(normalized difference vegetation index,NDVI)特征及其地形梯度特征具有一定理论意义。以贵州省境内包溪河流域为例,基于Landsat8遥感影像和FLAASH大气校正,探究NDVI分等级和分地形计算误差。得出以下结论: ①未进行大气校正所得NDVI在空间上总体被低估; ②未经大气校正所得NDVI不能反映各阶段趋势和比例关系,存在严重偏差,NDVI>0.6时绝对误差超过20%; ③大气校正与否影响各海拔NDVI变化趋势和数值高低,海拔在1 000 m以下时绝对误差上升,其后波动下降; ④大气校正与否对各坡度NDVI趋势和数值均有影响,随坡度上升,绝对误差先上升后下降,坡度在[45°,50°)时绝对误差最大; ⑤大气校正与否对各坡向NDVI趋势均有影响,西坡绝对误差最大,东坡最小。

关键词： NDVI ; 地形梯度特征 ; 大气校正 ; 包溪河流域 ; 山地流域

本文引用格式

程东亚, 李旭东. 大气校正前后山地流域NDVI变化特征对比. 国土资源遥感杂志[J], 2020, 32(1): 90-97 doi:10.6046/gtzyyg.2020.01.13

CHENG Dongya, LI Xudong. Comparison of change characteristics of NDVI in mountain basin before and after atmospheric correction. REMOTE SENSING FOR LAND & RESOURCES[J], 2020, 32(1): 90-97 doi:10.6046/gtzyyg.2020.01.13

0 引言

大气校正是遥感影像处理、反演的重要环节[1,2]。归一化植被指数(normalized difference vegetation index,NDVI),是遥感反演的常用指数,在生态治理、资源保护和城市环境等领域取得广泛应用。周雪欣等[3]结合NDVI对北盘江生态进行了评价; 吴端耀等[4]利用NDVI探究了贵州省植被变化; 江颖慧等[5]开展了城市环境与NDVI的相关性研究。当前,NDVI在人类生产生活和科学研究中发挥着重要作用,且其研究不断深入。此背景下,NDVI方法机制等相关研究也不断加深。王涛等[6]分析了不同NDVI提取方法对秦岭山地植被变化的影响; 廖顺宝等[7]利用NDVI图谱库探究了土地覆盖分类精度; 刘建文等[8]采用三次样条函数法、双逻辑斯蒂函数法和奇异谱分析法开展了NDVI时间序列重建。还有学者利用NDVI提取了森林冰冻受灾范围精度[9]、利用改进的NDVI密度分割方法提取了冬小麦信息[10]等。综合来看,目前对NDVI应用、方法改进和大气校正等研究已有一定基础,但涉及大气校正前后的地形梯度误差特征的研究相对较少。山地流域植被变化与其他区域不同,大气校正对其NDVI的影响尚不明确,不同地形梯度和空间效应上的影响及其程度等问题需进一步研究。由于尺度太大不易看出NDVI的空间变化特征,尺度太小不具有山地流域代表性。因此,以尺度合适、地形复杂、人口较少的包溪河流域为例,探究大气校正前后在空间尺度、地形尺度上NDVI的误差程度。

1 研究区概况及数据源

1.1 研究区概况

包溪河流域为乌江与沅江的分界,位于贵州省东部,主要分布于铜仁市和黔南州,面积约196.87 km2。该流域属亚热带季风气候区,降水丰沛,气候温和。该区海拔落差大,西部和南部地区最高海拔达1 500 m左右; 北部和中部河谷地带海拔相对较低,都在800 m以下(图1)。区内坡度总体东北部较低,西南部和南部部分地区坡度相对较高。由于该流域为山区,平地分布相对较少。

图1   包溪河流域海拔

Fig.1   Altitude of the Baoxi River Basin

1.2 数据源及其预处理

数字高程模型(digital elevation model,DEM)和遥感影像来源于中国科学院计算机网络信息中心地理空间数据云平台(http://www.gscloud.cn)。NDVI验证数据[11]来源于中国科学院资源环境科学数据中心(http://www.resdc.cn),空间分辨率为1 km。DEM空间分辨率为30 m,遥感影像为Landsat8产品,空间分辨率为30 m。遥感影像数据标识为LC81260412016242LGN00,数据类型为OLI_TIRS,中心经纬度为 27.43°N,108.46°E ,获取日期为2016年8月29日,云量为4%。影像质量好,可较好反映植被生长期的NDVI。

由于NDVI数据、地形数据和验证NDVI数据投影坐标系不一致,直接采用ArcGIS软件中动态投影(或公共坐标系)运算,但可能会产生微小差异。基础数据处理后,统一