两种NDVI时间序列数据拟合方法比较.pdfVIP

第25卷 第 11期 农 业 工 程 学 报 VO1．25 No．11 2009年 11月 TransactionsoftheCSAE NOV．2009 183 两种NDVI时间序列数据拟合方法比较 吴文斌 1,2,3，杨 鹏 ，唐华俊 ，，周清波 1,2，ShibasakiRyosuke。，张 莉 ， (1．中国农业科学院农业资源与农业区划研究所，北京 100081； 2．农业部资源遥感与数字农业重点 开放实验室，北京100081： 3．日本东京大学空间情报科学研究中心，日本东京 153—8505) 摘 要：归一化植被指数 (NDVI)时间序列数据拟合 目的是降低时序数据的噪声水平，重建高质量的NDVI时序数据， 有利于多种参数反演和信息提取。针对国际上普遍应用的两种NDVI时间序列数据拟合方法，即Savitzky·Golay滤波法 和非对称性高斯函数拟合法，该文在介绍两种方法基本概念的基础上，利用直接比较法和间接比较法在中国对两种拟合 方法进行了比较分析。结果表明，Savi~ky—Golay滤波法和非对称性高斯函数拟合法的拟合效果总体上一致，但二者之 间还是存在区域差异性，这种区域差异与两种方法的自身特点和中国区域 自然条件紧密相关。不同数据拟合方法的比较 研究可以弄清每种方法的优缺点和区域适宜性，有助于研究人员针对不同研究 目的和研究区域选择适宜的NDVI数据拟 合方法，减少遥感数据处理中的误差，提高研究精度。 关键词：时间序列分析，滤波，函数，拟合，Savitzky-Golay滤波法，非对称性高斯函数拟合法 doi：10．3969~．issn．1002～6819．2009．11．033 中图分类号：TP79：Q948．1 文献标识码：A 文章编号：1002—6819(2009)一1卜Ol83一O6 吴文斌，杨 鹏，唐华俊，等．两种NDVI时间序列数据拟合方法比较 [J]．农业工程学报，2009，25(I1)：183一I88． WuWenbin，YangPeng，TangHuajtm，eta1．ComparisonoftwofittingmethodsofNDVItimeseriesdatasets[J]．Transactionsof theCSAE，2009，25(11)：183—188．(inChinesewithEnglishabstract) 动周期和 NDVI再增长的百分比阂值主要取决于分析者 0 引 言 的经验和要求，人为影响因素大，使得大多数阈值去除 卫星传感器所获取的归一化植被指数 (NDVI)时间序 法无法得到普遍应用。因此，基于滤波的平滑方法 (如 列数据能够精确地反映地表植被的活力状态和植被季相 Savitzky．Golay滤波法[、小波变换 【】和傅立叶变换 【“1) 变化特征 j，已经在全球、洲际以及地区范围的植被覆盖 和非线性拟合法(如Logistic函数拟合法 j【_2J和非对称性高 制图、陆地生态系统动态监测和模拟方面发挥了重要作 斯函数拟合法f])成为 NDVI时序数据平滑处理的主要 用2【。目前常用的NDVI时间序列数据多是利用国际通 方法，已经广泛应用于全球不同区域。为了弄清不同数 用的最大值合成法 (MvC)生成，虽然最大值合成法能 据拟合方法的优缺点和区域适宜性，部分国外学者对不 消除部分云、气溶胶、太阳高度角和地物双向性反射等 同的拟合方法进行了比较研究。Sakamoto等选择 日本为 的干扰 ，但云和大气等的噪声影响仍然很大l5】，使得 研究区域，对小波变换和傅立叶变换方法进行了比较分 NDVI时序数据呈锯齿状的不规则波动，曲线季节变化趋 析，结果表明，小波变换总体上要优于傅立叶变换 l【引。 势不明显，较大地影响了陆地生态系统监测、植被参数