Monitoring of rice cropping intensity in the upper Mekong Delta,Vietnam using time-series MODIS data

Information on rice growing areas is important for policymakers to devise agricultural plans. This research explores the monitoring of rice cropping intensity in the upper Mekong Delta, Vietnam (from 2001 to 2007) using time-series MODIS NDVI 250-m data. Data processing includes three steps: (1) noise is filtered from the time-series NDVI data using empirical mode decomposition (EMD); (2) endmembers are extracted from the filtered time-series data and trained in a linear mixture model (LMM) for classification of rice cropping systems; and (3) classification results are verified by comparing them with the ground-truth and statistical data. The results indicate that EMD is a good filter for noise removal from the time-series data. The classification results confirm the validity of LMM, giving an overall accuracy of 90.1% and a Kappa coefficient of 0.7. The lowest producer and user accuracies were associated with single crop rain-fed rice class due to the mixed pixel problems. A strong yearly correlation at the district level was revealed in the MODIS-derived areas (R² ? 0.9). Investigation of interannual changes in rice cropping intensity from 2001 to 2007 showed a remarkable conversion from double to triple crop irrigated rice from 2001 to 2003, especially in the Thoai Son and Phu Tan districts. A big conversion from triple crop rice back to double crop rice cultivation was also observed in Phu Tan from 2005 to 2006. These changes were verified by visual interpretation of Landsat images and examination of NDVI profiles.     

Study on extraction of crop information using time-series MODIS data in the Chao Phraya Basin of Thailand

In order to acquire the crop-related information in Chao Phraya Basin, time-series MODIS data were used in this paper. Although the spatial resolution of MODIS data is not very high, it is still useful for detecting very large-scale phenomenon, such as changes in seasonal vegetation patterns. After the data processing a general crop-related LULC (land use and land cover) map, cropping intensity map and cropping patterns map were produced. Analysis of these maps showed that the main land use type in the study area was farmland, most of which was dominated by rice. Rice fields mostly concentrated in the flood plains and double or triple rice-cropping system was commonly employed in this area. Maize, cassava, sugarcane and other upland crops were mainly distributed in the high alluvial terraces. Because these area often have water shortage problem particularly in the dry season which can support only one crop in a year, the cropping intensity was very low. However, some upland areas can be cultivated twice a year with crops which have short growing seasons. The crop information extracted from MODIS data sets were assessed by CBERS data, statistic data and so on. It was shown that MODIS derived crop information coincided well with the statistic data at the provincial level. At the same time, crop information extracted by MODIS data sets and CBERS were compared with each other which also showed similar spatial patterns.     

Object-based rice mapping using time-series and phenological data

Remote sensing techniques are often used in mapping rice, but high quality time-series remote sensing data are difficult to obtain due to the cloudy weather of rice growing areas and long satellite revisit interval. As such, rice mapping is usually based on mono-temporal Landsat TM/ETM+ data, which have large uncertainties due to the spectral similarity of different vegetation types. Moreover, conventional pixel-based classification method is unable to meet the required accuracy for rice mapping. Therefore, this study proposes a new strategy for mapping rice in cloud-prone areas using fused data of Landsat-8 OLI time-series and phenological parameters, based on the object-based method. We determine the critical growth stages of paddy rice from observed phenological data and MODIS-NDVI time-series data. The spatial and temporal adaptive reflectance fusion model (STARFM) is used to blend the MODIS and Landsat data to obtain a multi-temporal Landsat-like dataset for classification. Finally, an object-oriented algorithm is used to extract rice paddies from the Landsat-like, time-series dataset. The validation experiments show that the proposed method can provide high accuracy rice mapping, with an overall accuracy of 92.38% and a kappa coefficient of 0.85.     

A study of forest vegetation dynamics in the south of the Krasnoyarskii Krai in spring

