Next Article in Journal
A Population Spatialization Model at the Building Scale Using Random Forest
Next Article in Special Issue
Soybean EOS Spatiotemporal Characteristics and Their Climate Drivers in Global Major Regions
Previous Article in Journal
Nitrous Oxide Profiling from Infrared Radiances (NOPIR): Algorithm Description, Application to 10 Years of IASI Observations and Quality Assessment
Previous Article in Special Issue
Spatial Transferability of Random Forest Models for Crop Type Classification Using Sentinel-1 and Sentinel-2
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Developing High-Resolution Crop Maps for Major Crops in the European Union Based on Transductive Transfer Learning and Limited Ground Data

Academy of Disaster Reduction and Emergency Management Ministry of Emergency Management & Ministry of Education, School of National Safety and Emergency Management, Beijing Normal University, Beijing 100875, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(8), 1809; https://doi.org/10.3390/rs14081809
Submission received: 2 March 2022 / Revised: 1 April 2022 / Accepted: 6 April 2022 / Published: 8 April 2022

Abstract

:
Precise and timely information on crop spatial distribution over large areas is paramount to agricultural monitoring, food security, and policy development. Currently, automatically classifying crop types at a large scale is challenging due to the scarcity of ground data. Although previous studies have indicated that transductive transfer learning (TTL) is a promising method to address this problem, it performs poorly within regions where crop compositions and phenology differ largely. Here we transferred random forest classifiers trained in limited regions with diversified growing conditions and land covers to the rest of the study area where ground data are scarce, with more than 130,000 Sentinel-2 images processed using the Google Earth Engine (GEE) platform. We established the 10 m crop maps for four major crops (i.e., maize, rapeseed, winter, and spring Triticeae crops) across 10 European Union (EU) countries from 2018 to 2019. The final crop maps had a high accuracy with overall accuracy generally greater than 0.89, with user’s accuracy and producer’s accuracy ranging from 0.72 to 0.98. Moreover, the resulting maps were consistent with the NUTS-2 level official statistics, with R2 consistently greater than 0.9. We further analyzed the crop rotation patterns and found that the rotation intervals across these EU countries were generally at least one year. Maize was dominantly rotated with winter Triticeae crops or converted to other land covers in the following year. Rapeseed was generally grown in rotation with winter Triticeae crops, whereas the rotation patterns of winter and spring Triticeae crops were more diversified. Red Edge Position (REP) and Normalized Difference Yellow Index (NDYI) played significant roles in crop classification across the EU. This study highlights the potential of the developed TTL method for crop classification over large spatial extents where labeled data are limited and the differences in crop compositions and phenology are relatively large.

Graphical Abstract

1. Introduction

Unprecedented global population growth led to a continual increase in food demand and agricultural production intensification [1,2]. However, land degradation, shortages of cropland and water resources, and more frequent extreme weather events caused by global warming have hampered the increase in food production and threatened food security [3,4,5]. Accurate and near-real-time monitoring of crop spatial distribution over large areas has important implications for crop growth monitoring, agricultural land use and agronomic management optimization, and food security early warning [6].
With an enormous increase in the spatial and temporal information derived from satellites, remote sensing has been widely used in agricultural applications for decades [7,8,9,10]. Previous studies have primarily used the Normalized Difference Vegetation Index (NDVI) or Enhanced Vegetation Index (EVI) time series derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) to map crop types over large extents [11,12,13]. However, the coarse spatial resolution with a large number of mixed pixels could lead to uncertainties in the results, particularly in areas with complex terrain and various land cover types. Although finer-resolution crop maps have been successfully made using Landsat data [14], combining two contemporaneous satellites (i.e., Lansat 7 and 8) nominally provides only an eight-day revisit cycle and the images can be polluted by cloud, which cannot capture the rapid variations that occurred in crop phenology over a short time [15]. The Sentinel-2 constellation with a 10 m spatial resolution and a five-day revisit time has shown a great advantage in crop classification at the parcel level over large areas [16].
Many previous studies have conducted crop mapping using spectral features derived from coarse- to high-resolution satellite imagery based on supervised or unsupervised classification methods [17,18,19]. Additionally, the emerging deep learning approach has also been successfully applied in crop classification, but such studies focus on mapping crop types at a small scale [20,21]. Nevertheless, the largest constraint to applications of the prevalent supervised classification methods including random forest (RF), support vector machines (SVM), and decision tree (DT) is currently scarce training samples rather than abundant imagery [22]. Field surveys are the main ways to collect ground samples for land cover; however, they are time-consuming, extremely expensive, are updated infrequently, and have limited spatial coverage. Recent studies attempt to use transfer learning (TL) or unsupervised classification approaches to classify crop types in areas where training data is lacking [23,24]. TL has been recognized as a trustworthy tool compared to the unsupervised classification methods, which are readily vulnerable to outliers and noisy data. Pan and Yang (2010) categorized TL into three settings [25], one of which was transductive transfer learning (TTL), characterized by labeled data available only in the source domain. For instance, Wang et al. (2019) trained a robust supervised classification model in a certain region or year with plentiful georeferenced ground data and then applied it to other areas or years where labeled data are unavailable, reasonably assuming that crop compositions and phenological characteristics are similar across these regions or years [24]. Hao et al. (2020) successfully transferred an RF model trained on Cropland Data Layer (CDL) data to identify crop types across three test sites with an average overall accuracy of 93% [26]. However, studies thus far have focused on transferring a pre-trained classification model to small areas, while less work has been devoted to exploring the viability of such a model transferred to large spatial extents. Since the variations in crop phenology caused by different environmental conditions from one region to another will degrade the classification performance, it is challenging but necessary to investigate whether TTL can work well in crop mapping on a national or sub-continental scale. As a promising study area for such effort, some countries in the European Union (EU) have similar planting structures with Triticeae crops (e.g., wheat, barley, oats, rye, maize, and rapeseed) primarily cultivated, hence we aim to develop an automatic TTL approach to develop crop maps for these regions where ground data are limited.
Spatially explicit information on major crop types in the EU is paramount for many studies such as improving crop yield prediction, assessing climate change impact, investigating the land-use change and its drivers, and evaluating agricultural disaster losses [27,28,29,30]. In Europe, detailed information on the spatial distribution of agricultural parcels was contained in the Land Parcel Information System (LPIS) and the Geospatial Aid Application (GSAA) data, which were mostly unavailable to date as they contained confidential information [31]. Existing crop maps that are publicly available include: (1) Global crop area maps with a coarse resolution of 10 km such as SPAM [32] and MIRCA2000 [33]; (2) National-scale crop type maps at moderate- to high-resolution including the CDL and Corn-Soy Data Layer in the United States [34,35], the Crop Inventory (CI) dataset generated by Agriculture and Agri-Food Canada (AAFC), the ChinaCropArea1km [36], the land cover map for Germany [17,37], maps of the major crops in Northeast China [38], and rice paddy maps in Bangladesh [39]. Therefore, datasets of high-resolution crop maps for the EU are particularly needed to facilitate the development of related scientific research on food security and policy making.
The objective of this study is to (1) propose an automatic approach to identifying staple crops with limited ground samples by transferring an RF classifier, based on Sentinel-2 imagery through the Google Earth Engine (GEE) cloud computing platform; (2) demonstrate the reliability of this approach in mapping maize, rapeseed, winter Triticeae crops (winter wheat, barley, and rye), and spring Triticeae crops (spring wheat, barley, and oat) from 2018 to 2019 across the major cultivation areas in the EU.

2. Data

2.1. Study Area

