An Improved Approach Considering Intraclass Variability for Mapping Winter Wheat Using Multitemporal MODIS EVI Images

Winter wheat is one of the major cereal crops in the world. Monitoring and mapping its spatial distribution has significant implications for agriculture management, water resources utilization, and food security. Generally, winter wheat has distinguished phenological stages during the growing season, which form a unique EVI (Enhanced Vegetation Index) time series curve and differ considerably from other crop types and natural vegetation. Since early 2000, the MODIS EVI product has become the primary dataset for satellite-based crop monitoring at large scales due to its high temporal resolution, huge observation scope, and timely availability. However, the intraclass variability of winter wheat caused by field conditions and agricultural practices might lower the mapping accuracy, which has received little attention in previous studies. Here, we present a winter wheat mapping approach that integrates the variables derived from the MODIS EVI time series taking into account intraclass variability. We applied this approach to two winter wheat concentration areas, the state of Kansas in the U.S. and the North China Plain region (NCP). The results were evaluated against crop-specific maps or statistical data at the state/regional level, county level, and site level. Compared with statistical data, the accuracies in Kansas and the NCP were 95.1% and 92.9% at the state/regional level with R2 (Coefficient of Determination) values of 0.96 and 0.71 at the county level, respectively. Overall accuracies in confusion matrix were evaluated by validation samples in both Kansas (90.3%) and the NCP (85.0%) at the site level. Comparisons with methods without considering intraclass variability demonstrated that winter wheat mapping accuracies were improved by 17% in Kansas and 15% in the NCP using the improved approach. Further analysis indicated that our approach performed better in areas with lower landscape fragmentation, which may partly explain the relatively higher accuracy of winter wheat mapping in Kansas. This study provides a new perspective for generating multiple subclasses as training inputs to decrease the intraclass differences for crop type detection based on the MODIS EVI time series. This approach provides a flexible framework with few variables and fewer training samples that could facilitate its application to multiple-crop-type mapping at large scales. Remote Sens. 2019, 11, 1191; doi:10.3390/rs11101191 www.mdpi.com/journal/remotesensing Remote Sens. 2019, 11, 1191 2 of 24

[1]  D. Tilman,et al.  Global food demand and the sustainable intensification of agriculture , 2011, Proceedings of the National Academy of Sciences.

[2]  Laurence Hubert-Moy,et al.  A Comparison of Parametric Classification Procedures of Remotely Sensed Data Applied on Different Landscape Units , 2001 .

[3]  Aigong Xu,et al.  Winter wheat mapping using temporal signatures of MODIS vegetation index data , 2012 .

[4]  John A. Richards,et al.  Remote Sensing Digital Image Analysis , 1986 .

[5]  J. Campbell Introduction to remote sensing , 1987 .

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

[7]  R. Lunetta,et al.  Land-cover characterization and change detection using multitemporal MODIS NDVI data , 2005, International Workshop on the Analysis of Multi-Temporal Remote Sensing Images, 2005..

[8]  C. Field,et al.  Global scale climate–crop yield relationships and the impacts of recent warming , 2007, Environmental Research Letters.

[9]  Zhu Ji-guan Agricultural Efficiency and Its Decomposition Based on DEA in the Huang-Huai-Hai Plain , 2013 .

[10]  P. C. Doraiswamy,et al.  Crop classification in the U.S. Corn Belt using MODIS imagery , 2007, 2007 IEEE International Geoscience and Remote Sensing Symposium.

[11]  Zhengwei Yang,et al.  Making Cropland Data Layer Data Accessible and Actionable in GIS Education , 2014 .

[12]  Huajun Tang,et al.  Regional yield estimation for winter wheat with MODIS-NDVI data in Shandong, China , 2008, Int. J. Appl. Earth Obs. Geoinformation.

[13]  Russell G. Congalton,et al.  Assessing the accuracy of remotely sensed data : principles and practices , 1998 .

[14]  Russell G. Congalton,et al.  MODIS phenology-derived, multi-year distribution of conterminous U.S. crop types , 2017 .