Remote sensing applications have greatly enhanced ability to monitor and manage in the areas of forestry. Accurate measurements of regional and global scale vegetation dynamics (phenology) are required to improve models and understanding of inter-annual variability in terrestrial ecosystem carbon exchange and climate–biosphere interactions. Study of vegetation phenology is required for understanding of variability in ecosystem. In this paper, monitoring of vegetation dynamics using time series of satellite data is presented. Vegetation variability (vegetation rate) in different topoclimatic areas is investigated. Original software using IDL interactive language for processing of satellite long-term data series was developed. To investigate growth dynamics vegetation rate inferred from remote sensing was used. All estimations based on annual time series of Moderate Resolution Imaging Spectroradiometer (MODIS) imagery. Vegetation rate for Enhanced Vegetation Index (EVI) and Normalized Difference Vegetation Index (NDVI) was calculated using MODIS data. The time series covers spring seasons of each of 9 years, from 2000 to 2008. Comparison of EVI and NDVI derived growth rates has shown that NDVI derived rates reveal spatial structure better. Using long-term data of vegetation rates variance was estimated that helps to reveal areas with anomalous growth rate. Such estimation shows sensitivity degree of different areas to different topoclimatic conditions. Woods of heights depend on spatial topoclimatic variability unlike woods of lowlands. Principal components analysis shows vegetation with different rate conditions. Also it reveals vegetation of same type in areas with different conditions. It was demonstrated that using of methods for estimating the dynamic state of vegetation based on remote sensing data enables successful monitoring of vegetation phenology.     

Spatio-temporal dynamics in the phenology of croplands across the Indo-Gangetic Plains

Spatio-temporal dynamics in land surface phenology parameters observed over croplands can inform on crop-climate interactions and, elucidate local to regional scale vulnerabilities either due to climate change or prevailing sub-optimal agricultural practices. Here, we observe spatio-temporal trends in land surface phenology parameters (cropping intensity, length of growing season and productivity) for kharif and rabi cropping seasons from satellite data across the Indo-Gangetic Plains from 1982 to 2006. The productivity of the Indo-Gangetic Plains croplands is of regional importance and is a vital component of Indian national food security efforts. Aside from local and intra-state heterogeneity in observed trends there was a clear west-to-east gradient in cropping intensity. Key observed trends include increasing cropping intensity in the eastern IGP, increasing number of growing days per year in Bihar, Uttar Pradesh and Haryana and increasing productivity in both cropping seasons across the IGP. This information is a crucial input to integrated assessments of the croplands to ensure management of the agricultural system shifts towards a trajectory of climate-resilience and environmental sustainability.     

Monitoring of long-term land subsidence from 2003 to 2017 in coastal area of Semarang,Indonesia by SBAS DInSAR analyses using Envisat-ASAR,ALOS-PALSAR,and Sentinel-1A SAR data

Land subsidence is a critical issue that large cities located in coastal areas, such as Semarang, Indonesia, must address. The monitoring of land subsidence is vital for predicting and mitigating the disasters that such subsidence may cause. Therefore, an economical and effective monitoring method, which can continuously provide accurate measurements over extensive areas, is highly required. Differential Interferometry Synthetic Aperture Radar (DInSAR) has the potential to be a powerful technique that can meet the above demands. Actually, DInSAR has been applied to monitor the subsidence in Semarang, but it was for a limited period before 2012.In order to clarify the transition of the long-term subsidence behavior in Semarang, the Small Baseline Subset (SBAS) method, which is one type of time-series DInSAR, is employed in this research. The sets of data of Envisat-ASAR (2003–2007), ALOS-PALSAR (2007–2011), and Sentinel-1A (2015–2017) are employed for the analyses. Then, the validity of the SBAS results is discussed from the viewpoints of both spatial distribution and temporal transition using GPS displacement measurement results and the geological conditions of the ground.On the other hand, as the lifespan of SAR satellites is commonly designed to be around 5–7?years, an appropriate method, which can connect the subsidence provided independently by the unlinked time-series data sets of the three different SAR satellite data, is required. This study uses the Hyperbolic Method (HM) to connect the above unlinked SBAS results. The HM is often used to fit the monitored subsidence in practice as a geotechnical engineering tool. Using this method, 14?years of the temporal behavior of the subsidence in Semarang is evaluated.It is found that the transition of the subsidence is different depending on the location, and that the subsidence rate is still increasing in the north and northeast parts of the coastal area. This study shows that SBAS DInSAR can be a useful tool for long-term continuous subsidence monitoring.     

A comparative analysis of SLR,MLR, ANN,XGBoost and CNN for crop height estimation of sunflower using Sentinel-1 and Sentinel-2