The study area included 10 countries, i.e., England, the Netherlands, Germany, Denmark, France, Italy, Poland, Hungary, Slovakia, and the Czech Republic, covering approximately 195.8 Mha of land area across the EU (Figure 1). Triticeae crops, maize, and rapeseed are staple crops in the EU, accounting for 65.9% of the harvested area for primary crops [40]. We chose these countries for two reasons: (1) To cover the major planting areas of the EU. In particular, the harvested area of the above crops in these countries makes up a total of 62.7% across the EU, and in each country accounts for beyond 2%, except for the Netherlands, which is selected due to its detailed crop field data; (2) To span a widespread growing environment where the scalability and reliability of our methods are investigated, consisting of dry, temperate, continental, and polar and alpine climates [41].

2.2. Data

2.2.1. Satellite Data

The Sentinel-2A/B (S2) Multi-Spectral Instrument (MSI) Level-1C top-of-atmosphere (TOA) reflectance data were collected to identify crop types during 2018 and2019. We selected blue (B2), green (B3), red (B4), red edge 1 (B5), red edge 2 (B6), red edge 3 (B7), near-infrared (B8), short-wave infrared 1 (B11), and short-wave infrared 2 (B12) spectral bands at a 10–20 m spatial resolution. In total, more than 130,000 S2 images were processed by the GEE platform (Table S1).

2.2.2. Reference Data

