Next Article in Journal
Thermal Contribution of the Local Climate Zone and Its Spatial Distribution Effect on Land Surface Temperature in Different Macroclimate Cities
Next Article in Special Issue
Assessment Analysis of Flood Susceptibility in Tropical Desert Area: A Case Study of Yemen
Previous Article in Journal
Global Multi-Scale Information Fusion for Multi-Class Object Counting in Remote Sensing Images
Previous Article in Special Issue
Evaluation of Drought Vulnerability of Maize and Influencing Factors in Songliao Plain Based on the SE-DEA-Tobit Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rapid and Automated Mapping of Crop Type in Jilin Province Using Historical Crop Labels and the Google Earth Engine

1
School of Environment, Northeast Normal University, Changchun 130024, China
2
Key Laboratory for Vegetation Ecology, Ministry of Education, Changchun 130024, China
3
State Environmental Protection Key Laboratory of Wetland Ecology and Vegetation Restoration, Northeast Normal University, Changchun 130024, China
4
Jilin Province Science and Technology Innovation Center of Agro-Meteorological Disaster Risk Assessment and Prevention, Northeast Normal University, Changchun 130024, China
5
College of Geographical Science, Inner Mongolia Normal University, Hohhot 010022, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(16), 4028; https://doi.org/10.3390/rs14164028
Submission received: 9 July 2022 / Revised: 7 August 2022 / Accepted: 16 August 2022 / Published: 18 August 2022

Abstract

:
In the context of climate change, the remote sensing identification of crops is extremely important for the rapid development of agricultural economy and the detailed assessment of the agro-meteorological disasters. The Jilin Province is the main grain production area in China, with a reputation of being a “golden corn belt”. The main crops in the Jilin Province are rice, corn, and soybean. A large amount of remote sensing data and programming codes from the Google Earth engine (GEE) platform allow for large-area farmland recognition. However, the substantial amount of crop sample information hinders the mapping of crop types over large farmland areas. To save costs and quickly and accurately map the crop types in a study area, multi-source remote sensing data and historical crop labels based on the GEE platform were used in this study, together with the random forest classification method and optimal feature selection to classify farming areas in the Jilin Province. The research steps were as follows: (1) select samples based on the historical crop layer of the farmland; and (2) obtain the classification characteristics of rice, corn, and soybean using multi-source remote sensing data, calculating the feature importance scores. Using different experimental combinations, an optimal classification method was then selected to classify crops in the Jilin Province. The results indicated variable impacts of vegetation indices (of different periods) on crop classification. The normalized difference vegetation index (NDVI), green normalized difference vegetation index (GNDVI), and green chlorophyll vegetation index (GCVI) in June exerted a significant impact on the classification of rice, corn, and soybean, respectively. The overall accuracy of crop classification during different periods based on historical crop labels reached 0.70, which is acceptable in crop classification research. The study results demonstrated that the proposed method has promising potential for mapping large-scale crop areas.

Graphical Abstract

1. Introduction

Market fluctuations, land degradation, pests, and climate change are causing increasing difficulties in agricultural practices and in maintaining population growth and health [1,2]. Frequent extreme events, such as high and low temperatures, floods, droughts, and frosts, urgently require the effective planning and management of agricultural land. The timely, accurate, macro, and dynamic extraction of crop planting areas can provide an important scientific decision-making base for government departments at all levels, as well as for farmers and promotes the scientific management of agricultural production, which is crucial for food security and economic development. Accurate informational extraction of crop planting areas can provide strong support for disseminating early warnings and predicting agrometeorological disasters under climate change conditions, and formulating crop planting plans [3,4,5].
Previous studies on the extraction of crop planting areas are based on statistical investigation methods, which require substantial human, material, and financial resources, and are vulnerable to human factors. Remote sensing technology provides an irreplaceable tool for the detailed monitoring of crops over a large area. The continuous updating of remote sensing platforms and the free opening of remote sensing data has ushered in an unprecedented era of remote sensing big data. In particular, the 2008 US National Address Survey (USGS) started allowing free access to all historical archived and real-time Landsat data. In addition, the launch of, and data sharing from, the China Resources Satellite and ESA Sentinel-2 satellite further supplemented the sources of medium spatial resolution data.
Currently, remote sensing technology is being increasingly applied to crop identification, yield estimation, and health monitoring, as well as cultivated land rotation, to provide reliable data to support agricultural management [6,7,8,9,10,11,12]. Over the past ten years, several studies have explored the use of remote sensing data to spatially extract crop planting areas and have presented distribution maps based on various remote sensing products at different spatiotemporal resolutions [13,14,15]. Based on the releases of Landsat 5 TM and 7 ETM+ data, the National Bureau of Agricultural Statistics (NASS) of the U.S. Department of Agriculture (USDA) has been annually releasing a farmland data layer (CDL) at a spatial resolution of 30/56 m since 2008 [13,14]. Combining detailed plot boundary information provided by the general land unit data of the Farm Service Agency (FSA) and a large amount of ground survey data, NASS has identified approximately 110 different types of crops in the United States. Agriculture and Agri-Food Canada (AAFC) has been releasing AAFC crop-type distribution data since 2011, using a decision tree algorithm and various optical (Landsat, resourceat-1, DMC, spot) and radar (radarsat-2) data, achieving an overall nationwide accuracy of 85% [15]. The remote-sensing-data products for crop spatial distribution in China are still relatively backward, and a regularly updated crop spatial distribution database is lacking, which undoubtedly creates certain difficulties for research in agrometeorology, agricultural management, agricultural sustainable development, etc. The Institute of Agricultural Resources and Agricultural Regionalization of the Chinese Academy of Agricultural Sciences, together with the Institute of Geographical Sciences and Resources of the Chinese Academy of Sciences, utilized data from various sources, including the Chinese Agricultural Resources and Regionalization database, 52 scientific research projects, and remote sensing data, and extracted the primary planting distribution areas of five major crops (rice, corn, barley, soybean, and cotton) in China from 2000 to 2015 at a spatial resolution of 1 km by combining area extraction and basic background data [16,17]. A researcher from the Institute of Geographical Sciences and Resources of the Chinese Academy of Sciences, in cooperation with the China Agricultural University, Northeast Agricultural University, and the University of Oklahoma, U.S., used available Sentinel-2 data to perform remote-sensing monitoring of crop planting areas and distributed a crop distribution product at 10 m resolution for northeast China from 2017 to 2019, thus providing the most detailed distribution information on the three major crops (rice, corn, and soybean) in the region [18].
Most crop classification research uses the vegetation features over the entire growing season as classification features [19,20,21]. However, doing this requires downloading and processing large amounts of remote sensing data. Therefore, researchers typically choose small-scale research areas to reduce the difficulties associated with large image processing. The development and adoption of remote sensing cloud computing has provided an effective solution for processing even large-scale remote sensing data. The Google Earth Engine (GEE) is a cloud-based platform. A large number of publicly available geospatial datasets are stored in this data directory, including observations from various satellites and aerial imaging systems in optical and non-optical wavelengths, environmental variables, weather and climate prediction and post-prediction, land cover, terrain, and socio-economic datasets [22]. It has large satellite images and rapid cloud-based geospatial processing capabilities. An increasing number of studies on crop classification based on the GEE platform are now available. For example, based on multi-source data from the GEE platform, corn, soybean, rice, garlic, winter wheat, sugarcane, and potato have been classified using supervised and unsupervised classification methods [23,24,25,26,27].
The importance of a ground classification feature directly affects classification accuracy. Calculating the importance of each classification feature and their optimal selection can improve the classification speed and reduce data redundancy [28,29]. Normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), soil-adjusted vegetation index (SAVI), and normalized difference water index (NDWI) are widely used in land cover classifications. These general vegetation indices are also applied to crop classification. To distinguish between different crop types, researchers have explored multiple classification features to distinguish subtle differences in the phenological characteristics of different crops, such as those available in the modified normalized difference water index (MNDWI), normalized difference snow index (NDSI), and green chlorophyll vegetation index (GCVI). Vegetation changes with phenology, and the vegetation index characteristics also have different effects on classification. Although the land surface water index (LSWI) has an important impact on the early classification of rice in northeast China, it has an insignificant effect on the classification of corn and soybean. The normalized difference soil vegetation index (NDSVI) has a significant impact on the late classification of rice; however, its impact on the entire growing season of corn and soybean is negligible [30,31,32,33,34,35].
Currently, most studies of data classify the crops of a particular year using crop training data from the same year. Although this method can improve the accuracy of classification, it requires substantial manpower and material resources. Therefore, to improve the utilization of historical crop data and the quick retrieval of crop classification data, the importance of different crop classification characteristics and the coupling of historical crop labels within different years were analyzed in this study to improve the accuracy of crop classification; the most robust classification model was then used to classify rice, corn, and soybean cultivation in the study area.