Sustainable monitoring and determining the biophysical characteristics of crops is of global importance due to the increase in demand for food. In this context, remote sensing data provide valuable information on crops. This study investigates the relationship between the variables determined from both Synthetic Aperture Radar (SAR) and optical images and crop height. For this purpose, backscatter (σVH, σVV, σVH / σVV) and coherence (?VH, ?VV) of multi-temporal dual-polarized Sentinel-1 and vegetation indices of multi-temporal Sentinel-2 data are analyzed. Two indices, namely, Normalized Difference Vegetation Index (NDVI) and NDVI with the red-edge band (NDVIred), are interpreted to identify the contribution of the red-edge band over the near-infrared band. The Zile District of Tokat province in Turkey where dominantly sunflower cultivation is carried out, was selected as the study area. In the analysis of the data, Simple Linear Regression (SLR), Multiple Linear Regression (MLR), Artificial Neural Network (ANN), EXtreme Gradient Boosting (XGBoost), and Convolutional Neural Network (CNN) were used. In the results of the study, ANN showed the lowest RMSE = 3.083 cm (RMSE%= 11.284) in the stem elongation period. The CNN followed the lowest RMSE for the Inflorescence development and flowering stages 19.223 cm (RMSE%=15.458) and 8.731 cm (RMSE%=5.821), respectively. In the ripening period, XGBoost achieved the lowest RMSE = 8.731 cm (RMSE%=6.091). All the best models in four methods were created using common variables of σVH, σVV, ?VH, ?VV and NDVIred, except ANN which exclude coherence variables. The results concluded that NDVIred contributed more than NDVI which is widely interpreted in previous studies.     

Mapping recent agricultural developments in China from satellite data

Landsat data have been employed to study and map agricultural developments in three regions of China: 1) Pearl River delta; 2) Nen River basin; and 3) Xinjiang Autonomous Region. Manual interpretation procedures used in conjunction with multi-date Landsat images and collateral information permitted rice yields to be estimated for the Pearl River delta in 1978. A combination of manual and computer-assisted analyses of Landsat data of Northeast China revealed that more than 15,000 km² of agricultural land in a 184,500 km² study area had been reclaimed from rangeland and marshland. These analyses also indicated a shift in cropping practices, with the foodcrops wheat and corn replacing cash crops such as soybeans. In the arid west, Landsat image data provided valuable input to a geographic information system (GIS). It appears the GIS approach will prove useful for evaluating agricultural land potential in the remote areas of China.     

Impact of rain on wave retrieval from Sentinel-1 synthetic aperture radar images in tropical cyclones

The main objective of our work was to investigate the impact of rain on wave observations from C-band (~5.3 GHz) synthetic aperture radar (SAR) in tropical cyclones. In this study, 10 Sentinel-1 SAR images were available from the Satellite Hurricane Observation Campaign, which were taken under cyclonic conditions during the 2016 hurricane season. The third-generation wave model, known as Simulating WAves Nearshore (SWAN) (version 41.31), was used to simulate the wave fields corresponding to these Sentinel-1 SAR images. In addition, rainfall data from the Tropical Rainfall Measuring Mission satellite passing over the spatial coverage of the Sentinel-1 SAR images were collected. The simulated results were validated against significant wave heights (SWHs) from the Jason-2 altimeter and European Centre for Medium-Range Weather Forecasts data, revealing a root mean square error (RMSE) of ~0.5 m with a 0.25 scatter index. Winds retrieved from the VH-polarized Sentinel-1 SAR images using the Sentinel-1 Extra Wide-swath Mode Wind Speed Retrieval Model after Noise Removal were taken as prior information for wave retrieval. It was discovered that rain did indeed affect the SAR wave retrieval, as evidenced by the 3.21-m RMSE of SWHs between the SAR images and the SWAN model, which was obtained for the ~1000 match-ups with raindrops. The raindrops dampened the wave retrieval when the rain rate was < ~5 mm/hr; however, they enhanced wave retrieval for higher rain rates. It was also found that the portion of the rain-induced ring wave with a wave number > 0.05 rad/m (~125 m wavelength) was clearly observed in the SAR-derived wave spectra.     

Estimating canopy LAI and chlorophyll of tropical forest plantation (North India) using Sentinel-2 data

