植被遥感时间序列数据重建方法简述及示例分析
A review of reconstruction methods for remote-sensing-based time series data of vegetation and some examples
通讯作者: 侯美亭(1982-),男,高级工程师,博士,主要研究方向为生态遥感。Email:houmt@outlook.com。
责任编辑: 李瑜
收稿日期: 2021-03-15 修回日期: 2021-05-28
基金资助: |
|
Received: 2021-03-15 Revised: 2021-05-28
作者简介 About authors
李伟光(1981-),男,高级工程师,主要研究方向为遥感在气象灾害监测中的应用。Email:
随着遥感数据的不断积累,植被遥感产品逐渐形成了完善的时间序列数据,这些数据对阐明生态系统动态变化及分析有关的驱动因素具有重要价值。然而,云遮挡、仪器误差等因素严重制约着植被遥感产品的观测质量,往往造成连续观测数据的缺失。对存在数据缺失的序列进行时空重建是准确提取序列变化特征的重要前提,时空重建就是充分利用遥感数据的时空相关性对数据缺失进行插值以及平滑滤波,以重建完整时间序列。本研究主要以植被遥感时间序列数据为例,对以往常用的时间序列重建方法进行了简要回顾,时间序列重建大致包括插值、平滑2大步骤,插值又分为基于时间的、基于空间的、时空相结合的插值3大类型。然后,以模拟的归一化植被指数(normalized difference vegetation index,NDVI)时间序列、GIMMS NDVI真实时间序列为例,并随机制造不同比例的数据缺失,对比了线性插值、奇异谱分析、Whittaker和时间序列谐波分析等4种常用的不同类型的数据重建方法的时间序列重建效果。结果显示,4种方法的表现各有优劣,但Whittaker方法显示了整体更好的性能。由于随区域的不同,插值方法的表现可能有所差异,故不同数据重建方法还有待进一步推广验证。
关键词:
Remote-sensing-based time series data of vegetation have been increasingly available with the accumulation of remote sensing data. These data are vital for ascertaining the changes in an ecosystem and analyzing relevant driving factors. However, some factors (e.g., cloud cover and instrument errors) restrict the observation quality of the vegetation products of remote sensing, creating data gaps in continuous and high-quality observation data. The data gaps can be filled based on the spatio-temporal dependence of the earth’s surface characteristics, which is called the spatio-temporal reconstruction of time series data. High-quality spatio-temporal reconstruction of time series data is an important prerequisite for the accurate extraction of changes in time series data. Taking the remote-sensing-based time series data of vegetation indices as examples, this study briefly reviewed the widely used reconstruction methods of time series data firstly. These methods generally include two steps: interpolation and smoothing. The interpolation can be divided into three major types, namely time-based, space-based, and spatio-temporal interpolation. Then, taking the simulated normalized vegetation index (NDVI) time series and actual GIMMS NDVI time series as examples, different proportions of data gaps in the two time series were created. Then, this study compared the effects of four types of data reconstruction methods (i.e., linear interpolation, singular spectrum analysis (SSA), Whittaker, and time series harmonic analysis (HANTS)) on the reconstruction results of the two time series. The results show that the four methods have their own advantages and disadvantages, and the Whittaker method showed relatively good performance overall. However, the performance of interpolation methods might vary within different regions, and thereby the data reconstruction methods need to be further verified.
Keywords:
本文引用格式
李伟光, 侯美亭.
LI Weiguang, HOU Meiting.
0 引言
卫星时间序列数据已被证明对阐明生态系统动态变化和驱动因素特别有价值[5,6]。Cord等[7]汇总了卫星遥感应用于生态系统监测的服务类别,指出食物供应和气候服务(温室气体调节)是卫星遥感应用最为广泛的2大生态服务领域,土地利用/土地覆盖(land use and land cover,LULC)和归一化植被指数(normalized difference vegetation index,NDVI)是最为常用的2类卫星遥感产品。NDVI与绿叶密度、净初级生产力和CO2通量高度相关,是地表生物量在景观水平上的一个代表[8]。尽管随遥感技术的发展,空间、光谱和时间分辨率不断增强,对植被空间异质性和连续性的理解不断加深,但基于遥感对离散的植被单元进行分类一直是LULC分类方法的重要基础[9]。遥感反演的植被产品是遥感应用于生态监测领域重要的数据之一。
伴随着遥感数据的不断积累,植被遥感产品形成了完善的时间序列数据,利用这些数据分析生态系统变化,例如评估生态系统过程中的扰动、检测植物物候学的时空变化、分析气候变化对植被的影响、监测LULC变化、植被时空异质性评价等,基本上构成了近几十年来最为活跃的生态遥感研究主题。植被覆盖和生产力的年际变化被认为是理解全球气候变化对生态系统动力学影响的关键[10]。目前,研究者已经开发出许多软件或程序来处理遥感时间序列数据,例如TIMESAT[11],BFAST[12],Timestats[13],SPIRITS[14],DBEST[15],phenor[16],FORCE[17],BEAST[18]和DATimeS[19]等。
1 植被遥感时间序列数据重建方法
卫星数据时间序列可以看作是一种3-D数组,对于这种时间序列,最常见的问题就是数据缺失,也可称为数据异常。数据缺失问题是遥感数据预处理中面临的最主要问题之一。缺失数据对定量遥感研究的影响会很严重,导致参数估计有偏差、信息丢失、可信度下降等。导致数据缺失/异常的原因可以归纳为气象条件(云、雪、灰尘和气溶胶)、仪器误差/传感器退化、数据传输过程中图像数据丢失、时间分辨率低等。其中,传感器退化引起的数据异常往往导致时间序列存在某种趋势,对于这种因素往往无法通过常见的数据重建方法进行校正,需要数据生产方对传感器进行评估和校准,比如MODIS C6版本数据(较C5 版本)对Terra/MODIS退化的校正[26]。本研究中的数据重建方法主要针对传感器退化因素以外引起的数据缺失和异常。
数据缺失/异常产生的噪声包含了不真实的时间变异,对于这类噪声,以往研究大致有以下2类处理方法: ①利用平滑方法,对包含噪声的原始时间序列直接进行平滑处理。例如,Mondal等[27]使用基于离散傅里叶变换(discrete Fourier transform,DFT)的平滑方法来处理原始的MOD13Q1数据,减少NDVI时间序列中的噪声,相比其他平滑技术,DFT需要的模型参数更少。②只使用原始时间序列中的高质量数据,通过数据插值和平滑方法,重建时间序列。原始数据的质量可基于质量标记产品来判定,或者基于相应数据点的标准偏差σ,一般认为位于平均值±3σ以外的数据点是低质量的[28,29,30]。例如,Jeganathan等[29]将任一像元1 a内特定时段的GIMMS NDVI值与其历年相应时段的多年平均值进行比较,然后移除NDVI值在平均值±3σ以外的那些数据点,经过这些处理后,剩余的数据点可认为是“高质量”的,再进一步利用插值方法来填补上述低质量数据消除过程中产生的数据空白。与第一类方法相比,第二类方法在数据重建时首先直接去除了异常值,相对更为可靠。本研究将着重概述第二类方法有关的研究进展。对于该类遥感时间序列数据重构方法,简而言之,包括插值和平滑2步,即通过利用遥感数据的时空相关性对已有的存在数据缺失的高质量数据序列进行插值、平滑滤波,以达到重建时间序列的目的。
1.1 数据插值
基于时间的插值方法主要是考虑到所插值的数据是与生物过程时间动态直接联系在一起的[32],应用尤为广泛。比较简单的基于时间的数据填补方法,可以由数据缺失像元的历史同期平均或最近的有效数据(特别是对于受到积雪影响的植被数据)替代[33],或者使用2个相邻有效值的移动平均值替换可能由云和大气气溶胶等因素引起的局部不规则变化和短间隙[34,35]。更为常用的时间插值方法为线性插值[36]和样条插值[37],其中,样条法也可以用于空间插值[38]和数据序列平滑[39]。与线性插值方法相比,时间序列谐波分析(harmonic analysis of time series,HANTS)、奇异谱分析(singular spectrum analysis,SSA)是更为复杂的时间插值方法。HANTS算法是快速傅里叶变换(fast Fourier transform,FFT)的一种改进,采用基于谐波分量的最小二乘曲线拟合方法,在植被遥感时间序列插值中展示了良好的性能[40,41]。利用SSA进行插值是一种近年来发展的数据插值方法,其原理在于使用迭代方法对缺失数据进行估计,然后使用这些估计来计算滞后协方差矩阵及其经验正交函数分解分量,再使用交叉验证来优化窗口宽度和SSA主模态的数量,以填补数据空白。尽管SSA插值方法在植被遥感数据中应用不多[42],但在海温、降水和土壤呼吸等多种类型的地学数据插值中显示了其价值[43,44]。然而,类似SSA和HANTS这样的方法是基于数据序列周期性的,如果数据缺失持续的时间长,就不能成功识别周期及振荡,从而导致利用这些方法得到的结果缺乏可信度[45]。
基于空间的插值方法主要是利用存在数据缺失的像元周边的(不存在数据异常的)像元信息进行插值,除样条法以外,常用的方法还包括最邻近插值、反距离加权插值、地统计方法(如克里格法)等。需要注意的是,数据缺失像元通常存在空间聚集的情况,在进行空间插值时,有必要评价空间自相关效应[46]。而且,由于混合像素的存在,植被指数在短距离内变化可能很大,导致基于空间的插值可能无法代表真实景观的复杂性。
尽管目前有很多种时间和空间插值方法,但这2类方法都不能在所有情况下超越对方,而时间和空间信息相结合的插值方法正在显现出更多的优势。Borak等[47]以MODIS LAI时间序列为例,对几种时间和空间插值技术进行了评估,发现插值技术的精度取决于下垫面覆盖,时间插值在非森林覆盖类型上执行得更有效,空间插值在森林区域执行得更好,而空间和时间方法的结合提供了比任何单一的时间和空间方法更好的插值能力。Pede等[38]发现,对于受到云覆盖影响的MODIS 8 d LST数据,空间插值方法、空间和时间相结合的插值方法比时间插值方法显示出了更大的优势,而且空间和时间相结合的方法受地形、气候等环境因素的影响更小。近年来更多研究使用来自邻近像素的历史数据和时间曲线来结合空间和时间信息进行插值,也显现出了比仅用时间或空间插值更大的优势[48,49]。Yan等[50]提出了一种应用于高空间分辨率、低时间分辨率的Landsat数据进行间隙填补的时空相似替代算法,通过合并空间和时间上相邻的相似像素来分割图像时间序列,然后将这些片段分组到不同的类中,需要填充的像素位置与属于同一类的其他像素进行比较,以识别替代像素,不过其算法不考虑地表变化趋势,因此较适用于在地表没有发生变化期间获取的图像。
最近还涌现了一些基于机器学习的插值算法,如神经网络、随机森林和决策树等,这些方法不同于传统的参数方法(如线性方法),信息来自于训练数据,不需要假设它们的统计分布或变量之间的相互关系,但这些方法的效果如何,还有待深入评价[19]。总体上,插值的效果不仅依赖于插值算法,更依赖于原始数据的质量,当原始数据序列具有较大的高质量数据比例时(比如大于85%),普通的线性插值一般就可以满足需要[51]。而且,对于常用的GIMMS NDVI数据,Burrell等[52]认为当原始时间序列的有效值大于85%时,质量即满足需求,可以直接使用原始序列。另外,在一些情形下(如物候特征提取),当某像元缺值数据超过一定数量时,比如连续缺值超过1个月,Pastor-Guzman等[53]选择不使用相应位置的数据,而通过插值方法获得相应缺失的数据。
1.2 时间序列平滑
由于插值算法本身局限性以及原始时间序列数据缺失程度的影响,插值后的时间序列数据可能仍然包含噪声,这时候需要进行平滑,平滑的主要目的是降噪,以保证重建时间序列的质量。
简单的平滑方法包括最大值合成(maximum value composite,MVC)方法、通过使用预定义的阈值移除异常值的最佳指数斜率提取算法(best index slope extraction,BISE)[54]、迭代插值方法(iterative interpolation for data reconstruction,IDR)[32]。Holben[55]是最早提出使用MVC在时间序列中抑制NDVI噪声的研究之一,而且MVC产品还减少了NDVI时间序列中云污染像素造成的数据间隙。然而,即使在MVC过程之后,特别是当云层覆盖持续时间长于合成期(如在雨季)时,一些残留的噪声可能仍然会出现在数据中。这时就需要使用更加稳健的平滑方法对序列进行降噪。
曲线拟合方法是目前使用较多的平滑方法[36],如SG(Savitzky-Golay)、TSGF(temporal smoothing and gap filling)、非对称高斯函数(asymmetric Gaussian,AG)、双逻辑函数(double logistic,DL)和傅里叶变换等应用于局部时间窗口的曲线拟合方法[56],以及可应用于整个时间序列的Whittaker滤波方法(一种样条平滑方法)等。其中,AG,DL和SG方法被集成到时间序列分析软件(TIMESAT)中,使得这3种方法更易于实现[57]。曲线拟合方法除了应用于平滑降噪外,也被用于提取植被物候变化参数,因为拟合最主要的依据就是植被时间变化特征。特别是在局部窗口内拟合卫星数据时间序列,多采用高阶谐波方法或更多参数,能够更好地捕捉植被动态,尤其适于不规则、不对称的植被指数曲线,不过它们对局部波动和数据噪声也更为敏感,某些情况下,可以对经过局部拟合的且仍然包含有不规则数据的序列,再应用局部中值滤波去除局部不规则值[34]。相比基于整个时间序列或低阶谐波的方法往往得到更平滑的曲线,但也可能偏离实际植被生长轨迹,导致过度校正,抑制了实际植被波动。
每种时间序列平滑算法都有其独特的优点和缺点,在选择平滑算法时,不仅应考虑时间序列的应用目的、数据中噪声的性质、时间序列曲线的典型形状和使用方法的方便性等[56],还应考虑原始时间序列的数据质量[58]。如果数据质量足够高,可以选择保持季节内变化的滤波方法,如SG滤波[22, 59]和Whittaker滤波[31, 60]; 在高质量数据较少的情况下,比如观测频率相对较低的Landsat和Sentinel-2数据,函数拟合(如DL)可能是实现稳健性的唯一选择[61]。另外,好的插值算法也能间接提高平滑方法的质量。比如,Padhee等[62]基于局地(尤其是同一土地覆盖类型)物候时空相似性的考虑,使用参考物候曲线来预测NDVI时间序列中的缺失值的方法对MODIS NDVI时间序列进行预插值,改善了云高覆盖地区时间序列中存在较大数据间隙而导致的HANTS方法的不足。
尽管目前没有公认的标准来比较不同的时间序列重建方法[63],但是可以通过评估重建后的时间序列对原始无云数据的保真度和识别受云污染值的能力来判断方法合适与否。由于地面真实数据通常不可用,Julien等[32]以GIMMS NDVI为例,提出了以下判断标准: 重建后的值与原始观测值的距离、以及与上包络线的接近度。与原始观测值的距离旨在量化对未污染数据的保真度。对于上包络近似准则,它基于大气污染有降低NDVI值的趋势这一考虑,这一标准在整个时间序列中被定义为NDVI低于原始观测值的重建观测值的频率(范围在0~1)。对于该准则,较低的值对应于更接近原始时间序列的上包络线。这2个标准对于评价时间序列重建方法具有较好的效果,特别是当大气污染或云对遥感反演的生物物理参数有负面影响时。
2 常用的时间序列数据重建方法比较示例
这里选取线性插值、SSA,Whittaker和HANTS等4类常用的不同类型的时间序列数据重建方法,通过在模拟的NDVI时间序列、原始NDVI完整时间序列中随机制造不同比例的数据缺失,使用上述4类方法进行数据重建,然后计算重建时间序列与原始时间序列的均方根误差(root mean square error,RMSE),以衡量重建序列同原始序列之间的偏差,从而判断不同方法的优劣。
RMSE亦称标准误差,值越小,拟合效果越好。计算方法为:
式中: y为原始序列; z为拟合序列; N为序列长度。
2.1 模拟NDVI序列
模拟NDVI序列包含短周期的波动和长期趋势[64],按照下式生成为:
图1
图1
模拟NDVI序列数据重建效果对比
Fig.1
Comparison of reconstruction results of simulated NDVI time series
表1 基于模拟NDVI序列的不同数据重建方法的RMSE对比
Tab.1
RMSEall | RMSEgap | |||||||
---|---|---|---|---|---|---|---|---|
数据缺失比例 /% | 线性插值 | SSA | Whittaker | HANTS | 线性插值 | SSA | Whittaker | HANTS |
10 | <0.000 1 | 0.096 6 | <0.000 1 | 0.045 6 | 0.012 6 | 0.105 4 | 0.000 8 | 0.061 2 |
20 | <0.000 1 | 0.004 8 | <0.000 1 | 0.045 4 | 0.014 7 | 0.004 0 | 0.000 8 | 0.067 7 |
30 | <0.000 1 | 0.077 5 | <0.000 1 | 0.043 6 | 0.015 0 | 0.085 3 | 0.001 0 | 0.073 1 |
40 | <0.000 1 | 0.102 3 | <0.000 1 | 0.023 8 | 0.024 1 | 0.110 0 | 0.002 4 | 0.049 2 |
2.2 原始NDVI时间序列
图2
图2
NDVI序列不同数据重建方法效果对比
Fig.2
Comparison of different reconstruction methods for NDVI time series
图3
图3
NDVI序列数据重建方法效果局部放大对比
Fig.3
Local detailed comparison of different reconstruction methods for NDVI time series
表2 基于GIMMS NDVI时间序列的不同数据重建方法的RMSE
Tab.2
RMSEall | RMSEgap | |||||||
---|---|---|---|---|---|---|---|---|
数据缺失比例 /% | 线性插值 | SSA | Whittaker | HANTS | 线性插值 | SSA | Whittaker | HANTS |
10 | <0.000 1 | 0.034 1 | 0.007 8 | 0.044 4 | 0.048 2 | 0.052 3 | 0.044 0 | 0.055 5 |
20 | <0.000 1 | 0.035 8 | 0.026 1 | 0.042 3 | 0.056 5 | 0.043 0 | 0.056 6 | 0.061 9 |
30 | <0.000 1 | 0.034 0 | 0.000 9 | 0.042 3 | 0.058 4 | 0.046 1 | 0.060 4 | 0.057 3 |
40 | <0.000 1 | 0.031 2 | <0.000 1 | 0.040 4 | 0.061 0 | 0.045 3 | 0.066 0 | 0.061 2 |
图2和表2显示了针对GIMMS NDVI序列不同数据缺失比例下的4种数据重建方法的性能,图3截取了其中一段展示其局部效果。可以发现,线性方法保持了模拟序列时的表现,随着数据缺失比例增加,RMSEgap逐渐上升。而处理模拟NDVI序列时表现最好的Whittaker方法依然保持了良好的性能,而且其只需要通过设置唯一参数即可平衡拟合的保真度和粗糙度,在保证拟合精度的同时,也提高了曲线的平滑度,不过,该方法的RMSEgap有时稍高于SSA。相比,在该示例中,SSA成为最好的数据重建方法,而且随缺失数据的增多,误差变化不大。HANTS方法对原始序列的保真和断点位置的插值基本稳定,但该方法参数设置较为复杂,需要多次试验才能确定合适的参数,提高平滑度的同时也牺牲了拟合的精度,而且其RMSE在不同缺失比例下都要高于SSA,这也与Ghafarian Malamiri等[42]的研究相一致。
3 结论
完整的高质量植被遥感时间序列是利用遥感准确提取植被生态系统变化特征的重要前提。然而,云覆盖等因素严重制约着植被遥感产品的观测质量,造成时间序列中存在不同程度的数据缺失或低质量数据。本研究回顾了以往常用的时间序列重建方法,这些方法可以概括为2大类型: ①利用平滑方法对包含噪声的原始时间序列直接进行平滑处理; ②仅使用原始时间序列中的高质量数据,通过数据插值和平滑方法,重建完整时间序列。后者是目前数据重建的主流方法,且已经发展了相应的多种多样的数据插值和平滑方法。
在回顾以往数据重建方法的基础上,本研究以在模拟的NDVI时间序列、原始NDVI完整时间序列中随机制造不同比例的数据缺失为例,具体对比了线性插值、SSA方法、Whittaker方法和HANTS方法等4类方法的数据重建性能,研究发现,这些方法在NDVI数据缺失插补中均取得了较好的精度。然而,它们有各自的优点和缺点,在不同的情况下填补数据空缺时需要加以考虑。其中,线性插值是一种广泛使用的方法,然而当缺失数据过多时,这种方法并不适用; 相比Whittaker 方法、SSA方法和HANTS方法显示出了更好的性能,尽管一些情形下,SSA及HANTS对应的RMSE高于Whittaker方法,但这3种方法的整体差异并不大; 不过,HANTS方法需要设置的参数相对较多,SSA方法在某些情形下的重建序列具有较大的振幅,而Whittaker方法在不同的数据缺失情形下都显示出了较好的效果,且需要的参数较少,是一种值得推广的插值方法。当然,除NDVI以外,EVI,EVI2以及近年来的SIF等也是较为常用的植被遥感产品。由于不同植被遥感时间序列最重要的共同特征之一是具有显著的季节动态,而常见的植被遥感时间序列重建方法一般具有很好的普适性(比如Whittaker 方法),也多基于序列周期性(比如HANTS方法),因此它们也适用于NDVI以外的其他植被遥感时间序列。
总体上,本研究主要侧重于回顾近年来的植被遥感时间序列重建方法,但仅采用了模拟数据和一个样点的遥感数据进行时间序列重建方法比较。而随区域的不同,插值方法的表现可能有所差异,故研究结果还有待进一步推广试验。
(责任编辑: 李瑜)
参考文献
What does remote sensing do for Ecology?
[J]. ,DOI:10.2307/1941546 URL [本文引用: 1]
North American vegetation patterns observed with the NOAA-7 advanced very high resolution radiometer
[J]. ,DOI:10.1007/BF00033449 URL [本文引用: 1]
Overview of solar-induced chlorophyll fluorescence (SIF) from the orbiting carbon observatory-2:Retrieval,cross-mission comparison,and global monitoring for GPP
[J]. ,DOI:10.1016/j.rse.2018.02.016 URL [本文引用: 1]
Agree on biodiversity metrics to track from space
[J]. ,DOI:10.1038/523403a URL [本文引用: 1]
Comparative study of three satellite image time-series decomposition methods for vegetation change detection
[J]. ,DOI:10.1080/22797254.2018.1465360 URL [本文引用: 1]
中国西南喀斯特地区植被变化时空特征及其成因
[J]. ,
Spatial and temporal patterns of changing vegetation and the influence of environmental factors in the Karst region of Southwest China
[J].
Priorities to advance monitoring of ecosystem services using earth observation
[J]. ,DOI:10.1016/j.tree.2017.03.003 URL [本文引用: 1]
Calibration of remotely sensed,coarse resolution NDVI to CO2 fluxes in a sagebrush-steppe ecosystem
[J]. ,DOI:10.1016/S0034-4257(03)00004-X URL [本文引用: 1]
Remote sensing of plant functional types
[J]. ,DOI:10.1111/nph.2010.186.issue-4 URL [本文引用: 1]
Satellite observed widespread decline in Mongolian grasslands largely due to overgrazing
[J]. ,DOI:10.1111/gcb.2014.20.issue-2 URL [本文引用: 1]
TIMESAT:A program for analyzing time-series of satellite sensor data
[J]. ,DOI:10.1016/j.cageo.2004.05.006 URL [本文引用: 1]
Detecting trend and seasonal changes in satellite image time series
[J]. ,DOI:10.1016/j.rse.2009.08.014 URL [本文引用: 1]
TimeStats:A software tool for the retrieval of temporal patterns from global satellite archives
[J]. ,DOI:10.1109/JSTARS.4609443 URL [本文引用: 1]
Image time series processing for agriculture monitoring
[J]. ,DOI:10.1016/j.envsoft.2013.10.021 URL [本文引用: 1]
Detecting changes in vegetation trends using time series segmentation
[J]. ,DOI:10.1016/j.rse.2014.09.010 URL [本文引用: 1]
An integrated phenology modelling framework in R
[J]. ,DOI:10.1111/mee3.2018.9.issue-5 URL [本文引用: 1]
FORCE-Landsat + Sentinel-2 analysis ready data and beyond
[J]. ,DOI:10.3390/rs11091124 URL [本文引用: 1]
Detecting change-point,trend,and seasonality in satellite time series data to track abrupt changes and nonlinear dynamics:A Bayesian ensemble algorithm
[J]. ,DOI:10.1016/j.rse.2019.04.034 URL [本文引用: 1]
DATimeS:A machine learning time series GUI toolbox for gap-filling and vegetation phenology trends detection
[J]. ,DOI:10.1016/j.envsoft.2020.104666 URL [本文引用: 3]
Noise reduction of NDVI time series:An empirical comparison of selected techniques
[J]. ,DOI:10.1016/j.rse.2008.09.003 URL [本文引用: 1]
Performance of smoothing methods for reconstructing NDVI time-series and estimating vegetation phenology from MODIS data
[J]. ,DOI:10.3390/rs9121271 URL [本文引用: 1]
Comparison of eight techniques for reconstructing multi-satellite sensor time-series NDVI data sets in the Heihe River basin,China
[J]. ,DOI:10.3390/rs6032024 URL [本文引用: 2]
Global evaluation of gap-filling approaches for seasonal NDVI with considering vegetation growth trajectory,protection of key point,noise resistance and curve stability
[J]. ,DOI:10.1016/j.rse.2016.11.023 URL [本文引用: 1]
长时间序列NDVI数据重建方法比较研究进展
[J]. ,
Advance in method comparison of reconstructing remote sensing time series data sets
[J].
遥感时间序列数据滤波重建算法发展综述
[J]. ,
Review on methods of remote sensing time-series data reconstruction
[J].
Reanalysis of global terrestrial vegetation trends from MODIS products:Browning or greening?
[J]. ,DOI:10.1016/j.rse.2016.12.018 URL [本文引用: 1]
Mountain agriculture extraction from time-series MODIS NDVI using dynamic time warping technique
[J]. ,DOI:10.1080/01431161.2018.1444289 URL [本文引用: 1]
Global relationships among traditional reflectance vegetation indices (NDVI and NDII),evapo-transpiration (ET),and soil moisture variability on weekly timescales
[J]. ,DOI:10.1016/j.rse.2018.10.020 URL [本文引用: 1]
Remotely sensed trends in the phenology of northern high latitude terrestrial vegetation,controlling for land cover change and vegetation type
[J]. ,DOI:10.1016/j.rse.2013.11.020 URL [本文引用: 2]
Monitoring grassland seasonal carbon dynamics,by integrating MODIS NDVI,proximal optical sampling,and eddy covariance measurements
[J]. ,DOI:10.3390/rs8030260 URL [本文引用: 1]
A perfect smoother
[J]. ,The well-known and popular Savitzky-Golay filter has several disadvantages. A very attractive alternative is a smoother based on penalized least squares, extending ideas presented by Whittaker 80 years ago. This smoother is extremely fast, gives continuous control over smoothness, interpolates automatically, and allows fast leave-one-out cross-validation. It can be programmed in a few lines of Matlab code. Theory, implementation, and applications are presented.
Comparison of cloud-reconstruction metho-ds for time series of composite NDVI data
[J]. ,DOI:10.1016/j.rse.2009.11.001 URL [本文引用: 3]
Increasing altitudinal gradient of spring vegetation phenology during the last decade on the Qinghai-Tibetan Plateau
[J]. ,DOI:10.1016/j.agrformet.2014.01.003 URL [本文引用: 1]
Reconstruction of a complete global time series of daily vegetation index trajectory from long-term AVHRR data
[J]. ,DOI:10.1016/j.rse.2014.10.012 URL [本文引用: 2]
Land surface phenology from MODIS:Characterization of the Collection 5 global land cover dynamics product
[J]. ,DOI:10.1016/j.rse.2010.04.005 URL [本文引用: 1]
A comparison of methods for smoothing and gap filling time series of remote sensing observations:Application to MODIS LAI products
[J]. ,DOI:10.5194/bg-10-4055-2013 URL [本文引用: 2]
Extracting the spatiotemporal pattern of cropping systems from NDVI time series using a combination of the spline and HANTS algorithms:A case study for Shandong Province
[J]. ,DOI:10.1080/07038992.2017.1252906 URL [本文引用: 1]
An empirical comparison of interpolation methods for MODIS 8-day land surface temperature composites across the conterminous Unites States
[J]. ,DOI:10.1016/j.isprsjprs.2018.06.003 URL [本文引用: 2]
Evaluation of the potential of MODIS satellite data to predict vegetation phenology in different biomes:An investigation using ground-based NDVI measurements
[J]. ,DOI:10.1016/j.rse.2013.01.010 URL [本文引用: 1]
Reconstruction of global MODIS NDVI time series:Performance of harmonic analysis of time series (HANTS)
[J]. ,DOI:10.1016/j.rse.2015.03.018 URL [本文引用: 1]
Optimizing and comparing gap-filling techniques using simulated NDVI time series from remotely sensed global data
[J]. ,DOI:10.1016/j.jag.2018.11.008 URL [本文引用: 1]
Comparison of harmonic analysis of time series (HANTS) and multi-singular spectrum analysis (M-SSA) in reconstruction of long-gap missing data in NDVI time series
[J]. ,DOI:10.3390/rs12172747 URL [本文引用: 2]
Spatio-temporal filling of missing points in geophysical data sets
[J]. ,DOI:10.5194/npg-13-151-2006 URL [本文引用: 1]
Meissner H.Gap-filling continuously-measured soil respiration data:A highlight of time-series-based methods
[J]. ,DOI:10.1016/j.agrformet.2020.107912 URL [本文引用: 1]
The operational method of filling information gaps in satellite imagery using numerical models
[J]. ,DOI:10.1016/j.jag.2018.09.002 URL [本文引用: 1]
Spatial heterogeneity of greening and browning between and within bioclimatic zones in northern west Siberia
[J]. ,DOI:10.1088/1748-9326/11/11/115002 URL [本文引用: 1]
Effective interpolation of incomplete satellite-derived leaf-area index time series for the continental United States
[J]. ,DOI:10.1016/j.agrformet.2008.08.017 URL [本文引用: 1]
Interpolation of the mean anomalies for cloud filling in land surface temperature and normalized difference vegetation index
[C]. ,
A simple method to improve the quality of NDVI time-series data by integrating spatiotemporal information with the Savitzky-Golay filter
[J]. ,DOI:10.1016/j.rse.2018.08.022 URL [本文引用: 1]
Spatially and temporally complete Landsat reflectance time series modelling:The fill-and-fit approach
[J]. ,DOI:10.1016/j.rse.2020.111718 URL [本文引用: 1]
Large-scale detection of vegetation dynamics and their potential drivers using MODIS images and BFAST:A case study in Quebec,Canada
[J]. ,DOI:10.1016/j.rse.2017.11.017 URL [本文引用: 1]
The impact of dataset selection on land degradation assessment
[J]. ,DOI:10.1016/j.isprsjprs.2018.08.017 URL [本文引用: 1]
Remote sensing of mangrove forest phenology and its environmental drivers
[J]. ,DOI:10.1016/j.rse.2017.11.009 URL [本文引用: 1]
Influence of heterogeneous landscapes on computed green-up dates based on daily AVHRR NDVI observations
[J]. ,DOI:10.1016/j.rse.2009.07.020 URL [本文引用: 1]
Characteristics of maximum-value composite images from temporal AVHRR data
[J]. ,DOI:10.1080/01431168608948945 URL [本文引用: 1]
A review of vegetation phenological metrics extraction using time-series,multispectral satellite data
[J]. ,DOI:10.1016/j.rse.2019.111511 URL [本文引用: 2]
A method for robust estimation of vegetation seasonality from Landsat and Sentinel-2 time series data
[J]. ,DOI:10.3390/rs10040635 URL [本文引用: 1]
A simple method for reconstructing a high-quality NDVI time-series data set based on the Savitzky-Golay filter
[J]. ,DOI:10.1016/j.rse.2004.03.014 URL [本文引用: 1]
A time series for monitoring vegetation activity and phenology at 10-daily time steps covering large parts of South America
[J]. ,DOI:10.1080/17538947.2010.505664 URL [本文引用: 1]
Detecting interannual variation in deciduous broadleaf forest phenology using Landsat TM/ETM+ data
[J]. ,DOI:10.1016/j.rse.2013.01.011 URL [本文引用: 1]
Spatio-temporal reconstruction of MODIS NDVI by regional land surface phenology and harmonic analysis of time-series
[J]. ,DOI:10.1080/15481603.2019.1646977 [本文引用: 1]
Satellite-derived Normalized Difference Vegetation Index (NDVI) is frequently obstructed by adverse atmospheric components resulting in data gaps in time series. The harmonic analysis of time series (HANTS) algorithm is a widely used technique to reconstruct missing NDVI time series. However, due to restriction of HANTS to act within temporal dimension, its direct application is bound to endure practical problems in spatiotemporal reconstruction due to large data gaps. This study proposes Moving Offset Method (MOM), a novel prefilling method applied on NDVI time series prior to application of HANTS. MOM restores the missing NDVI time series by assuming that it tends to follow a reference pattern of land surface phenology (NDVIref). The NDVIref is prepared by using a recursive search and fill algorithm (SFA) for data availability without null values. It restores null values in NDVIref at a pixel by using coefficients of linear regression with NDVIref at another pixel having identical conditions. Finally, the prefilling is prior to application of HANTS. The proposed approach is demonstrated by using MODIS 16-daily time series data for Northeast India and Bhutan region which is covered with frequent seasonal clouds. Besides direct application of HANTS, it is also compared with similar approaches which includes prefilling by inverse distance weighted (IDW) and cubic spline, prior to application of HANTS. The fitting indicators, overall reconstruction error (ORE) and normalized noise related error (NNRE) are found to be best for proposed approach in spatiotemporal comparison. Also, restoration of seasonality trait the NDVI time series better in the proposed approach. This approach is concluded to be an enhancement for HANTS that could be helpful in improving quality of NDVI reconstruction for regions with frequent seasonal obstructions around the globe.
An overview of global leaf area index (LAI):Methods,products,validation,and applications
[J]. ,DOI:10.1029/2018RG000608 URL [本文引用: 1]
Phenological change detection while accounting for abrupt and gradual trends in satellite image time series
[J]. ,DOI:10.1016/j.rse.2010.08.003 URL [本文引用: 1]
/
〈 | 〉 |