A multi-temporal phenology based classification approach for Crop Monitoring in Kenya

The SBAM (Satellite Based Agricultural Monitoring) project, funded by the Italian Space Agency aims at: developing a validated satellite imagery based method for estimating and updating the agricultural areas in the region of Central-Africa; implementing an automated process chain capable of providing periodical agricultural land cover maps of the area of interest and, possibly, an estimate of the crop yield. The project aims at filling the gap existing in the availability of high spatial resolution maps of the agricultural areas of Kenya. A high spatial resolution land cover map of Central-Eastern Africa including Kenya was compiled in the year 2000 in the framework of the Africover project using Landsat images acquired, mostly, in 1995. We investigated the use of phenological information in supporting the use of remotely sensed images for crop classification and monitoring based on Landsat 8 and, in the near future, Sentinel 2 imagery. Phenological information on crop condition was collected using time series of NDVI (Normalized Difference Vegetation Index) based on Landsat 8 images. Kenyan countryside is mainly characterized by a high number of fragmented small and medium size farmlands that dramatically increase the difficulty in classification; 30 m spatial resolution images are not enough for a proper classification of such areas. So, a pan-sharpening FIHS (Fast Intensity Hue Saturation) technique was implemented to increase image resolution from 30 m to 15 m. Ground test sites were selected, searching for agricultural vegetated areas from which phenological information was extracted. Therefore, the classification of agricultural areas is based on crop phenology, vegetation index behaviour retrieved from a time series of satellite images and on AEZ (Agro Ecological Zones) information made available by FAO (FAO, 1996) for the area of interest. This paper presents the results of the proposed classification procedure in comparison with land cover maps produced in the past years by other projects. The results refer to the Nakuru County and they were validated using field campaigns data. It showed a satisfactory overall accuracy of 92.66 % which is a significant improvement with respect to previous land cover maps

A Review of Vegetation Phenological Metrics Extraction Using Time-Series, Multispectral Satellite Data

Vegetation dynamics and phenology play an important role in inter-annual vegetation changes in terrestrial ecosystems and are key indicators of climate-vegetation interactions, land use/land cover changes, and variation in year-to-year vegetation productivity. Satellite remote sensing data have been widely used for vegetation phenology monitoring over large geographic domains using various types of observations and methods over the past several decades. The goal of this paper is to present a detailed review of existing methods for phenology detection and emerging new techniques based on the analysis of time-series, multispectral remote sensing imagery. This paper summarizes the objective and applications of detecting general vegetation phenology stages (e.g., green onset, time or peak greenness, and growing season length) often termed “land surface phenology,” as well as more advanced methods that estimate species-specific phenological stages (e.g., silking stage of maize). Common data-processing methods, such as data smoothing, applied to prepare the time-series remote sensing observations to be applied to phenological detection methods are presented. Specific land surface phenology detection methods as well as species-specific phenology detection methods based on multispectral satellite data are then discussed. The impact of different error sources in the data on remote-sensing based phenology detection are also discussed in detail, as well as ways to reduce these uncertainties and errors. Joint analysis of multiscale observations ranging from satellite to more recent ground-based sensors is helpful for us to understand satellite-based phenology detection mechanism and extent phenology detection to regional scale in the future. Finally, emerging opportunities to further advance remote sensing of phenology is presented that includes observations from Cubesats, near-surface observations such as PhenoCams, and image data fusion techniques to improve the spatial resolution of time-series image data sets needed for phenological characterization

Classification of C3 and C4 Vegetation Types Using MODIS and ETM+ Blended High Spatio-Temporal Resolution Data

The distribution of C3 and C4 vegetation plays an important role in the global carbon cycle and climate change. Knowledge of the distribution of C3 and C4 vegetation at a high spatial resolution over local or regional scales helps us to understand their ecological functions and climate dependencies. In this study, we classified C3 and C4 vegetation at a high resolution for spatially heterogeneous landscapes. First, we generated a high spatial and temporal land surface reflectance dataset by blending MODIS (Moderate Resolution Imaging Spectroradiometer) and ETM+ (Enhanced Thematic Mapper Plus) data. The blended data exhibited a high correlation (R2 = 0.88) with the satellite derived ETM+ data. The time-series NDVI (Normalized Difference Vegetation Index) data were then generated using the blended high spatio-temporal resolution data to capture the phenological differences between the C3 and C4 vegetation. The time-series NDVI revealed that the C3 vegetation turns green earlier in spring than the C4 vegetation, and senesces later in autumn than the C4 vegetation. C4 vegetation has a higher NDVI value than the C3 vegetation during summer time. Based on the distinguished characteristics, the time-series NDVI was used to extract the C3 and C4 classification features. Five features were selected from the 18 classification features according to the ground investigation data, and subsequently used for the C3 and C4 classification. The overall accuracy of the C3 and C4 vegetation classification was 85.75% with a kappa of 0.725 in our study area

