Mapping Maize Cultivated Area Combining MODIS EVI Time Series and the Spatial Variations of Phenology over Huanghuaihai Plain

: Crop phenology is a significant factor that affects the precision of crop area extraction by using the multi-temporal vegetation indices (VIs) approach. Considering the phenological differences of maize among the different regions, the summer maize cultivated area was estimated by using enhanced vegetation index (EVI) time series images from the Moderate Resolution Imaging Spectroradiometer (MODIS) over the Huanghuaihai Plain in China. By analyzing the temporal shift in summer maize calendars, linear regression equations for simulating the summer maize phenology were obtained. The simulated maize phenology was used to correct the MODIS EVI time series curve of summer maize. Combining the mean absolute distance (MAD) and p-tile algorithm, the cultivated areas of summer maize were distinguished over the Hunaghuaihai Plain. The accuracy of the extraction results in each province was above 85%. Comparing the maize area of two groups from MODIS-estimated and statistical data, the validation results showed that the R 2 reached 0.81 at the city level and 0.69 at the county level. It demonstrated that the approach in this study has the ability to effectively map the summer maize area over a large scale and provides a novel idea for estimating the planting area of other crops.


Introduction
Maize, as the primary staple, plays a vital role in agricultural production in China. It is very significant to estimate the cultivated area of maize for national food security and sustainable economic development [1,2]. Huanghuaihai Plain covers approximately 4.5% of the whole country, whereas the maize planting area accounts for about 30% of the total maize acreage, which is the main production areas of summer maize in China. The distribution of maize in this region affects grain policy making, adjustment of the cropping system structure, and academic studies related to maize in China [3,4]. Therefore, there is an increasing need for objective, timely, and accurate estimation of the maize area over the Huanghuaihai Plain.
As a kind of earth observation technology, remote sensing can effectively acquire the spatial distribution and spectral information of ground objects, with wide geographic coverage and the merit of large information capacity, high accuracy, and speed, which provides opportunities for crop identification [5][6][7][8][9]. Different crops have various spectral characteristics in remotely sensed images [10,11]. The multi-temporal vegetation indices (VIs) reflect the changes in the spectral features of vegetation over time, which can discriminate different crop types [12][13][14][15][16]. Many previous research studies employed the multi-temporal VIs approaches to map the maize cultivated area. Wardlow and Egbert used temporal normalized difference vegetation index (NDVI) data from 250 m moderate resolution imaging spectroradiometer (MODIS) images to classify the maize over the U.S. Central Great Plains [17]. Gu et al. reconstructed a mid-scale time series NDVI dataset derived from the fusion of MODIS and thematic mapper (TM) images based on the wavelet transform to estimate the maize area in Yuanyang County of Henan Province [18]. Maguranyanga and Murwira combined the temporal NDVI series and ancillary field data to map the maize area in the large-scale commercial farms of Zimbabwe though the maximum entropy method [19]. Tang et al. adopted the identification method for spring maize based on spectral and phenological features derived from the MODIS land surface reflectance time series data in Northeast China [20]. In contrast, although these studies were valuable for the recognition of maize, there are still some challenges to be faced in mapping the spatial distribution of maize on a large scale due to the spatial variations of maize phenology.
The spectral features of summer maize on remote sensing images are different at the same time, as the maize phenology varies in different regions. In order to weaken the influence of phenology differences on the accuracy of crop area extraction, some researchers have considered the influence of latitude on summer maize phenology, when using the multi-temporal VIs approach to map the maize area on a large scale. For example, Zhang et al. used the lag values of phenology change with latitude to modify the MODIS-enhanced vegetation index (EVI) time series curve of maize and mapped the spring maize in Northeast China [21]. Wang et al. also considered the differences in maize phenology with latitude to discriminate the summer maize based on the multi-temporal MODIS EVI images over the Huanghuaihai Plain [22]. In both studies, the influence of latitude on crop phenology was taken into account when using multi-temporal MODIS images mapping the summer maize area on a large scale, which had improved the accuracy of maize area extraction.
In fact, the phenology of summer maize is not only affected by latitude but also sensitive to the change in environmental conditions, as well as changes with the environmental variation. The maize phenology is sensitive to changes in environmental conditions, resulting in the maize phenology that would be altered by environmental variability [23][24][25][26][27][28]. In the past several years, the apparent warming climate caused early flowing and maturity and consequently shortened the reproductive growth stage of maize [29]. The mean temperature as one of the thermal conditions is closely related to the sowing and seedling emergence time [30][31][32]. The accumulated temperature is routinely used to determine crop planting schedules, so the accumulated temperature shift will lead to the change in maize growth stages [33]. Moisture is an essential requirement for maize growth, so precipitation has effects on the growth of maize. It was reported that the correlations between precipitation and the length of various maize growing phases were positive [26,34]. Furthermore, geographical position is a vital factor for maize growth [35,36]. Liu et al. found that for every 1° increase in latitude, northward, the growth durations of sowing to emergence and emergence to silking were significantly increased by 0.7 and 1.25 d, respectively [37]. This means the maize phenology has been changed in different regions. It can be clearly seen that the phenology of maize varied with the changes in annual average temperature, accumulated temperature, precipitation, and geographical location. On a large scale, with a great difference in the environment of different regions, the phenologies of maize are different, which lead to the difference in the spectral characteristics of summer maize in the remote sensing image. There is, therefore, the phenological differences of maize caused by various environmental factors being combined with the multi-temporal MODIS EVI images to map the cultivated area of summer maize for better extraction accuracy on a large scale.
In the present study, multiple linear regression equations were constructed to simulate the summer maize phenology through analyzing the relationship between environmental factors and summer maize phenology over the Huanghuaihai Plain. The simulated phenology of summer maize was obtained through environmental factors data (longitude, latitude, annual average temperature, annual precipitation, and ≥10 °C accumulated temperature), which acted as an independent variable to modify the summer maize MODIS EVI time series curve and generate the standard maize EVI time series curve over Huanghuaihai Plain. Combined with the mean absolute distance (MAD) and the Ptile algorithm, the summer maize planting area was discriminated successfully in the study area. It provides a novel idea for improving the accuracy of maize area extraction using multi-temporal remote sensed data on a large scale.

