Crop classification using temporal stacks of multispectral satellite imagery

The increase in performance, availability, and coverage of multispectral satellite sensor constellations has led to a drastic increase in data volume and data rate. Multi-decadal remote sensing datasets at the petabyte scale are now available in commercial clouds, with new satellite constellations generating petabytes/year of daily high-resolution global coverage imagery. The data analysis capability, however, has lagged behind storage and compute developments, and has traditionally focused on individual scene processing. We present results from an ongoing effort to develop satellite imagery analysis tools that aggregate temporal, spatial, and spectral information and can scale with the high-rate and dimensionality of imagery being collected. We investigate and compare the performance of pixel-level crop identification using tree-based classifiers and its dependence on both temporal and spectral features. Classification performance is assessed using as ground-truth Cropland Data Layer (CDL) crop masks generated by the US Department of Agriculture (USDA). The CDL maps contain 30m spatial resolution, pixel-level labels for around 200 categories of land cover, but are however only available post-growing season. The analysis focuses on McCook county in South Dakota and shows crop classification using a temporal stack of Landsat 8 (L8) imagery over the growing season, from April through October. Specifically, we consider the temporal L8 stack depth, as well as different normalized band difference indices, and evaluate their contribution to crop identification. We also show an extension of our algorithm to map corn and soy crops in the state of Mato Grosso, Brazil.

[1]  Leo Breiman,et al.  Classification and Regression Trees , 1984 .

[2]  R. Myneni,et al.  The interpretation of spectral vegetation indexes , 1995 .

[3]  B. Gao NDWI—A normalized difference water index for remote sensing of vegetation liquid water from space , 1996 .

[4]  S. K. McFeeters The use of the Normalized Difference Water Index (NDWI) in the delineation of open water features , 1996 .

[5]  Leo Breiman,et al.  Random Forests , 2001, Machine Learning.

[6]  Garrett T. Kenyon,et al.  Large-scale functional models of visual cortex for remote sensing , 2009, 2009 IEEE Applied Imagery Pattern Recognition Workshop (AIPR 2009).

[7]  Olivier Arino,et al.  Global Land Cover Map for 2009 (GlobCover 2009) , 2012 .

[8]  Daniela I. Moody,et al.  Undercomplete learned dictionaries for land cover classification in multispectral imagery of Arctic landscapes using CoSA: clustering of sparse approximations , 2013, Defense, Security, and Sensing.

[9]  Daniela I. Moody,et al.  Land cover classification in multispectral imagery using clustering of sparse approximations over learned feature dictionaries , 2014 .

[10]  Mark Johnson,et al.  Seeing the Earth in the Cloud: Processing one petabyte of satellite imagery in one day , 2015, 2015 IEEE Applied Imagery Pattern Recognition Workshop (AIPR).

[11]  Michael S. Warren,et al.  Land cover classification in fused multisensor multispectral satellite imagery , 2016, 2016 IEEE Southwest Symposium on Image Analysis and Interpretation (SSIAI).