A high-performance and in-season classification system of field-level crop types using time-series Landsat data and a machine learning approach

Abstract Accurate and timely spatial classification of crop types based on remote sensing data is important for both scientific and practical purposes. Spatially explicit crop-type information can be used to estimate crop areas for a variety of monitoring and decision-making applications such as crop insurance, land rental, supply-chain logistics, and financial market forecasting. However, there is no publically available spatially explicit in-season crop-type classification information for the U.S. Corn Belt (a landscape predominated by corn and soybean). Instead, researchers and decision-makers have to wait until four to six months after harvest to have such information from the previous year. The state-of-the-art research on crop-type classification has been shifted from relying on only spectral features of single static images to combining together spectral and time-series information. While Landsat data have a desirable spatial resolution for field-level crop-type classification, the ability to extract temporal phenology information based on Landsat data remains a challenge due to low temporal revisiting frequency and inevitable cloud contamination. To address this challenge and generate accurate, cost-effective, and in-season crop-type classification, this research uses the USDA's Common Land Units (CLUs) to aggregate spectral information for each field based on a time-series Landsat image data stack to largely overcome the cloud contamination issue while exploiting a machine learning model based on Deep Neural Network (DNN) and high-performance computing for intelligent and scalable computation of classification processes. Experiments were designed to evaluate what information is most useful for training the machine learning model for crop-type classification, and how various spatial and temporal factors affect the crop-type classification performance in order to derive timely crop type information. All experiments were conducted over Champaign County located in central Illinois, and a total of 1322 Landsat multi-temporal scenes including all the six optical spectral bands spanning from 2000 to 2015 were used. Computational experiments show the inclusion of temporal phenology information and evenly distributed spatial training samples in the study domain improves classification performance. The shortwave infrared bands show notably better performance than the widely used visible and near-infrared bands for classifying corn and soybean. In comparison with USDA's Crop Data Layer (CDL), this study found a relatively high Overall Accuracy (i.e. the number of the corrected classified fields divided by the number of the total fields) of 96% for classifying corn and soybean across all CLU fields in the Champaign County from 2000 to 2015. Furthermore, our approach achieved 95% Overall Accuracy by late July of the concurrent year for classifying corn and soybean. The findings suggest the methodology presented in this paper is promising for accurate, cost-effective, and in-season classification of field-level crop types, which may be scaled up to large geographic extents such as the U.S. Corn Belt.

[1]  Martha C. Anderson,et al.  Toward mapping crop progress at field scales through fusion of Landsat and MODIS imagery , 2017 .

[2]  C. Field,et al.  A reanalysis using improved leaf models and a new canopy integration scheme , 1992 .

[3]  Zhiming Feng,et al.  Cross-Comparison of Vegetation Indices Derived from Landsat-7 Enhanced Thematic Mapper Plus (ETM+) and Landsat-8 Operational Land Imager (OLI) Sensors , 2013, Remote. Sens..

[4]  D. Lobell,et al.  Cropland distributions from temporal unmixing of MODIS data , 2004 .

[5]  ImageNet Classification with Deep Convolutional Neural , 2013 .

[6]  Changsheng Li,et al.  Mapping paddy rice agriculture in South and Southeast Asia using multi-temporal MODIS images , 2006 .

[7]  B. Wardlow,et al.  Analysis of time-series MODIS 250 m vegetation index data for crop classification in the U.S. Central Great Plains , 2007 .

[8]  Thomas R. Loveland,et al.  A review of large area monitoring of land cover change using Landsat data , 2012 .

[9]  Geoffrey E. Hinton,et al.  Learning representations by back-propagating errors , 1986, Nature.

[10]  Chengquan Huang,et al.  Global, 30-m resolution continuous fields of tree cover: Landsat-based rescaling of MODIS vegetation continuous fields with lidar-based estimates of error , 2013, Int. J. Digit. Earth.

[11]  Martha C. Anderson,et al.  Landsat-8: Science and Product Vision for Terrestrial Global Change Research , 2014 .

[12]  J. H. Everitt,et al.  Evaluating high resolution SPOT 5 satellite imagery to estimate crop yield , 2009, Precision Agriculture.

[13]  David M. Johnson,et al.  A multi-resolution approach to national-scale cultivated area estimation of soybean , 2017 .

[14]  Saskia Foerster,et al.  Crop type mapping using spectral-temporal profiles and phenological information , 2012 .

[15]  A. Viña,et al.  Remote estimation of leaf area index and green leaf biomass in maize canopies , 2003 .

[16]  Tara N. Sainath,et al.  Deep Neural Networks for Acoustic Modeling in Speech Recognition: The Shared Views of Four Research Groups , 2012, IEEE Signal Processing Magazine.

[17]  Douglas K. Bolton,et al.  Forecasting crop yield using remotely sensed vegetation indices and crop phenology metrics , 2013 .

[18]  Chengquan Huang,et al.  Global characterization and monitoring of forest cover using Landsat data: opportunities and challenges , 2012, Int. J. Digit. Earth.