Study area and Reference Area
Huanghuaihai Plain (Figure 1) (32°-40° N, 114°-121° E), which includes Beijing, Tianjin, Hebei, Henan, and Shandong, is the main planting area of summer maize in China. The region belongs to the monsoon climate of medium latitudes with four clearly distinct seasons [38]. The annual precipitation is 500-900 mm. The precipitation is mostly concentrated in summer, accounting for about 70% of the annual precipitation. The mean annual temperature is 8-15 °C [39]. Summer maize phenology is significantly different due to the difference in climate among regions over Huanghuaihai Plain.
Ruzhou City (Figure 1), with a relatively single planting structure, was selected as the reference area, which is located in the middle of Henan Provence (33°56′-34°20′ N,112°31′-113°07′ E). Summer maize is the main summer crop in the region, accounting for about 70% of the summer crop planting area. Therefore, choosing Ruzhou City as the reference area can reduce the event of mixed classification.

Landsat7 ETM+ Image and Pre-Processing
The spatial resolution of Landsat7 ETM+ (https://glovis.usgs.gov) is 30 m, enabling highprecision extraction of crop area on a small scale [40,41]. According to the weather conditions, image quality, and phonological characteristics of summer maize, the ETM+ images on June 9th (sowing), September 13th (milk-ripe), and October 31st (harvested stage), 2012, were selected to obtain the summer maize cultivated area in the reference area. Striping, radiation calibration, atmospheric correction based on the FLAASH model, and clipping were performed for 3 images before classification.

MODIS EVI Data and Pre-Processing
MOD13Q1 data were derived from the National Aeronautics and Space Administration (NASA) (https://search.earthdata.nasa.gov). MOD13Q1 data provide an EVI vegetation layer every 16 days at 250 m spatial resolution. The summer maize growth period (2012-06-01-2012-10-31) contained MOD13Q1 data for 9 phases, covering h26v04, h26v05, h27v04, and h27v05. First, the images in the same phase were spliced, and the spliced images were then subjected to projection conversion processing. The projected coordinate system was UTM_Zone_51N, and the geographic coordinate system was CGS_WGS_1984. Finally, the boundary data of Huanghuaihai and Ruzhou were used to mask the processed images.

Land Use Data
The land use data (Globeland30, http://www.globallandcover.com/GLC30Download/index.aspx) of spatial 30 m resolution in 2010 were employed to modify the extracted area of summer maize in the reference area. The cultivated land (code 10) dataset in the land use data was overlaid with the spatial distribution of summer maize from ETM+ images to exclude the influence of other land use types on the extraction results in the reference area.

Phenological and Meteorological Data
Summer maize phenological and meteorological data of Huanghuaihai Plain in 2012 were derived from the National Meteorological Information Center (http://data.cma.cn). In this study, the observation phenology data from 32 agricultural meteorological stations are gathered, including latitude, longitude, maize fertility name, and fertility date that was converted into the day of the year (DOY) at each site. The location of agro-meteorological stations is shown in Figure 2, and the phenological information is listed in Table 1. The grid dataset of the daily value of surface temperature and precipitation in China with the spatial resolution of 0.5° was employed to obtain the meteorological elements in the study area. The grid dataset was generated by the thin plate spline (TPS) method based on the data of precipitation and air temperature of the high-density site on the ground in China in the special national information materials, as well as the GTOPO30 data. The grid dataset was defined as the CGS_WGS_1984 geographic coordinate system and was projected into the UTM_Zone_51N coordinate system. We resampled the grid data with 250 m spatial resolution by the bilinear algorithm. Based on the processed grid dataset, the annual average temperature, annual precipitation, and ≥10 °C accumulated temperature in Huanghuaihai Plain were calculated from the grid data with 250 m spatial resolution. In addition, the annual average temperature, annual precipitation, and ≥10 °C accumulated temperature at 32 agro-meteorological stations were obtained ( Table 1).

Statistical Data
The statistics of the summer maize cultivated area are from the statistical yearbooks [42][43][44][45][46][47][48][49][50][51]. The statistics of the total surface area of summer maize on 43 city-levels were collected, including Beijing, Tianjin, 7 cities in Hebei, 16 cities in Henan, and 18 cities in Shandong. Meanwhile, the statistics of the total surface area of summer maize on 245 county-levels were gathered, including 9 counties in Beijing, 6 counties in Tianjin, 51 counties in Hebei, 99 counties in Henan, and 80 counties in Shandong.

Methodology
The overall process for the recognition of the summer maize cultivated area in Huanghuaihai Plain is displayed in Figure 3. The main process includes: (1) The support vector machine (SVM) algorithm was adopted to distinguish the spatial distribution of summer maize in the reference area. The summer maize MODIS EVI time series curve was acquired through masking multiple periods of MODIS EVI images and the distribution of summer maize in the reference area, which was smoothed with the Savizky-Golay method [52][53][54][55]. (2) Through the stepwise regression method to analyze the relationship between summer maize phenological observation data and environmental factors from 32 agrometeorological stations, multiple linear regression equations were constructed to simulate the summer maize phenology. The simulated phenology acted as the corrected parameter for the MODIS EVI time series curve to obtain the equations of the standard summer maize EVI time series curve in the study area. (3) By calculating the mean absolute distances (MAD) map between standard summer maize EVI time series and actual MODIS EVI time series images, and setting the appropriate thresholds with the p-tile method, the summer maize planting area was estimated in Huanghuaihai Plain.

Method of Maize Area Identification in the Reference Area
By surveying the phenological characteristics of the main crops in the reference area (Figure 4), it can be seen that summer maize was sowed in early and middle June, which presented as bare soil on remote sensing images. The sowing dates of peanut and cotton were between late April and early May. The planting area of peanut and cotton showed vegetation information on remote sensing images in mid-June. It is the best time to distinguish summer maize from other summer crops in mid-June. From July to September, the summer maize was in the period from seven-leaf to maturation, which showed vegetation information on the remote sensing images. In late October, summer maize had been harvested and winter wheat was sowing or had seedling emergence in the reference area, so the summer maize planting area displayed the bare soil spectrum feature on the remote sensing images. Hence, the images with good quality on 9 June, 13 September, and 31 October, 2012, represented different spectrum features at different growth stages of maize, which were used to extract the summer maize acreage in the reference area. The specific steps were as follows: (1). Combined with Google Earth, the training samples of road, vegetation, habitation, water body, and bare soil were selected from three periods of ETM+ data by visual interpretation. The SVM algorithm was used to classify the three images.

Multiple Linear Stepwise Regression
In processes of stepwise regression analysis, contributions of environmental factors to variances in the regression equation are taken as criteria [8,56]. Significant factors that affect summer maize phenology are introduced into the regression equation one by one, and factors with no significant influence on summer maize phenology are removed one by one. Stepwise regression analysis can not only establish "optimal" regression observations, but also overcome the multicollinearity to a certain extent [57][58][59]. Given the widely reported sensitivity of maize growth to surroundings, it seems important to consider the effects of various factors on the maize phenology in different regions. Therefore, in this research, we employed multiple linear stepwise regression to fit the relationship between summer maize phenology and environmental factors (i.e., longitude, latitude, annual average temperature, annual precipitation, and ≥10 °C accumulated temperature).

The Mean Absolute Distance
The similar measurement method is one of the main approaches to classify crop planting information based on the multi-temporal vegetation indices (VIs) [60][61][62]. The basic idea of the similar measurement method is to establish VI curves of different features, and then compare the similarity of the time series VI curves between the pixel to be divided and the referenced pixel, which can determine its type. In the current research, MAD [18] was selected as the indicator of similar measurement. By comparing the similarity between the actual MODIS EVI pixel and the standard summer maize EVI pixel, the MAD map was calculated in the study area. The specific formula of MAD is as follows: where ij d is the value of the MAD; mk is the number of the pixels of the template image in the k sequence, taking the pixel of each raster of the standard summer maize EVI as the template, and mk is 1; xik is the actual value of MODIS EVI in the k sequence; xjk is the standard EVI value of summer maize in the sequence; n is the total number of the time series (n = 9).

The p-Tile Algorithm
The p-tile algorithm [21,63], a classical way to select the threshold, was proposed by Doyle in 1962. The research objects were divided into the target mode Po and the background mode Pb, and the quantile p-value was obtained by the ratio of Po and Pb, which is equal to the prior probability. The suitable threshold value can be found according to the prior probability, such that the target mode was greater than the threshold value and the background mode was less than the threshold value. The specific formulae are as follows: where Po is the target mode, Pb is the background mode, f(x,y) is the research objects (MAD map), and T is the threshold. The ratio of the statistical summer maize acreage and the total area of every province was calculated as a prior probability in Huanghuaihai Plain. The histogram of the MAD map of each province was accumulated until the proportion of accumulated value was greater or equal to the prior probability, and this is when the mad value was the optimal threshold T.

Extraction and Reconstruction of Summer Maize MODIS EVI Time Series Curve in the Reference area
Based on ETM+ data, the spatial distribution of summer maize (38.5 × 10 3 ha) in Ruzhou is shown in Figure 5. The result of summer maize extraction in Ruzhou City was verified by the statistics of the total surface area of summer maize (43.4 × 10 3 ha), and the accuracy was higher than 88%. It was revealed that the way of intersection operation based on different spectral characteristics in maize critical periods can acquire the maize cultivated area with high precision in the relatively single area. The summer maize MODIS EVI time series images of the reference area were obtained by masking the spatial distribution area of summer maize and multi-temporal MODIS EVI images from the reference area. The average value of summer maize MODIS EVI images in each period was calculated, which was regarded as the summer maize EVI value to obtain the original MODIS EVI time series curve of summer maize in the reference area. Smoothing the original MODIS EVI time series curve of summer maize with Savitzky-Golay filtering, the smoothed MODIS EVI time series curve of summer maize was gained in the reference area ( Figure 6). The Gaussian function was used to fit the smoothed MODIS EVI time series curve of summer maize in the reference area. The formula is as follows: where i is the day of the year (DOY), EVIi is the simulated EVI value of summer maize in the reference area on DOY i, and the parameters of the simulated EVI function are shown in Table 2. It is thus clear from Figure 6 that the value of summer maize MODIS EVI increased first and then decreased with time in the reference area, consisting of the growth trajectory of summer maize. From DOY 157 to 225, summer maize was in the periods of sowing to tasseling, which was the rapid growth stage, and the EVI value gradually increased. It entered the period of tasseling around DOY 225, which was the most vigorous growth stage, and the value of EVI also reached the peak value. From DOY 225 to 258, summer maize was in the periods of tasseling to maturation, and gradually matured with the decrease in EVI value.

Estimation of Summer Maize Phenology
The summer maize phenology is divided into eight key stages, which are sowing, seedling emergence, trefoil, seven-leaf, jointing, tasseling, milk-ripe, and maturation. The phenology of summer maize in different regions varies due to the different natural environments. In this study, employing the multiple linear stepwise regression method, the phenological observation data and the corresponding meteorological data of 32 agrometeorological stations were analyzed, and the multiple linear regression equations between summer maize phenology and natural environmental factors were constructed. The specific formulas are as follows:  (13)) at 32 agro-meteorological stations, respectively. The correlation coefficient (R) values were between 0.559 and 0.790 (Figure 7), which showed that the correlation between the simulated phenology and the actual phenology of summer maize was significant. It was illustrated that the summer maize phenology simulated by multiple linear regression equations were reliable in the current study. According to the multiple linear regression equations, the simulated results of summer maize phenology in the study area were calculated and are shown in Figure 8. Generally, the eight growth stages of summer maize in Hunaghuaihai Plain showed a trend of delaying from the southwest to northeast. The differences in seedling emergence, seven-leaf, and joint stages were the smallest, and the DOY discrepancy among the regions was about 16 days in the study area, whereas the difference in the maturation phase of summer maize was largest, up to about 32 days. The occurrence of this event may be related to the climate with high temperature in summer and a larger temperature difference between North and South in winter in Huanghaihai Plain [64,65]. Based on the summer maize phenology of Ruzhou City (Table 1) and the simulated phenology of summer maize in the study area, the lag at each pixel from the sowing stage to maturation stage of maize was roughly calculated using the following equations. The corresponding values of i and s are shown in Table 3.
where i is the day of the year (DOY); s represents the time series of summer maize phenology; Lagi are the day's delay from DOY 161 to 289, and Ys and Yrs represent the simulated summer maize phenology of the study area and the phenology of summer maize in the reference area on DOY i, respectively. The growth stages of summer maize used in Equation (14) are the periods of seedling emergence, seven-leaf, jointing, tasseling, milk-ripe, and maturation. The phenology of summer maize is different in various regions on a large scale. The features of the summer maize MODIS EVI time series curve in the reference area cannot represent the MODIS EVI time series curve of the entire study area. It was necessary to modify the summer maize MODIS EVI time series according to the various analog phenologies of summer maize in Huanghuaihai Plain. Hence, in the current research, employing the lag value of summer maize phenology to correct Equation (2), the standard EVI time series equation of summer maize in the study area was obtained. The formula is as follows: ( 145,177 289) a e a e i (15) where SEVIi is the standard summer maize EVI value of the study area on DOY i, and the parameters of the standard EVI function are shown in Table 2.

Validation of Standard Summer Maize EVI Time Series Curve in the Study Area
Ten counties were randomly selected in the study area to validate the standard summer maize EVI time series curve. According to the method in Section 3.1, the planting area of summer maize in 10 validated regions was mapped by Landsat-7 ETM+ images. The maize distribution region extracted by Landsat was used to mask the MODIS EVI time images, and the actual MODIS EVI time series values of summer maize in 10 validated areas were obtained by calculation. The actual MODIS EVI time series values of summer maize were compared with the standard EVI time series values of summer maize by Equation (15). The results showed that the correction coefficient between the actual summer maize MODIS EVI time series curve and the standard summer maize EVI time series curve was all greater than 0.9 with p < 0.01 (Table 4). It indicated that the standard curve simulated by Equation (15) was consistent with the actual curve. Table 4. Correlation analysis between the actual MODIS EVI curve and the simulation standard EVI curve of maize.

The MAD between Standard Summer Maize EVI and Actual MODIS EVI
The MAD map of the study area is shown in Figure 9, where the minimum value is 0 and the maximum is 6. The lower the MAD value, the greater the similarity between the actual EVI value and the simulated EVI value, which indicated that the summer maize planting area is more likely. On the contrary, the summer maize planting area is less likely. According to the p-tile algorithm, the optimal thresholds of five provinces were determined: 0.85 for Beijing, 0.89 for Tianjin, 0.71 for Hebei, 0.59 for Henan, and 0.74 for Shandong. The value of MAD that was smaller than the optimal threshold was considered to be the summer maize area by comparing every raster value of the MAD map with the optimal threshold value, which was used to extract the spatial distribution of summer maize in Huanghuaihai Plain.

Extraction Results of Summer Maize in Huanghuaihai Plain
The spatial distribution of the summer maize planting area of the study area is shown in Figure  10. It can be seen from Figure 10 that the distribution of summer maize is sparse in Beijing, Tianjin, and the western and southern parts of Henan Province, and uniform in other regions of Huanghuaihai Plain. The event is owed to the agricultural planting structure, social development, and ecological environment in China. There is less area for agricultural planting in Beijing and Tian as the developed cities. Rice and wheat are the main crops in the south of Henan Province, and woodland covers most of the western Henan Province. Hebei, Shandong, and the central and northern parts of Henan Province are flat with good photoperiod and temperature conditions, which are suitable for maize growth [66,67]. The spatial distribution of summer maize in the current study was consistent with the previous extracted results [68] in northern and southern Hebei Province, northwest of Shandong Province, and northern Henan Province, which can manifest the extracted results with a certain reliability.

Validation of the Extracted Summer Maize Area in Huanghuaihai Plain
The cultivated areas of summer maize in Beijing, Tianjin, Hebei, Henan, and Shandong extracted by MODIS EVI time series images were 129.7 × 10 3 , 175.4 × 10 3 , 2995.4 × 10 3 , 3322.7 × 10 3 , and 3428.2 × 10 3 ha, respectively. The planting area of summer maize was largest in Henan Province and the smallest in Beijing. The statistics and the extraction of summer maize were used to analyze the overestimation and underestimation of area extraction. If the comparison results were positive, it indicated that the extraction area was larger than the statistical area, which was overestimation. Otherwise, if the comparison was negative, it was underestimated. As can be seen from Table 5, the extraction of summer maize was overestimated in Hebei, Henan, and Shandong, and underestimated in Beijing and Tianjin. The extraction accuracy of Tianjin Province was the highest, reaching 98.04%, and Shandong Province was the lowest, 86.41%. The accuracy of the extraction results of five Provinces was over 85% in the study. Note: High (+) or low (−)/% = (statistical area-extracted area)/statistical area×100, and accuracy (%) is 100-Abs (high or low).
The statistical data of the maize area from 43 cities and 245 counties in Huanghuaihai Plain were employed to verify the extraction results. Comparing the maize acreage of two groups from MODISestimated and statistical data, the validation results showed that the R 2 reached 0.81 with an Root Mean Square Error (RMSE) of 143.8 ha at the city level and the R 2 reached 0.69 with an RMSE of 18.7 ha at the county level ( Figure 11). Those demonstrated that there was a good correlation between the maize area estimated by MODIS satellite data and statistical data. Compared to the previous research that used the multi-temporal remote sensing images to map the summer maize area, the extraction accuracy was significantly improved in the current research [20,21,69]. Six county-level regions (Changyuan, Gaotang, Linzhang, Lushan, Yongnian, Xinji) were selected from the study area as the sampling validation areas. Based on Landsat ETM+ images with a spatial resolution of 30 m, the distribution regions of summer maize in the six regions in 2012 were obtained by the method used to estimate summer maize area in the reference area in Section 3.1. As shown in Figure 12, the summer maize distribution areas mapped by ETM+ images were roughly consistent with those drawn by MODIS images. The extraction results were compared to the summer maize distribution areas drawn based on the MODIS data. The two kinds of maize distribution areas were superposed and calculated to obtain the area of the public part. The public area is 29.4 ha in Changyuan, 37.7 ha in Gaotang, 36.6 ha in Linzhang, 24.7 ha in Lushan, 33.9 ha in Yongnian, and 43.2 ha in Xinji ( Table 6). The proportions of the public part to the extraction results by Landsat images, the extraction results by MODIS images, and the statistical data were used to quantitatively analyze the relative accuracy of summer maize extracted with MODIS data. The proportion of the public part area to the extraction results based on MODIS and Landsat images was compared. The results showed that the gap was largest in Lushan county, with 90.7% and 62.5%, respectively; the gap was smallest in Changyuan, 67.7% and 68.6%, respectively ( Table 6). The mean proportion of the public part area to the statistical data was 82.3%, which meets the request for precision. These suggest that the method proposed in this paper is reliable.

Extraction Results in Different Years
The area extraction algorithm in the present research was employed to map the summer maize cultivated area in 2013 over Huanghuaihai Plain. The spatial distribution characteristics of summer maize in 2013 were basically consistent with those in 2012. It can be seen from Table 7 that the summer maize acreage from MODIS images in Beijing, Tianjin, Hebei, Henan, and Shandong were 105.1 × 10 3 , 178.2 × 10 3 , 2630.4 × 10 3 , 3104.3 × 10 3 , and 2707.0 × 10 3 ha, respectively. Compared to the statistical area, the precision of summer maize extraction results in five provinces was all above 85%. Meanwhile, the extraction results of the summer maize area from 43 cities and 245 countries of the study area in 2013 were compared to the statistical data, which showed that R 2 was 0.81 with an RMSE of 132.6 × 10 3 ha at the city level and R 2 was 0.68 with RMSE of 17.3 × 10 3 ha at the county level ( Figure 13). It can be found that the precision of extraction results in different years had a high accuracy by using the standard summer maize EVI time series to map the summer maize, indicating that the method has a certain universality over a large scale.  Note: High (+) or low (−)/% = (statistical area-extracted area)/statistical area×100, and accuracy (%) is 100-Abs (high or low).

Conclusions
The present study demonstrates a promising approach for mapping a summer maize cultivated area over a large scale, which considered the phenological differences of maize among the different regions. In practicality, the relationship between summer maize calendars and environmental factors was constructed and environmental factor data were employed to simulate the summer maize phenology in the whole study area. Simulated phenology of summer maize acted as a parameter to obtain a standard EVI time series curve of summer maize. A MAD map was calculated through comparison between time series images of standard summer maize EVI and the actual MODIS-EVI. Appropriate thresholds were set in five different provinces to map spatial distributions of summer maize in the study area. This research has demonstrated that, in terms of precision analysis, the accuracy of the extraction result in each province was above 85%. Statistical data from the Statistical Yearbook were gathered for accuracy evaluation, and the results showed that the R 2 reached 0.81 with an RMSE of 143.8 ha at the city level and the R 2 reached 0.69 with an RMSE of 18.7 ha at the county level. Those indicated a robust potential for identifying areas where the summer maize was cultivated over a large scale. Of course, there are still many challenges in the estimation of the crop planting area while considering the complex agricultural structure, including crop types, management practices, and crop interplanting operation. Further studies can also attempt to use more high-spatial-resolution images to map the crop cultivated area and improve the accuracy of extraction results on a large scale.