[15]  Eloise Kendy,et al.  Groundwater Exploitation and Its Impact on the Environment in the North China Plain , 2001 .

[16]  Mark A. Friedl,et al.  Mapping Crop Cycles in China Using MODIS-EVI Time Series , 2014, Remote. Sens..

[17]  Yang Song,et al.  Mapping Winter Wheat Planting Area and Monitoring Its Phenology Using Sentinel-1 Backscatter Time Series , 2019, Remote. Sens..

[18]  Armando Apan,et al.  Early-season crop area estimates for winter crops in NE Australia using MODIS satellite imagery , 2010 .

[19]  Francisca López-Granados,et al.  Airborne multi-spectral imagery for mapping cruciferous weeds in cereal and legume crops , 2011, Precision Agriculture.

[20]  Naoto Yokoya,et al.  Hyperspectral Pansharpening: A Review , 2015, IEEE Geoscience and Remote Sensing Magazine.

[21]  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 .

[22]  Chongcheng Chen,et al.  Winter wheat mapping combining variations before and after estimated heading dates , 2017 .

[23]  Giles M. Foody,et al.  A relative evaluation of multiclass image classification by support vector machines , 2004, IEEE Transactions on Geoscience and Remote Sensing.

[24]  C. Justice,et al.  A generalized regression-based model for forecasting winter wheat yields in Kansas and Ukraine using MODIS data , 2010 .

[25]  Russell G. Congalton,et al.  A review of assessing the accuracy of classifications of remotely sensed data , 1991 .

[26]  S. Piao,et al.  Spring vegetation green-up date in China inferred from SPOT NDVI data: A multiple model analysis , 2012 .

[27]  D. J. Booth,et al.  A comparison of classification algorithms in terms of speed and accuracy after the application of a post-classification modal filter , 1989 .

[28]  Ramanathan Sugumaran,et al.  Classification of Iowa wetlands using an airborne hyperspectral image: a comparison of the spectral angle mapper classifier and an object-oriented approach , 2005 .

[29]  Michael Frankfurter,et al.  Numerical Recipes In C The Art Of Scientific Computing , 2016 .

[30]  S. Robinson,et al.  Food Security: The Challenge of Feeding 9 Billion People , 2010, Science.

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

[32]  Christopher Conrad,et al.  Temporal segmentation of MODIS time series for improving crop classification in Central Asian irrigation systems , 2011 .

[33]  Xuehong Bai,et al.  Multiple cropping intensity in China derived from agro-meteorological observations and MODIS data , 2014, Chinese Geographical Science.

[34]  Yanjun Shen,et al.  Effects of irrigation on water balance, yield and WUE of winter wheat in the North China Plain , 2006 .

[35]  Giles M. Foody,et al.  Toward intelligent training of supervised image classifications: directing training data acquisition for SVM classification , 2004 .

[36]  Miao Zhang,et al.  Mapping Winter Wheat Biomass and Yield Using Time Series Data Blended from PROBA-V 100- and 300-m S1 Products , 2016, Remote. Sens..

[37]  Yang Liu,et al.  The relationship between threshold-based and inflexion-based approaches for extraction of land surface phenology , 2017 .

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

[39]  Xin Huang,et al.  A Phenology-Based Method to Map Cropping Patterns under a Wheat-Maize Rotation Using Remotely Sensed Time-Series Data , 2018, Remote. Sens..

[40]  J. I. Ortiz-Monasterio,et al.  Extreme heat effects on wheat senescence in India , 2012 .

[41]  Giorgos Mallinis,et al.  A Hidden Markov Models Approach for Crop Classification: Linking Crop Phenology to Time Series of Multi-Sensor Remote Sensing Data , 2015, Remote. Sens..

[42]  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 .

[43]  B. Datt,et al.  On the relationship between training sample size and data dimensionality: Monte Carlo analysis of broadband multi-temporal classification , 2005 .

[44]  Yu Wang,et al.  Mapping winter wheat using phenological feature of peak before winter on the North China Plain based on time-series MODIS data , 2017 .