With the free and full access to images from Sentinel-2 satellite, the interest to use this data for quantitative retrieval of vegetation parameters is ever-increasing. LAI and chlorophyll are two key variables which are desired for studying productivity, nutrient and stress status of vegetation. Studies carried out on croplands using simulated Sentinel-2 MSI and parametric approach have identified vegetation indices (VIs) with high sensitivity to LAI and chlorophyll. To test how Sentinel-2 red-edge based VIs perform for retrieval of LAI and Chlorophyll of tropical mixed forest canopies, this study has been performed. The field measurements of LAI and chlorophyll content were recorded in a total of 28 ESUs (Elementary Sampling Units) in Bhakra range in the Tarai Central Forest Division, Uttarakhand (India). The in-situ measurements were statistically correlated with Sentinel-2VIs and strength of correlation was validated using Predicted Residual Error Sum of Squares (PRESS) statistic. Field LAI corrected for foliage clumpiness effect improved correlation of VIs with LAI. Among all VIs tested, Normalized Difference Index (NDI) offered highest positive correlation (R² = 0.79, p < 0.05) with LAI while Red-Edge Chlorophyll Index (RECI) (R² = 0.83, RMSE = 0.24 g/m², p < 0.05) and Simple Ratio (SR) 740/705 (R² = 0.79, RMSE = 0.27 g/m², p < 0.05) were the most closely related to chlorophyll content. VIs with red-edge and NIR combinations offered best results.     

Multi-seasonal evaluation of hybrid poplar (P. Deltoides) plantations using Worldview-3 imagery and State-Of-The-Art ensemble learning algorithms

Forest resources are the primary components of the ecosystem environment. Poplars (Populus sp.), a member of the fast-growing trees, are one of the most productive forest tree species for industrial production thanks to their desirable traits comprising rapid growth, hybridization ability, and ease of propagation. Determining poplar cultivated areas and mapping their geographical distributions is critical for planners and decision-makers to increase the ecological and economic benefits of poplar trees. Due to the biodiversity of each geographical region and seasonal vegetation variations, classification based on remotely sensed imagery is essential for cropland monitoring. The main goal of this study is to evaluate the potential of high-resolution multi-temporal (growing season and end of the growing season) Worldview-3 imagery in mapping poplar plantations in the Akyaz? district of Sakarya, Turkey. For this purpose, pixel- and object-based image analysis with up-to-date ensemble learning algorithms, namely random forest (RF), categorical boosting (CB), and extreme gradient boosting tree (XGB), were employed for mapping poplar fields. Results indicated that the object-based classification approach provided statistically significant improvements in map-level (about 4%) and class-level accuracy (e.g., approximately 7% and %2 for poplar and young poplar classes, respectively) than pixel-based classification. While the CB performed superior classification performance for the object-based approach (92.56%), the highest classification performance was obtained with the XGB algorithm for the pixel-based approach (90.42%) for the end of the growing season data. McNemar’s statistical test also confirmed that the performances of CB and XGB algorithms were statistically similar in pixel-based classification. Finally, analysis of multi-season images revealed that sensitivity of the vegetation phenology and seasonal effects considerably affect the separability of poplar tree species.     

Validation of non-linear split window algorithm for land surface temperature estimation using Sentinel-3 satellite imagery: Case study; Tehran Province,Iran

In recent years, land surface temperature (LST) has become critical in environmental studies and earth science. Remote sensing technology enables spatiotemporal monitoring of this parameter on large scales. This parameter can be estimated by satellite images with at least one thermal band. Sentinel-3 SLSTR data provide LST products with a spatial resolution of 1 km. In this research, direct and indirect validation procedures were employed to evaluate the Sentinel-3 SLSTR LST products over the study area in different seasons from 2018 to 2019. The validation method was based on the absolute (direct) evaluation of this product with field data and comparison (indirect) evaluation with the MODIS LST product and the estimated LST using the non-linear split-window (NSW) algorithm. Also, two emissivity estimation methods, (1) NDVI thresholding method (NDVI-THM) and (2) classification-based emissivity method (CBEM), were used to estimate the LST using the NSW method according to the two thermal bands of Sentinel-3 images. Then, the accuracy of these methods in estimating LST was evaluated using field data and temporal changes of vegetation, which the NDVI-THM method generated better results. For indirect evaluation between the Sentinel-3 LST product, MODIS LST product, and LST estimated using NSW, four filters based on spatial and temporal separates between pairs of pixels and pixel quality were used to ensure the accuracy and consistency of the compared pairs of a pixel. In general, the accuracy results of the LST products of MODIS and Sentinel-3, and LST estimated using NSW showed a similar trend for LST changes during the seasons. With respect to the two absolute and comparative validations for the Sentinel-3 LST products, summer with the highest values of bias (?1.24 K), standard deviation (StDv = 2.66 K), and RMSE (2.43 K), and winter with the lowest ones (bias of 0.14 K, StDv of 1.13 K, and RMSE of 1.12 K) provided the worst and best results for the seasons in the period of 2018–2019, respectively. According to both absolute and comparative evaluation results, the Sentinel-3 SLSTR LST products provided reliable results for all seasons on a large temporal and spatial scale over our studied area.     