[19]  David B. Lobell,et al.  The use of satellite data for crop yield gap analysis , 2013 .

[20]  Fan Ruixia A Multi-resolution Approach for Image Tracking , 2002 .

[21]  J. Townshend,et al.  Global land cover classi(cid:142) cation at 1 km spatial resolution using a classi(cid:142) cation tree approach , 2004 .

[22]  Jason Weston,et al.  A unified architecture for natural language processing: deep neural networks with multitask learning , 2008, ICML '08.

[23]  Neil Flood,et al.  Continuity of Reflectance Data between Landsat-7 ETM+ and Landsat-8 OLI, for Both Top-of-Atmosphere and Surface Reflectance: A Study in the Australian Landscape , 2014, Remote. Sens..

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

[25]  Guigang Zhang,et al.  Deep Learning , 2016, Int. J. Semantic Comput..

[26]  Stefano Ermon,et al.  Deep Gaussian Process for Crop Yield Prediction Based on Remote Sensing Data , 2017, AAAI.

[27]  Matthew C. Hansen,et al.  Corn and Soybean Mapping in the United States Using MODIS Time‐Series Data Sets , 2007 .

[28]  Tim R. McVicar,et al.  Determining temporal windows for crop discrimination with remote sensing: a case study in south-eastern Australia , 2004 .

[29]  Feng Gao,et al.  Landsat Ecosystem Disturbance Adaptive Processing System LEDAPS algorithm description , 2013 .

[30]  Qiming Qin,et al.  Estimating crop water stress with ETM+ NIR and SWIR data , 2008 .

[31]  M. Bauer,et al.  Land cover classification and change analysis of the Twin Cities (Minnesota) Metropolitan Area by multitemporal Landsat remote sensing , 2005 .

[32]  Mathew R. Schwaller,et al.  On the blending of the Landsat and MODIS surface reflectance: predicting daily Landsat surface reflectance , 2006, IEEE Transactions on Geoscience and Remote Sensing.

[33]  Limin Yang,et al.  Development of a 2001 National land-cover database for the United States , 2004 .

[34]  Chenghai Yang,et al.  Original paper: Evaluating high resolution SPOT 5 satellite imagery for crop identification , 2011 .

[35]  Limin Yang,et al.  COMPLETION OF THE 1990S NATIONAL LAND COVER DATA SET FOR THE CONTERMINOUS UNITED STATES FROM LANDSAT THEMATIC MAPPER DATA AND ANCILLARY DATA SOURCES , 2001 .

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

[37]  Jiyuan Liu,et al.  Characterization of forest types in Northeastern China, using multi-temporal SPOT-4 VEGETATION sensor data , 2002 .

[38]  Jürgen Schmidhuber,et al.  Deep learning in neural networks: An overview , 2014, Neural Networks.

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

[40]  D. Roy,et al.  Monitoring conterminous United States (CONUS) land cover change with Web-Enabled Landsat Data (WELD) , 2014 .

[41]  Wang Chong-chan Land cover classification using RADARSAT-2 full polarmetric SAR data , 2015 .

[42]  Peijuan Wang,et al.  Integrating remote sensing data from multiple optical sensors for ecological and crop condition monitoring , 2013, Optics & Photonics - Optical Engineering + Applications.

[43]  Peijuan Wang,et al.  Fusing Landsat and MODIS Data for Vegetation Monitoring , 2015, IEEE Geoscience and Remote Sensing Magazine.

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

[45]  David P. Roy,et al.  Continuous fields of land cover for the conterminous United States using Landsat data: first results from the Web-Enabled Landsat Data (WELD) project , 2011 .

[46]  Steven E. Franklin,et al.  A comparison of pixel-based and object-based image analysis with selected machine learning algorithms for the classification of agricultural landscapes using SPOT-5 HRG imagery , 2012 .

[47]  J. Legarsky,et al.  Land-cover classification using radarsat and landsat imagery for St. Louis, Missouri , 2007 .

[48]  Conghe Song,et al.  Long-term land cover dynamics by multi-temporal classification across the Landsat-5 record , 2013 .

[49]  B. Wardlow,et al.  Large-area crop mapping using time-series MODIS 250 m NDVI data: An assessment for the U.S. Central Great Plains , 2008 .

[50]  D. Lobell,et al.  A scalable satellite-based crop yield mapper , 2015 .

[51]  Matthew C. Hansen,et al.  National-scale soybean mapping and area estimation in the United States using medium resolution satellite imagery and field survey , 2017 .

[52]  Zhengwei Yang,et al.  Monitoring US agriculture: the US Department of Agriculture, National Agricultural Statistics Service, Cropland Data Layer Program , 2011 .

[53]  H. Tian,et al.  Spatial and temporal patterns of China's cropland during 1990¿2000: An analysis based on Landsat TM data , 2005 .

[54]  I. Papila,et al.  A COMPARISON OF PIXEL-BASED AND OBJECT-BASED CLASSIFICATION METHODS , A CASE STUDY : ISTANBUL , TURKEY , 2011 .