The accurate and timely spatial distribution information of various crop types is vital for food security. In this study, the 2-band enhanced vegetation index (EVI2) data from Moderate Resolution Imaging Spectroradiometer (MODIS) were combined with sparse representation approach to identify the distribution of various crop types in Shandong Province, China. Three groups of input variables derived from EVI2 including annual time series EVI2 (TS), harmonic features (HF), and combined vector formed by harmonic features and texture features (HFT) were used. The online dictionary learning and orthogonal matching pursuit algorithms were applied to generate the dictionary and solve the sparse coefficients of the identified samples, respectively. Then, the label of the identified samples can be obtained according to the minimum residuals between the dictionary and the sparse coefficients of the identified samples. At the provincial level, the validation based on the statistical data showed that three groups of input variables presented lower than ±25% at relative errors, and input variables of HFT performed better than the other two. At the municipal level, the results achieved by using input variables of HFT also agreed well with the statistics with the coefficient of determination R 2 > 0.85 for wheat and maize, as well as R 2 > 0.71 for peanut and cotton during 2014-2016. These results demonstrate that the combination of the input variables of HFT derived from time series MODIS EVI2 data and sparse representation approach can be used for crop identification in the study area. INDEX TERMS Crop area identification, harmonic features, sparse representation, texture features, time series EVI2.