[45]  Huadong Guo,et al.  Detecting winter wheat phenology with SPOT-VEGETATION data in the North China Plain , 2014 .

[46]  Yang Shao,et al.  An evaluation of time-series smoothing algorithms for land-cover classifications using MODIS-NDVI multi-temporal data , 2016 .

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

[48]  Limin Wang,et al.  Mapping crop phenology using NDVI time-series derived from HJ-1 A/B data , 2015, Int. J. Appl. Earth Obs. Geoinformation.

[49]  Huib Hengsdijk,et al.  Quantifying production potentials of winter wheat in the North China Plain , 2006 .

[50]  Felix Rembold,et al.  Estimation of Inter-annual Crop Area Variation by the Application of Spectral Angle Mapping to Low Resolution Multitemporal NDVI Images , 2006 .

[51]  G. M. Paulsen,et al.  Genetic Gain in Yield Attributes of Winter Wheat in the Great Plains , 2001 .

[52]  A. Savitzky,et al.  Smoothing and Differentiation of Data by Simplified Least Squares Procedures. , 1964 .

[53]  Bingfang Wu,et al.  Crop classification using HJ satellite multispectral data in the North China Plain , 2013 .

[54]  Giles M. Foody,et al.  Training set size requirements for the classification of a specific class , 2006 .

[55]  Alan H. Strahler,et al.  Global land cover mapping from MODIS: algorithms and early results , 2002 .

[56]  Xiaoguang Jiang,et al.  Comparison of artificial neural networks and support vector machine classifiers for land cover classification in Northern China using a SPOT-5 HRG image , 2012 .

[57]  Chi-Farn Chen,et al.  A Phenology-Based Classification of Time-Series MODIS Data for Rice Crop Monitoring in Mekong Delta, Vietnam , 2013, Remote. Sens..

[58]  N. Fuhs,et al.  Evaluation of several schemes for classification of remotely sensed data , 1980 .

[59]  H. Grau,et al.  Deforestation and fragmentation of Chaco dry forest in NW Argentina (1972–2007) , 2009 .

[60]  Mark Monmonier,et al.  Measures of Pattern Complexity for Choroplethic Maps , 1974 .

[61]  Alan H. Strahler,et al.  The Moderate Resolution Imaging Spectroradiometer (MODIS): land remote sensing for global change research , 1998, IEEE Trans. Geosci. Remote. Sens..

[62]  A. Huete,et al.  A 30-m landsat-derived cropland extent product of Australia and China using random forest machine learning algorithm on Google Earth Engine cloud computing platform , 2018, ISPRS Journal of Photogrammetry and Remote Sensing.

[63]  Rajendra Prasad,et al.  Comparison of support vector machine, artificial neural network, and spectral angle mapper algorithms for crop classification using LISS IV data , 2015 .

[64]  Li Zhang,et al.  Monitoring winter wheat growth in North China by combining a crop model and remote sensing data , 2008, Int. J. Appl. Earth Obs. Geoinformation.

[65]  Timothy A. Warner,et al.  Kernel-based extreme learning machine for remote-sensing image classification , 2013 .

[66]  Qiang Yu,et al.  Modeling Wheat and Maize Productivity as Affected by Climate Variation and Irrigation Supply in North China Plain , 2010 .

[67]  J. I. Ortiz-Monasterio,et al.  Climate change: Can wheat beat the heat? , 2008 .

[68]  C. Lo,et al.  Using a time series of satellite imagery to detect land use and land cover changes in the Atlanta, Georgia metropolitan area , 2002 .

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

[70]  Damien Sulla-Menashe,et al.  Winter wheat area estimation from MODIS-EVI time series data using the Crop Proportion Phenology Index , 2012 .

[71]  Yang Shao,et al.  Monitoring agricultural cropping patterns across the Laurentian Great Lakes Basin using MODIS-NDVI data , 2010, Int. J. Appl. Earth Obs. Geoinformation.

[72]  Clement Atzberger,et al.  Advances in Remote Sensing of Agriculture: Context Description, Existing Operational Monitoring Systems and Major Information Needs , 2013, Remote. Sens..