Canadian Landsat studies for monitoring agricultural intensification and urbanization: A summary

This report summarizes a major paper reviewing the progression of work using Landsat for land use and agricultural monitoring in Canada since 1973. From the launch of Landsat to the present, the focus has moved from using simple visual techniques for interpretation of photographic products to the use of conceptually simple methods which now employ the power of special-purpose image analysis hardware and the standard geometrically corrected products available in Canada. Outlined in the major paper are studies on urban growth, large area land use mapping, crop studies (related to both crop area estimates and erosion potential), assessment of changes in vegetation vigour and clearing of land in areas which were devoted to both dry land farming and forestry. This report emphasizes developments in land use studies.     

Land use/land cover classification using time series Landsat 8 images in a heavily urbanized area

It is of great significance to timely, accurately, and effectively monitor land use/cover in city regions for the reasonable development and utilization of urban land resources. The remotely sensed dynamic monitoring of Land use/land cover (LULC) in rapidly developing city regions has increasingly depended on remote-sensing data at high temporal and spatial resolutions. However, due to the influence of revisiting periods and weather, it is difficult to acquire enough time-series images with high quality at both high temporal and spatial resolution from the same sensor. In this paper we used the temporal-spatial fusion model ESTARFM (Enhanced Spatial and Temporal Adaptive Reflectance Fusion Model) to blend Landsat8 and MODIS data and obtain time-series Landsat8 images. Then, land cover information is extracted using an object-based classification method. In this study, the proposed method is validated by a case study of the Changsha City. The results show that the overall accuracy and Kappa coefficient were 94.38% and 0.88, respectively, and the user/producer accuracies of vegetation types were all over 85%. Our approach provides an accurate and efficient technical method for the effective extraction of land use/cover information in the highly heterogeneous regions.     

POBI interpolation algorithm for CYGNSS near real time flood detection research: A case study of extreme precipitation events in Henan,China in 2021

In late 2016, NASA launched the first constellation of the global navigation satellite system reflectometry (GNSS-R) small satellites called the Cyclone Global Navigation Satellite System (CYGNSS). The stable data quality and continuous free availability of CYGNSS scientific data provided a new method for flood monitoring. However, owing to the pseudorandom distribution of CYGNSS data, researchers must always choose between high temporal resolution and high spatial resolution during the performance of flood monitoring based on CYGNSS data. For floods caused by extreme precipitation with sudden and short durations, the current flood mapping based on CYGNSS data cannot be updated in near real time. However, the near real time update of the flood distribution range is meaningful for postdisaster emergency response and rapid rescue. This study aimed to address this problem using a newly proposed spatial interpolation method based on previously observed behaviour (POBI). First, a method for calculating the surface reflectivity of the CYGNSS was introduced, followed by the principle of the POBI spatial interpolation method. The applicability of the POBI method in Henan Province, China, was then analysed, and by using the flood in Henan Province, China, in July 2021 as an example, the feasibility of CYGNSS near real time flood mapping based on the POBI method was evaluated. Based on the results, near real time and 3 km flood distribution monitoring results can be obtained using the proposed new method. The results were evaluated using MODIS (Moderate Resolution Imaging Spectroradiometer) images and compared with the observations of SMAP (Soil Moisture Active Passive) and GPM (Global Precipitation Measurement) in the same period. The results show that the flooded areas obtained by CYGNSS correspond to the inundated areas in MODIS images and are also in high agreement with the SMAP. In addition, CYGNSS allows for finer mapping and quantification of inundation areas and flood duration. Moreover, we also discussed the potential of CYGNSS to detect floods in shorter periods of time (a few hours) and did a preliminary evaluation using precipitation data from meteorological stations. The results are also highly consistent.     

Validity assessment of SAMOSA retracking for fully-focused SAR altimeter waveforms