Reference data were acquired from existing nationwide crop field datasets or land cover maps. The first data for England, provided by the Rural Payment Agency (https://ckan.publishing.service.gov.uk/dataset, accessed on 7 April 2022) and named crop map of England (CROME), were generated by an RF classifier combining Sentinel-1 and S2 images with an overall accuracy (OA) of 95.4% (81%) and a Kappa coefficient of 0.95 (0.8) for 2019 (2018). CROME included over 20 main crop types, grassland, and non-agricultural land coves (e.g., woodland and water), and the user’s accuracy (UA) and producer’s accuracy (PA) of major crop types were generally greater than 80%. The second data, 10 m land cover maps for France (https://www.theia-land.fr/en/product/land-cover-map, accessed on 7 April 2022), classified 23 land cover classes since 2018 based on S2 images by an RF classifier with an OA of 89% and F1 scores for crop classes exceeding 90% [42]. The third dataset was the crop field dataset retrieved from the Base Registration Crop Parcels (BRP) in the Netherlands, providing the cultivated crop information at the parcel level (i.e., polygon) in Agricultural Area Netherlands (www.PDOK.nl, accessed on 7 April 2022). In addition, the Land Use/Cover Area frame Survey (LUCAS) in situ data was used to directly validate the classification results for all EU countries in 2018, conducted every three years since 2006 [43,44]. The previous LUCAS core observation (i.e., a circle with a 1.5 m radius and area of 7.07 m2) that corresponds to a fraction of an S1/2 pixel cannot be directly useful for such decametric sensors. The up-to-date LUCAS 2018 data delineated a polygon of 0.52 ha with homogeneous land cover by a new LUCAS Copernicus module, making it suitable for verifying high-resolution land cover maps.
The land cover maps for England and France produced by the RF classifier necessitated quality control to ensure the high-quality training and validation of the resultant crop maps. They were both used for sampling as more widespread ground samples covering diverse growing conditions over space perhaps made large-area crop mapping using TTL successful. The classification results in CROME were represented as hexagonal cells of 0.41 ha (approximate to a 64 m grid); however, the predicted confidence information for each cell was lacking. Thus, we first generated a 12 km-by-12 km grid across England and then selected parcels larger than 1 ha for each land cover, which corresponds to 3 cells with the same label. Finally, we randomly sampled 250 points in each grid cell and filtered out the points falling outside the parcels. Such sampling method was repeated for land cover maps of France; the differences were that (1) a high confidence mask with a threshold of 95% was used to restrict the sampling extent; (2) the “straw cereals” class was excluded from sampling as it contained both winter and spring Triticeae crops and was not directly usable; and (3) a 50 km grid cell was generated across France. In total, 56,877 and 55,144 training samples were acquired for 2018 and 2019, respectively, containing five land cover classes including other land covers (e.g., sugar beet, potato, sunflower, grassland), maize, rapeseed, winter Triticeae crops (i.e., winter wheat, barley, and rye) and spring Triticeae crops (i.e., spring wheat, barley, and oat) (Table 1).

2.2.3. Agricultural Statistical Data

The cultivation area of each crop type at a national (i.e., Nomenclature of Territorial Units for Statistics-1 regions, NUTS-1 regions) and subnational level (i.e., NUTS-2 regions) from 2018 to 2019 was acquired from Eurostat (https://ec.europa.eu/eurostat/databrowser, accessed on 7 April 2022). Such data were used as additional validation of the resultant maps, particularly in areas where ground samples were not available.

3. Methods

The flowchart for TTL-based crop mapping is presented schematically in Figure 2, consisting of image preprocessing, feature selection, transferring supervised classification, and accuracy assessment.

3.1. Image Processing

We first used the Landsat simple cloud score algorithm to mask out low-quality pixels due to cloud contamination [45]. Such a method assumed that: (1) clouds are reasonably bright in the blue and cirrus bands; (2) clouds are reasonably bright in all visible bands; (3) clouds are moist; (4) clouds are not snow. It calculated cloud score based on four spectral bands (Aerosols, Blue, Green, and Red band) and two spectral indices (Normalized Difference Moisture Index (NDMI) and Normalized Difference Snow Index (NDSI)). Clouds were detected more thoroughly using this method than using the QA60 band which could only mask out dense and cirrus clouds [18]. Moreover, cloud shadows were masked using the Temporal Dark Outlier Mask (TDOM) method [46]. Then, we reconstructed the S2 time series by a moving median composite method as the number of clear observations showed obvious spatiotemporal heterogeneity [18]. In particular, we first generated a time series of 15-day temporal interval composites by calculating the median value of observations at each interval. The reason for choosing the 15-day composite was that there was no substantial shift in plant growth within two weeks and the overall classification performance was similar between 10-day and 15-day composites [18,47]. Second, a gap-filling method was conducted on the composite time series in which a given observation was replaced by the median value of three adjacent observations (i.e., previous, current, and subsequent observations). Such a method is readily realized in the GEE platform and can generate homogenous results with high computational efficiency [48]. It is also not affected by outliers in time series compared to the mean value composite method.

3.2. Feature Selection

Six spectral bands and three spectral indices were selected for crop mapping, including B5, B6, B7, B8, B11, B12, NDVI, Normalized Difference Yellow Index (NDYI), and Red Edge Position (REP) (Figure 3). Previous studies have demonstrated that these spectral bands and indices could provide valuable information to distinguish different land covers such as SWIR and red-edge bands that could discriminate between maize and soybean [18]. Moreover, NDVI was a commonly used indicator that could represent seasonal patterns of phenology, contributing to the discrimination between winter and spring Triticeae crops nearly in April, as well as crops and non-agricultural land covers such as grassland and forest (Figure 3a). NDYI was useful for differentiating rapeseed and other non-flowering crops as it could capture the increasing yellowness with flowering [49] (Figure 3b). REP also contributed to differentiating maize from other land covers (Figure 3c). Other visible bands were excluded owing to their remarkably strong scattering caused by the atmosphere [50]. The three spectral indices mentioned above were calculated as follows:
NDVI = ρ n i r ρ r e d ρ n i r + ρ r e d
NDYI = ρ g r e e n ρ b l u e ρ g r e e n + ρ b l u e
REP = 705 + 35 × 0.5 × ( ρ r e 3 + ρ r e d ) ρ r e 1 ρ r e 2 ρ r e 1
where ρ b l u e , ρ g r e e n , ρ r e d , ρ r e 1 , ρ r e 2 , ρ r e 3 , ρ n i r represent the TOA reflectance values of B2 (496.6 nm (S2A)/492.1 nm (S2B)), B3 (560 nm (S2A)/559 nm (S2B)), B4 (664.5 nm (S2A)/665 nm (S2B)), B5 (703.9 nm (S2A)/703.8 nm(S2B)), B6 (740.2 nm (S2A)/739.1 nm (S2B)), B7 (782.5 nm (S2A)/779.7 nm (S2B)), and B8 (835.1 nm (S2A)/833 nm (S2B)) in the S2 MSI sensor.

3.3. Transferring Supervised Classification

We used RF classifier to map crop extent, which has been widely applied in land cover classification across different geographical areas [24,37,48,51]. Moreover, previous studies have demonstrated that RF was immune to high data dimensionality and multicollinearity and robust to noisy training data [7,52]. RF is an ensemble classifier consisting of multiple decision trees in which two-thirds of the training samples are trained and the remaining one-third is used for cross-validation [53]. The final classification output is selected as the class with the maximum votes against all decision trees. The number of decision trees (numberOfTrees) was set at 100 according to the trade-off between the classification accuracy and computational efficiency, while other parameters were set to the default in the GEE platform.
We first trained the RF classifier on training samples derived from the land cover maps of England and France (mentioned in Section 2.2.2) using a total of 126 spectral features (14 composite periods with 9 features in each period) (Table 2). Our supervised models classified land covers into five classes, i.e., “others”, “maize”, “rapeseed”, “winter Triticeae crops”, and “spring Triticeae crops”. We tried a range of time windows (e.g., 1 January to 31 December or 1 February to 31 December) and determined to select 1 April to 31 October considering the trade-off between the number of good satellite observations and classification accuracy. The model was then transferred to other countries where labeled data were lacking, to test its spatial scalability and robustness in large-area crop mapping with limited training data.

3.4. Accuracy Assessment

First, point-based validation was implemented in two ways. One was performed based on crop field datasets or land cover maps in England, France, and the Netherlands for two years, and the other was based on LUCAS in situ data covering the whole studied areas only for 2018 due to data unavailability in 2019 (Table 2). In particular, we first adopted a stratified random sampling design, 5000 samples were randomly distributed across each country and the number of each class’s samples was determined according to the areal ratio in each country [54]. We used a sampling unit consisting of 3 × 3 pixels to mitigate the adverse impact of geolocation inaccuracy and then labeled these samples according to reference data. The samples located in homogeneous parcels and larger than 1 ha or corresponding confidence greater than 95% were retained to ensure high-quality validation. Overall, 4953 (4968), 4991 (4984), and 2620 (2634) validation samples were reserved in 2018 (2019) for England, France, and the Netherlands, respectively. Additionally, the LUCAS Copernicus data were linked with LUCAS core points to augment ground truth information based on the conditions that the LUCAS theoretical geolocation fell within the LUCAS polygon and the reported Copernicus level-2 land cover was consistent. As a result, 28,882 LUCAS Copernicus points were retained for validation. Finally, the accuracy was summarized by the confusion matrix, including four metrics (i.e., PA, UA, OA, and F1 score) which were calculated as:
PA i = N i R i × 100
UA i = N i C i × 100
OA = N c n × 100
F 1 i = 2 × PA × UA PA + UA
where N c is the total number of correctly classified validation samples, n represents total number of validation samples, N i denotes the number of correctly identified validation samples of class i, R i donates the number of validation samples of class i, C i implies the number of validation samples classified as class i.
In addition, area comparison using statistical data was a commonly used method for accuracy assessments on large scales [12,35,37]. The area estimates for each class from the resultant crop maps were compared to the NUTS-2 level statistics, except for spring Triticeae crops that were only recorded at the NUTS-1 level. The coefficient of determination (R2) and the root mean square error (RMSE) were calculated to evaluate the consistency in the area between classification maps and statistical data, which were defined as:
R 2 = 1 i = 1 n ( O i E i ) 2 i = 1 n ( O i O ¯ ) 2
RMSE = i = 1 n ( O i E i ) 2 n
where O i and E i are the statistical and estimated area of a NUTS-2 (NUTS-1) geopolitical entity i, respectively. n represents the total number of countries.

4. Results

4.1. Accuracy Assessment of the Developed TTL-Based Method

To examine whether RF models trained in limited regions with plentiful training data could be applied to classify crops over large areas where ground data were lacking, we conducted rigorous accuracy assessments combining both point-based validation and subnational-scale area comparison (Figure 4 and Figure 5, Table 3 and Table 4, Tables S2–S8). In total, the OA was generally greater than 0.89, whereas PA and UA varied between different countries or crops. The crop maps for three countries with abundant validation samples showed higher OA ranging from 0.93 to 0.96. The best-performing crops were winter Triticeae crops in England and maize in the Netherlands with PA and UA above 0.9. Winter Triticeae crops generally showed higher accuracy with PA and UA > 0.88 than spring Triticeae crops (PA and UA > 0.7). Rapeseed achieved moderately high accuracy (PA and UA > 0.84). Since the land cover map for France and the LUCAS data did not distinguish between winter and spring Triticeae crops, we defined them as a single class named Triticeae crops. For France, Triticeae crops were mapped with high accuracy with an F1 score of 0.9, followed by rapeseed and maize (F1 ≥ 0.86). Overall, validating the crop map for 2018 using the LUCAS data covering the whole study area indicated that the classification performance was the best for maize with an F1 score of 0.81, followed by rapeseed and Triticeae crops (F1 ≥ 0.77).
The comparisons showed that the estimated and statistical areas for the winter Triticeae crops, maize, and rapeseed were evenly distributed around the 1:1 line, whereas those for the spring Triticeae crops appeared to be more scattered, and the corresponding mapped area exhibited a slight underestimation, especially in 2018 compared with the other three crops (Figure 5). The area estimates were consistent with the statistics with R2 consistently greater than 0.9 and RMSE less than 126.39 Kha, indicating that our TTL-based method succeeded in large-area crop mapping where ground data were limited.

4.2. Crop Type Classification Results

The resultant crop maps for 10 EU countries in 2018 and 2019 at 10 m spatial resolutions were presented in Figure 6 and Figure S1. Overall the outlines of individual agricultural parcels were distinctly identified in the current maps, indicating that both inter-field heterogeneity and intra-field homogeneity were well characterized. The agricultural lands in the EU were dominated by medium-sized fields (2–25 ha) and the crop maps captured such diversified agricultural parcels fairly well. Although the field size was relatively small (even < 0.2 ha) in the Netherlands, a 10 m pixel was enough and the classification results clearly delineated these small parcels.
We further calculated the important values of the features in the RF models to investigate the critical factors for crop mapping across EU countries (Figure S2). NDVI during August was consistently ranked as the top feature, followed by REP, red edge 1, NDYI, and swir1 in unfixed order, highlighting the important roles of the REP, red edge bands, and NDYI in classifying the major crop types across EU countries.
A visual comparison of the current crop maps and the reference dataset for France and the Netherlands revealed that crop spatial distributions were correctly captured across numerous varieties of agricultural landscapes (Figure 7 and Figure S3). Our products outperformed existing land cover maps for England and France and filled the gaps in high-resolution crop maps for EU countries. One drawback of CROME was that the classification results were represented as hexagonal cells of 0.41 ha (approximate to a 64 m grid), which could contain several land cover types. By contrast, our results successfully made clear distinctions among the different crops. Furthermore, we distinguished winter and spring Triticeae crops even though the land cover maps for France incorporated them as one class.
Winter Triticeae crops were the most widely cultivated crops across the whole study area (accounting for 44.3% of the total four crop areas), followed by maize (26.4%), rapeseed (15.6%), and spring Triticeae crops (13.7%). They showed similarity in the spatial pattern along latitudes and were mainly distributed with latitudes ranging from 47–56°N (Figure S4a–d). More specifically, winter Triticeae crops were also cultivated south of 47°N, whereas spring Triticeae crops were rarely planted, and maize was mostly distributed with latitudes ranging from 47 to 49 and 51 to 54°N (roughly accounting for 54% of the total maize area). However, these crops differed largely in the spatial distribution along longitudes. Winter Triticeae crops were widely distributed with longitudes ranging from 2°W to 5°E and 8 to 21°E, whereas the distribution of spring Triticeae crops was more concentrated with longitudes ranging from 8 to 12°E and 15 to 23°E (Figure S4e,h). The distribution of maize was relatively widespread with longitudes ranging from 3°W to 2°E, 5 to 14°E, and 16 to 23°E. Rapeseed was distributed at longitudes of 2°W–5°E and 9–20°E (Figure S4f,g).

4.3. Crop Rotation Analysis

The rotation patterns of each crop were further analyzed based on high-accuracy crop maps from 2018 to 2019 (Figure 8, Figure 9 and Figure 10). The rotation pattern in a certain crop pixel was considered as no rotation (≥1-year break) if the crop was identified in the pixel in both 2018 and 2019 (either 2018 or 2019). Overall, the crop rotation intervals were dominated by at least a 1-year break in 85.9% of pixels, differing among different countries or crops (Figure 8). Rapeseed was the most remarkable crop with all the countries consistently characterized by a rotation break of ≥1 year in more than 99% of pixels, followed by spring Triticeae crops (90.9%), winter Triticeae crops (84.1%), and maize (78.7%) (Figure 8). Denmark had the lowest proportion of pixels (<80%) with a rotation for more than one year. Moreover, a relatively less proportion of rotation was found in France, Italy, the Netherlands, and Poland for maize, and Germany and England for winter Triticeae crops. The spatial details of the rotation patterns further demonstrated the differences in rotation patterns among different crops (Figure 9). Interestingly, the crop parcels from the previous year were distributed in adjacent fields in the following year (Figure 9).
We also investigated the land cover type used for crop rotation in the following year by calculating the proportions of pixels for a certain crop in 2018 converted to other land cover types in 2019 (Figure 10). Maize was generally grown in rotation with winter Triticeae crops (50.7% of the total pixels had undergone this rotation) or converted to other land covers (40.7%), whereas the rotation of rapeseed and winter Triticeae crops was the dominant pattern for rapeseed (86.5%). The rotation patterns of the winter and spring Triticeae crops were more complicated. The former was mainly rotated with other land covers (44.1%), rapeseed (23.8%), maize (20.5%), and a little with the latter (11.6%). On the contrary, spring Triticeae crops were primarily rotated with winter Triticeae crops (44.6%), followed by other land covers (31.7%) and rapeseed (14.9%). Additionally, the land cover type used for crop rotation (except for rapeseed) varied between countries. In Czechia and Denmark, maize was also rotated with spring Triticeae crops (>27%). Unlike the predominant rotation of maize and winter Triticeae crops in Germany, France, and England (>66%), maize was grown in rotation with other land covers in Hungary, Italy, and the Netherlands (>69%). Moreover, the winter Triticeae crops were largely rotated with the spring Triticeae crops in Denmark (>45%). In addition, the rotation patterns of the three crops other than rapeseed in the Netherlands were consistent, that is, dominantly rotated with other land covers (>77%).

5. Discussions

5.1. Advantages of the TTL-Based Method for Large-Area Crop Mapping

Although a large number of satellites imageries are increasingly accessible with unprecedented temporal, spatial, and spectral frequency, the availability of ground data did not keep up, hampering the accurate training and validation of commonly used supervised classification methods [22,24,26,55]. Since large-area crop mapping is paramount for ensuring food security and supervised classification methods are predominant tools for realizing such a mission, it is necessary to investigate the feasibility of transferring a supervised model trained in limited regions with abundant ground data to large spatial extents with scarce data. The 10 m crop maps for the 10 EU countries exhibited high accuracy and clearly identified the field boundaries across various agricultural landscapes, suggesting that our TTL-based method is a reliable and promising method to apply to crop classification over large areas with limited ground data. Nevertheless, some points need to be noted: (1) the training samples should cover different growing conditions as reasonably as possible because the phenological differences are relatively large across different regions, and (2) the training samples are required to include plentiful land cover types, especially those that were absent in sampling regions but widespread in other areas. This is well supported by the findings of Kluger et al. (2021) [55], in that it is better to collect training samples used for TTL-based crop classification in a region with diverse crop-type composition than to collect training samples in a region with low crop-type diversity. Previous studies used growing degree days (GDD) to represent the differences in crop phenology and indicated that lower similarity of GDD in two regions could result in poorer performance of an RF model transferred from one to the other [24]. Given our study area, maize phenology is advanced with latitude decreasing (Figure S5). Thus, how to reasonably cope with such intra-class phenology differences is essential to achieving sound TTL-based crop mapping over large areas. To avoid omission of maize fields that differed from training regions in phenology, we randomly sampled a small number of maize samples from land cover products from France as France was characterized by diversified climatic (e.g., temperate, continental, and polar and alpine climates) characteristics, which to a large extent contain the climate types of the rest of study area (except for dry climate) [41]. This procedure allowed for the capture of phenology differences in other regions such as Italy and Germany, whereas merely using training samples generated in England cannot work. Moreover, the CROME data for England did not contain some types of land cover that frequently occur in other countries such as soybean, sunflower, rice, and mineral surfaces, which are readily misclassified as desired crop types. Therefore, we also collected samples for these classes to augment the diversity of the training sets.
In short, the TTL-based method makes ample use of ground data acquired from limited regions to classify crop types over large spatial extents where labeled data are not available. Compared with previous studies, we first attempt to reasonably enrich the diversity of the training sets to overcome the disadvantages of the TTL method that performs poorly in settings where the differences in crop composition and phenology are relatively large.

5.2. Detailed Investigation of Crop Rotation Patterns

The rigorous validation by both point-based validation and subnational-scale area comparison demonstrated high accuracy and reliability of the resultant 10 m crop maps during 2018 and2019. Therefore, it is possible to further investigate the rotation patterns of the four crops, including the rotation interval and land cover types. We found that the crop rotation interval of these EU countries was characterized by at least 1 year, which was consistent with actual situations [56]. Crops are traditionally grown in rotation in a large number of European countries, and it is recognized that (1) most crops benefit from the nutrients released by the residues of the preceding crop [57] and (2) crop rotation contributes to boosting yield, ameliorating soil fertility, and reducing pests and diseases, and consequently improving the resilience of the crop production system under climate change [58,59,60]. In addition, we revealed that maize was generally rotated with winter Triticeae crops or converted to other land covers. Rapeseed was dominantly grown in rotation with winter Triticeae crops, whereas the rotation patterns of winter and spring Triticeae crops were diversified. These findings were supported by some field experiments in Europe [56,61,62,63]. We provided the refined spatial characteristics of crop rotation patterns across the EU at the parcel level. Such information is novel and essential for many purposes such as agro-ecosystem modeling, land use, and cropping system optimization.

5.3. Uncertainty

Although the resultant crop maps have been proven to be reliable, there are still some uncertainties. First, the largest constraint is the availability of S2 data which is primarily restricted by clouds. Even though the moving median composite method partly fills the gaps due to cloud contamination, it continues to malfunction, especially in regions (e.g., Northern Europe) with dense cloud coverage, resulting in a high omission error. This issue is expected to be solved using Harmonized Landsat and Sentinel-2 (HLS) data, which have been applied to crop mapping [26,37,64]. However, such data have not been made available in the GEE platform, hampering rapid and highly efficient crop mapping over large areas. Second, the TOA reflectance product was used rather than surface reflectance (SR) owing to the limited temporal availability of the latter on the GEE platform for our study area. Although previous studies have substantiated that TOA reflectance products are reliable for crop classification [18,49], it might dampen the classification accuracy because of the pronounced scattering and absorption of visible radiation by atmospheric contamination [65]. Here we not only excluded shorter-wavelength bands (e.g., red and blue bands) that were most susceptible to atmospheric effects from feature selection but also used indices based on spectral band ratios to reduce such effects [50,66]. With more S2 SR data available for the study area since 2019, it is promising to improve the classification accuracy in subsequent years. Third, we did not collect training samples for winter and spring Triticeae crops covering different environmental conditions and such ground data were lacking in regions other than England. It not only led to an underestimation of area compared with the official statistics but also restricted the scalability of our TTL-based method. For instance, winter wheat phenology is advanced with latitude decreasing, especially in Spain where the anthesis date can be approximately 40 days earlier than that in England [67,68], hampering the application of the TTL method in such regions with crops intensively cultivated due to the scarcity of corresponding ground data. Fourth, we grouped winter and spring Triticeae crops as single classes, respectively. Nevertheless, there are slight differences in phenology among wheat, barley, rye, and oat. Further researches should attempt to apply the TTL method to refine crop classification (both in crop category and spatial resolution) over larger spatial extents using increasing and precise ground data.

6. Conclusions

We produced the crop maps at 10 m spatial resolution for four main crops across 10 EU countries from 2018 to 2019 through the GEE platform, by transferring an RF model trained in limited areas with plentiful ground data to other larger extents with scarce ground samples. Compared with previous efforts to classify crops using the TTL method, we reasonably augmented the diversity of training data by additional sampling covering different growing conditions and land cover types, addressing the poor performance of the TTL method when the differences in crop composition and phenology are relatively large. We also demonstrated the important roles of REP and NDYI in crop mapping across the EU. The overall accuracy of the final maps was generally greater than 0.89, with PA and UA ranging from 0.72 to 0.98. In comparison with the NUTS-2 level official statistics, the resultant maps proved to be highly accurate with R2 consistently greater than 0.9. These rigorous validations indicated that the TTL method was robust in classifying crops over large spatial extents where ground data are limited. The crop rotation analyses showed that the rotation interval across these EU countries was predominantly at least one year. We found that maize was generally rotated with winter Triticeae crops or converted to other land covers in the following year. Rapeseed was grown in rotation with winter Triticeae crops, whereas the rotation patterns of winter and spring Triticeae crops were complicated. The TTL method used in this study can be readily replicated to classify crop types in other large regions, given the strong computing power and ever-increasing satellite observations in the GEE platform. The datasets (available at https://doi.org/10.17632/4frtdtxsk9.2, accessed on 1 March 2022) are crucial for many purposes including agro-ecosystem modeling, crop growth monitoring, land use and cropping system optimization, and policy making.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/rs14081809/s1, Figure S1: Crop maps at 10 m resolution of ten EU countries for 2018. Figure S2: Feature importance values for the tope of 50 variables from random forest models in 2018 (a) and 2019 (b). Figure S3: Visual comparison between our final crop maps (a1–a8) and existing reference datasets in 2018 for France (b1–b6) and Netherlands (b7,b8). For France, winter and spring Triticeae crops was concluded as a single category (green color). Figure S4: Distribution of crop area along latitude gradients (a–d) and longitude gradients (e–h). Figure S5: Temporal profiles of NDVI for maize of different countries. Lines depicts the mean values. Shaded area depicts error bars with one positive/negative standard deviation of maize. Table S1: The number of Sentinel-2 images processed in this study for each country in 2018 and 2019. Table S2: Confusion matrix of the point-based validation for England in 2019. Table S3: Confusion matrix of the point-based validation for England in 2018. Table S4: Confusion matrix of the point-based validation for Netherlands in 2019. Table S5: Confusion matrix of the point-based validation for Netherlands in 2018. Table S6: Confusion matrix of the point-based validation for France in 2019. Table S7: Confusion matrix of the point-based validation for France in 2018. Table S8: Confusion matrix of the point-based validation for the whole study area using LUCAS data in 2018.

Author Contributions

Conceptualization, Z.Z.; Data curation, J.C.; Formal analysis, Y.L.; Methodology, Y.L.; Validation, L.Z.; Visualization, J.H.; Writing—original draft, Y.L.; Writing—review & editing, J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the National Natural Science Foundation of China (42061144003, 41977405).

Data Availability Statement

The data presented in this study are openly available at https://doi.org/10.17632/4frtdtxsk9.2.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Godfray, H.C.J.; Beddington, J.R.; Crute, I.R.; Haddad, L.; Lawrence, D.; Muir, J.F.; Pretty, J.; Robinson, S.; Thomas, S.M.; Toulmin, C. Food Security: The Challenge of Feeding 9 Billion People. Science 2010, 327, 812–818. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Meyfroidt, P. Trade-offs between environment and livelihoods: Bridging the global land use and food security discussions. Glob. Food Secur. 2018, 16, 9–16. [Google Scholar] [CrossRef]
  3. Bindraban, P.S.; van der Velde, M.; Ye, L.M.; van den Berg, M.; Materechera, S.; Kiba, D.I.; Tamene, L.; Ragnarsdottir, K.V.; Jongschaap, R.; Hoogmoed, M.; et al. Assessing the impact of soil degradation on food production. Curr. Opin. Environ. Sustain. 2012, 4, 478–488. [Google Scholar] [CrossRef]
  4. Fitton, N.; Alexander, P.; Arnell, N.; Bajzelj, B.; Calvin, K.; Doelman, J.; Gerber, J.S.; Havlik, P.; Hasegawa, T.; Herrero, M.; et al. The vulnerabilities of agricultural land and food production to future water scarcity. Glob. Environ. Chang. 2019, 58, 101944. [Google Scholar] [CrossRef]
  5. Tirado, M.; Clarke, R.; Jaykus, L.; McQuatters-Gollop, A.; Frank, J. Climate change and food safety: A review. Food Res. Int. 2010, 43, 1745–1765. [Google Scholar] [CrossRef]
  6. Bajzelj, B.; Richards, K.S.; Allwood, J.M.; Smith, P.; Dennis, J.S.; Curmi, E.; Gilligan, C.A. Importance of food-demand management for climate mitigation. Nat. Clim. Chang. 2014, 4, 924–929. [Google Scholar] [CrossRef] [Green Version]
  7. Belgiu, M.; Dragut, L. Random forest in remote sensing: A review of applications and future directions. ISPRS J. Photogramm. Remote Sens. 2016, 114, 24–31. [Google Scholar] [CrossRef]
  8. Friedl, M.A.; McIver, D.K.; Hodges, J.C.F.; Zhang, X.Y.; Muchoney, D.; Strahler, A.H.; Woodcock, C.E.; Gopal, S.; Schneider, A.; Cooper, A.; et al. Global land cover mapping from MODIS: Algorithms and early results. Remote Sens. Environ. 2002, 83, 287–302. [Google Scholar] [CrossRef]
  9. Fritz, S.; See, L.; McCallum, I.; You, L.Z.; Bun, A.; Moltchanova, E.; Duerauer, M.; Albrecht, F.; Schill, C.; Perger, C.; et al. Mapping global cropland and field size. Glob. Chang. Biol. 2015, 21, 1980–1992. [Google Scholar] [CrossRef]
  10. Moran, M.; Inoue, Y.; Barnes, E. Opportunities and limitations for image-based remote sensing in precision crop management. Remote Sens. Environ. 1997, 61, 319–346. [Google Scholar] [CrossRef]
  11. Gumma, M.K.; Thenkabail, P.S.; Teluguntla, P.; Rao, M.N.; Mohammed, I.A.; Whitbread, A.M. Mapping rice-fallow cropland areas for short-season grain legumes intensification in South Asia using MODIS 250 m time-series data. Int. J. Digit. Earth 2016, 9, 981–1003. [Google Scholar] [CrossRef] [Green Version]
  12. Konduri, V.S.; Kumar, J.; Hargrove, W.W.; Hoffman, F.M.; Ganguly, A.R. Mapping crops within the growing season across the United States. Remote Sens. Environ. 2020, 251, 112048. [Google Scholar] [CrossRef]
  13. Ozdogan, M. The spatial distribution of crop types from MODIS data: Temporal unmixing using Independent Component Analysis. Remote Sens. Environ. 2010, 114, 1190–1204. [Google Scholar] [CrossRef]
  14. Qiu, B.W.; Lu, D.F.; Tang, Z.H.; Chen, C.C.; Zou, F.L. Automatic and adaptive paddy rice mapping using Landsat images: Case study in Songnen Plain in Northeast China. Sci. Total Environ. 2017, 598, 581–592. [Google Scholar] [CrossRef]
  15. Roy, D.P.; Yan, L. Robust Landsat-based crop time series modelling. Remote Sens. Environ. 2020, 238, 110810. [Google Scholar] [CrossRef]
  16. Defourny, P.; Bontemps, S.; Bellemans, N.; Cara, C.; Dedieu, G.; Guzzonato, E.; Hagolle, O.; Inglada, J.; Nicola, L.; Rabaute, T.; et al. Near real-time agriculture monitoring at national scale at parcel resolution: Performance assessment of the Sen2-Agri automated system in various cropping systems around the world. Remote Sens. Environ. 2019, 221, 551–568. [Google Scholar] [CrossRef]
  17. Preidl, S.; Lange, M.; Doktor, D. Introducing APiC for regionalised land cover mapping on the national scale using Sentinel-2A imagery. Remote Sens. Environ. 2020, 240, 111673. [Google Scholar] [CrossRef]
  18. You, N.S.; Dong, J.W. Examining earliest identifiable timing of crops using all available Sentinel 1/2 imagery and Google Earth Engine. ISPRS J. Photogramm. Remote Sens. 2020, 161, 109–123. [Google Scholar] [CrossRef]
  19. Zhang, G.; Xiao, X.; Biradar, C.M.; Dong, J.; Qin, Y.; Menarguez, M.A.; Zhou, Y.; Zhang, Y.; Jin, C.; Wang, J.; et al. Spatiotemporal patterns of paddy rice croplands in China and India from 2000 to 2015. Sci. Total Environ. 2017, 579, 82–92. [Google Scholar] [CrossRef]
  20. Zhong, L.H.; Hu, L.; Zhou, H.; Tao, X. Deep learning based winter wheat mapping using statistical data as ground references in Kansas and northern Texas, US. Remote Sens. Environ. 2019, 233, 111411. [Google Scholar] [CrossRef]
  21. Zhong, L.; Hu, L.; Zhou, H. Deep learning based multi-temporal crop classification. Remote Sens. Environ. 2019, 221, 430–443. [Google Scholar] [CrossRef]
  22. Burke, M.; Driscoll, A.; Lobell, D.B.; Ermon, S. Using satellite imagery to understand and promote sustainable development. Science 2021, 371, 1219. [Google Scholar] [CrossRef] [PubMed]
  23. Nowakowski, A.; Mrziglod, J.; Spiller, D.; Bonifacio, R.; Ferrari, I.; Mathieu, P.P.; Garcia-Herranz, M.; Kim, D.-H. Crop type mapping by using transfer learning. Int. J. Appl. Earth Obs. Geoinf. 2021, 98, 102313. [Google Scholar] [CrossRef]
  24. Wang, S.; Azzari, G.; Lobell, D.B. Crop type mapping without field-level labels: Random forest transfer and unsupervised clustering techniques. Remote Sens. Environ. 2019, 222, 303–317. [Google Scholar] [CrossRef]
  25. Pan, S.J.; Yang, Q. A Survey on Transfer Learning. IEEE Trans. Knowl. Data Eng. 2010, 22, 1345–1359. [Google Scholar] [CrossRef]
  26. Hao, P.; Di, L.; Zhang, C.; Guo, L. Transfer Learning for Crop classification with Cropland Data Layer data (CDL) as training samples. Sci. Total Environ. 2020, 733, 138869. [Google Scholar] [CrossRef] [PubMed]
  27. Cao, J.; Zhang, Z.; Luo, Y.C.; Zhang, L.L.; Zhang, J.; Li, Z.Y.; Tao, F.L. Wheat yield predictions at a county and field scale with deep learning, machine learning, and google earth engine. Eur. J. Agron. 2021, 123, 126204. [Google Scholar] [CrossRef]
  28. Skakun, S.; Kussul, N.; Shelestov, A.; Kussul, O. The use of satellite data for agriculture drought risk quantification in Ukraine. Geomat. Nat. Hazards Risk 2016, 7, 901–917. [Google Scholar] [CrossRef]
  29. Yang, J.; Gong, P.; Fu, R.; Zhang, M.; Chen, J.; Liang, S.; Xu, B.; Shi, J.; Dickinson, R.E. The role of satellite remote sensing in climate change studies. Nat. Clim. Chang. 2013, 3, 875–883. [Google Scholar] [CrossRef]
  30. Zhong, L.; Hu, L.; Yu, L.; Gong, P.; Biging, G.S. Automated mapping of soybean and corn using phenology. ISPRS J. Photogramm. Remote Sens. 2016, 119, 151–164. [Google Scholar] [CrossRef] [Green Version]
  31. Tóth, K.; Kučas, A. Spatial information in European agricultural data management. Requirements and interoperability supported by a domain model. Land Use Policy 2016, 57, 64–79. [Google Scholar] [CrossRef]
  32. Yu, Q.Y.; You, L.Z.; Wood-Sichra, U.; Ru, Y.T.; Joglekar, A.K.B.; Fritz, S.; Xiong, W.; Lu, M.; Wu, W.B.; Yang, P. A cultivated planet in 2010-Part 2: The global gridded agricultural-production maps. Earth Syst. Sci. Data 2020, 12, 3545–3572. [Google Scholar] [CrossRef]
  33. Portmann, F.T.; Siebert, S.; Döll, P. MIRCA2000-Global monthly irrigated and rainfed crop areas around the year 2000: A new high-resolution data set for agricultural and hydrological modeling. Glob. Biogeochem. Cycles 2010, 24, 782. [Google Scholar] [CrossRef]
  34. Johnson, D.M.; Mueller, R. The 2009 Cropland Data Layer. Photogramm. Eng. Remote Sens. 2010, 76, 1201–1205. [Google Scholar]
  35. 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]
  36. 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]
  37. 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]
  38. 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]
  39. Singha, M.; Dong, J.; Zhang, G.; Xiao, X. High resolution paddy rice maps in cloud-prone Bangladesh and Northeast India using Sentinel-1 data. Sci. Data 2019, 6, 26. [Google Scholar] [CrossRef]
  40. Faostat, F. FAOSTAT Statistical Database; FAO (Food and Agriculture Organization of the United Nations): Rome, Italy, 2019. [Google Scholar]
  41. Beck, H.E.; Zimmermann, N.E.; McVicar, T.R.; Vergopolan, N.; Berg, A.; Wood, E.F. Present and future Koppen-Geiger climate classification maps at 1-km resolution. Sci. Data 2018, 5, 180214. [Google Scholar] [CrossRef] [Green Version]
  42. Inglada, J.; Vincent, A.; Arias, M.; Tardy, B.; Morin, D.; Rodes, I. Operational High Resolution Land Cover Map Production at the Country Scale Using Satellite Image Time Series. Remote Sens. 2017, 9, 95. [Google Scholar] [CrossRef] [Green Version]
  43. D’Andrimont, R.; Verhegghen, A.; Meroni, M.; Lemoine, G.; Strobl, P.; Eiselt, B.; Yordanov, M.; Martinez-Sanchez, L.; van der Velde, M. LUCAS Copernicus 2018: Earth-observation-relevant in situ data on land cover and use throughout the European Union. Earth Syst. Sci. Data 2021, 13, 1119–1133. [Google Scholar] [CrossRef]
  44. D’Andrimont, R.; Yordanov, M.; Martinez-Sanchez, L.; Eiselt, B.; Palmieri, A.; Dominici, P.; Gallego, J.; Reuter, H.I.; Joebges, C.; Lemoine, G.; et al. Harmonised LUCAS in-situ land cover and use database for field surveys from 2006 to 2018 in the European Union. Sci. Data 2020, 7, 352. [Google Scholar] [CrossRef] [PubMed]
  45. Oreopoulos, L.; Wilson, M.J.; Várnai, T. Implementation on Landsat Data of a Simple Cloud-Mask Algorithm Developed for MODIS Land Bands. IEEE Geosci. Remote Sens. Lett. 2011, 8, 597–601. [Google Scholar] [CrossRef] [Green Version]
  46. Housman, I.W.; Chastain, R.A.; Finco, M.V. An Evaluation of Forest Health Insect and Disease Survey Data and Satellite-Based Remote Sensing Forest Change Detection Methods: Case Studies in the United States. Remote Sens. 2018, 10, 1184. [Google Scholar] [CrossRef] [Green Version]
  47. Xu, X.; Conrad, C.; Doktor, D. Optimising Phenological Metrics Extraction for Different Crop Types in Germany Using the Moderate Resolution Imaging Spectrometer (MODIS). Remote Sens. 2017, 9, 254. [Google Scholar] [CrossRef] [Green Version]
  48. Teluguntla, P.; Thenkabail, P.S.; Oliphant, A.; Xiong, J.; Gumma, M.K.; Congalton, R.G.; Yadav, K.; Huete, A. A 30-m landsat-derived cropland extent product of Australia and China using random forest machine learning algorithm on Google Earth Engine cloud computing platform. ISPRS J. Photogramm. Remote Sens. 2018, 144, 325–340. [Google Scholar] [CrossRef]
  49. Han, J.; Zhang, Z.; Luo, Y.; Cao, J.; Zhang, L.; Zhang, J.; Li, Z. The RapeseedMap10 database: Annual maps of rapeseed at a spatial resolution of 10 m based on multi-source data. Earth Syst. Sci. Data 2021, 13, 2857–2874. [Google Scholar] [CrossRef]
  50. Song, X.-P.; Potapov, P.V.; Krylov, A.; King, L.; Di Bella, C.M.; Hudson, A.; Khan, A.; Adusei, B.; Stehman, S.V.; Hansen, M.C. National-scale soybean mapping and area estimation in the United States using medium resolution satellite imagery and field survey. Remote Sens. Environ. 2017, 190, 383–395. [Google Scholar] [CrossRef]
  51. Phalke, A.R.; Ozdogan, 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. Remote Sens. 2020, 167, 104–122. [Google Scholar] [CrossRef]
  52. Inglada, J.; Arias, M.; Tardy, B.; Hagolle, O.; Valero, S.; Morin, D.; Dedieu, G.; Sepulcre, G.; Bontemps, S.; Defourny, P.; et al. Assessment of an Operational System for Crop Type Map Production Using High Temporal and Spatial Resolution Satellite Optical Imagery. Remote Sens. 2015, 7, 12356–12379. [Google Scholar] [CrossRef] [Green Version]
  53. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  54. Olofsson, P.; Foody, G.M.; Herold, M.; Stehman, S.V.; Woodcock, C.E.; Wulder, M.A. Good practices for estimating area and assessing accuracy of land change. Remote Sens. Environ. 2014, 148, 42–57. [Google Scholar] [CrossRef]
  55. Kluger, D.M.; Wang, S.; Lobell, D.B. Two shifts for crop mapping: Leveraging aggregate crop statistics to improve satellite-based maps in new regions. Remote Sens. Environ. 2021, 262, 112488. [Google Scholar] [CrossRef]
  56. Kollas, C.; Kersebaum, K.C.; Nendel, C.; Manevski, K.; Müller, C.; Palosuo, T.; Armas-Herrera, C.M.; Beaudoin, N.; Bindi, M.; Charfeddine, M.; et al. Crop rotation modelling—A European model intercomparison. Eur. J. Agron. 2015, 70, 98–111. [Google Scholar] [CrossRef]
  57. Askegaard, M.; Eriksen, J. Residual effect and leaching of N and K in cropping systems with clover and ryegrass catch crops on a coarse sand. Agric. Ecosyst. Environ. 2008, 123, 99–108. [Google Scholar] [CrossRef]
  58. Huynh, H.T.; Hufnagel, J.; Wurbs, A.; Bellingrath-Kimura, S.D. Influences of soil tillage, irrigation and crop rotation on maize biomass yield in a 9-year field study in Muncheberg, Germany. Field Crop Res. 2019, 241, 107565. [Google Scholar] [CrossRef]
  59. Loges, R.; Bunne, I.; Reinsch, T.; Malisch, C.; Kluß, C.; Herrmann, A.; Taube, F. Forage production in rotational systems generates similar yields compared to maize monocultures but improves soil carbon stocks. Eur. J. Agron. 2018, 97, 11–19. [Google Scholar] [CrossRef]
  60. Teixeira, E.I.; de Ruiter, J.; Ausseil, A.-G.; Daigneault, A.; Johnstone, P.; Holmes, A.; Tait, A.; Ewert, F. Adapting crop rotations to climate change in regional impact modelling assessments. Sci. Total Environ. 2018, 616–617, 785–795. [Google Scholar] [CrossRef] [PubMed]
  61. Fischer, J.; Bohm, H.; Hess, J. Maize-bean intercropping yields in Northern Germany are comparable to those of pure silage maize. Eur. J. Agron. 2020, 112, 125947. [Google Scholar] [CrossRef]
  62. Poyda, A.; Wizemann, H.-D.; Ingwersen, J.; Eshonkulov, R.; Högy, P.; Demyan, M.S.; Kremer, P.; Wulfmeyer, V.; Streck, T. Carbon fluxes and budgets of intensive crop rotations in two regional climates of southwest Germany. Agric. Ecosyst. Environ. 2019, 276, 31–46. [Google Scholar] [CrossRef]
  63. Vasileiadis, V.P.; Sattin, M.; Otto, S.; Veres, A.; Palinkas, Z.; Ban, R.; Pons, X.; Kudsk, P.; van der Weide, R.; Czembor, E.; et al. Crop protection in European maize-based cropping systems: Current practices and recommendations for innovative Integrated Pest Management. Agric. Syst. 2011, 104, 533–540. [Google Scholar] [CrossRef]
  64. Claverie, M.; Ju, J.; Masek, J.G.; Dungan, J.L.; Vermote, E.F.; Roger, J.C.; Skakun, S.V.; Justice, C. The Harmonized Landsat and Sentinel-2 surface reflectance data set. Remote Sens. Environ. 2018, 219, 145–161. [Google Scholar] [CrossRef]
  65. Ju, J.; Roy, D.P.; Vermote, E.; Masek, J.; Kovalskyy, V. Continental-scale validation of MODIS-based and LEDAPS Landsat ETM+ atmospheric correction methods. Remote Sens. Environ. 2012, 122, 175–184. [Google Scholar] [CrossRef] [Green Version]
  66. Gao, F.; Jin, Y.; Schaaf, C.B.; Strahler, A.H. Bidirectional NDVI and atmospherically resistant BRDF inversion for vegetation canopy. IEEE Trans. Geosci. Remote Sens. 2002, 40, 1269–1278. [Google Scholar] [CrossRef]
  67. Harkness, C.; Semenov, M.A.; Areal, F.; Senapati, N.; Trnka, M.; Balek, J.; Bishop, J. Adverse weather conditions for UK wheat production under climate change. Agric. For. Meteorol. 2020, 282–283, 107862. [Google Scholar] [CrossRef]
  68. Oteros, J.; García-Mozo, H.; Botey, R.; Mestre, A.; Galán, C. Variations in cereal crop phenology in Spain over the last twenty-six years (1986–2012). Clim. Chang. 2015, 130, 545–558. [Google Scholar] [CrossRef]