2. Materials and Methods

2.1. Study Area

The Jilin Province (121°38′–131°19′ E, 40°50′–46°19′ N) is located in central northeast China, with visible differences in geomorphic morphology across the entire province. The terrain generally inclines from the southeast to the northwest. There are several mountains and hills in the east and plains in the middle and west (Figure 1). The province has a temperate continental monsoon climate, with clear distinctions among the four seasons. The spring is dry with several peaks and summer is hot with reorganized precipitation. Autumn is cool and suitable, whereas winter is cold and long. The climate is suitable for crop growth [36,37]. Crops are planted once a year.
As an important grain production base in China, the Jilin Province is important for ensuring national food security and the steady development of agriculture. The Jilin Province is one of the three golden corn belts worldwide. The Jilin Province had a grain sown area of 5681.8 thousand hectares in 2020, accounting for approximately 4.86% of the national total (Statistical Yearbook, 2021); the grain output from the region was 38.0317 million tons, accounting for approximately 5.83% of the national output. Rice, corn, and soybean are the three most cultivated food crops, accounting for 89% of the total cultivated land in the Jilin Province. The sowing time of corn and soybean varies with the year and location, and largely depends on temperature conditions, soil moisture, and farm management decisions. The rice-growing season spans from early April to mid-September and generally lasts 120–195 days. The growth period of corn is generally 90–130 days from the end of April to the end of September, while that of soybean is from mid-May to mid-September, generally spanning 100–120 days. The data used in this study are shown in Figure 2.

2.2. Data Sources

2.2.1. Remote Sensing Data in GEE

The GEE is an open-source intelligent cloud platform, specifically used for the interpretation and operation of satellite images and other spatial data. It fully uses Google’s powerful storage capacity and advanced computing power and can quickly process large geospatial datasets. The GEE platform provides a centralized remote sensing image warehouse. These images were geo-referenced, atmospherically corrected, and converted to surface reflectance for this study. All image processing and analysis were performed using the GEE application programming interface [22]. The data were analyzed on the GEE platform. Therefore, resolution matching between different data sources was not required during the calculation process. Vector clipping of administrative regions was performed through online coding, batch cloud removal, and image cloud content screening.
The GEE cloud computing platform contains Landsat and MODIS images, including a thematic map (TM), enhanced thematic map (ETM+), operational land imager (OLI) and MODIS series images [22]. During image analysis, Landsat images in GEE were selected as a secondary product. This product was processed using radiometric correction and system-level geometric correction. It has satisfactory radiometric characteristics and is suitable for time-series analysis, with a spatial resolution of 30 m, a temporal resolution of 16 days, and 11 bands. MODIS/terra surface reflectance products (mod09q1) include two reflectance waves with a spatial resolution of 250 m and a mass layer, which have been corrected for atmospheric gases and aerosols. Surface reflectance was obtained after calibration between different Landsat instruments and atmospheric correction using auxiliary data from MODIS, such as water vapor, ozone, and aerosol optical depth. Surface reflectance provides a high-quality data source for agricultural remote sensing and crop-type mapping. Several studies have used Landsat and MODIS images to map land-use types and crop planting areas [27,38].

2.2.2. Cropland Data Layer

The cropland data layer in northeast China was obtained from You, N et al. [18]. Using Sentinel-2 data, a random forest (RF) algorithm, and a complex feature selection method; a 10 m-resolution GEE-based map was drawn for rice, corn, and soybean crop types in northeast China during 2017–2019, with an overall accuracy (OA) of 0.81–0.86. The correlation coefficient between the crop type map and the statistical yearbook data for 2017 and 2018 was higher than 0.83, and the R2 values of the corn area and statistical yearbook data were 0.98 and 0.99, respectively. The R2 of the soybean area and statistical data were 0.83 and 0.94, respectively. This dataset was used as a test sample to obtain the best classification combination by coupling samples from different years and classifying rice, corn, and soybean crops in the study area. As the accuracy of the test samples directly affects the results of data classification, the high classification accuracy of this dataset provided reliable data support for this study.

2.2.3. Validation Data

A total of 1269 farmland surface samples were collected from 2019, including 308 for rice, 477 for corn, and 484 for soybean (Figure 1). The location and crop type of each sample were recorded in the fields along a route using mobile GIS equipment (The seventh generation of iPad produced by Apple Company of the United States was equipped with a GIS software ovital map). Other land cover types (such as construction land, forests, and water bodies) were also recorded. The samples were used for accuracy evaluation.
In addition, the Jilin Provincial Bureau of Statistics releases an annual report on the planting area of major crops in each city and county. We obtained the statistical yearbook of planting area at the county and municipal levels from 2015 to 2020 and recorded the planting areas of rice, corn, and soybean in 2015 and 2020. The dataset was used to compare the results of crop type maps at the prefecture level.

2.2.4. Auxiliary Data