We demonstrate in this work how we can take advantage of known unfocused SAR (UF-SAR) retracking methods (e.g. the physical SAMOSA model) for retracking of fully-focused SAR (FF-SAR) waveforms. Our insights are an important step towards consistent observations of sea surface height, significant wave height and backscatter coefficient (wind speed) with both UF-SAR and FF-SAR. This is of particular interest for SAR altimetry in the coastal zone, since coastal clutter may be filtered out more efficiently in the high-resolution FF-SAR waveform data, which has the potential to improve data quality. We implemented a multi-mission FF-SAR altimetry processor for Sentinel-3 (S3) and Sentinel-6 Michael Freilich (S6), using a back-projection algorithm, and analysed ocean waveform statistics compared to multilooked UF-SAR. We find for Sentinel-3 that the averaged power waveforms of UF-SAR and FF-SAR over ocean are virtually identical, while for Sentinel-6 the FF-SAR power waveforms better resemble the UF-SAR zero-Doppler beam. We can explain and model the similarities and differences in the data via theoretical considerations of the waveform integrals. These findings suggest to use the existing UF-SAR SAMOSA model for retracking S3 FF-SAR waveforms but the SAMOSA zero-Doppler beam model for S6 FF-SAR waveforms, instead. Testing the outlined approach over short track segments, we obtain range biases between UF-SAR and FF-SAR lower than 2 mm and significant wave height biases lower than 5 cm.     

基于DTW的HJ-1B卫星地表温度产品质量检验

HJ-1B卫星热红外数据应用广泛,但其地表温度反演产品的质量检验工作尚不完善。以黑河流域为研究区,利用普适性单通道算法得到HJ 1B地表温度,基于7个地面站点（下垫面为荒漠、沙漠、植被、农作物、雪地和湿地）同步观测资料和MODIS地表温度产品（MOD11A1）,引入动态时间规整方法（DTW）对站点处HJ 1B地表温度进行验证。试验结果表明：HJ 1B地表温度反演产品与地面观测值的偏差值在沙漠和荒漠站点为1K以内,均方根误差在05K左右;对于植被和农作物站点的偏差在2K以内,均方根误差为1~2K;基于DTW的验证对时序不匹配的数据评价结果与现有指标表现一致。HJ-1B地表温度反演产品与地面站点的相对偏差均低于其与MODIS地表温度反演产品的相对偏差。     

Detection of surface temperature anomaly of the Sea of Marmara

Monitoring sea surface temperature (SST) over a long-term and detecting the anomalies highly contribute to understanding the prevailing water quality of the sea. Earth observation satellite images are the key data sources that offer the long-term SST detection in a cost and time effective way. Since the Sea of Marmara in Türkiye is surrounded by the highly populated provinces, the water quality of the sea has gained importance for scientific and public communities over the years. This article emphasizes on the significance of detecting SST trend and corresponding anomalies of the Sea of Marmara over the past 32 years. To address the SST variations of the Sea of Marmara in time, a comprehensive set of both field and satellite data regarding SSTs were obtained within the context of this study. The SST trend and its anomalies between the years 1990 and 2021 were detected by applying Seasonal-Trend decomposition procedure based on LOESS (STL) method to NOAA OISST V2 data. On the other hand, spatial SST distribution was detected with Landsat-8, Sentinel-3 and NOAA OISST V2 satellite data. SST results were verified with the in-situ data within the scope of accuracy assessment. The results showed that SST time-series data performed an increasing trend and had anomalies mostly during the spring months in the recent years.     

A simplified data assimilation method for reconstructing time-series MODIS NDVI data

The Normalized Difference Vegetation Index (NDVI) is an important vegetation index, widely applied in research on global environmental and climatic change. However, noise induced by cloud contamination and atmospheric variability impedes the analysis and application of NDVI data. In this work, a simplified data assimilation method is proposed to reconstruct high-quality time-series MODIS NDVI data. We extracted 16-Day L3 Global 1 km SIN Grid NDVI data sets for western China from MODIS vegetation index (VI) products (MOD13A2) for the period 2003–2006. NDVI data in the first three years (2003–2005) were used to generate the background field of NDVI based on a simple three-point smoothing technique, which captures annual features of vegetation change. NDVI data for 2006 were used to test our method. For every time step, the quality assurance (QA) flags of the MODIS VI products were adopted to empirically determine the weight between the background field and NDVI observations. Ultimately, more reliable NDVI data can be produced. The results indicate that the newly developed method is robust and effective in reconstructing high-quality MODIS NDVI time-series.