Using Landsat 8 image time series for crop mapping in a region of Cerrado, Brazil.

Abstract: The objective of this research is to classify agricultural land use in a region of the Cerrado (Brazilian Savanna) biome using a time series of Enhanced Vegetation Index (EVI) from Landsat 8 OLI. Phenological metrics extracted from EVI time series, a Random Forest algorithm and data mining techniques are used in the process of classification. The area of study is a region in the Cerrado in a region of the municipality of Casa Branca, São Paulo state, Brazil. The results are encouraging and demonstrate the potential of phenological parameters obtained from time series of OLI vegetation indices for agricultural land use classification

Reconstruction of Daily 30 m Data from HJ CCD, GF-1 WFV, Landsat, and MODIS Data for Crop Monitoring

With the recent launch of new satellites and the developments of spatiotemporal data fusion methods, we are entering an era of high spatiotemporal resolution remote-sensing analysis. This study proposed a method to reconstruct daily 30 m remote-sensing data for monitoring crop types and phenology in two study areas located in Xinjiang Province, China. First, the Spatial and Temporal Data Fusion Approach (STDFA) was used to reconstruct the time series high spatiotemporal resolution data from the Huanjing satellite charge coupled device (HJ CCD), Gaofen satellite no. 1 wide field-of-view camera (GF-1 WFV), Landsat, and Moderate Resolution Imaging Spectroradiometer (MODIS) data. Then, the reconstructed time series were applied to extract crop phenology using a Hybrid Piecewise Logistic Model (HPLM). In addition, the onset date of greenness increase (OGI) and greenness decrease (OGD) were also calculated using the simulated phenology. Finally, crop types were mapped using the phenology information. The results show that the reconstructed high spatiotemporal data had a high quality with a proportion of good observations (PGQ) higher than 0.95 and the HPLM approach can simulate time series Normalized Different Vegetation Index (NDVI) very well with R2 ranging from 0.635 to 0.952 in Luntai and 0.719 to 0.991 in Bole, respectively. The reconstructed high spatiotemporal data were able to extract crop phenology in single crop fields, which provided a very detailed pattern relative to that from time series MODIS data. Moreover, the crop types can be classified using the reconstructed time series high spatiotemporal data with overall accuracy equal to 0.91 in Luntai and 0.95 in Bole, which is 0.028 and 0.046 higher than those obtained by using multi-temporal Landsat NDVI data

Crop monitoring and yield estimation using polarimetric SAR and optical satellite data in southwestern Ontario

Optical satellite data have been proven as an efficient source to extract crop information and monitor crop growth conditions over large areas. In local- to subfield-scale crop monitoring studies, both high spatial resolution and high temporal resolution of the image data are important. However, the acquisition of optical data is limited by the constant contamination of clouds in cloudy areas. This thesis explores the potential of polarimetric Synthetic Aperture Radar (SAR) satellite data and the spatio-temporal data fusion approach in crop monitoring and yield estimation applications in southwestern Ontario. Firstly, the sensitivity of 16 parameters derived from C-band Radarsat-2 polarimetric SAR data to crop height and fractional vegetation cover (FVC) was investigated. The results show that the SAR backscatters are affected by many factors unrelated to the crop canopy such as the incidence angle and the soil background and the degree of sensitivity varies with the crop types, growing stages, and the polarimetric SAR parameters. Secondly, the Minimum Noise Fraction (MNF) transformation, for the first time, was applied to multitemporal Radarsat-2 polarimetric SAR data in cropland area mapping based on the random forest classifier. An overall classification accuracy of 95.89% was achieved using the MNF transformation of the multi-temporal coherency matrix acquired from July to November. Then, a spatio-temporal data fusion method was developed to generate Normalized Difference Vegetation Index (NDVI) time series with both high spatial and high temporal resolution in heterogeneous regions using Landsat and MODIS imagery. The proposed method outperforms two other widely used methods. Finally, an improved crop phenology detection method was proposed, and the phenology information was then forced into the Simple Algorithm for Yield Estimation (SAFY) model to estimate crop biomass and yield. Compared with the SAFY model without forcing the remotely sensed phenology and a simple light use efficiency (LUE) model, the SAFY incorporating the remotely sensed phenology can improve the accuracy of biomass estimation by about 4% in relative Root Mean Square Error (RRMSE). The studies in this thesis improve the ability to monitor crop growth status and production at subfield scale

Temporal optimisation of image acquisition for land cover classification with random forest and MODIS time-series

