一种基于图论重构MODIS EVI时间序列数据集的新方法

摘要/Abstract

摘要： MODIS的增强型植被指数（EVI）时间序列数据早已广泛应用于植被观测、生态环境和全球气象变化等研究领域,但即使EVI时间序列数据已经经过严格的预处理,其中仍然存在着一些噪声。因此,本文开发了一种简单有效的方法来重构EVI时间序列数据,减少EVI 时间序列数据中的噪声,尤其是一些由大气云层和冰雪覆盖产生的噪声。新方法的理论来源于图论,利用拉普拉斯矩阵的关系对EVI中选定的邻域窗口的像元权重进行赋值,得到中心像元的拟合。新方法已应用于2016—2018 年的MODIS MOD13A1 产品,并与S-G滤波法、谐波函数法、双逻辑斯蒂拟合法和非对称高斯函数法进行了比较。结果表明,在荒漠、草原和林地中,新方法留一验证测试的绝对差值最小,相较于其他方法效果较优;在拟合不同植被类型的EVI时间序列数据时,图论邻点方法呈现出更好的细节拟合曲线;其在5 类植被类型中的RMSE 值分别为200.59、46.58、63.48、165.47 和40.95,在5 种方法中均为最小值,在获取高保真和高质量的EVI时间序列数据方面优势更明显有效。本文的方法研究可以给植被遥感时序数据的去噪和生态环境的研究提供有益借鉴。

关键词: 图论邻点方法, EVI, 时序数据, 植被遥感, MODIS, 曲线拟合, 去噪, 重构方法

Abstract: The MODIS Enhanced Vegetation Index (EVI) time- series data has been widely used in many research fields such as vegetation observation, ecological environment, and global meteorological changes. However, even though the EVI time series data has undergone strict preprocessing, there are still some noises in it. Therefore, this paper develops a simple and effective method to reconstruct EVI time-series data and eliminate the noise in EVI time- series data, especially some noise caused by atmospheric clouds and snow cover. The theory of the new method is derived from graph theory, using the relationship of the Laplacian matrix to assign the weight of the pixel of the selected neighborhood window in EVI to get the fitting of the center pixel. The new method has been applied to MODIS MOD13A1 products from 2016 to 2018 and compared with the S-G filtering method, Harmonic Analysis of Time Series method, Double Logistic function method, and Asymmetric Gaussian model function method. The results show that in the desert, grassland, and woodland, the absolute difference of the leave-one verification test of the new method is the smallest, which is better than other methods; when fitting EVI time-series data of different vegetation types, the graph theory neighbor method presents a better detailed fitting curve; the RMSE values of the new method in the five vegetation types are 200.59, 46.58, 63.48, 165.47, and 40.95 respectively, which are the smallest values among the five methods and are more effective in obtaining high- fidelity and high- quality EVI time- series data. The method research in this article can provide a useful reference for the denoising of vegetation remote sensing time- series data and the study of the ecological environment.

Key words: graph theory neighbor point method, EVI, time- series data, vegetation remote sensing, MODIS, curve fitting, denoising, reconstruction method

谌稳, 孙立群, 李晴岚, 陈晨, 李家叶. 一种基于图论重构MODIS EVI时间序列数据集的新方法[J]. 地球信息科学学报, , (): 1-2.

CHENWen, SUN Liqun, LI Qinglan, CHEN Chen, LI Jiaye. A New Method to Reconstruct MODIS EVI Time Series Data Set based on Graph Theory[J]. Journal of Geo-information Science, , (): 1-2.

参考文献

