谐波改进的植被指数时间序列重建算法  

Improved algorithm for reconstructing vegetation index image time series based on Fourier Harmonic Analysis

作　　者：张霞[1] 李儒[2] 岳跃民[3] 刘波[1] 刘海霞[1] 

机构地区：[1]遥感科学国家重点实验室,中国科学院遥感应用研究所,北京100101 [2]香港中文大学太空与地球信息科学研究所,香港 [3]中国科学院亚热带农业生态研究所,湖南长沙410125

出　　处：《遥感学报》2010年第3期437-447,共11页NATIONAL REMOTE SENSING BULLETIN

基　　金：国家自然科学基金项目(编号:40971205);国家科技支撑课题(编号:2007BAH15B01);遥感科学国家重点实验室开放基金(编号:03Q00300649)~~

摘　　要：提出一种基于傅里叶谐波分析的改进算法,引入异常值检测算法,检测拟合过程中的异常值,增加数据拟合的真实性;迭代前动态估算出待处理序列点的峰值个数(即频数),解决整个区域预设单一频数的不合理性;引入拟合影响因子,自动控制迭代终止条件,避免传统方法中人为设置阈值导致的不确定性。利用2003年华北平原MODIS_EVI时间序列图像验证表明,较之HANTS算法,改进算法能够有效修正噪声污染像元值,修正后的EVI时序曲线更能反映地物内在的物候变化规律,并能够更好地保真原始曲线上的特征(点),如作物EVI最大值、最小值出现的时间和大小关系。This paper,based on Fourier-transform-based Harmonic Analysis Algorithm,proposes an improved algorithm to overcome the drawback of artificially setting key parameters and to reconstruct high-quality time-series data.Firstly,outlier detection algorithm,instead of threshold setting method,is used to find unreasonable data points before curve fitting.Then,regarding the inherent phenological regulation for each land cover,Numbers of Frequency（NOF） is calculated pixel by pixel by polynomial fitting,which is more reasonable than setting a unique global NOF for the whole scene which contains complex land cover types.Fitting-effect Index,instead of manually setting one fitting tolerance,is employed to decide automatically when to terminate the iteration.The improved method is validated by the MODIS_EVI time series of Huabei plain of 2003.The widely used Harmonic Analysis of Time Series（HANTS） is chosen as a comparison.The result shows that both of the two methods can reflect the phenological regulations of land covers,the reconstructed EVI temporal profile of single-season cropland（e.g.cotton） takes on one peak pattern,and those of double-season cropland and inland water take on double peak pattern and low-steady curve.But the improved method performs better in tracking the change tendency of the original curve.Moreover,the peak time and peak value of croplands are mostly consistent with the original curve,which will be useful for VI-based crop yield prediction.

关 键 词：植被指数图像时间序列 滤波 傅里叶谐波 异常值检测 拟合影响因子 

分 类 号：TP751[自动化与计算机技术—检测技术与自动化装置]