Using Enhanced Gap-Filling and Whittaker Smoothing to Reconstruct High Spatiotemporal Resolution NDVI Time Series Based on Landsat 8, Sentinel-2, and MODIS Imagery

Normalized difference vegetation index (NDVI) time series data, derived from optical images, play a crucial role for crop mapping and growth monitoring. Nevertheless, optical images frequently exhibit spatial and temporal discontinuities due to cloudy and rainy weather conditions. Existing algorithms for reconstructing NDVI time series using multi-source remote sensing data still face several challenges. In this study, we proposed a novel method, an enhanced gap-filling and Whittaker smoothing (EGF-WS), to reconstruct NDVI time series (EGF-NDVI) using Google Earth Engine. In EGF-WS, NDVI calculated from MODIS, Landsat-8, and Sentinel-2 satellites were combined to generate high-resolution and continuous NDVI time series data. The MODIS NDVI was employed as reference data to fill missing pixels in the Sentinel-Landsat NDVI (SL-NDVI) using the gap-filling method. Subsequently, the filled NDVI was smoothed using a Whittaker smoothing filter to reduce residual noise in the SL-NDVI time series. With reference to the all-round performance assessment (APA) metrics, the performance of EGF-WS was compared with the conventional gap-filling and Savitzky-Golay filter approach (GF-SG) in Fusui County of Guangxi Zhuang Autonomous Region. The experimental results have demonstrated that the EGF-WS can capture more accurate spatial details compared with GF-SG. Moreover, EGF-NDVI of Fusui County exhibited a low root mean square error (RMSE) and a high coefficient of determination (R-2). In conclusion, EGF-WS holds significant promise in providing NDVI time series images with a spatial resolution of 10 m and a temporal resolution of 8 days, thereby benefiting crop mapping, land use change monitoring, and various ecosystems, among other applications.