We used the national land cover database product to mask the non-cultivated pixels. The land cover database classifies various land cover types, such as water, development, and forests, including a single “cultivated crop” category that combines all crop types. The ground resolution of the grid was 30 m (http://www.globallandcover.com/, accessed on 7 February 2022).

2.3. Methods

This study combined Landsat and MODIS time-series remote sensing images and agricultural survey statistical data to perform crop type mapping and area estimation. The study was divided into two parts (Figure 3). In the first, we selected the pure pixel label of the crop type based on farmland data from 2017 to 2019. Simultaneously, training samples from different year combinations were generated. In the second part, the Landsat and MODIS data were preprocessed using the GEE platform, and the vegetation index was calculated. Considering the best spectral time characteristics as independent variables, an optimal RF model was obtained by selecting sample data from different periods. The sub-pixel crop components for each crop type in the entire study area were predicted. Finally, field survey samples and statistical data were used to evaluate the accuracy of the classification results.

2.3.1. Pure Pixel Sample Selection

To select pure pixels, the 25 × 25 moving window method was used to filter the 10 m-resolution farmland data from 2017 to 2019 and then systematically select specific pixels that were most likely to contain only a single crop type, making them eligible for use in model training. Pixels at a 250 m resolution are competent for model training only if (1) other pixels of the same crop type are surrounded by a 25 × 25 pixel moving window and (2) if 100% of the complete 10 m farmland data layer—called “pure pixels” is used. Through this screening process, the specific spectral characteristics of the different crop labels can be obtained. From 2017 to 2019, 6000 pixels were resampled annually, including 2000 pixels for rice, corn, and soybean, respectively, converted into point features, and spliced to create a training dataset. One drawback of this process is that “other” crop classes are oversampled, by 25 times, compared with any known crop type, which can introduce modeling bias, as shown in Figure 4.

2.3.2. Feature Selection

Different vegetation indices have different effects on crop classification. Vegetation feature extraction requires time and computational resources. Feature importance analysis can help select robust features, reduce the number of features to be extracted, and improve the classification algorithms. To explore the impact of different vegetation indexes on crop classification in the study area, the monthly average values of NDVI, EVI, GCVI, red-edge chlorophyll index (RECI), optimized soil-adjusted vegetation index (OSAVI), green normalized difference vegetation index (GNDVI), modified soil adjusted vegetation index 2 (MSAVI2), LSWI, NDWI, and NDSVI, for the growing season (April–October) were selected, with the importance of the characteristics calculated, and rice, corn, and soybean in the study area, classified. The calculation formulae and characteristic descriptions of the vegetation indices are presented in Table 1.

2.3.3. RF Classification Method and Feature Importance Selection

RF, which was proposed in 2001, combines bagging integrated learning theory and random sub-control methods, and provides unique advantages for processing high-dimensional data and building decision trees [12,39]. Previous studies have demonstrated that RF has stronger robustness, accuracy, stability, and accuracy than several traditional classifiers such as maximum likelihood, single decision tree, and single-layer neural networks [40]. The RF algorithm on the GEE platform has been successfully applied to land cover change detection, agricultural land monitoring, and crop-type classification. While training the crop classifiers, we adjusted two parameters of RF on GEE: (1) number of trees: the number of trees determines the number of binary cart trees used to build RF models. Notably, the accuracy improves slightly with an increasing number of trees, and the calculation cost increases linearly. The number of trees used in this study was set as 300. (2) Minleafpopulation: the minimum number of samples required for a leaf node. We set the minleafpopulation to 10 to limit the depth of each tree and avoid overfitting. The other four parameters, including variablespersplit (number of variables per split, characteristics of the square root of the number by default), bagfraction (part of the bag/input tree, 0.5 by default), outofbag mode (whether the classifier should run in out-of-bag (OBB) mode), and seed (random seed), were set as default. In the study, pure pixel samples from 2017 to 2019 were used as training samples, and the field survey data in 2019 were used as validation samples for random forest classification.
Two methods are generally used to calculate the importance score of RF characteristics, based on the reduction of the Gini washing impurity and the error rate of OBB data (Yin et al., 2020). In this study, the reduced value of Gini rinsing impurity for each feature, before and after node branching, was calculated and normalized as the feature importance score. Assuming that there are x classification features, we calculated the feature X j . The specific calculation process is as follows:
  • Gini coefficient calculation:
    G i n i m = k = 1 K k k P m k P m k = 1 k = 1 K P m k 2
    Here, K represents K categories and P m k   the proportion of category K in node m ;
  • Characteristic X j : The importance of J at node m , that is, the change in value of the Gini coefficient impurity before and after node m branching. The greater the change in value, the greater the characteristic X j can quickly divide the samples into different sets with higher purity, that is, the stronger the classification ability;
    F I M j m G i n i = G i n i m G i n i m l G i n i m r
    Here G i n i m _ l and G i n i m _ r , respectively, represent the purity of G i n i   coefficients of the left and right child nodes after branching;
  • If feature X j occurs M times in the decision tree t , then the feature importance of feature   X j in decision tree t is:
    F I M t j G i n i = m = 1 M F I M j m G i n i
  • Assuming that there are t trees in the RF, the characteristic importance of characteristic   X j in the RF is calculated as:
    F I M j G i n i = i = 1 n F I M i j G i n i
  • Finally, the importance of all features is normalized as the final feature importance score:
    F I M j G i n i = F I M j G i n i x = 1 X F I M x G i n i

2.3.4. Experimental Design

The primary problems described earlier in this paper are as follows: (1) what is the amount of precision that can be achieved in crop-type classification by integrating time-series Landsat data and machine learning methods? (2) Using historical crop labels to classify crop types in different years, how does the selection of training data from different years affect the classification accuracy?
To solve the above problems, the training and test data were grouped by different years to test the impact of the time sampling. Accuracy was assessed using independent test data. We determined the starting or ending year in the experiment. Specifically, (1) the crop types in 2019 were predicted using data from a single year. Because the sample data in this study pertained to only 2017, 2018, and 2019, the pure pixel sample data for the same years were used to predict the crop types in 2019, and the accuracy was verified using field survey data from 2019. (2) Using multi-year data to predict crop types in 2019, data from 2017 and 2018 were used to predict the crop types in 2019, and the pure pixel sample data for 2017–2019 were used to predict the crop types in 2019.

2.3.5. Accuracy Verification

For evaluating the accuracy of the experimental data, the confusion matrix was calculated using the field survey data for the verification and the historical experiment label as the sample data [41]. The formula can be described as follows:
Producer s   accuracy = n u m b e r   o f   c o r r e c t l y   c l a s s i f i e d   p i x e l s   o f   a   c r o p   t y p e s t o t a l   n u m b e r   o f   p i x e l s   i n   t h e   v a l i d a t i o n   d a t e × 100
User s   accuracy = n u m b e r   o f   c o r r e c t l y   c l a s s i f i e d   p i x e l s   o f   a   c r o p   t y p e s t o t a l   n u m b e r   o f   t h e   c r o p   t y p e   p i x e l s   i n   t h e   c l a s s i f i e d   m a p × 100
A comparative analysis of the classification results and statistical data for 2015 and 2020 was used to reveal the relative error.

3. Results

3.1. Feature Importance Analysis

Feature importance is shown in Figure 5. The importance scores of different characteristic variables reveal large differences in the importance of different characteristics. For the classification of rice, the NDVI value in June was the most central to the classification, followed by the NDWI value in July, along with the EVI and GCVI values for August. For the classification of corn, the GNDVI value in June was the most significant, followed by the GNDVI and GCVI values in June. The GNDVI value in May was also highly important, whereas the NDWI, RECI, and OSAVI feature values were of low importance to the entire classification of corn. Regarding the classification characteristics of soybean, the GNDVI, GCVI, and NDWI values in June were highly important, whereas the OSAVI value had an insignificant impact. From the importance maps of the three classification characteristics, it can be concluded that the classification importance of crops varies in different months. The characteristics observed for June are the most important for classifying rice, corn, and soybean. The feature values in July had an impact on the classification of rice and corn, and the feature values in April and October had a negligible impact on the classification of the three crops. Farming in the Jilin Province starts in late April and early May or in late September and early October. The importance of classification features thus directly affects the accuracy of the crop classification.

3.2. Impact of Time Sampling on Classification

To fully utilize the existing crop type labels and reduce the crop classification cost, the crop types in the study area were identified according to the experimental scheme described in Section 2.3.4. Using the current year data to train the model and using field sampling data for verification, we verified the accuracy of the classification results and evaluated them. The experimental results are presented in Figure 6 and Table 2. When only the historical crop labels of 2017 and 2018 were used for classification, the classification user accuracy for these two years was calculated as 0.75 and 0.76 for rice, 0.64 and 0.68 for corn, and the lowest 0.56 and 0.59 for soybean, respectively. The classification producer’s accuracy was calculated as 0.62 and 0.63 for rice, 0.59 and 0.56 for corn, and the lowest 0.40 and 0.42 for soybean, respectively. When the coupling crop labels of 2017 and 2018 were used to classify the crops in 2019, the user’s accuracy improved significantly (0.84, 0.74, and 0.69 for rice, corn, and soybean, respectively), and the producer’s accuracy was 0.79, 0.69, and 0.65 for rice, corn, and soybean, respectively. When only 2019 crop labels were used, the classification accuracies for rice, corn, and soybean improved further to 0.90, 0.82, and 0.70, respectively. When crop labels for 2017–2019 were used to classify the 2019 crop types, classification accuracies of 0.91, 0.83, and 0.69 were obtained for rice, corn, and soybean, respectively. The results indicate that the accuracy is low when only one-year historical crop notes are used for classification, and coupling crop labels for several years can effectively improve the results. The classification accuracy was higher when the crop label training model of the current year was used to classify crops, especially for rice and corn. Although using historical and current crop labels as training samples slightly increased the classification accuracy, the improvement was not significant. However, when only historical crop labels are used to classify crops in the following year, coupling multi-year historical labels can improve the results.

3.3. Spatial Distribution of Crops in Jilin Province

The spatial distribution of rice, corn, and soybean in the Jilin Province in 2015 and 2020 is illustrated in Figure 7 using the sample data and classification characteristic index of the multi-year historical crop memo combination obtained using the experimental method. Because the study area is a part of the golden corn belt region, corn is abundantly distributed across it, whereas soybean and rice are concentrated in the central and western regions. In 2015, soybean planting was mainly distributed in the north-central and northwest of the study area with a small and scattered distribution; corn planting was wide and concentrated, and rice planting was mainly on the river banks. Comparatively, in 2020, the cultivated land area decreased, and the planting areas of soybean and rice increased. Soybean was planted mainly in the middle of the Jilin Province. Compared with 2015, its distribution was concentrated, and the planting area of rice increased significantly. In 2015, rice was mainly distributed in Changchun, Baicheng, and Jilin, with small and scattered planting areas in other areas. Three urban areas with the highest corn planting were the cities of Changchun, Siping, and Songyuan. Baishan City had the smallest corn planting area, and the lowest proportion of cultivated land area in the study area. Soybean planting was low in the entire study area. In 2020, rice planting in the study area was mainly limited to the cities of Baicheng, Changchun, and Jilin. Planting areas in other locations were scattered and small. Corn was the largest crop in the study area, with the planting area mainly distributed east of the study area. The cities of Changchun, Songyuan, and Jilin had the largest planting areas. The corn planting area in Changchun City increased compared with that in 2015, whereas that in SiPing City decreased; however, the corn planting area displayed an overall increasing trend. Soybean had less planting area in the study area than that of rice and corn. The largest soybean planting area was observed in the Yanbian Korean Autonomous Prefecture. The planting areas in other locations were smaller and scattered. The relative errors between the results of the crop area extracted using remote sensing and statistical data in 2015 and 2020 are presented in Table 3 and Table 4, and the scatter diagrams of the extracted and statistical areas are presented in Figure 8 and Figure 9.
As illustrated in Figure 8 and Figure 9 and Table 3 and Table 4, the relative error between the point far from the 1:1 line and the statistical value is large. According to the 2015 data, the relative error was the lowest for rice (0.02–0.37), between 0.08 and 0.39 for corn, and higher for soybean (0.02–0.043). Comparatively, based on the extracted and statistical data in 2020, the relative error was the smallest for rice (0.04–0.33), between 0.17 and 0.36 for corn, and the highest for soybean (0.06–0.55). The classification data for 2015 and 2020 revealed the highest classification accuracy for rice, followed by that for corn, and the lowest for soybean.

4. Discussion

4.1. Select Multiple Vegetation Index

Using Landsat and MODIS data, 70 eigenvalues of 10 vegetation indices were obtained and calculated using different spectral channels, which fully captured the subtle differences in the spectral characteristics of rice, corn, and soybean during the growing season. If we want to use time-series remote-sensing imagery, multi-source data are able to complement each other compared to using one type of remote-sensing data alone. When the cloudiness of a certain image cannot meet the research needs, other images can be used as a substitute to continue the research [8,10]. The results indicated that using the comprehensive vegetation index and individual vegetation index or spectral characteristics provides greater conduciveness to crop classification, which is consistent with previous studies [28,38]. In addition, information on key time periods is extremely important for crop classification. The eigenvalue score was the highest in June, and the spectral characteristics were the most sensitive to crop classification. In this study, the importance score of progressive forest features was used to explain the classification contribution of different features during different periods. Our research demonstrates that selecting the best features based on the importance of RF-derived features can improve mapping accuracy and computational efficiency.

4.2. Availability of Remote Sensing Data

The GEE platform has substantial satellite data storage, computing power, implementation of various machine learning algorithms, and a flexible user interface, which is the core of this research. Not only can we sample and export contaminated pixels easily from the carefully planned image set, but we can also reuse the light decision tree trained offline to classify clouds, haze, and shadows. The scalability of GEE was also used to calculate different vegetation indices in pixel growth seasons, and the RF model was applied.
The availability analysis of Landsat data in the study area for 2015 and 2020 is presented in Figure 10. The available image map for 2015 was higher than that of 2020, which is consistent with the crop classification accuracies for 2015 and 2020. If the availability of remote sensing images is high, the accuracy of crop classification can be improved. Sentinel-2 images have a time resolution of 10 days and a spatial resolution of 10 m, with high data availability. In order to improve the accuracy of the classification, researchers used Sentinel-2 to classify crops [7,8]. However, the possible year of Sentinel data is short. In order to further apply crop classification data in the later stage, Landsat and MODIS data were selected in this study. A small number of Landsat images were available for Baishan City and Yanbian Korean Autonomous Prefecture, in the east of the study area, along with a small number of cloudless remote sensing images in June and August for this area. Although June and August are two important months for crop classification, a reduced number of available images lead to an overestimated extraction area for rice and beans in this area, increasing the extraction error.

4.3. Spectral Characteristics of Different Crops

In remote sensing, identical ground object types display the same or similar spectral and spatial information characteristics under equivalent conditions (texture, terrain, illumination, vegetation coverage, etc.); however, spectral confusion can occur over a certain period of time, such as that among different crop types like cotton, corn, soybean, and millet. The similarity of spectral features poses a challenge to the accuracy of crop classification. As the planting and harvest times, as well as the growth rate of corn and soybean are similar, it is additionally difficult to accurately classify them. The spectral characteristics of corn and soybean are similar, and the spectral curves overlap in the growing season, which is consistent with previous studies [38]. The spectral reflectance of corn and soybean can overlap, as illustrated in Figure 11. In the figure, the time-series feature values of EVI and LSWI are similar for soybean and corn, and the LSWI of rice is significantly different from that of corn and soybean, which also confirms that the classification accuracy is higher for rice than that for corn and soybean.
The growth rates of different crops are different, and the vegetation characteristics in different months have different effects on crop classification. The farming system in the study area had one crop per year. The transplanting time of rice is in the middle of May and the harvest time is at the end of September. The planting time of corn and soybean is similar. It is planted in the middle of April and harvested at the end of September. Research shows that the vegetation characteristics in June have the greatest impact on crop classification. In June, different vegetation entered the rapid growth season, and corn entered the elongating stage in mid-June. Leaf area and leaf number increased rapidly. The leaf area and leaf number for soybean and rice were not as high as for corn. The difference in chlorophyll content led to NDVI being the most important for the vegetation classification in June. The second is NDWI in July. Rice enters the heading stage in July, and irrigation is needed to ensure rice yield. Therefore, the index has a high score of classification importance for rice. There are many available classified vegetation indexes for rice, and the importance scores of different vegetation indexes from May to September are higher, so the classification result of rice is the highest. The second is corn. The important months of corn classification are May, June, July, and September. The minimum available month for classification correlates to soybeans, with the most obvious being June. The vegetation index scores of other months are relatively low, which also leads to the lowest classification accuracy of soybeans. The historical crop labels from 2017 to 2019 are used in this study, and the accuracy of soybean classification in this data set is the lowest, which is also the direct reason for the low accuracy of soybean classification and the fewer vegetation features available for classification.

5. Conclusions

Against the background of climate change, the accurate mapping of crop types in large-scale farmlands is highly significant in ensuring food security and improving the agricultural economy. In this study, the use of time-series data from remote sensing images (Landsat and MODIS) and soybean, corn, and rice farmland maps of northeast China from 2017 to 2019, used as the training sample, were utilized to draw a planting area map of the Jilin Province using the RF method. An experimental analysis was performed to determine the most useful featured information for classifying rice, corn, and soybean and determine ways to improve crop classification based on historical crop labels. The research method mainly included two parts: (1) pure pixel extraction, based on a historical farmland data layer, and (2) the determination of important featured information for rice, corn, and soybean classification using field survey and farmland data layers and the RF feature. The crop types were then mapped based on remote sensing data from different years and training space samples. This study demonstrates that a combination of multi-year training samples can effectively improve the crop classification of rice, corn, and soybean. The relative errors in the area and statistical data for rice, corn, and soybean in the Jilin Province from 2015–2020 ranged between 0.02 and 0.56. The results indicate that this method can be used for large-scale crop classification.

Author Contributions

Conceptualization, F.Z. and Z.D.; methodology, F.Z.; software, F.Z.; validation, F.Z. and J.Z.; formal analysis, F.Z.; investigation, F.Z. and Y.B. (Yongbin Bao); resources, S.G. and A.H.; data curation, F.Z.; writing—original draft preparation, F.Z.; writing—review and editing, Y.B. (Yulong Bao); visualization, F.Z.; supervision, J.Z.; project administration, J.Z.; funding acquisition, J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

National Key R&D Program of China (2019YFD1002201); The National Natural Science Foundation of China (U21A2040); The National Natural Science Foundation of China (41877520); The National Natural Science Foundation of China (42077443); Industrial technology research and development Project of Development and Reform Commission of the Jilin Province (2021C044-5); The Key Scientific and Technology Research and Development Program of the Jilin Province (20200403065SF); The Key Scientific and Technology Research and Development Program of the Jilin Province (20210203153SF); The Science and Technology Development Planning of the Jilin Province (20190303018SF); The Construction Project of Science and Technology Innovation Center (20210502008ZP); The Changchun Key Scientific and Technology Program (19SS007).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nicholson, C.C.; Emery, B.F.; Niles, M.T. Global relationships between crop diversity and nutritional stability. Nat. Commun. 2021, 12, 5310. [Google Scholar] [CrossRef] [PubMed]
  2. Becker-Reshef, I.; Justice, C.; Sullivan, M.; Vermote, E.; Tucker, C.; Anyamba, A.; Small, J.; Pak, E.; Masuoka, E.; Schmaltz, J.; et al. Monitoring global croplands with coarse resolution earth observations: The global agriculture monitoring (glam) project. Remote Sens. 2010, 2, 1589–1609. [Google Scholar] [CrossRef]
  3. Ferreira, T.R.; Silva, B.B.D.; Moura, M.S.B.D.; Verhoef, A.; Nóbrega, R.L.B. The use of remote sensing for reliable estimation of net radiation and its components: A case study for contrasting land covers in an agricultural hotspot of the brazilian semiarid region. Agr. For. Meteorol. 2020, 291, 108052. [Google Scholar] [CrossRef]
  4. Fick, S.E.; Hijmans, R.J. Worldclim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  5. Ajadi, O.A.; Barr, J.; Liang, S.-Z.; Ferreira, R.; Kumpatla, S.P.; Patel, R.; Swatantran, A. Large-scale crop type and crop area mapping across brazil using synthetic aperture radar and optical imagery. Int. J. Appl. Earth Obs. 2021, 97, 102294. [Google Scholar] [CrossRef]
  6. Ibrahim, E.S.; Rufin, P.; Nill, L.; Kamali, B.; Nendel, C.; Hostert, P. Mapping crop types and cropping systems in nigeria with Sentinel-2 imagery. Remote Sens. 2021, 13, 3523. [Google Scholar] [CrossRef]
  7. Griffiths, P.; Nendel, C.; Hostert, P. Intra-annual reflectance composites from Sentinel-2 and Landsat for national-scale crop and land cover mapping. Remote Sens. Environ. 2019, 220, 135–151. [Google Scholar] [CrossRef]
  8. Prins, A.J.; Van Niekerk, A. Crop type mapping using LiDAR, Sentinel-2 and aerial imagery with machine learning algorithms. Geo-Spat. Inf. Sci. 2021, 24, 215–227. [Google Scholar] [CrossRef]
  9. Kpienbaareh, D.; Sun, X.; Wang, J.; Luginaah, I.; Bezner Kerr, R.; Lupafya, E.; Dakishoni, L. Crop Type and Land Cover Mapping in Northern Malawi Using the Integration of Sentinel-1, Sentinel-2, and PlanetScope Satellite Data. Remote Sens. 2021, 13, 700. [Google Scholar] [CrossRef]
  10. Blickensdörfer, L.; Schwieder, M.; Pflugmacher, D.; Nendel, C.; Erasmi, S.; Hostert, P. Mapping of Crop Types and Crop Sequences with Combined Time Series of Sentinel-1, Sentinel-2 and Landsat 8 Data for Germany. Remote Sens. Environ. 2012, 269, 112831. [Google Scholar] [CrossRef]
  11. Rustowicz, R.; Cheong, R.; Wang, L.; Ermon, S.; Burke, M.; Lobell, D. Semantic Segmentation of Crop Type in Africa: A Novel Dataset and Analysis of Deep Learning Methods. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR) Workshops, Long Beach, CA, USA, 16–20 June 2019. [Google Scholar]
  12. Mazarire, T.T.; Ratshiedana, P.E.; Nyamugama, A.; Adam, E.; Chirima, G. Exploring machine learning algorithms for mapping crop types in a heterogeneous agriculture landscape using Sentinel-2 data. A case study of Free State Province, South Africa. S. Afr. J. Geomat. 2020, 9, 333–347. [Google Scholar] [CrossRef]
  13. Boryan, C.; Yang, Z.W.; Di, L.P. Deriving 2011 cultivated land cover data sets using usda national agricultural statistics service historic cropland data layers. In Proceedings of the 2012 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Munich, Germany, 22–27 July 2012; pp. 6297–6300. [Google Scholar] [CrossRef]
  14. Boryan, C.G.; Yang, Z.W. Implementation of a new automatic stratification method using geospatial cropland data layers in Nass area frame construction. In Proceedings of the 2014 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Quebec City, QC, Canada, 13–18 July 2014; pp. 2110–2113. [Google Scholar] [CrossRef]
  15. Fisette, T.; Davidson, A.; Daneshfar, B.; Rollin, P.; Aly, Z.; Campbell, L. Annual space-based crop inventory for Canada. In Proceedings of the 2014 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Quebec City, QC, Canada, 13–18 July 2014; pp. 5095–5098. [Google Scholar] [CrossRef]
  16. Luo, Y.; Zhang, Z.; Chen, Y.; Li, Z.; Tao, F. ChinaCropPhen1km: A high-resolution crop phenological dataset for three staple crops in China during 2000–2015 based on leaf area index (LAI) products. Earth Syst. Sci. Data 2020, 12, 197–214. [Google Scholar] [CrossRef]
  17. Luo, Y.; Zhang, Z.; Li, Z.; Chen, Y.; Zhang, L.; Cao, J.; Tao, F. Identifying the spatiotemporal changes of annual harvesting areas for three staple crops in China by integrating multi-data sources. Environ. Res. Lett. 2020, 15, 074003. [Google Scholar] [CrossRef]
  18. You, N.; Dong, J.; Huang, J.; Du, G.; Zhang, G.; He, Y.; Yang, T.; Di, Y.; Xiao, X. The 10-m crop type maps in Northeast China during 2017–2019. Sci. Data 2021, 8, 41. [Google Scholar] [CrossRef] [PubMed]
  19. Peña-Barragán, J.M.; Ngugi, M.K.; Plant, R.E.; Six, J. Object-based crop identification using multiple vegetation indices, textural features and crop phenology. Remote Sens. Environ. 2011, 115, 1301–1316. [Google Scholar] [CrossRef]
  20. Adrian, J.; Sagan, V.; Maimaitijiang, M. Sentinel sar-optical fusion for crop type mapping using deep learning and google earth engine. ISPRS. J. Photogramm. 2021, 175, 215–235. [Google Scholar] [CrossRef]
  21. Belgiu, M.; Csillik, O. Sentinel-2 cropland mapping using pixel-based and object-based time-weighted dynamic time warping analysis. Remote Sens. Environ. 2018, 204, 509–523. [Google Scholar] [CrossRef]
  22. Gorelick, N.; Hancher, M.; Dixon, M.; Ilyushchenko, S.; Thau, D.; Moore, R. Google earth engine: Planetary-scale geospatial analysis for everyone. Remote Sens. Environ. 2017, 202, 18–27. [Google Scholar] [CrossRef]
  23. Phalke, A.R.; Özdoğan, M.; Thenkabail, P.S.; Erickson, T.; Gorelick, N.; Yadav, K.; Congalton, R.G. Mapping croplands of europe, middle east, russia, and central asia using landsat, random forest, and google earth engine. ISPRS J. Photogramm. 2020, 167, 104–122. [Google Scholar] [CrossRef]
  24. You, N.; Dong, J. Examining earliest identifiable timing of crops using all available Sentinel 1/2 imagery and google earth engine. ISPRS J. Photogramm. 2020, 161, 109–123. [Google Scholar] [CrossRef]
  25. Wang, J.; Xiao, X.; Liu, L.; Wu, X.; Qin, Y.; Steiner, J.L.; Dong, J. Mapping sugarcane plantation dynamics in Guangxi, China, by time series Sentinel-1, Sentinel-2 and landsat images. Remote Sens. Environ. 2020, 247, 111951. [Google Scholar] [CrossRef]
  26. Wuyun, D.; Sun, L.; Sun, Z.; Chen, Z.; Hou, A.; Teixeira Crusiol, L.G.; Reymondin, L.; Chen, R.; Zhao, H. Mapping fallow fields using Sentinel-1 and Sentinel-2 archives over farming-pastoral ecotone of Northern China with Google Earth Engine. GIsci. Remote Sens. 2022, 59, 333–353. [Google Scholar] [CrossRef]
  27. Tian, H.; Pei, J.; Huang, J.; Li, X.; Wang, J.; Zhou, B.; Qin, Y.; Wang, L. Garlic and winter wheat identification based on active and passive satellite imagery and the Google Earth Engine in Northern China. Remote Sens. 2020, 12, 3539. [Google Scholar] [CrossRef]
  28. Yin, L.; You, N.; Zhang, G.; Huang, J.; Dong, J. Optimizing feature selection of individual crop types for improved crop mapping. Remote Sens. 2020, 12, 162. [Google Scholar] [CrossRef]
  29. Lin, C.; Zhong, L.; Song, X.-P.; Dong, J.; Lobell, D.B.; Jin, Z. Early- and in-season crop type mapping without current-year ground truth: Generating labels from historical information via a topology-based approach. Remote Sens. Environ. 2022, 274, 112994. [Google Scholar] [CrossRef]
  30. Lambert, M.-J.; Traoré, P.C.S.; Blaes, X.; Baret, P.; Defourny, P. Estimating smallholder crops production at village level from Sentinel-2 time series in Mali’s cotton belt. Remote Sens. Environ. 2018, 216, 647–657. [Google Scholar] [CrossRef]
  31. Kwak, G.-H.; Park, N.-W. Impact of texture information on crop classification with machine learning and uav images. Appl. Sci. 2019, 9, 643. [Google Scholar] [CrossRef]
  32. Boschetti, M.; Busetto, L.; Manfron, G.; Laborte, A.; Asilo, S.; Pazhanivelan, S.; Nelson, A. Phenorice: A method for automatic extraction of spatio-temporal information on rice crops using satellite data time series. Remote Sens. Environ. 2017, 194, 347–365. [Google Scholar] [CrossRef]
  33. Wang, S.; Di Tommaso, S.; Deines, J.M.; Lobell, D.B. Mapping twenty years of corn and soybean across the us midwest using the landsat archive. Sci. Data 2020, 7, 307. [Google Scholar] [CrossRef] [PubMed]
  34. Xu, J.; Zhu, Y.; Zhong, R.; Lin, Z.; Xu, J.; Jiang, H.; Huang, J.; Li, H.; Lin, T. Deepcropmapping: A multi-temporal deep learning approach with improved spatial generalizability for dynamic corn and soybean mapping. Remote Sens. Environ. 2020, 247, 111946. [Google Scholar] [CrossRef]
  35. Veloso, A.; Mermoz, S.; Bouvet, A.; Le Toan, T.; Planells, M.; Dejoux, J.-F.; Ceschia, E. Understanding the temporal behavior of crops using Sentinel-1 and Sentinel-2-like data for agricultural applications. Remote Sens. Environ. 2017, 199, 415–426. [Google Scholar] [CrossRef]
  36. Li, L. Study on Climatic Suitability and Frost Damage of Soybean in Jilin Province; Nanjing University of Information Engineering: Nanjing, China, 2021. [Google Scholar] [CrossRef]
  37. Ma, Y.; Guga, S.; Xu, J.; Liu, X.; Tong, Z.; Zhang, J. Assessment of maize drought risk in midwestern jilin province: A comparative analysis of topsis and vikor models. Remote Sens. 2022, 14, 2399. [Google Scholar] [CrossRef]
  38. Cai, Y.P.; Guan, K.Y.; Peng, J.; Wang, S.W.; Christopher, S.; Brian, W.; Li, Z. A High-Performance and in-Season Classification System of Field-Level Crop Types Using Time-Series Landsat Data and a Machine Learning Approach. Remote Sens. Environ. 2018, 210, 35–47. [Google Scholar] [CrossRef]
  39. Rodriguez-Galiano, V.F.; Chica-Olmo, M.; Abarca-Hernandez, F.; Atkinson, P.M.; Jeganathan, C. Random forest classification of mediterranean land cover using multi-seasonal imagery and multi-seasonal texture. Remote Sens. Environ. 2012, 121, 93–107. [Google Scholar] [CrossRef]
  40. Zhong, L.; Gong, P.; Biging, G.S. Efficient corn and soybean mapping with temporal extendability: A multi-year experiment using landsat imagery. Remote Sens. Environ. 2014, 140, 1–13. [Google Scholar] [CrossRef]
  41. Dahal, D.; Wylie, B.; Howard, D. Rapid crop cover mapping for the conterminous united states. Sci. Rep. 2018, 8, 8631. [Google Scholar] [CrossRef]
