Estimating crop area using seasonal time series of enhanced vegetation index from MODIS satellite imagery

Cereal grain is one of the main export commodities of Australian agriculture. Over the past decade, crop yield forecasts for wheat and sorghum have shown appreciable utility for industry planning at shire, state, and national scales. There is now an increasing drive from industry for more accurate and cost-effective crop production forecasts. In order to generate production estimates, accurate crop area estimates are needed by the end of the cropping season. Multivariate methods for analysing remotely sensed Enhanced Vegetation Index (EVI) from 16-day Moderate Resolution Imaging Spectroradiometer (MODIS) satellite imagery within the cropping period (i.e. April-November) were investigated to estimate crop area for wheat, barley, chickpea, and total winter cropped area for a case study region in NE Australia. Each pixel classification method was trained on ground truth data collected from the study region. Three approaches to pixel classification were examined: (i) cluster analysis of trajectories of EVI values from consecutive multi-date imagery during the crop growth period; (ii) harmonic analysis of the time series (HANTS) of the EVI values; and (iii) principal component analysis (PCA) of the time series of EVI values. Images classified using these three approaches were compared with each other, and with a classification based on the single MODIS image taken at peak EVI. Imagery for the 2003 and 2004 seasons was used to assess the ability of the methods to determine wheat, barley, chickpea, and total cropped area estimates. The accuracy at pixel scale was determined by the percent correct classification metric by contrasting all pixel scale samples with independent pixel observations. At a shire level, aggregated total crop area estimates were compared with surveyed estimates. All multi-temporal methods showed significant overall capability to estimate total winter crop area. There was high accuracy at pixel scale (>98% correct classification) for identifying overall winter cropping. However, discrimination among crops was less accurate. Although the use of single-date EVI data produced high accuracy for estimates of wheat area at shire scale, the result contradicted the poor pixel-scale accuracy associated with this approach, due to fortuitous compensating errors. Further studies are needed to extrapolate the multi-temporal approaches to other geographical areas and to improve the lead time for deriving cropped-area estimates before harvest.

[1]  A. Huete,et al.  Overview of the radiometric and biophysical performance of the MODIS vegetation indices , 2002 .

[2]  D. Stephens.  Objective criteria for estimating the severity of drought in the wheat cropping areas of Australia , 1998 .

[3]  A. Strahler,et al.  Monitoring vegetation phenology using MODIS , 2003 .

[4]  A. Huete,et al.  A comparison of vegetation indices over a global set of TM images for EOS-MODIS , 1997 .

[5]  Graeme L. Hammer,et al.  Infusing the use of seasonal climate forecasting into crop management practice in North East Australia using discussion support software , 2002 .

[6]  Changsheng Li,et al.  Mapping paddy rice agriculture in southern China using multi-temporal MODIS images , 2005 .

[7]  Holger Meinke,et al.  Improving wheat simulation capabilities in Australia from a cropping systems perspective: water and nitrogen effects on spring wheat in a semi-arid environment , 1997 .

[8]  Graeme L. Hammer,et al.  A simple regional-scale model for forecasting sorghum yield across North-Eastern Australia , 2005 .

[9]  J. C. Price.  Comparing MODIS and ETM+ data for regional and global land classification , 2003 .

[10]  Graeme L. Hammer,et al.  Spatial and temporal patterns in Australian wheat yield and their relationship with ENSO , 2002 .

[11]  A. Strahler,et al.  Application of the MODIS global supervised classification model to vegetation and land cover mapping of Central America , 2000 .

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

[13]  Holger Meinke,et al.  Improving wheat simulation capabilities in Australia from a cropping systems perspective III. The integrated wheat model (I_WHEAT) , 1998 .

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

[15]  D. Legates,et al.  Crop identification using harmonic analysis of time-series AVHRR NDVI data , 2002 .

[16]  C. Justice,et al.  Development of vegetation and soil indices for MODIS-EOS , 1994 .

[17]  J. Townshend,et al.  Detection of land cover changes using MODIS 250 m data , 2002 .

[18]  G. Hammer,et al.  Prediction of global rainfall probabilities using phases of the Southern Oscillation Index , 1996, Nature.

[19]  Graeme L. Hammer,et al.  Advances in application of climate prediction in agriculture , 2001 .