On the Generation of Gapless and Seamless Daily Surface Reflectance Data

The land surface reflectance data are indispensable to generate many other land products. Global land surface reflectance data have been routinely produced from remote sensing sensors aboard different satellite platforms. However, the original data, especially the daily data, suffer from a large number of spatial gaps, which result from atmospheric contamination and instrument deficiencies. This seriously limits their further applications. Many composite products with less spatial gaps have been generated to solve the above problem, but they easily sacrifice their temporal resolutions of original data. Even worse, they cannot be directly implemented in realistic applications because of the noise and composite seams. This paper proposes a temporal–spatial reconstruction method (TSRM) to generate daily gapless and seamless land surface reflectance data. The TSRM integrates both temporal and spatial information for recovering different land cover types using three processing steps. First, spatial gaps are coarsely filled with multiyear weighted average (Step1). After that, all the gaps that are not filled in the first step are interpolated by using harmonic analysis of time series with true value constraint (Step2). Finally, the reconstructed results in the last step are seamlessly processed using the Poisson image editing method, and the seamless daily reflectance data set is generated (Step3). The Moderate Resolution Imaging Spectroradiometer reflectance data set (MOD09GA and MYD09GA) on two testing areas is selected to verify the performance of the proposed TSRM. Experimental results show that the TSRM has good performance with regard to maintaining the temporal and spatial integrity of the daily land surface reflectance data. Results on different testing sites also demonstrate that the TSRM preserves spectral integrity with clear seasonal trends for each spectral band.

[1]  Patrick Pérez,et al.  Poisson image editing , 2003, ACM Trans. Graph..

[2]  C. Justice,et al.  The generation of global fields of terrestrial biophysical parameters from the NDVI , 1994 .

[3]  Bo Huang,et al.  Comparison of three time-series NDVI reconstruction methods based on TIMESAT , 2012, IGARSS.

[4]  Gang Yang,et al.  A Moving Weighted Harmonic Analysis Method for Reconstructing High-Quality SPOT VEGETATION NDVI Time-Series Data , 2015, IEEE Transactions on Geoscience and Remote Sensing.

[5]  C. Tucker,et al.  Climate-Driven Increases in Global Terrestrial Net Primary Production from 1982 to 1999 , 2003, Science.

[6]  José A. Sobrino,et al.  Comparison of cloud-reconstruction methods for time series of composite NDVI data , 2010 .

[7]  P. Beck,et al.  Improved monitoring of vegetation dynamics at very high latitudes: A new method using MODIS NDVI , 2006 .

[8]  Liangpei Zhang,et al.  Cloud removal for remotely sensed images by similar pixel replacement guided with a spatio-temporal MRF model , 2014 .

[9]  J. Townshend,et al.  Spatially and temporally continuous LAI data sets based on an integrated filtering method: Examples from North America , 2008 .

[10]  Hao He,et al.  A Changing-Weight Filter Method for Reconstructing a High-Quality NDVI Time Series to Preserve the Integrity of Vegetation Phenology , 2012, IEEE Transactions on Geoscience and Remote Sensing.

[11]  Steven Platnick,et al.  Spatially complete global spectral surface albedos: value-added datasets derived from Terra MODIS land products , 2005, IEEE Transactions on Geoscience and Remote Sensing.

[12]  Massimo Menenti,et al.  Mapping vegetation-soil-climate complexes in southern Africa using temporal Fourier analysis of NOAA-AVHRR NDVI data , 2000 .

[13]  José A. Sobrino,et al.  Changes in land surface temperatures and NDVI values over Europe between 1982 and 1999 , 2006 .

[14]  Xiaoliang Lu,et al.  Removal of Noise by Wavelet Method to Generate High Quality Temporal Data of Terrestrial MODIS Products , 2007 .

[15]  Alan H. Strahler,et al.  Bidirectional NDVI and atmospherically resistant BRDF inversion for vegetation canopy , 2002, IEEE Trans. Geosci. Remote. Sens..

[16]  F. L. Gadallah,et al.  Destriping multisensor imagery with moment matching , 2000 .

[17]  R. Pachauri.  Managing the Risks of Extreme Events and Disasters , 2012 .

[18]  Farid Melgani,et al.  Missing-Area Reconstruction in Multispectral Images Under a Compressive Sensing Perspective , 2013, IEEE Transactions on Geoscience and Remote Sensing.

[19]  Liangpei Zhang,et al.  Reconstructing MODIS LST Based on Multitemporal Classification and Robust Regression , 2015, IEEE Geoscience and Remote Sensing Letters.

[20]  David P. Roy,et al.  Generation of Temporally Complete Daily Nadir MODIS Reflectance Time Series , 2010 .

[21]  David P. Roy,et al.  The impact of misregistration upon composited wide field of view satellite data and implications for change detection , 2000, IEEE Trans. Geosci. Remote. Sens..

[22]  C. O. Justice,et al.  Improvements in the global biospheric record from the Advanced Very High Resolution Radiometer (AVHRR) , 2000 .

[23]  Maosheng Zhao,et al.  Drought-Induced Reduction in Global Terrestrial Net Primary Production from 2000 Through 2009 , 2010, Science.