[1] Marcos B, Gonçalves J, Alcaraz-Segura D, et al.Improving the detection of wildfire disturbances in space and time based on indicators extracted from MODIS data: A case study in northern Portugal[J]. International Journal of Applied Earth Observation and Geoinformation, 2019, 78:77-85.
[2] 王正兴, 刘闯,Huete Alfredo.植被指数研究进展: 从AVHRR-NDVI 到MODIS-EVI[J].生态学报,2003,23(5):979-987.
[Wang Z X, Liu C, Huete Alfredo.Research progress of vegetation Index: from AVHRR-NDVI to MODIS-EVI[J]. Acta Ecologica Sinica, 2003,23(5):979-987.]
[3] 孙立双,马运涛,毕天平,等.辽宁地区不同地表覆盖类型EVI 和NDVI 特征[J].沈阳建筑大学学报(自然科学版),2013,29(6):1024-1029.
[Sun L S, Ma Y T, Bi T P, et al.EVI and NDVI characteristics of different surface cover types in Liaoning area[J]. Journal of Shenyang Construction University (Natural Science Edition), 2013,29(6):1024-1029.]
[4] 周惠慧,王楠,黄瑶,等.不同时间间隔下的遥感时间序列重构模型比较分析[J].地球信息科学学报,2016,18(10):1410-1417.
[Zhou H H, Wang N, Huang Y, et al.Comparison and analysis of remote sensing time series reconstruction models under different time intervals[J]. Journal of Geo-information Science, 2016,18(10):1410-1417.]
[5] 刘建文,周玉科.站点尺度的青藏高原时序NDVI 重构方法比较与应用[J].地理科学进展,2018,37(3):427-437.
[Liu J W, Zhou Y K.Comparison and application of sitescale NDVI reconstruction methods for Qinghai- Tibet Plateau time series[J]. Progress in Geography, 2018,37(3): 427-437.]
[6] 贾若楠,杜鑫,李强子,等.近15 年锡林郭勒盟植被变化时空特征及其对气候的响应[J].中国水土保持科学,2016,14(5):47-56.
[Jia R N, Du X, Li Q Z, et al.Spatio-temporal characteristics of vegetation change and its response to climate in Xilingol League in recent 15 years[J]. Science of Soil andWater Conservation, 2016,14(5):47-56.]
[7] 刘倩楠,岳彩荣,欧阳志云,等.基于MODIS—NDVI 时序数据的重庆市植被变化研究[J].测绘与空间地理信息,2012,35(3):99-102.
[Liu Q N, Yue C R, Ouyang Z Y, et al.Study on vegetation change in Chongqing based on MODIS-NDVI time series data[J]. Geomatics & Spatial Information Technology, 2012,35(3):99-102.]
[8] Holben B N.Characteristics of maximum-value composite images from temporal AVHRR data[J]. International journal of remote sensing, 1986,7(11):1417-1434.
[9] Viovy N, Arino O, Belward A S.The Best Index Slope Extraction (BISE): A method for reducing noise in NDVI time- series[J]. International Journal of remote sensing, 1992,13(8):1585-1590.
[10] Lovell J L, Graetz R D.Filtering pathfinder AVHRR land NDVI data for Australia[J]. International Journal of Remote Sensing, 2001,22(13):2649-2654.
[11] Chen Jin, Jönsson P, Tamura M, et al.A simple method for reconstructing a high- quality NDVI time- series data set based on the Savitzky-Golay filter[J]. Remote sensing of Environment, 2004,91(3-4):332-344.
[12] Ma MingGuo, Veroustraete F. Reconstructing pathfinder AVHRR land NDVI time-series data for the Northwest of China[J]. Advances in Space Research, 2006,37(4):835-840.
[13] Ryo Michishita, Zhenyu Jin, Jin Chen, et al.Empirical comparison of noise reduction techniques for NDVI timeseries based on a new measure[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2014,91(5):17-28.
[14] 蒋雪冰,胡月明,刘振华,等.基于线性内插的扩展卡尔曼滤波法NDVI 时间序列重构研究[J].科技通报,2017,33(2):137-142.
[Jiang X B, Hu M H, Liu Z H, et al.Research on NDVI time series reconstruction based on Extended Kalman Filter based on linear interpolation[J]. Bulletin of Science and Technology, 2017,33(2):137-142.]
[15] Roerink G J, Menenti M, Verhoef W.Reconstructing cloudfree NDVI composites using Fourier analysis of time series[J]. International Journal of Remote Sensing, 2000,21(9):1911-1917.
[16] Beck P S A, Atzberger C, Høgda K A, et al. Improved monitoring of vegetation dynamics at very high latitudes: A new method using MODIS NDVI[J]. Remote sensing of Environment, 2006,100(3):321-334.
[17] Jonsson P, Eklundh L.Seasonality extraction by function fitting to time- series of satellite sensor data[J]. IEEE transactions on Geoscience and Remote Sensing, 2002,40(8):1824-1832.
[18] Sellers P J, Randall D A, Collatz G J, et al.A revised land surface parameterization (SiB2) for atmospheric GCMs. 1. Model formulation[J]. Journal of Climate, 1996,9(4):676-705.
[19] Zhang X, Li R, Yue Y M, et al.Improved algorithm for reconstructing vegetation index image time series based on Fourier Harmonic Analysis[J]. Journal of Remote Sensing, 2010,14(3):437-447.
[20] 何月,樊高峰,张小伟,等.浙江省植被物候变化及其对气候变化的响应[J].自然资源学报,2013,28(2):220-233.
[He Y, Pan G F, Zhang X W, et al.Phenological change of vegetation in Zhejiang province and its response to climate change[J]. Journal of Natural Resources, 2013,28(2): 220-233.]
[21] 张晗,任志远.多种时序NDVI重建方法比较与应用分析[J]. 中国农业科学,2014,47(15):2998-3008.
[Zhang H, Ren Z Y.Comparison and application analysis of several temporal NDVI reconstruction methods[J]. Scientia Agricultura Sinica, 2014,47(15):2998-3008.]
[22] Zhou J X,Chen J, Chen X H, et al.Sensitivity of six typical spatiotemporal fusion methods to different influential factors: A comparative study for a normalized difference vegetation index time series reconstruction[J]. Remote Sensing of Environment, 2021,252:112-130.
[23] Zurita M R, Clevers J G, Schaepman M E.Unmixingbased Landsat TM and MERIS FR data fusion[J]. IEEE geoscience and remote sensing letters, 2008,5(3):453-457.
[24] Rao Y, Zhu X, Chen J, et al.An improved method for producing high spatial- resolution NDVI time series datasets with multi-temporal MODIS NDVI data and Landsat TM/ ETM+images[J]. Remote Sensing, 2015,7(6):7865-7891.
[25] Gao F, Hilker T, Zhu X, et al.Fusing Landsat and MODIS data for vegetation monitoring[J]. IEEE Geoscience and Remote Sensing Magazine, 2015,3(3):47-60.
[26] Wang Q, Peter MA.Spatio-temporal fusion for daily Sentinel- 2 images[J]. Remote Sensing of Environment, 2018, 204:31-42.
[27] Huang B, Song H.Spatiotemporal reflectance fusion via sparse representation[J]. IEEE Trans. Geoscience and Remote Sensing, 2012,50(10):3707-3716.
[28] Zhu X, Helmer E H, Gao F, et al.A flexible spatiotemporal method for fusing satellite images with different resolutions[J]. Remote Sensing of Environment, 2016,172: 165-177.
[29] Atzberger C, Eilers P H C. A time series for monitoring vegetation activity and phenology at 10- daily time steps covering large parts of South America[J]. International Journal of Digital Earth, 2011,4(5):365-386.
[30] Julien Y, Sobrino J A.Comparison of cloud- reconstruction methods for time series of composite NDVI data[J]. Remote Sensing of Environment, 2010,114(3):618-625.
[31] 耿丽英,马明国.长时间序列NDVI 数据重建方法比较研究进展[J].遥感技术与应用,2014,29(2):362-368.
[Geng L Y, Ma M G.Progress in comparative study of NDVI data reconstruction methods for long time series[J]. Remote Sensing Technology andApplication, 2014,29(2):362-368.]
[32] 王乾坤,于信芳,舒清态,等. MODIS EVI 时序数据重建方法及拟合分析[J].地球信息科学学报,2015,17(6):732-741.
[Wang Q K, Yu X F, Shu Q T, et al.Analysis of time- series data reconstruction method and fitting based on EVI of MODIS[J]. Journal of Geo- information Science, 2015,17(6):732-741.]
[33] Jönsson P, Eklundh L.TIMESAT: A program for analyzing time-series of satellite sensor data[J]. Computers and Geosciences, 2004,30(8):33-845.