The analysis and classification of land cover is one of the principal applications in terrestrial remote sensing. Due to the seasonal variability of different vegetation types and land surface characteristics, the ability to discriminate land cover types changes over time. Multi-temporal classification can help to improve the classification accuracies, but different constraints, such as financial restrictions or atmospheric conditions, may impede their application. The optimisation of image acquisition timing and frequencies can help to increase the effectiveness of the classification process. For this purpose, the Feature Importance (FI) measure of the state-of-the art machine learning method Random Forest was used to determine the optimal image acquisition periods for a general (Grassland, Forest, Water, Settlement, Peatland) and Grassland specific (Improved Grassland, Semi-Improved Grassland) land cover classification in central Ireland based on a 9-year time-series of MODIS Terra 16 day composite data (MOD13Q1). Feature Importances for each acquisition period of the Enhanced Vegetation Index (EVI) and Normalised Difference Vegetation Index (NDVI) were calculated for both classification scenarios. In the general land cover classification, the months December and January showed the highest, and July and August the lowest separability for both VIs over the entire nine-year period. This temporal separability was reflected in the classification accuracies, where the optimal choice of image dates outperformed the worst image date by 13% using NDVI and 5% using EVI on a mono-temporal analysis. With the addition of the next best image periods to the data input the classification accuracies converged quickly to their limit at around 8–10 images. The binary classification schemes, using two classes only, showed a stronger seasonal dependency with a higher intra-annual, but lower inter-annual variation. Nonetheless anomalous weather conditions, such as the cold winter of 2009/2010 can alter the temporal separability pattern significantly. Due to the extensive use of the NDVI for land cover discrimination, the findings of this study should be transferrable to data from other optical sensors with a higher spatial resolution. However, the high impact of outliers from the general climatic pattern highlights the limitation of spatial transferability to locations with different climatic and land cover conditions. The use of high-temporal, moderate resolution data such as MODIS in conjunction with machine-learning techniques proved to be a good base for the prediction of image acquisition timing for optimal land cover classification results

Mapping of peanut crops in Queensland, Australia using time-series PROBA-V 100-m normalized difference vegetation index imagery

Mapping of peanut crops is essential in supporting peanut production, yield prediction, and commodity forecasting. While ground-based surveys can be used over small areas, the development of remote-sensing technologies could provide rapid and inexpensive crop area estimates with high accuracy over large regions. Some of these recent earth observation satellite systems, such as the Project for On-Board Autonomy Vegetation (PROBA-V), have the advantage of increased spatial and temporal resolution. With a study area located in the South Burnett region, Queensland, Australia, the primary aim of this study was to assess the ability of time-series PROBA-V 100-m normalized difference vegetation index (NDVI) for peanut crop mapping. Two datasets, i.e., PROBA-V NDVI time-series imagery and the corresponding phenological parameters generated from TIMESAT data analysis technique, were classified using maximum likelihood classification, spectral angle mapper, and minimum distance classification algorithms. The results show that among all methods used, the application of MLC in PROBA-V NDVI time series produced very good overall accuracy, i.e., 92.75%, with producer and user accuracy of each class ≥78.79  %  . For all algorithms tested, the mapping of peanut cropping areas produced satisfactory classification results, i.e., 75.95% to 100%. Our study confirmed that the use of finer resolution 100 m of PROBA-V imagery (i.e., relative to MODIS 250-m data) has contributed to the success of mapping peanut and other crops in the study area

Assessment of Time-Series MODIS Data for Cropland Mapping in the U.S. Central Great Plains

The goal of this study was to further investigate the potential of MODIS NDVI 250-m data for crop spectral characterization, discrimination, and mapping in the Great Plains of the USA using various exploratory approaches. GIS operations, and reference data refinement using clustering and visual assessment of each crop's NDVI cluster profiles in Nebraska, demonstrated that it is possible to devise an alternative reference data set and refinement plan that redresses the unexpected loss of training and validation data. A pixel-level analysis of the time-series MODIS 250-m NDVI for 1,288 field sites representing each of the eight cover types under investigation across Nebraska found that each crop type had a distinctive MODIS 250-m NDVI profile corresponding to the crop calendar. A visual and statistical comparison of the average NDVI profiles showed that the crop types were separable at different times of the growing season based on their phenology-driven spectral-temporal differences. In Kansas, an initial investigation revealed that there was near-complete agreement between the winter wheat crop profiles but that there were some minor differences in the crop profiles for alfalfa and summer crops between 2001 and 2005. However, the profiles of summer crops - corn, grain sorghum, and soybeans - displayed a shift to the right by at least 1 composite date, indicative of possible late crop planting and emergence. Alfalfa and summer crops, seem to suggest that time series NDVI response curves for crops over a growing period for one year of valid ground reference data may not be used to map crops for a different year without taking into account the climatic and/or environmental conditions of each year