[24]  John R. G. Townshend,et al.  Global data sets for land applications from the Advanced Very High Resolution Radiometer: an introduction , 1994 .

[25]  José A. Sobrino,et al.  Global land surface phenology trends from GIMMS database , 2009 .

[26]  Gang Yang,et al.  Recovering Quantitative Remote Sensing Products Contaminated by Thick Clouds and Shadows Using Multitemporal Dictionary Learning , 2014, IEEE Transactions on Geoscience and Remote Sensing.

[27]  D. Roy,et al.  The MODIS Land product quality assessment approach , 2002 .

[28]  W. Verhoef,et al.  Reconstructing cloudfree NDVI composites using Fourier analysis of time series , 2000 .

[29]  Philip Lewis,et al.  Temporal Constraints on Linear BRDF Model Parameters , 2010, IEEE Transactions on Geoscience and Remote Sensing.

[30]  C. Tucker.  Red and photographic infrared linear combinations for monitoring vegetation , 1979 .

[31]  S. Bruin,et al.  Analysis of monotonic greening and browning trends from global NDVI time-series , 2011 .

[32]  David P. Roy,et al.  The Global Impact of Clouds on the Production of MODIS Bidirectional Reflectance Model-Based Composites for Terrestrial Monitoring , 2006, IEEE Geoscience and Remote Sensing Letters.

[33]  W. Verhoef,et al.  A colour composite of NOAA-AVHRR-NDVI based on time series analysis (1981-1992) , 1996 .

[34]  Chao-Hung Lin,et al.  Cloud Removal From Multitemporal Satellite Images Using Information Cloning , 2013, IEEE Transactions on Geoscience and Remote Sensing.

[35]  A. Belward,et al.  The Best Index Slope Extraction ( BISE): A method for reducing noise in NDVI time-series , 1992 .

[36]  Alan H. Strahler,et al.  The Moderate Resolution Imaging Spectroradiometer (MODIS): land remote sensing for global change research , 1998, IEEE Trans. Geosci. Remote. Sens..

[37]  Gang Yang,et al.  Missing Information Reconstruction of Remote Sensing Data: A Technical Review , 2015, IEEE Geoscience and Remote Sensing Magazine.

[38]  D. Roy,et al.  An overview of MODIS Land data processing and product status , 2002 .

[39]  Alan H. Strahler,et al.  An algorithm for the retrieval of albedo from space using semiempirical BRDF models , 2000, IEEE Trans. Geosci. Remote. Sens..

[40]  Massimo Menenti,et al.  Reconstruction of global MODIS NDVI time series: performance of harmonic analysis of time series (HANTS). , 2015 .

[41]  Josef Cihlar,et al.  Evaluation of compositing algorithms for AVHRR data over land , 1994, IEEE Trans. Geosci. Remote. Sens..

[42]  A. Savitzky,et al.  Smoothing and Differentiation of Data by Simplified Least Squares Procedures. , 1964 .

[43]  Zhongbo Su,et al.  Reconstruction of a cloud - free vegetation index time series for the Tibetan plateau , 2004 .

[44]  Tinghua Ai,et al.  A spatial and temporal reflectance fusion model considering sensor observation differences , 2013 .

[45]  Per Jönsson,et al.  Seasonality extraction by function fitting to time-series of satellite sensor data , 2002, IEEE Trans. Geosci. Remote. Sens..

[46]  Frédéric Baret,et al.  The CACAO Method for Smoothing, Gap Filling, and Characterizing Seasonal Anomalies in Satellite Time Series , 2013, IEEE Transactions on Geoscience and Remote Sensing.

[47]  Chao Zeng,et al.  Recovering missing pixels for Landsat ETM + SLC-off imagery using multi-temporal regression analysis and a regularization method , 2013 .

[48]  Christopher Justice,et al.  Towards a Generalized Approach for Correction of the BRDF Effect in MODIS Directional Reflectances , 2009, IEEE Transactions on Geoscience and Remote Sensing.

[49]  J. Townshend,et al.  Developing a spatially continuous 1 km surface albedo data set over North America from Terra MODIS products , 2007 .

[50]  F. Veroustraete,et al.  Reconstructing pathfinder AVHRR land NDVI time-series data for the Northwest of China , 2006 .

[51]  J. Faundeen,et al.  The 1 km AVHRR global land data set: first stages in implementation , 1994 .

[52]  B. Holben.  Characteristics of maximum-value composite images from temporal AVHRR data , 1986 .

[53]  H. Mooney,et al.  Modeling the Exchanges of Energy, Water, and Carbon Between Continents and the Atmosphere , 1997, Science.

[54]  John F. Hermance,et al.  Stabilizing high‐order, non‐classical harmonic analysis of NDVI data for average annual models by damping model roughness , 2007 .

[55]  Wout Verhoef,et al.  Mapping agroecological zones and time lag in vegetation growth by means of Fourier analysis of time series of NDVI images , 1993 .

[56]  Jin Chen,et al.  A simple method for reconstructing a high-quality NDVI time-series data set based on the Savitzky-Golay filter , 2004 .

[57]  Aaron Moody,et al.  Land-Surface Phenologies from AVHRR Using the Discrete Fourier Transform , 2001 .