站点尺度的青藏高原时序NDVI重构方法比较与应用-地理科学进展.pdf

第37 卷 第3 期 地理科学进展 Vol.37, No.3 2018 年3 月 ProgressinGeography Mar. 2018 站点尺度的青藏高原时序NDVI 重构方法 比较与应用 1，2 3* 刘建文 ，周玉科 (1. 福州大学空间信息工程研究中心，福州350003 ；2. 数据挖掘与信息共享教育部重点实验室，福州350003 ； 3. 中国科学院地理科学与资源研究所生态系统网络观测与模拟院重点实验室，北京 100101) 摘 要：基于卫星遥感的植被指数时序数据广泛应用于植被覆盖监测、生物量反演等多个研究领域，但由于传感器 本身、大气条件、环境特征等因素引起的噪声会影响数据的应用效果，因此开展植被指数时序数据重构研究具有实 际意义。本文基于2000-2015 年MODIS 归一化差异植被指数(NDVI)数据，采用三次样条函数法、双逻辑斯蒂函数 法和奇异谱分析法3 种常用方法，对青藏高原106个气象站点所在的典型覆被NDVI 时序数据进行重构，并以植被 物候信息提取作为应用，比较分析了3 种算法的保真性、细节拟合能力及物候特征提取效果。研究表明，D-L 及 Spline 函数分别对受冰雪及云层影响较大(荒漠、灌木、林地)及较小的覆被类型(草原、农作物)表现出较好的细节拟 合能力；SSA 方法拟合能力较差，易出现NDVI 重构曲线整体“下移”的现象，造成峰值拟合结果偏低，并且表现出 NDVI 绝对值越小拟合效果越差的现象。从保持原始数据真值的能力来看，受噪声点影响较大的覆被类型(林地、 灌木、草原)Spline 函数略优于D-L 函数法；而林地类型中SSA 方法表现优于D-L 函数法。从物候信息提取结果来 看，D-L 函数法所提取的生长季稍有提前，Spline 函数及SSA 方法分别表现出生长季开始点及结束点滞后的现象， 灌木、林地类型表现出明显的年际波动变化的特征，荒漠类型由于NDVI 绝对值偏低，3 种方法物候提取结果一致 性表现出锯齿状不规则波动。此外，D-L 方法生长季开始期(SOS)和生长季结束期(EOS)在各覆被区均小于其他方 法，波动较大；SSA 方法提取的EOS 在大部分覆被地区大于其他方法；Spline 提取结果的年际波动与SSA 高度相 似。该研究可为高原植被不同覆被类型下NDVI 时序数据噪声去除的方法选择提供借鉴。 关键词：时序数据重构；植被物候；双逻辑斯蒂函数；三次样条函数；奇异谱分析法 1 引言 al, 2002; 聂勇等, 2012) 。理想的植被指数能增强对 植物信息的敏感性，最大限度地反映植被真实的生 基于卫星遥感获取的植被指数能够较好地反 长状况及物候规律。但由于传感器性能各异、光照 映植被的生长态势，已经广泛应用于植被变化监 与观测角度不同、数据传输误差、地形与大气条件 测、全球陆地表面建模、生物量计算等多个领域(Yu 因时而异，观测值通常带有不可预测的噪声，从而 et al, 2004; Julien et al, 2006; 李斌等, 2017) 。气象 影响植被指数的应用。因此，研究NDVI 时间序列 观测数据与植被指数的关系也得到广泛的研究，基 数据的重构方法，尽量真实地刻画植被生长情况， 于多个站点的区域平均植被指数时序数据可有效 是植被生态遥感领域的重要研究课题，也是植被综 表征相应覆被类型区内的植被生长状况(Chen et 合对比研究的热点。 收稿日期：2017-04-12；修订日期：2017-11-08。 基金项目：国家自然科学基金项目 ；国家重点研发计划