Figure 1. The study area covers 10 countries across the European Union.
Figure 1. The study area covers 10 countries across the European Union.
Remotesensing 14 01809 g001
Figure 2. Flow chart of the methodology for the large-area crop mapping using transductive transfer learning approach.
Figure 2. Flow chart of the methodology for the large-area crop mapping using transductive transfer learning approach.
Remotesensing 14 01809 g002
Figure 3. Temporal profiles of (ai) for major crop and non-agricultural land cover types in the European Union. Lines depict the mean values of different crop and non-agricultural land cover types. Shaded areas depict error bars with one positive/negative standard deviation of maize, rapeseed, and winter and spring Triticeae crops.
Figure 3. Temporal profiles of (ai) for major crop and non-agricultural land cover types in the European Union. Lines depict the mean values of different crop and non-agricultural land cover types. Shaded areas depict error bars with one positive/negative standard deviation of maize, rapeseed, and winter and spring Triticeae crops.
Remotesensing 14 01809 g003
Figure 4. Validation results of the resultant crop maps. ENG: England, NLD: The Netherlands, FRA: France, PA: producer’s accuracy, UA: user’s accuracy.
Figure 4. Validation results of the resultant crop maps. ENG: England, NLD: The Netherlands, FRA: France, PA: producer’s accuracy, UA: user’s accuracy.
Remotesensing 14 01809 g004
Figure 5. Comparisons between the mapped and statistical areas at NUTS-2 level (except for spring Triticeae crops at country level) in 2018 and 2019 for (a) winter Triticeae crops, (b) maize, (c) rapeseed, and (d) spring Triticeae crops.
Figure 5. Comparisons between the mapped and statistical areas at NUTS-2 level (except for spring Triticeae crops at country level) in 2018 and 2019 for (a) winter Triticeae crops, (b) maize, (c) rapeseed, and (d) spring Triticeae crops.
Remotesensing 14 01809 g005
Figure 6. Crop maps at 10 m resolution for 8 of the 10 EU countries in 2019.
Figure 6. Crop maps at 10 m resolution for 8 of the 10 EU countries in 2019.
Remotesensing 14 01809 g006
Figure 7. Visual comparison between (a1a8) the resultant crop maps and existing reference datasets in 2019 for (b1b6) France and (b7,b8) the Netherlands. For France, winter and spring Triticeae crops were grouped as a single category (green color).
Figure 7. Visual comparison between (a1a8) the resultant crop maps and existing reference datasets in 2019 for (b1b6) France and (b7,b8) the Netherlands. For France, winter and spring Triticeae crops were grouped as a single category (green color).
Remotesensing 14 01809 g007
Figure 8. Proportions of areas that are grown in rotation intervals of 0 and ≥1-year break across the 10 European Union (EU) countries.
Figure 8. Proportions of areas that are grown in rotation intervals of 0 and ≥1-year break across the 10 European Union (EU) countries.
Remotesensing 14 01809 g008
Figure 9. Spatial details of rotation patterns in France (a,c), Germany (b), and England (d) from 2018 to 2019 for maize, winter Triticeae crops, rapeseed, and spring Triticeae crops.
Figure 9. Spatial details of rotation patterns in France (a,c), Germany (b), and England (d) from 2018 to 2019 for maize, winter Triticeae crops, rapeseed, and spring Triticeae crops.
Remotesensing 14 01809 g009
Figure 10. Area proportions of (ad) converted to other land cover types.
Figure 10. Area proportions of (ad) converted to other land cover types.
Remotesensing 14 01809 g010
Table 1. The number of training samples in 2018 and 2019.
Table 1. The number of training samples in 2018 and 2019.
Class2018Total2019Total
Other land coversSugar beet47127,73646927,014
Potato27092437
Sunflower925873
Grassland87338120
Forest41483874
Water150146
Building223265
Vineyards709716
Fallow573328
Legumes36023499
Rice11291241
Others43645046
Winter Triticeae cropsWinter wheat746212,521741512,216
Winter barley40373818
Winter rye1022983
Maize-6749674965266526
Rapeseed-4394439445784578
Spring Triticeae cropsSpring wheat1324547710824810
Spring barley33723039
Spring oat781689
Table 2. Summary of training and point-based validation.
Table 2. Summary of training and point-based validation.
YearTrainingNumber of SamplesValidationNumber of Samples
2018England, France56,87710 countries43,882
2019England, France55,144England, France, The Netherlands15,000
Table 3. Classification accuracy of the point-based validation for England and the Netherlands.
Table 3. Classification accuracy of the point-based validation for England and the Netherlands.
ClassesENG 2019ENG 2018NLD 2019NLD 2018
PAUAF1PAUAF1PAUAF1PAUAF1
Others0.980.980.980.970.970.970.980.980.980.980.970.98
Winter Triticeae crops0.950.920.930.930.910.920.930.920.930.880.900.89
Maize0.810.870.840.730.940.820.960.980.970.910.970.94
Rapeseed0.920.930.920.910.940.920.880.880.880.880.920.90
Spring Triticeae crops0.920.880.900.840.900.80.910.830.870.860.730.79
OA0.970.950.960.95
Table 4. Classification accuracy of the point-based validation for France and the whole study area using LUCAS data.
Table 4. Classification accuracy of the point-based validation for France and the whole study area using LUCAS data.
ClassesFRA 2019FRA 2018LUCAS 2018
PAUAF1PAUAF1PAUAF1
Others0.970.940.950.980.950.960.960.910.93
Triticeae crops0.890.900.890.890.930.910.720.840.77
Maize0.770.970.860.780.960.860.780.890.83
Rapeseed0.840.880.860.890.910.900.740.840.79
OA0.930.940.89
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Luo, Y.; Zhang, Z.; Zhang, L.; Han, J.; Cao, J.; Zhang, J. Developing High-Resolution Crop Maps for Major Crops in the European Union Based on Transductive Transfer Learning and Limited Ground Data. Remote Sens. 2022, 14, 1809. https://doi.org/10.3390/rs14081809

AMA Style

Luo Y, Zhang Z, Zhang L, Han J, Cao J, Zhang J. Developing High-Resolution Crop Maps for Major Crops in the European Union Based on Transductive Transfer Learning and Limited Ground Data. Remote Sensing. 2022; 14(8):1809. https://doi.org/10.3390/rs14081809

Chicago/Turabian Style

Luo, Yuchuan, Zhao Zhang, Liangliang Zhang, Jichong Han, Juan Cao, and Jing Zhang. 2022. "Developing High-Resolution Crop Maps for Major Crops in the European Union Based on Transductive Transfer Learning and Limited Ground Data" Remote Sensing 14, no. 8: 1809. https://doi.org/10.3390/rs14081809

APA Style

Luo, Y., Zhang, Z., Zhang, L., Han, J., Cao, J., & Zhang, J. (2022). Developing High-Resolution Crop Maps for Major Crops in the European Union Based on Transductive Transfer Learning and Limited Ground Data. Remote Sensing, 14(8), 1809. https://doi.org/10.3390/rs14081809

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