Figure 1. Location map of the study area.
Figure 1. Location map of the study area.
Remotesensing 14 04028 g001
Figure 2. Data framework.
Figure 2. Data framework.
Remotesensing 14 04028 g002
Figure 3. The workflow for mapping crop-type fields by time series Landsat and Modis data.
Figure 3. The workflow for mapping crop-type fields by time series Landsat and Modis data.
Remotesensing 14 04028 g003
Figure 4. Pure pixel crop label generation.
Figure 4. Pure pixel crop label generation.
Remotesensing 14 04028 g004
Figure 5. Feature importance.
Figure 5. Feature importance.
Remotesensing 14 04028 g005
Figure 6. Classification performance changes with time; (a) 2017 data vs. 2019 classification results; (b) 2018 data on 2019 crop classification results; (c) 2019 data on 2019 crop classification results; (d) coupling 2017 and 2018 data to 2019 crop classification results; (e) classification results of 2017, 2018 and 2019 data for 2019.
Figure 6. Classification performance changes with time; (a) 2017 data vs. 2019 classification results; (b) 2018 data on 2019 crop classification results; (c) 2019 data on 2019 crop classification results; (d) coupling 2017 and 2018 data to 2019 crop classification results; (e) classification results of 2017, 2018 and 2019 data for 2019.
Remotesensing 14 04028 g006
Figure 7. Distribution map of the main crops in the study area.
Figure 7. Distribution map of the main crops in the study area.
Remotesensing 14 04028 g007
Figure 8. Comparison of the estimated area and the statistical data for rice, corn, and soybean in the study area in 2015.
Figure 8. Comparison of the estimated area and the statistical data for rice, corn, and soybean in the study area in 2015.
Remotesensing 14 04028 g008
Figure 9. Comparison of the estimated area and the statistical data for rice, corn and soybean in the study area in 2020.
Figure 9. Comparison of the estimated area and the statistical data for rice, corn and soybean in the study area in 2020.
Remotesensing 14 04028 g009
Figure 10. Available quantity of Landsat images. (a): Available quantity of Landsat images in 2015; (b): Available quantity of Landsat images in 2020.
Figure 10. Available quantity of Landsat images. (a): Available quantity of Landsat images in 2015; (b): Available quantity of Landsat images in 2020.
Remotesensing 14 04028 g010
Figure 11. Time series vegetation index characteristics of rice, corn and soybean.
Figure 11. Time series vegetation index characteristics of rice, corn and soybean.
Remotesensing 14 04028 g011
Table 1. The calculation formulae and characteristic descriptions of the vegetation indices.
Table 1. The calculation formulae and characteristic descriptions of the vegetation indices.
FeaturesCalculation FormulaDescription of Features
NDVI(NIR/RED)/(NIR + RED)Vegetation has strong reflection in the near-infrared band, with high reflectivity, while it has strong absorption in the red band, with low reflectivity, which can reflect the health of vegetation and the growth of vegetation. It is widely used in agriculture, forestry, ecological environment, and other fields, but it is very sensitive to soil brightness and atmospheric impact.
RECI(NIR/RED) − 1Because the chlorophyll content directly depends on the nitrogen content in the plant, this index is sensitive to the chlorophyll content in the leaves nourished by nitrogen, so it is helpful to detect the yellow or deciduous areas of the plant.
MSAVI(2× NIR + 1 − sqrt((2 × NIR + 1)2 – 8 × (NIR − RED)))/2It is applicable when NDVI cannot provide accurate values, especially in areas with a high proportion of bare soil, sparse vegetation, or low chlorophyll content in plants. This index is useful at the beginning of the crop growth season when seedlings begin to grow.
GNDVI(NIR − GREEN)/(NIR + GREEN)This index is a modification of NDVI, which is sensitive to withered or aging crops, monitors the nitrogen content in leaves, and is sensitive to dense canopy or mature vegetation.
GDSVI(SWIR1 − RED)/(SWIR1 + RED)This index is sensitive to the aging and yellowing of vegetation and can be used to distinguish the characteristics of different crops in different growth seasons.
NDWI(GREEN − NIR)/(GREEN + NIR)The difference ratio between green light band and near-infrared band is used to enhance the water information and weaken the information of vegetation, soil, buildings, and other ground features. This index has great advantages in pure water extraction and is widely used in farmland inundation and wetland feature extraction.
OSAVI(NIR − RED)/(NIR + RED + 0.16)It is applicable when the canopy coverage is low and has better sensitivity to the canopy coverage of more than 50%.
EVI2.5 × (NIR − RED)/((NIR) + (6 × RED) − (7.5× BULE) + 1)Since NDVI is easily disturbed by soil background and atmosphere, EVI adapts to atmospheric and soil noise, especially in the vegetation area of cats, reducing saturation. Used to analyze dense vegetation areas with large amounts of chlorophyll.
GCVINIR/GREEN − 1It is used to estimate the chlorophyll content in various plants, which reflects the physiological state of vegetation.
LSWI(NIR − SWIR1)/(NIR + SWIR1)Surface water index, which can characterize the change of soil moisture.
Table 2. Classification producer’s accuracy and user’s accuracy of samples from different years.
Table 2. Classification producer’s accuracy and user’s accuracy of samples from different years.
SampleRiceCornSoybean
Producer’s
Accuracy
User’s
Accuracy
Producer’s
Accuracy
User’s
Accuracy
Producer’s
Accuracy
User’s
Accuracy
20170.62590.75290.59240.64820.40120.5676
20180.63150.76190.56150.68580.42620.5937
20190.82540.90410.78540.85020.6980.7528
2017 + 20180.79260.842560.69520.748280.65420.7014
2017 + 2018 + 20190.83410.915480.76210.829470.6340.69291
Table 3. Relative error of crop extraction results in 2015 (104 ha).
Table 3. Relative error of crop extraction results in 2015 (104 ha).
City Rice CornSoybean
Statistical AreaEstimated AreaRelative ErrorStatistical AreaEstimated AreaRelative ErrorStatistical AreaEstimated AreaRelative Error
ChangChun16.8316.47−0.02105.11113.840.0816.4716.13−0.02
JiLin13.3312.62−0.0547.4256.320.1615.629.34−0.40
SiPing4.933.82−0.2380.9375.68−0.073.825.290.39
LiaoYuan1.991.32−0.3421.5617.34−0.241.320.75−0.43
TongHua9.067.13−0.2119.0424.080.216.134.22−0.31
BaiShan0.100.05−0.503.435.600.393.792.46−0.35
SongYuan11.0313.780.2582.3864.13−0.2813.7818.550.35
BaiCheng14.8915.130.0252.9042.23−0.2518.1310.60−0.42
YanBian3.684.450.2120.8224.870.1627.4519.26−0.30
Table 4. Relative error of crop extraction results in 2020 (104 ha).
Table 4. Relative error of crop extraction results in 2020 (104 ha).
CityRiceCornSoybean
Statistical AreaEstimated AreaRelative ErrorStatistical AreaEstimated AreaRelative ErrorStatistical AreaEstimated AreaRelative Error
ChanChun18.3915.12−0.18133.48110.77−0.172.703.680.36
JiLin13.3414.810.1148.6160.490.243.164.320.37
SiPing4.384.11−0.0651.2768.570.341.631.980.21
LiaoYuan1.881.48−0.2119.9112.95−0.350.540.840.55
TongHua5.815.950.0214.1211.21−0.210.881.250.42
BaiShan0.040.050.333.695.020.362.943.130.06
SongYuan12.9112.43−0.0479.0953.25−0.334.796.100.27
BaiCheng20.3422.610.1151.6241.67−0.193.224.670.45
YanBian3.754.940.3220.0127.020.3512.0713.470.12
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhi, F.; Dong, Z.; Guga, S.; Bao, Y.; Han, A.; Zhang, J.; Bao, Y. Rapid and Automated Mapping of Crop Type in Jilin Province Using Historical Crop Labels and the Google Earth Engine. Remote Sens. 2022, 14, 4028. https://doi.org/10.3390/rs14164028

AMA Style

Zhi F, Dong Z, Guga S, Bao Y, Han A, Zhang J, Bao Y. Rapid and Automated Mapping of Crop Type in Jilin Province Using Historical Crop Labels and the Google Earth Engine. Remote Sensing. 2022; 14(16):4028. https://doi.org/10.3390/rs14164028

Chicago/Turabian Style

Zhi, Feng, Zhenhua Dong, Suri Guga, Yongbin Bao, Aru Han, Jiquan Zhang, and Yulong Bao. 2022. "Rapid and Automated Mapping of Crop Type in Jilin Province Using Historical Crop Labels and the Google Earth Engine" Remote Sensing 14, no. 16: 4028. https://doi.org/10.3390/rs14164028

APA Style

Zhi, F., Dong, Z., Guga, S., Bao, Y., Han, A., Zhang, J., & Bao, Y. (2022). Rapid and Automated Mapping of Crop Type in Jilin Province Using Historical Crop Labels and the Google Earth Engine. Remote Sensing, 14(16), 4028. https://doi.org/10.3390/rs14164028

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop