Next Article in Journal
Identifying and Monitoring Gardens in Urban Areas Using Aerial and Satellite Imagery
Previous Article in Journal
Using the InVEST-PLUS Model to Predict and Analyze the Pattern of Ecosystem Carbon storage in Liaoning Province, China
Previous Article in Special Issue
Wavelet Analysis of GPR Data for Belowground Mass Assessment of Sorghum Hybrid for Soil Carbon Sequestration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Resolution National-Scale Mapping of Paddy Rice Based on Sentinel-1/2 Data

1
College of Environmental and Resource Sciences, Zhejiang University, Hangzhou 310058, China
2
Zhejiang Ecological Civilization Academy, Huzhou 313300, China
3
Land Satellite Remote Sensing Application Center, MNR, Beijing 100048, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(16), 4055; https://doi.org/10.3390/rs15164055
Submission received: 12 May 2023 / Revised: 2 July 2023 / Accepted: 2 August 2023 / Published: 16 August 2023
(This article belongs to the Special Issue Digital Farming with Remote Sensing)

Abstract

:
Rice has always been one of the major food sources for human beings, and the monitoring and planning of cultivation areas to maintain food security and achieve sustainable development is critical for this crop. Traditional manual ground survey methods have been recognized as being laborious, while remote-sensing technology can perform the accurate mapping of paddy rice due to its unique data acquisition capabilities. The recently emerged Google Earth Engine (GEE) cloud-computing platform was found to be capable of storing and computing the resources required for the rapid processing of massive quantities of remote-sensing data, thereby revolutionizing traditional analysis patterns and offering unique advantages for large-scale crop mapping. Since the phenology of paddy rice depends on local climatic conditions, and considering the vast expanse of China with its outstanding geospatial heterogeneity, a zoning strategy was proposed in this study to separate the monsoon climate zone of China into two regions based on the Qinling Mountain–Huaihe River Line (Q-H Line), while discrepant basic data and algorithms have been adopted to separately map mid-season rice nationwide. For the northern regions, optical indices have been calculated based on Sentinel-2 images, growth spectral profiles have been constructed to identify phenological periods, and rice was mapped using One-Class Support Vector Machine (OCSVM); for the southern regions, microwave sequences have been constructed based on Sentinel-1 images, and rice was mapped using Random Forest (RF). By applying this methodological system, mid-season rice at 10 m spatial resolution was mapped on the GEE for the entire Chinese monsoon region in 2021. According to the accuracy evaluation coefficients and publicly released local statistical yearbook data, the relative error of the mapped areas in each province was limited to 10%, and the overall accuracy exceeded 85%. The results could indicate that mid-season rice can be mapped more accurately and efficiently on a China-wide scale with relatively few samples based on the proposed zoning strategy and mapping methods. By adjusting the parameters, the time interval for mapping could also be further extended. The powerful cloud-computing competence of the GEE platform was used to map rice on a large spatial scale, and the results can help governments to ascertain the distribution of mid-season rice across the country in a short-term period, which would be well suited to meeting the increasingly efficient and fine-grained decision-making and management requirements.

1. Introduction

As one of the staple foods of human beings, rice accounts for about 12% of the world’s arable land and feeds nearly half of the world’s population [1]. Thus, the visualization of rice cultivation areas has become a critical task. The traditional mapping of crop cultivation extent requires manual field surveys and step-by-step reporting; these processes are usually time consuming and pose difficulties in obtaining the exact distribution of crops [2,3]. Remote sensing can play a unique role in the monitoring, evaluation, and management of agricultural production activities in open areas due to its strengths in relation to its expansive coverage, short revisit periods, fast information acquisition, and relatively low cost [4]. In the 21st century, remote-sensing monitoring has become an effective approach for mapping paddy rice due to the combination of high-resolution data and the cross-fertilization of other novel technologies [5].
An in-depth understanding of rice phenology is a non-negligible prerequisite for deciphering the extent of rice cultivation. Although the phenological characteristics exhibited by rice at different growth stages can vary according to geographical environment in which the rice is grown, the genetically determined characteristics of a species cause the growth stages of rice in different locations to follow certain patterns [6,7,8]. On the other hand, among all food crops, rice is the only one that requires a substantial amount of water at the growth stage, while transplanting is needed [9,10]; therefore, it has become critical to identify rice fields from satellite images enriched with multiple ground objects by examining the dynamic evolution of water–land–vegetation components during an entire growth cycle based on the physiological characteristics of rice during different growth stages.
During the past two decades, remote-sensing-based rice mapping has undergone an evolution from manual interpretation to intelligent monitoring. The earlier visual interpretation methods have been found to be inefficient and costly, while the subsequent machine learning algorithms (e.g., RF and convolutional neural networks (CNNs)) and object-based segmentation algorithms have provided significant improvements [11]. Currently, the recognized features for rice field identification can be classified into two main types: optical vegetation indices and microwave flood signals. Accordingly, two major methodological systems, optics-based mapping and microwave-based mapping, have been developed.
Regarding the optics-based methodology, MODIS, Landsat, and Sentinel-2 data are frequently-used sources [12,13,14,15,16,17]. To better determine the growth states of rice, researchers usually couple rice phenology with multiple vegetation indices to analyze the profile characteristics of the index time series and establish threshold formulas for pixel-by-pixel rice field extraction [18]. The advantages of this method are its simple principles and ease of use, while its prominent disadvantages are its vulnerability to weather factors (e.g., cloud occlusion) and mixed pixels; therefore, the reliability of the corresponding results depends on the number and quality of images at critical growing stages.
As for the microwave-based methodology, the popular satellite data sources employed are predominantly RadarSat-2, ENVISAT-1, and Sentinel-1 [19,20,21,22]. Synthetic aperture radar (SAR) emits microwave signals that can penetrate clouds, theoretically enabling the acquisition of images under various weather conditions. The growth habits of rice dictate that the main grown period is in the rainy season, so SAR data can overcome the disadvantage wherein optical data might be susceptible to cloudy weather [23]. In practice, the time series variation of radar backscatter coefficients during the rice growth phase is generally considered to be the basis for identifying paddy pixels [24,25]. However, insufficient temporal resolution and the speckle noise caused by signal fading have also become essential factors as they limit the application of SAR-based methods [26,27].
Beyond optics-based mapping and microwave-based mapping, the combination of optics and microwaves is a newly developed method that has gradually become a promising class of rice-mapping algorithm in recent years due to technological advances and data expansion. In 2019, Cai et al. mapped rice based on phenology data and time-series Sentinel-2 NDVI and Sentinel-1 SAR images based on an object-based classifier with optimal feature combinations acquired by determining the JBh distance [28]. In the same year, Mansaray et al. explored an approach combining four-source optical satellite (Landsat-8, Sentinel-2A, HJ-1, and Gaofen-1) and microwave satellite (Sentinel-1) images by constructing a spectral index dataset; the approach was used to detect rice fields [29]. In 2020, Fiorillo et al. plotted lowland rice crop areas in the Sédhiou region with an RF algorithm using composite Sentinel-1/2 imagery [30]. The above studies showed that the synergy of optics and microwave images can compensate for the loss of information due to cloudy weather during the critical seasons of paddy rice by incorporating more data and deeper features, thus achieving a generally better degree of mapping accuracy than that achieved using a single data source. In recent years, with the continuous updating of methodologies, the selection of suitable data combination forms and combination algorithms in order to reduce the impact of cloud pollution has become the key research direction.
Remarkable advances in precision agriculture have continuously driven rice mapping in the direction of adapting to the increasingly efficient decision-making needs of governments [31,32]. The advent of geographic cloud-computing platforms may be capable of providing neoteric solutions for large-scale, accurate, and rapid crop mapping. One of the most iconic methods in this regard is the use of the Google Earth Engine, which was developed by Google in collaboration with Carnegie Mellon University, the National Aeronautics and Space Administration (NASA), and the United States Geological Survey (USGS). To date, it is the world’s most advanced cloud-computing platform dedicated to processing remote-sensing imagery and other Earth observation data [33]. With millions of servers worldwide, the GEE is capable of analyzing trillions of images in parallel, leveraging Google’s computing infrastructure and open remote-sensing datasets to store more than 40 years of satellite data online, allowing researchers to diagnose near real-time variations in the Earth’s surface [34,35,36]. In contrast to local server operation, the GEE has distinctly improved the computational efficiency of geospatial big data while providing reliable support for processing extensive remote-sensing images. There have been numerous studies conducted regarding crop mapping using the GEE. In 2015, Lemoine et al. proposed the concept of using the GEE to monitor the acreage and growth status of crops globally [37]. In the same year, Lobell et al. selected soybean fields in the Mid-West U.S. as a study area and proposed a scalable crop-mapping method that did not require ground calibration data after fast pre-processing, cloud masking, or calculation of vegetation indices for all Landsat TM/ETM+ data from 2008 to 2013 with the aid of the GEE [38]. In 2019, Liu et al. used Sentinel-2B and Landsat-8 images to extract typical crop phenology and texture features in Tumut Right Banner County, the Inner Mongolia Autonomous Region, based on the GEE platform [39]. All of the above studies highlighted the advantages of GEE mapping quite adequately.
To date, national-scale rice-mapping studies have been conducted in China, Thailand, Malaysia, and other Southeast Asian countries [7,40,41,42], but their applications have been limited by the coarse spatial resolution (0.5°, 500 m, 250 m, etc.) of the products and the pixel-mixing effect, rendering them incapable of meeting the increasingly fine agricultural management needs in areas with complex topography and smallholder production patterns. Although the GEE may provide a feasible implementation solution for the rapid and accurate identification of rice cultivation ranges at large scales, rice mapping in a huge country like China has still been inadequately studied. Considering the geospatial heterogeneity of rice distribution at the national scale, a zoning strategy for differentiating the monsoon climate zone of China into two regions according to the Q-H Line was proposed. Furthermore, optical data and OCSVM for the northern region, in addition to microwave data and RF for the southern region, with the provincial administrative regions serving as units, have been utilized. As the area of mid-season rice is larger than that of other types of rice, it was selected as the mapping target in this study; thus, a map of the national mid-season rice cultivation extent in 2021, with a high spatial resolution of 10 m, was obtained. The purpose of this study was set to minimize the influence of cloud pollution under the conditions of a shortage of ground-truthing samples while achieving large-scale mapping of rice in the short term and to meet the increasing efficiency and increasingly fine requirements of government departments in administrative management and resource surveys.

2. Materials and Methods

2.1. Study Area

The People’s Republic of China (abbreviated as “China”) is located in eastern Asia and on the west coast of the Pacific Ocean (73°33′E–135°05′E, 3°51′N–53°33′N, Figure 1), possessing a land area of about 9.6 million km2 and a total sea area of about 4.73 million km2 [43]. As the country with the longest history of rice cultivation, China has the world’s largest rice acreage and production values and has also been the largest rice consumer [44,45].
Affected by alternating winter and summer winds over a considerable area, China is the country with the most extensive monsoon influence [46]. In geography, the northern boundary of summer wind influence is usually regarded as the dividing line between monsoon and non-monsoon areas. The specific division of monsoon and non-monsoon zones in China is defined by the “Greater Khingan Mountains—Yin Mountains—Helan Mountains—BaYanKaLa Mountains—Kailas Range” line [47]. The areas to the east and south of the line constitute the monsoon zone, and the areas to the west and north of the line constitute the non-monsoon zone. As shown in Table 1, there have been significant differences between the monsoon and non-monsoon climates in China. The high summer temperatures and superior thermal conditions in the monsoon zone of China provide a favorable environment for rice production [48]. Although the southern region (based on the Q-H Line) of China has served as the main area for rice production, it has also been cultivated in many provinces of the northern region. However, there is almost no distribution of rice fields in the non-monsoon region. In order to meet the demand for large-scale rice mapping, the monsoon region of China is addressed in this study.
The Qinling Mountain–Huaihe River Line (Q-H Line) is the geographic dividing line between the northern and southern regions of China [49]. To the north and south of this line there are distinct differences in natural conditions (especially precipitation and temperature), geographical landscapes, agricultural practices, and the living habits of people [50]. Accordingly, there have been notable differences in the growth phenology, field size, and distribution density of rice found on both sides of this line. As described in Section 1, based on the performed studies, it was shown that the spectral sequences of rice phenological periods constitute a key feature that distinguishes rice from other crops. In this regard, a zoning strategy dividing China into southern and northern regions based on the Q-H Line (Figure 1) was proposed in this study since the phenological periods of rice are fully dependent on regional climatic conditions and cropping practices.

2.2. Datasets

2.2.1. Sentinel-1/2 Data

The optical remote-sensing dataset introduced in this study is the Harmonized Sentinel-2 MSI: Multi-Spectral Instrument, Level-1C; and the microwave SAR dataset is Sentinel-1 SAR GRD: C-band Synthetic Aperture Radar, Ground Range Detected, log scaling. Both datasets were originally provided by the European Space Agency (ESA), are free to access on the GEE platform, and cover the whole year of 2021 and the entire life cycle of all mid-season rice in China. The basic parameters of these two datasets are shown in Table 2.
The Sentinel-1 SAR GRD product available on the GEE is a calibrated product that was pre-processed with the Sentinel-1 toolbox (S1TBX), which contains 4 data acquisition modes (SM, IW, EW, and WV) and two polarization bands (VV and VH). The Interferometric Wide (IW) mode and Vertical-transmit/Horizontal-receive (VH) band of this product have been adopted in this study since IW is the main operation mode of this satellite on land that can achieve stable long-term archive retrieval. It was shown that the VH band is more sensitive to rice growth compared to the VV band [23,51,52]. More details about the Sentinel-1 SAR product can be found in the User Guides on the ESA website (https://sentinel.esa.int/web/sentinel/user-guides/sentinel-1-sar/, accessed on 4 May 2023).
Among the publicly released multi-level products of Sentinel-2, the Level-1C product is the Top-of-Atmosphere (TOA) Reflectance product after atmospheric and geometric corrections, and it has eliminated the artifact interference compared to the Level-2A product [53]. In addition, the HARMONIZED dataset from the Sentinel-2 MSI Level-1C product of the GEE employed in this study was de-clouded, and the reflectance band-dependent offset was removed [54]. More details on the Sentinel-2 MSI product can be found in the User Guides on the ESA website (https://sentinel.esa.int/web/sentinel/user-guides/sentinel-2-msi, accessed on 4 May 2023).

2.2.2. Sample Data

In this study, a national-scale mapping sample database was established on the Google Earth platform based on visual interpretation methods concerning high-definition satellite images (spatial resolution of 0.54 m), including 7984 rice sample points (Figure 2) and 26,573 non-rice sample points (water bodies, buildings, bare land, forests, grasslands, corn, wheat, cotton, soybean, etc.).
The steps taken to establish the specific databases used in this study are as follows:
  • Topographic masking. The slope of entire China was calculated using Digital Elevation Model (DEM) data, and areas with slopes ≥ 15° were removed via the mask because rice would not be cultivatable within the vast majority of these areas. The DEM data introduced in this study were the NASA SRTM Digital Elevation 30 m data, which can provide worldwide digital elevation values with a spatial resolution of 1 arc-second (~30 m) [55]. The voids in this dataset were filled with the support of other open-source data (ASTER GDEM2, GMTED2010 and NED), which can be downloaded from the GEE platform.
  • Reference data overlaying. Maps of cropping patterns in China during the 2015–2021 period were overlaid on the analyzed area after removing the high-slope mask, and the single-season rice layers were selected and extracted from them. Maps of cropping patterns in China from 2015 to 2021 with a spatial resolution of 500 m have been adopted in this study [56]; they are available on the following open source website (https://figshare.com/articles/dataset/Maps_of_cropping_patterns_in_China_during_2015-2020/14936052, accessed on 4 May 2023).
  • Visual interpretation. The single-season rice layer obtained in step 2 was overlaid onto the Google Earth image layer, and the representative “pure” sample points were selected via manual visual interpretation based on prior knowledge (relating to aspects such as aerial photo comparison and the shape, color, and texture characteristics of the fields) (Figure 3). The correctness of the selected points was examined by dragging the time bar to access the historical images. The total number of sample points in each province and the proportions of their number distributions in each prefecture-level city have been determined according to the proportion of rice planted area in the official statistical yearbooks of each region; i.e., the larger the official statistical rice-sowing area, the more sample points are selected in the provinces and prefecture-level cities. The distribution of the number of rice and non-rice samples according to province can be found in the Supplementary Materials (Table S1).
  • Standardized pre-processing. To standardize the sample database for use in the machine learning model, the sample files (Shapefile format) were imported into the GEE platform and converted into Feature Collection format after point picking conducted on Google Earth. For the rice sample points (i.e., points of interest (POI)), the ‘Landcover’ attribute was added, and all of the points have been assigned a value of 1. For the non-rice sample points (i.e., non-points of interest (NPOI)), the ‘Landcover’ attribute has also been added, and all of the points have been assigned a value of 0. Subsequently, each data item in the sample database was buffered in a square area of N m × N m (the value of N depends on the average size of an independent contiguous paddy field in the provincial administrative region), and the final buffered data (in Feature Collection format) were selected as the sample set required for model training and accuracy verification.
  • Random sampling. The collected rice sample points and non-rice sample points were combined into one dataset and disrupted via random number alignment. A total of 70% of the sample points were randomly selected for model training, while 30% of the sample points were selected for accuracy validation.

2.3. Methodology

The overall technical route of this study (shown in Figure 4) involved 5 steps: zoning strategy development, sample database construction, mid-season rice mapping, post-processing, and accuracy evaluation.

2.3.1. Mapping Methodology for Northern Region

  • Optical feature extraction.
For rice feature extraction in the northern region, considering that the occurrence of rainy weather in the area to the north of the Q-H Line would be much less than that to the south and that the weather conditions would also be comparatively stable, the information received by the satellite from optical radiation was considered sufficient and capable of being applied to the optical-based mapping method. Aiming to further eliminate the interference of clouds to achieve better mapping accuracy, a de-clouding algorithm and time series median synthesis were implemented in this study. Specifically, first, a de-clouding function was applied to each Sentinel-2 image for the defined phenological period to remove cloudy pixels. Second, considering that clouds move over time, a median synthesis was performed for the images of the specified phenological period (the specific preprocessing code can be found at https://code.earthengine.google.com/c937f24454144ea8efbfc5e9c8abd396?noload=true, accessed on 2 July 2023). In addition, reliable studies have shown that the vegetation indices calculated based on each waveband of remote-sensing images are essential traits for distinguishing the spectra of target crops from those of other land cover classes [57,58,59,60]; therefore, the Bare Soil Index (BSI), Land Surface Water Index (LSWI), Normalized Difference Vegetation Index (NDVI), Plant Senescence Reflectance Index (PSRI), Enhanced Vegetation Index (EVI), and Green Chlorophyll Vegetation Index (GCVI) have been chosen and calculated for all valid pixels within each northern province throughout 2021 based on Sentinel-2 images after de-clouding pre-processing and median synthesis. The formulas for these indices are shown in Table 3 [61].
  • Phenological stage division.
For the phenological staging in the northern region, a pixel-based phenological feature composite method (Eppf-CM) proposed by Ni et al. was referenced to determine four major phenological stages of rice growth (Bare Soil Stage, Transplanting Stage, Growing Stage, and Mature Stage) [61]. In this study, it is shown that this method can enhance the spectral separability of paddy rice from other crops and make efficient use of information from multiple representative phenological periods of paddy rice to characterize the spectral profiles during the entire life cycle. By analyzing the time series profiles of the mean values of BSI, LSWI, NDVI, and PSRI for mid-season rice sample sites in the northern provinces, four different critical phenological periods were delineated. Specifically, for the bare soil stage, BSI detection was employed since it can effectively characterize the mineral content of soil, with a persistently high BSI and relatively low NDVI and LSWI. For the transplanting stage, LSWI was chosen since numerous studies have shown its excellent performance with respect to identifying flood signals with an LSWI rising rapidly and a BSI falling. For the growing stage, NDVI was chosen as it is a well-recognized vegetation sensitivity index, whose value increased and peaked during this period, with a decreasing BSI. For the mature stage, PSRI was adopted because of its high responsiveness to changes in carotenoid and chlorophyll content, with PSRI rising significantly and NDVI decreasing.
To efficiently evaluate the dynamics of the four optical indices at a given sample point to help better delineate the rice phenological periods, a script for viewing the time series profiles of BSI, LSWI, NDVI, and PSRI at the sample location for 2021 based on the Inspector function of the GEE was formulated. Users can simply run the script to load the images and select the corresponding rice sample points with Inspector to view the index profile from the window on the right side of the interface. The four optical index profiles for a rice sample point in Liaoning Province are shown in Figure 5. The code of this GEE script can be accessed via the following link (https://code.earthengine.google.com/c816a64f1307576bbf6de2b9b217f217, accessed on 4 May 2023).
  • Optical index sequence construction.
Previous studies have shown that EVI is insensitive to soil background and aerosol responses, which are less likely to be saturated in areas with high vegetation cover, and performs better than NDVI in environments with higher humidity [62,63]. Moreover, GCVI can reflect chlorophyll content in leaves through the reflection ratio of near-infrared and green light bands, thus allowing for the accurate determination of the growth dynamics of leaves during transplanting stages [64,65]. Therefore, in order to enhance the identification of the two time periods with prominent growth characteristics, the EVI and GCVI have been calculated separately according to the equation given in Table 3 for each image acquired based on the results of the weather period classification scheme mentioned above. The corresponding time series have also been obtained. Finally, a synthesis operation was performed for the respective phenological periods along with the four optical index time series already calculated in the above step. The six index sequences were then superimposed via band fusion to produce a new feature image for subsequent machine learning model training.
  • Machine learning algorithm—OCSVM.
The machine learning algorithm utilized for rice mapping in the northern region is One-Class Support Vector Machine (OCSVM), which is one of the built-in classifiers of the GEE [66]. This algorithm is a classical anomaly detection algorithm, and its principle is similar to that of a SVM: it assumes that there are two classes of points in the sample space, and it was designed to find a division hyperplane to separate these two classes of samples; the division hyperplane should choose the one with the best generalization ability [67]. Unlike SVM, which is concerned with the binary classification problem, OCSVM has only one classification of interest, which exactly meets the need for anomaly detection, and usually does not focus on anomalous data [68]. Compared with traditional anomaly detection methods, OCSVM can more adequately handle problems containing high-dimensional data, complex datasets, noisy data, and missing data [69].
In addition, it can be empirically gleaned that the rice fields in the northern region are highly contiguous, while the individual field sizes are larger than in the southern regions. Therefore, the sample buffer size (N value) was uniformly set to 10 m for each northern province. Furthermore, 10 m × 10 m buffers were generated for each sample/non-sample point, and the feature images output by the Eppf-CM were then randomly sampled in these buffers and inputted to the OCSVM model. Since these points have been randomly distributed over the entire region of the target ROI, it can be surmised that they adequately represent the variation in spectral information. All the model-training, pixel-based classification, and accuracy validation procedures were performed using the GEE.

2.3.2. Mapping Methodology for Southern Region

  • Microwave feature extraction.
Experiential knowledge has suggested that the range and trend of VH backscatter are promising features with which to distinguish rice from other types of land cover [70,71]. This was demonstrated by the low value of VH in rice pixels during the transplanting stage due to the inundation of water (Figure 6b). Subsequently, due to the growth of rice, the volume scattering of the plants gradually improved, and the backscattering increased significantly during the growing stage. Finally, the backscattering decreased due to the surface scattering from the bare land that remained after the rice had matured and been harvested [72]. Therefore, this “V”-shaped falling–rising feature of VH is the key factor for rice field identification in the cloudy southern region. To visually represent this feature, a script to view the VH backscattering time series profile at the target site over 2021 was developed on the GEE. Figure 6 shows the VH backscatter time series profiles for various land cover classes including mid-season rice in a region of Anhui Province, where the typical “V” shaped fluctuations can be seen in the red circle in Figure 6a; the figure also shows that the time coincided with the local rice-transplanting stage. The code for this GEE script can be accessed via the following link (https://code.earthengine.google.com/6f95a8a252182c4d736bcaae47bf2024, accessed on 4 May 2023).
  • Phenological stage division.
The complex topography, variable climate, and geographically distinctive water and heat conditions in the southern region resulted in the prominent spatial heterogeneity of rice phenology, and an analysis conducted based on characteristic time series profiles would no longer be effective. And there has no study been found which could identify rice phenology accurately and precisely for the entire southern region of China. Therefore, this study referred to the publicly available calendar of rice phenology information on the official website of the Ministry of Agriculture and Rural Affairs of the People’s Republic of China (Table 4). This calendar provides multi-stage phenological information on the growth of the main crops throughout China in a geographically zoned manner with a decadal (near 10 d) time scale.
  • VH sequence construction.
The prerequisite step for conducting rice extraction in the southern region is to construct SAR image sequences based on microwave images and phenological periods for pixel-by-pixel classification. The specific process of constructing SAR image sequences is as follows: Firstly, the Sentinel-1 images should be screened on the GEE according to the phenological calendar. Secondly, the median values of the images in each corresponding phenological period should be synthesized individually, while the mosaic operation should be applied to multiple images at the same location in the specified area. Finally, the synthesized images should be the SAR image units per phenological period, which are then used for the subsequent training of the machine learning model.
  • Machine learning algorithm—RF.
The machine learning algorithm that was utilized for rice mapping in the southern region is Random Forest (RF). RF is a classifier that contains multiple decision trees [73]. It can generate a new set of training samples by repeatedly randomly selecting n samples from the original training sample set through a bootstrap resampling technique; then, it generates m decision trees to form a random forest, and the classification result of the new data is determined by the number of votes formed by the classification trees [74,75]. In essence, RF is a modification of the decision tree algorithm wherein multiple decision trees are combined, and the establishment of each tree depends on independently drawn samples. Many relevant studies have shown that the RF algorithm presents high classification accuracy and can successfully handle high data dimensions, thus indicating its high computational speed and robustness toward overfitting [28,76,77].
In addition, our preliminary observations of the images have revealed that the distribution pattern of the rice fields in the southern region varies greatly, presenting both large fields with high contiguity and small, fragmented, and scattered fields; therefore, the sample buffer size (N value) for each southern province was optimized according to the actual situation of different places after several experiments. After N m × N m buffers were generated for each sample/non-sample point, the constructed SAR sequence images were randomly sampled in these buffers and inputted to the RF model (the number of trees in the RF model was set to 500, and the other parameters were set to the default values).

2.3.3. Post-Processing Method

Once the preliminary results (in GeoTIFF format) were obtained, two important challenges arose at the same time, which could not be ignored: the impulse noise inherent in the Sentinel series data and the hole phenomenon caused by the lack of consideration of spatial correlation in pixel-by-pixel classification methods [18,61]. Thus, the extracted maps have been imported to local PCs, and ArcGIS software was introduced to perform spatial filtering operations in order to produce more reliable classification results. These operations were incorporated to alleviate noise, repair holes, and minimize misclassification errors at neighboring edges of different features. Specifically, the extraction results have been optimized using the Majority Filter (MF) in the ArcMap10.7 toolbox.
The basic principle of MF is to replace the pixels in the raster based on the mode of neighboring values (Figure 7) [78]. The optional parameters of this tool include “number_neighbors” (which determines the number of neighboring pixels in the filter kernel, whose default value is “FOUR”, which means that the kernel of the filter will be the four direct (orthogonal) neighbors to the present cell) and “majority_definition” (which specifies the number of neighboring pixels that must have the same value, the default value of which is “MAJORITY”, which means that a majority of cells must have the same value and be contiguous). In this study, the two parameters have been set as default values.
As shown in Figure 8, the output raster was stabilized, some of the scattered misclassified rice field pixels were removed, and fine voids within some of the contiguous rice field patches were filled after running the plurality filter several times, while the result would then be selected as the final outcome of rice mapping.

2.3.4. Mapping Accuracy Validation

Due to the lack of publicly available China-wide high-resolution (≤10 m) mid-season rice-tagging data, it was not feasible to reliably evaluate the accuracy of the mapping results at the individual pixel scale. Therefore, the visual observation method, the confusion matrix-based classification coefficient method, and the statistical yearbook area comparison method have been carried out for multi-perspective accuracy verification of the results in this study.
  • Visual observation method. Visual observation is a direct method for evaluating mapping performance. It involves overlaying classification results onto high-definition remote-sensing reference images from Google Earth and matching them.
  • Classification coefficient method. By invoking the built-in confusion matrix calculation function of the GEE, various accuracy evaluation coefficients can be quantified, such as Producer’s Accuracy (PA), User ‘s Accuracy (CA), Overall Accuracy (OA), Kappa Coefficient (KC), etc., and the results can be quantitatively tested. The specific formulas for the metrics mentioned above can be found in the Supplementary Materials (Table S2).
  • Statistical yearbook comparison method. The statistical area comparison method is designed to evaluate the classification results from the perspective of quantitative statistics by counting the number of rice pixels within the mapping results of each province before converting them into rice field areas according to the raster size and then comparing them with the rice-sowing areas noted in the statistical yearbooks of each province to acquire the relative errors. The statistical yearbooks of each province can be found on the official websites of the statistical offices of the local governments.

3. Results

3.1. Mapping Results for Northern Region

For the northern region, three northeastern provinces (Heilongjiang, Jilin, and Liaoning), which have been highly representative and have the highest statistical area of rice cultivation among the northern provinces, have been selected for the presentation of the mid-season rice-mapping results and accuracy evaluation due to the space limitation of the paper.

3.1.1. Overall Distribution

The overall distribution of mid-season rice in the three northeastern provinces of China in 2021 is illustrated in Figure 9. Among these provinces, mid-season rice in Heilongjiang Province (Figure 9a) presented the widest distribution and the highest planting density and contiguity, which was spatially concentrated in the central and eastern flat terrain areas and broadly distributed in the west. The distribution of rice in the south turned out to be relatively small, and there are almost no large rice areas in the northern mountainous regions. In addition, most of the mid-season rice in Heilongjiang was observed to be distributed along rivers, which could be related to the water-loving characteristic of rice. Most of the mid-season rice in Jilin Province (Figure 9b) is distributed in the northwestern and central parts, and the field contiguity was determined to be the lowest among the three provinces, which is consistent with the topographic characteristics of Jilin Province (higher in the southeast and lower in the northwest). In Liaoning Province (Figure 9c), the topography was observed to be higher in the east and west and lower in the middle; therefore, most of the mid-season rice distribution is concentrated in the middle, which is the most morphologically concentrated but smallest area among the three provinces.

3.1.2. Local Visual Comparison

The partial results of rice mapping in the three northeastern provinces and their comparison with real high-definition images are shown in Figure 10. The visual matching of the three provinces led to a largely consistent discovery: the mid-season rice fields in the northeast have been found to be more regular in shape, larger in size, and higher in contiguity. The optical index mapping method based on Sentinel-2 data and the OCSVM algorithm showed excellent visual performance and could clearly distinguish mid-season rice from other land categories (buildings, water bodies, bare land, and other crops). The model extraction resulted in clear edge segmentation and the clear identification of narrow roads and bridges crisscrossing the field. However, many tiny holes were detected inside some of the individual contiguous paddy patches when zooming in for closer observation. These holes somehow matched the distribution pattern of noise in the Sentinel-2 data.

3.1.3. Accuracy Evaluation

The accuracy of the results of mid-season rice mapping in the three northeastern provinces based on classification coefficients and statistical yearbook areas is reported in Table 5. With a uniform buffer size of 10 m, the extraction results of the three northeastern provinces showed satisfactory classification accuracy. The OA in all three provinces exceeded 90%, while the KC values exceeded 85%. On the other hand, the relative errors with the rice sown area data in the official statistical yearbook were limited to 4%.

3.2. Mapping Results for Southern Region

For the southern region, seven provinces (Hunan, Hubei, Jiangxi, Anhui, Jiangsu, Zhejiang, and Guangxi) with highly representative and statistically large areas of rice cultivation have been selected for the presentation of mid-season rice-mapping results and accuracy evaluation due to the space limitation of the paper.

3.2.1. Overall Distribution

The general distribution of mid-season rice in seven provinces in southern China in 2021 is shown in Figure 11. The pattern of divergence in rice distribution with variations in topography, temperature, and hydrology is clearly reflected. It is noticeable that the distribution of mid-season rice in Hunan Province (Figure 11a) is widespread and dispersed with cultivation in all corners of the province, and the relatively most-concentrated area is located in and around the Dongting Lake basin. In Anhui Province (Figure 11b), mid-season rice varied greatly from north to south, with high planting density in the Huaihe Plain in the north where continuous rice fields have surrounded almost all cities and rivers, and scattered rice fields were found in the Jianghuai Tableland hills in the south. A relatively vast area of mid-season rice was observed in Hubei Province (Figure 11c), which has mostly clustered in the central lowland area surrounded by mountains. Jiangsu Province (Figure 11d) was identified as the most ideal area for the growth of rice as it is the province with the highest proportion of plain terrain (86.89%) nationwide and it has a sub-tropical monsoon climate; therefore, mid-season rice is widely distributed. Most of the mid-season rice in Jiangxi Province (Figure 11e) is concentrated in and around the Poyang Lake basin in the north, with a scattered presence in the mountainous hilly areas in the south. In Zhejiang Province (Figure 11f), the geomorphological characteristics can be summarized as “seven mountains, one water and two fields”; thus, the rice cultivation area was found to be relatively small, and most of it is clustered in the Hangjiahu Plain in the northeast. In the Guangxi Zhuang Autonomous Region (Figure 11g), the vast low mountains and hills precipitated a complex distribution of mid-season rice, which was mostly found in river valleys.

3.2.2. Local Visual Comparison

The partial results of rice mapping in the seven southern provinces and their comparison with the real high-definition images are shown in Figure 12. It can be found that compared with the northern region, the mid-season rice fields in the south have fragmented shapes and strong spatial heterogeneity of field sizes. The SAR sequence-mapping method based on Sentinel-1 data and the RF algorithm provided relatively effective differentiation of land classes, especially the rejection of water bodies. It is significant that the mapping procedure also showed many misclassifications, i.e., buildings, water bodies, trees, and other crops were misclassified as rice, or rice was identified as other land classes. Some irregularly randomly distributed multi-scale patches were also detected in the results. The uncertainties in classification caused by these phenomena were most evident in provinces with a complex topography and a small statistical area of rice, such as Zhejiang (Figure 12f) and Guangxi (Figure 12g).

3.2.3. Accuracy Evaluation

The accuracy of the results of mid-season rice mapping based on classification coefficients and statistical yearbook areas in the seven southern provinces are reported in Table 6. As mentioned above, the average size of rice fields varied substantially from province to province in the southern region; therefore, the relatively optimal buffer size (N value) was determined after multiple trials in each province in this study. The OA of all seven southern provinces exceeded 85%, while the KC values exceeded 70% and the F-scores were ≥ 75%; in the meantime, the relative errors in the rice sowing area data in the official statistical yearbook were limited to 10% and lower after majority filtering. With regard to this bias, some studies have found that a large number of rice fields might be abandoned while the paddy fields are converted to dry land (i.e., planted with other crops or vegetables) in areas with complex topography after field investigations due to the difficulty of implementing mechanized rice cultivation [11]. Regarding the official village-based statistics on rice cultivation area, many farmers might not have reported the actual rice cultivation area of their land, which could be the main factor leading to the deviation between the estimation of the statistical and remote-sensing data [79,80].

4. Discussion

4.1. Comparative Analysis of Zoning-Mapping Results

As can be observed in the graphs and tables in the Results section, the results of all three accuracy evaluation methods have indicated that the accuracy of rice mapping in the three northeastern provinces was better than that in the seven southern provinces. After a long-term survey of satellite high-resolution images, it was determined that the areas to the south and north of the Q-H Line varied in terms of topography, phenology, and agricultural activities to different degrees, leading to remarkable discrepancies in the adopted mapping methods and mapping effects. These differences are manifested in the following three aspects.
In terms of topography, the distribution of paddy fields in the lower terrain in the north was more regular and coherent than that in the south since there were fewer hilly mountains, which eased visual interpretation, and the samples were more reliable. On the contrary, the distribution of rice fields in several provinces in the southern region was more fragmented and hybrid in nature, and the mixing effect of pixels was prominent, which made the visual interpretation of rice sample points harder while rendering the samples less stable.
In terms of phenology, the three northeastern provinces are located in the middle temperate zone, and the rice maturation system is mono-annual; therefore, their phenological characteristics are obvious, and it is convenient to interpret the four key phenological stages from the spectral profiles of rice [81]. The local rice maturation system was found to be relatively diverse due to the water and heat conditions in the seven southern provinces. There was a mixture of mono-annual, bi-annual, and even tri-annual maturation, which resulted in flexible phenological characteristics and made it difficult to determine the reasonable boundaries of the phenological periods.
In terms of farming activities, the relatively simple rice maturation scheme in the northern region has allowed for only a single-phase planting process, making the inter-annual rice field distribution changes more stable and the samples more consistent over time. Even if the year of mapping and the year of reference image did not match, the actual location of the sample did not change drastically. In contrast, under the complex rice maturation scheme in the southern region, farmers tended to adopt complex planting measures such as crop rotation and intercropping [82], which led to the difficulty in visual the interpretation of the rice sample points and the lack of distinct V-shaped characteristics of the microwave sequences.
In addition, the identified rice fields were usually accompanied by a certain number of tiny holes within them, resulting in a slightly underestimated extracted area compared to the official statistical area in the rice-mapping results of the provinces in the northern region. Many irregular patches were observed in the extraction results of the southern region. This could be the result of the noise problem in the Sentinel series data source itself and the adopted pixel-by-pixel classification method’s failure to account for neighborhood information mining [83,84]. Both of these problems could be addressed by implementing morphological methods to filter out spurious pixels or repair voids, but there has not yet been a scientific solution that can completely solve them [84,85,86].

4.2. Features and Uncertainties of This Study

In response to the need for increasingly efficient and fine-grained arable land resource survey and management methods and considering the lack of national-scale rice mapping research, a rational and innovative zoning strategy based on the Q-H Line was proposed in this study in view of the geospatial heterogeneity of rice distribution and the corresponding phenology nationwide. On this basis, a high-resolution rapid mapping of mid-season rice in the entire monsoon region of China was conducted on the GEE cloud-computing platform, which possesses a massive amount of remote-sensing data storage and a built-in machine learning classifier. Furthermore, the large-scale mapping of rice in a short period of time under the condition of a relative scarcity of ground truth samples was achieved by minimizing the influence of cloud pollution, leading to the acquisition of encouraging results, which could provide valuable references for similar studies and applications in the future. In addition, the extension of the mapping years by adjusting the time span parameters was also achieved, which could provide a possible method of producing long time series and large-scale high-resolution rice datasets.
However, there are some obvious shortcomings in this study that should be overcome in future in-depth studies. Firstly, regarding the choice between accuracy and efficiency, which is often faced in remote-sensing image classification, efficiency was favored in this study, resulting in a lack of quantity and quality of the selected sample points. This limitation could have resulted in unstable mapping results and accuracy evaluation. Secondly, the classifiers adopted might be the most basic machine learning classifiers that have been applied for a long time, and the parameters of the algorithms might be too roughly designed and not sufficiently adapted to the local conditions. Thirdly, the identification of rice phenological stages might not have been precise enough, as the entire life cycle of rice was only divided into four stages in the northern region, while a simple statistically based phenological calendar was adopted in the southern region, which could also have negatively affected extraction accuracy. Fourthly, in the preprocessing steps, pixels with slopes greater than 15 degrees were removed in order to exclude terrains where rice may not be planted to enhance the efficiency of RF, while some of the paddy rice growing in the terraces were inevitably omitted, thus influencing the mapping accuracy. Moreover, the insufficient number of effective images, the low temporal resolution of the Sentinel-1 data, the inherent border noise of the Sentinel-1/2 data, and the single rice species selected for extraction would also deserve great attention.
Finally, another crucial factor of uncertainty affecting the mapping accuracy is the size of the buffer zone (i.e., the value of N) of the paddy field samples. In the Methodology section, it was stated that N needs to be initially adjusted according to the actual average size of the paddy fields in each region (whose optimal value can only be determined through undertaking several tests). However, in practice, this process may be time-consuming and require a combination of field research and empirical knowledge. If the selected N value is larger than the average size of the local paddy fields, the samples input to the classifier will be mixed with features of other land use types, and the sample purity will be reduced, resulting in a greater predicted area due to the misidentification of non-paddy fields. If the selected N value is smaller than the average size of local paddy fields, the sample features learned by the classifier will be insufficient, and the feature heterogeneity within the paddy fields will be ignored; thus, part of the paddy fields will be ignored, some paddy fields are not fully included in the classification results, resulting in a smaller predicted area. Therefore, in future studies, it’s suggested to develop an adaptive buffer size selection method, i.e., to find the optimal N values for various areas using a cyclic feedback tuning method.

4.3. Future Research Directions

In response to the shortcomings mentioned above and the common shortcomings of similar studies, futural exploratory research might be conducted around the following four aspects.
  • Fusion of multi-source data. In the future, researchers could try to fuse multi-source remote-sensing data (Sentinel-1/2, Landsat, MODIS, UAV, etc.) to improve the feature index system and enhance classification effects. Some experimental studies have been conducted, for which more reliable results have been obtained, which could provide new insights for mapping crop cultivation spatial distribution [87,88,89,90].
  • Synergy of cloud computing and deep learning. The deep learning algorithms that have emerged in the last decade have improved rice extraction on complex surfaces and in fragmented landscapes by building a moderate number of neuronal computation nodes and multi-layer operational hierarchies with higher classification accuracy compared to traditional machine learning algorithms [91,92]. However, more complex model structures would also require better hardware performance, longer training times, and a larger number of data labels [93]. Therefore, researchers could attempt to combine the advantages of the high accuracy of deep learning and the high efficiency of the GEE platform to construct local deep learning models on this cloud-computing platform.
  • Enhancement of post-processing. For the optimization of the results after GEE derivation, other morphological methods could be considered besides performing majority filtering; different processing window sizes (dynamic windows) could also be experimented with to further filter out noise and repair voids, thus compensating for the deficiencies of the underlying data and algorithms [94,95].
  • Development of sample-poor mapping technologies. A lack of accurate sample points might become the norm for future large-scale crop mapping efforts and a bottleneck for technological progress in related fields [76,96,97]. Therefore, the development of deep learning classification algorithms with stronger autonomous learning capacity via coupling deep learning frameworks (e.g., Tensor Flow) or large artificial intelligence models (e.g., the ChatGPT) with emerging technologies would be an important way in which to improve the accuracy of large-scale rice extraction models and their generalization capability [98,99,100].

5. Conclusions

In this study, a national-scale mid-season rice zonal mapping strategy was proposed under the conditions of the relative scarcity of ground-truthing sample data given the spatial heterogeneity of rice phenology across China. Secondly, based on Sentinel-1/2 data and machine learning algorithms, a 2021 Chinese monsoon zone mid-season rice cultivation extent dataset was created on the GEE platform, which contains maps of mid-season rice fields for 28 provincial administrative regions in China. The data file format is GeoTIFF, with a spatial resolution of 10 m, and the geo-reference system used was WGS84 (EPSG:4326). Finally, with the help of multiple methods, the accuracy evaluation results have shown that the relative errors of the mapped areas of mid-season rice in each province were limited to 10% or less compared with the statistical yearbook, while the overall accuracy obtained from the confusion matrix exceeded 85%. Restricted by the availability of samples and geographical complexity, the accuracy of mapping in the south turned out to be generally lower than that in the north. Additionally, the extension of the mapping years by adjusting the time period parameters can be achieved. The inspiring results obtained in this study provide a solution that can be referenced to produce long-time-series, large-scale and high-resolution spatial datasets of rice cultivation extent. In the future, efforts could be made toward the fusion of multi-source remote-sensing data, the fusion of cloud-computing platforms and deep learning models, the reinforcement of the post-processing of extraction results, and the development of intelligent mapping methods in the absence of abundant samples in order to realize further synergistic progress in accuracy and efficiency in research and applications related to large-scale crop mapping.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/rs15164055/s1, Table S1: Allocation of rice selection points by provincial administrative units within the whole China monsoon region.; Table S2: Formula for calculating classification coefficients based on confusion matrix.

Author Contributions

Conceptualization, C.H. and J.D.; methodology, C.H.; data curation, C.H. and J.Z.; investigation, J.Z.; validation, C.H.; writing—original draft preparation, C.H.; writing—review and editing, J.D. and P.L.; visualization, C.H.; supervision, S.Y. and A.L.; project administration, J.D.; funding acquisition, J.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Outsourcing Project of the Center for Remote Sensing Application of Land and Satellite, Ministry of Natural Resources, the Consulting Project of Chinese Academy of Engineering (2023-30-13), and the National Key Research and Development Program of China (2020YFC1807501).

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to projects that funded this study has not yet been fully recognized as complete, thus premature disclosure of the data would potentially create a privacy or moral hazard.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Elert, E. Rice by the Numbers: A Good Grain. Nature 2014, 514, S50–S51. [Google Scholar] [CrossRef] [PubMed]
  2. Zhao, R.; Li, Y.; Ma, M. Mapping Paddy Rice with Satellite Remote Sensing: A Review. Sustainability 2021, 13, 503. [Google Scholar] [CrossRef]
  3. Mosleh, M.; Hassan, Q.; Chowdhury, E. Application of Remote Sensors in Mapping Rice Area and Forecasting Its Production: A Review. Sensors 2015, 15, 769–791. [Google Scholar] [CrossRef] [PubMed]
  4. Mahlayeye, M.; Darvishzadeh, R.; Nelson, A. Cropping Patterns of Annual Crops: A Remote Sensing Review. Remote Sens. 2022, 14, 2404. [Google Scholar] [CrossRef]
  5. Orynbaikyzy, A.; Gessner, U.; Conrad, C. Crop Type Classification Using a Combination of Optical and Radar Remote Sensing Data: A Review. Int. J. Remote Sens. 2019, 40, 6553–6595. [Google Scholar] [CrossRef]
  6. Kuenzer, C.; Knauer, K. Remote Sensing of Rice Crop Areas. Int. J. Remote Sens. 2013, 34, 2101–2139. [Google Scholar] [CrossRef]
  7. Frolking, S.; Qiu, J.; Boles, S.; Xiao, X.; Liu, J.; Zhuang, Y.; Li, C.; Qin, X. Combining Remote Sensing and Ground Census Data to Develop New Maps of the Distribution of Rice Agriculture in China: Paddy Rice Cropland Maps for China. Glob. Biogeochem. Cycles 2002, 16, 38-1–38-10. [Google Scholar] [CrossRef]
  8. Peña-Arancibia, J.L.; Mahboob, M.G.; Islam, A.T.; Mainuddin, M.; Yu, Y.; Ahmad, M.D.; Ibn Murad, K.F.; Saha, K.K.; Hossain, A.; Moniruzzaman, M.; et al. The Green Revolution from Space: Mapping the Historic Dynamics of Main Rice Types in One of the World’s Food Bowls. Remote Sens. Appl. Soc. Environ. 2021, 21, 100460. [Google Scholar] [CrossRef]
  9. 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]
  10. Salmon, J.M.; Friedl, M.A.; Frolking, S.; Wisser, D.; Douglas, E.M. Global Rain-Fed, Irrigated, and Paddy Croplands: A New High Resolution Map Derived from Remote Sensing, Crop Inventories and Climate Data. Int. J. Appl. Earth Obs. Geoinf. 2015, 38, 321–334. [Google Scholar] [CrossRef]
  11. Xiao, X.; Boles, S.; Liu, J.; Zhuang, D.; Frolking, S.; Li, C.; Salas, W.; Moore, B. Mapping Paddy Rice Agriculture in Southern China Using Multi-Temporal MODIS Images. Remote Sens. Environ. 2005, 95, 480–492. [Google Scholar] [CrossRef]
  12. Gumma, M.K. Mapping Rice Areas of South Asia Using MODIS Multitemporal Data. J. Appl. Remote Sens 2011, 5, 053547. [Google Scholar] [CrossRef]
  13. Sun, H.; Huang, J.; Huete, A.R.; Peng, D.; Zhang, F. Mapping Paddy Rice with Multi-Date Moderate-Resolution Imaging Spectroradiometer (MODIS) Data in China. J. Zhejiang Univ. Sci. A 2009, 10, 1509–1522. [Google Scholar] [CrossRef]
  14. Shew, A.M.; Ghosh, A. Identifying Dry-Season Rice-Planting Patterns in Bangladesh Using the Landsat Archive. Remote Sens. 2019, 11, 1235. [Google Scholar] [CrossRef]
  15. Dong, J.; Xiao, X.; Kou, W.; Qin, Y.; Zhang, G.; Li, L.; Jin, C.; Zhou, Y.; Wang, J.; Biradar, C.; et al. Tracking the Dynamics of Paddy Rice Planting Area in 1986–2010 through Time Series Landsat Images and Phenology-Based Algorithms. Remote Sens. Environ. 2015, 160, 99–113. [Google Scholar] [CrossRef]
  16. Choudhary, K.; Shi, W.; Dong, Y.; Paringer, R. Random Forest for Rice Yield Mapping and Prediction Using Sentinel-2 Data with Google Earth Engine. Adv. Space Res. 2022, 70, 2443–2457. [Google Scholar] [CrossRef]
  17. Liu, L.; Xiao, X.; Qin, Y.; Wang, J.; Xu, X.; Hu, Y.; Qiao, Z. Mapping Cropping Intensity in China Using Time Series Landsat and Sentinel-2 Images and Google Earth Engine. Remote Sens. Environ. 2020, 239, 111624. [Google Scholar] [CrossRef]
  18. Pang, J.; Zhang, R.; Yu, B.; Liao, M.; Lv, J.; Xie, L.; Li, S.; Zhan, J. Pixel-Level Rice Planting Information Monitoring in Fujin City Based on Time-Series SAR Imagery. Int. J. Appl. Earth Obs. Geoinf. 2021, 104, 102551. [Google Scholar] [CrossRef]
  19. Nguyen, D.B.; Wagner, W. European Rice Cropland Mapping with Sentinel-1 Data: The Mediterranean Region Case Study. Water 2017, 9, 392. [Google Scholar] [CrossRef]
  20. Yang, H.; Pan, B.; Li, N.; Wang, W.; Zhang, J.; Zhang, X. A Systematic Method for Spatio-Temporal Phenology Estimation of Paddy Rice Using Time Series Sentinel-1 Images. Remote Sens. Environ. 2021, 259, 112394. [Google Scholar] [CrossRef]
  21. Ho Tong Minh, D.; Ndikumana, E.; Baghdadi, N.; Courault, D.; Hossard, L. Applying Deep Learning for Agricultural Classification Using Multitemporal SAR Sentinel-1 for Camargue, France. In Proceedings of the Image and Signal Processing for Remote Sensing XXIV; Bruzzone, L., Bovolo, F., Benediktsson, J.A., Eds.; SPIE: Berlin, Germany, 2018; p. 39. [Google Scholar]
  22. Zhang, Y.; Wang, C.; Wu, J.; Qi, J.; Salas, W.A. Mapping Paddy Rice with Multitemporal ALOS/PALSAR Imagery in Southeast China. Int. J. Remote Sens. 2009, 30, 6301–6315. [Google Scholar] [CrossRef]
  23. Xu, S.; Zhu, X.; Chen, J.; Zhu, X.; Duan, M.; Qiu, B.; Wan, L.; Tan, X.; Xu, Y.N.; Cao, R. A Robust Index to Extract Paddy Fields in Cloudy Regions from SAR Time Series. Remote Sens. Environ. 2023, 285, 113374. [Google Scholar] [CrossRef]
  24. Wei, P.; Chai, D.; Huang, R.; Peng, D.; Lin, T.; Sha, J.; Sun, W.; Huang, J. Rice Mapping Based on Sentinel-1 Images Using the Coupling of Prior Knowledge and Deep Semantic Segmentation Network: A Case Study in Northeast China from 2019 to 2021. Int. J. Appl. Earth Obs. Geoinf. 2022, 112, 102948. [Google Scholar] [CrossRef]
  25. Wei, P.; Chai, D.; Lin, T.; Tang, C.; Du, M.; Huang, J. Large-Scale Rice Mapping under Different Years Based on Time-Series Sentinel-1 Images Using Deep Semantic Segmentation Model. ISPRS J. Photogramm. Remote Sens. 2021, 174, 198–214. [Google Scholar] [CrossRef]
  26. Setiyono, T.; Quicho, E.; Gatti, L.; Campos-Taberner, M.; Busetto, L.; Collivignarelli, F.; García-Haro, F.; Boschetti, M.; Khan, N.; Holecz, F. Spatial Rice Yield Estimation Based on MODIS and Sentinel-1 SAR Data and ORYZA Crop Growth Model. Remote Sens. 2018, 10, 293. [Google Scholar] [CrossRef]
  27. Xiao, W.; Xu, S.; He, T. Mapping Paddy Rice with Sentinel-1/2 and Phenology-, Object-Based Algorithm—A Implementation in Hangjiahu Plain in China Using GEE Platform. Remote Sens. 2021, 13, 990. [Google Scholar] [CrossRef]
  28. Cai, Y.; Lin, H.; Zhang, M. Mapping Paddy Rice by the Object-Based Random Forest Method Using Time Series Sentinel-1/Sentinel-2 Data. Adv. Space Res. 2019, 64, 2233–2244. [Google Scholar] [CrossRef]
  29. Mansaray, L.R.; Yang, L.; Kabba, V.T.S.; Kanu, A.S.; Huang, J.; Wang, F. Optimising Rice Mapping in Cloud-Prone Environments by Combining Quad-Source Optical with Sentinel-1A Microwave Satellite Imagery. GISci. Remote Sens. 2019, 56, 1333–1354. [Google Scholar] [CrossRef]
  30. Fiorillo, E.; Di Giuseppe, E.; Fontanelli, G.; Maselli, F. Lowland Rice Mapping in Sédhiou Region (Senegal) Using Sentinel 1 and Sentinel 2 Data and Random Forest. Remote Sens. 2020, 12, 3403. [Google Scholar] [CrossRef]
  31. Hoffmann Silva Karp, F. Within-Field Yield Prediction for Sugarcane and Rice Focused on Precision Agriculture Applications. Master’s Thesis, Louisiana State University and Agricultural and Mechanical College, Baton Rouge, LA, USA, 2020. [Google Scholar]
  32. Wang, Y.-P.; Chen, S.-H.; Chang, K.-W.; Shen, Y. Identifying and Characterizing Yield Limiting Factors in Paddy Rice Using Remote Sensing Yield Maps. Precis. Agric. 2012, 13, 553–567. [Google Scholar] [CrossRef]
  33. Tamiminia, H.; Salehi, B.; Mahdianpari, M.; Quackenbush, L.; Adeli, S.; Brisco, B. Google Earth Engine for Geo-Big Data Applications: A Meta-Analysis and Systematic Review. ISPRS J. Photogramm. Remote Sens. 2020, 164, 152–170. [Google Scholar] [CrossRef]
  34. Zhao, Q.; Yu, L.; Li, X.; Peng, D.; Zhang, Y.; Gong, P. Progress and Trends in the Application of Google Earth and Google Earth Engine. Remote Sens. 2021, 13, 3778. [Google Scholar] [CrossRef]
  35. Mutanga, O.; Kumar, L. Google Earth Engine Applications. Remote Sens. 2019, 11, 591. [Google Scholar] [CrossRef]
  36. Yao, J.; Wu, J.; Xiao, C.; Zhang, Z.; Li, J. The Classification Method Study of Crops Remote Sensing with Deep Learning, Machine Learning, and Google Earth Engine. Remote Sens. 2022, 14, 2758. [Google Scholar] [CrossRef]
  37. Lemoine, G.; Leo, O. Crop Mapping Applications at Scale: Using Google Earth Engine to Enable Global Crop Area and Status Monitoring Using Free and Open Data Sources. In Proceedings of the 2015 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Milan, Italy, 26–31 July 2015; IEEE: Piscataway, NJ, USA; pp. 1496–1499. [Google Scholar]
  38. Lobell, D.B.; Thau, D.; Seifert, C.; Engle, E.; Little, B. A Scalable Satellite-Based Crop Yield Mapper. Remote Sens. Environ. 2015, 164, 324–333. [Google Scholar] [CrossRef]
  39. Liu, Y. Research on Remote Sensing Information Extraction of Typical Crops in Arid Areas Based on Cloud Computing Platform. Master’s Thesis, Yunnan Normal University, Yunnan, China, 2019. [Google Scholar]
  40. Xiao, X.; Boles, S.; Frolking, S.; Li, C.; Babu, J.Y.; Salas, W.; Moore, B. Mapping Paddy Rice Agriculture in South and Southeast Asia Using Multi-Temporal MODIS Images. Remote Sens. Environ. 2006, 100, 95–113. [Google Scholar] [CrossRef]
  41. Song, P.; Mansaray, L.R.; Huang, J.; Huang, W. Mapping Paddy Rice Agriculture over China Using AMSR-E Time Series Data. ISPRS J. Photogramm. Remote Sens. 2018, 144, 469–482. [Google Scholar] [CrossRef]
  42. Wei, S.; Zhang, H.; Wang, C.; Xu, L.; Wu, F.; Zhang, B. Large-Scale Rice Mapping of Thailand Using Sentinel-1 Multi-Temporal SAR Data. In Proceedings of the 2019 SAR in Big Data Era (BIGSARDATA), Beijing, China, 5–6 August 2019; IEEE: Piscataway, NJ, USA; pp. 1–4. [Google Scholar]
  43. Archive. Available online: http://english.www.gov.cn/archive/ (accessed on 8 May 2023).
  44. Peng, S.; Tang, Q.; Zou, Y. Current Status and Challenges of Rice Production in China. Plant Prod. Sci. 2009, 12, 3–8. [Google Scholar] [CrossRef]
  45. Ma, G.; Yuan, L. Hybrid Rice Achievements, Development and Prospect in China. J. Integr. Agric. 2015, 14, 197–205. [Google Scholar] [CrossRef]
  46. He, J.; Ju, J.; Wen, Z.; Lü, J.; Jin, Q. A Review of Recent Advances in Research on Asian Monsoon in China. Adv. Atmos. Sci. 2007, 24, 972–992. [Google Scholar] [CrossRef]
  47. Qian, W.; Ding, T.; Hu, H.; Lin, X.; Qin, A. An Overview of Dry-Wet Climate Variability among Monsoon-Westerly Regions and the Monsoon Northernmost Marginal Active Zone in China. Adv. Atmos. Sci. 2009, 26, 630–641. [Google Scholar] [CrossRef]
  48. Andriesse, W.; Fresco, L.O. A Characterization of Rice-Growing Environments in West Africa. Agric. Ecosyst. Environ. 1991, 33, 377–395. [Google Scholar] [CrossRef]
  49. Yuan, Q.; Wu, S.; Dai, E.; Zhao, D.; Zhang, X.; Ren, P. Spatio-Temporal Variation of the Wet-Dry Conditions from 1961 to 2015 in China. Sci. China Earth Sci. 2017, 60, 2041–2050. [Google Scholar] [CrossRef]
  50. Liu, J.; Yang, Q.; Liu, J.; Zhang, Y.; Jiang, X.; Yang, Y. Study on the Spatial Differentiation of the Populations on Both Sides of the “Qinling-Huaihe Line” in China. Sustainability 2020, 12, 4545. [Google Scholar] [CrossRef]
  51. Zhang, C.; Di, L.; Lin, L.; Li, H.; Guo, L.; Yang, Z.; Yu, E.G.; Di, Y.; Yang, A. Towards Automation of In-Season Crop Type Mapping Using Spatiotemporal Crop Information and Remote Sensing Data. Agric. Syst. 2022, 201, 103462. [Google Scholar] [CrossRef]
  52. Zhang, K.; Chen, Y.; Zhang, B.; Hu, J.; Wang, W. A Multitemporal Mountain Rice Identification and Extraction Method Based on the Optimal Feature Combination and Machine Learning. Remote Sens. 2022, 14, 5096. [Google Scholar] [CrossRef]
  53. Fatchurrachman; Rudiyanto; Soh, N.C.; Shah, R.M.; Giap, S.G.E.; Setiawan, B.I.; Minasny, B. High-Resolution Mapping of Paddy Rice Extent and Growth Stages across Peninsular Malaysia Using a Fusion of Sentinel-1 and 2 Time Series Data in Google Earth Engine. Remote Sens. 2022, 14, 1875. [Google Scholar] [CrossRef]
  54. Sentinel-2 Datasets in Earth Engine|Earth Engine Data Catalog|Google Developers. Available online: https://developers.google.com/earth-engine/datasets/catalog/sentinel-2 (accessed on 8 May 2023).
  55. NASA SRTM Digital Elevation 30m|Earth Engine Data Catalog. Available online: https://developers.google.com/earth-engine/datasets/catalog/USGS_SRTMGL1_003 (accessed on 8 May 2023).
  56. Qiu, B.; Hu, X.; Chen, C.; Tang, Z.; Yang, P.; Zhu, X.; Yan, C.; Jian, Z. Maps of Cropping Patterns in China during 2015–2021. Sci. Data 2022, 9, 479. [Google Scholar] [CrossRef]
  57. Xiao, X.; He, L.; Salas, W.; Li, C.; Moore, B.; Zhao, R.; Frolking, S.; Boles, S. Quantitative Relationships between Field-Measured Leaf Area Index and Vegetation Index Derived from VEGETATION Images for Paddy Rice Fields. Int. J. Remote Sens. 2002, 23, 3595–3604. [Google Scholar] [CrossRef]
  58. Wang, F.; Huang, J.; Tang, Y.; Wang, X. New Vegetation Index and Its Application in Estimating Leaf Area Index of Rice. Rice Sci. 2007, 14, 195–203. [Google Scholar] [CrossRef]
  59. Harrell, D.L.; Tubaña, B.S.; Walker, T.W.; Phillips, S.B. Estimating Rice Grain Yield Potential Using Normalized Difference Vegetation Index. Agron. J. 2011, 103, 1717–1723. [Google Scholar] [CrossRef]
  60. Zhang, G.; Xiao, X.; Dong, J.; Kou, W.; Jin, C.; Qin, Y.; Zhou, Y.; Wang, J.; Menarguez, M.A.; Biradar, C. Mapping Paddy Rice Planting Areas through Time Series Analysis of MODIS Land Surface Temperature and Vegetation Index Data. ISPRS J. Photogramm. Remote Sens. 2015, 106, 157–171. [Google Scholar] [CrossRef]
  61. Ni, R.; Tian, J.; Li, X.; Yin, D.; Li, J.; Gong, H.; Zhang, J.; Zhu, L.; Wu, D. An Enhanced Pixel-Based Phenological Feature for Accurate Paddy Rice Mapping with Sentinel-2 Imagery in Google Earth Engine. ISPRS J. Photogramm. Remote Sens. 2021, 178, 282–296. [Google Scholar] [CrossRef]
  62. Huete, A.; Didan, K.; Miura, T.; Rodriguez, E.P.; Gao, X.; Ferreira, L.G. Overview of the Radiometric and Biophysical Performance of the MODIS Vegetation Indices. Remote Sens. Environ. 2002, 83, 195–213. [Google Scholar] [CrossRef]
  63. Son, N.T.; Chen, C.F.; Chen, C.R.; Minh, V.Q.; Trung, N.H. A Comparative Analysis of Multitemporal MODIS EVI and NDVI Data for Large-Scale Rice Yield Estimation. Agric. For. Meteorol. 2014, 197, 52–64. [Google Scholar] [CrossRef]
  64. Guan, K.; Li, Z.; Rao, L.N.; Gao, F.; Xie, D.; Hien, N.T.; Zeng, Z. Mapping Paddy Rice Area and Yields Over Thai Binh Province in Viet Nam From MODIS, Landsat, and ALOS-2/PALSAR-2. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 2238–2252. [Google Scholar] [CrossRef]
  65. Feng, S.; Zhao, J.; Liu, T.; Zhang, H.; Zhang, Z.; Guo, X. Crop Type Identification and Mapping Using Machine Learning Algorithms and Sentinel-2 Time Series Data. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2019, 12, 3295–3306. [Google Scholar] [CrossRef]
  66. Yang, G.; Yu, W.; Yao, X.; Zheng, H.; Cao, Q.; Zhu, Y.; Cao, W.; Cheng, T. AGTOC: A Novel Approach to Winter Wheat Mapping by Automatic Generation of Training Samples and One-Class Classification on Google Earth Engine. Int. J. Appl. Earth Obs. Geoinf. 2021, 102, 102446. [Google Scholar] [CrossRef]
  67. Wang, Y.; Wong, J.; Miner, A.S. Anomaly Intrusion Detection Using One Class SVM. In Proceedings of the Fifth Annual IEEE SMC Information Assurance Workshop, West Point, NY, USA, 10–11 June 2004; pp. 358–364. [Google Scholar]
  68. Chen, Y.; Zhou, X.S.; Huang, T.S. One-Class SVM for Learning in Image Retrieval. In Proceedings of the 2001 International Conference on Image Processing (Cat. No.01CH37205), Thessaloniki, Greece, 7–10 October 2001; Volume 1, pp. 34–37. [Google Scholar]
  69. Bounsiar, A.; Madden, M.G. One-Class Support Vector Machines Revisited. In Proceedings of the 2014 International Conference on Information Science & Applications (ICISA), Seoul, Republic of Korea, 6–9 May 2014; pp. 1–4. [Google Scholar]
  70. Zhang, X.; Wu, B.; Ponce-Campos, G.; Zhang, M.; Chang, S.; Tian, F. Mapping Up-to-Date Paddy Rice Extent at 10 M Resolution in China through the Integration of Optical and Synthetic Aperture Radar Images. Remote Sens. 2018, 10, 1200. [Google Scholar] [CrossRef]
  71. Bazzi, H.; Baghdadi, N.; El Hajj, M.; Zribi, M.; Minh, D.H.T.; Ndikumana, E.; Courault, D.; Belhouchette, H. Mapping Paddy Rice Using Sentinel-1 SAR Time Series in Camargue, France. Remote Sens. 2019, 11, 887. [Google Scholar] [CrossRef]
  72. Nguyen, D.B.; Gruber, A.; Wagner, W. Mapping Rice Extent and Cropping Scheme in the Mekong Delta Using Sentinel-1A Data. Remote Sens. Lett. 2016, 7, 1209–1218. [Google Scholar] [CrossRef]
  73. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef]
  74. Wang, Y.; Zang, S.; Tian, Y. Mapping Paddy Rice with the Random Forest Algorithm Using MODIS and SMAP Time Series. Chaos Solitons Fractals 2020, 140, 110116. [Google Scholar] [CrossRef]
  75. Mansaray, L.R.; Wang, F.; Huang, J.; Yang, L.; Kanu, A.S. Accuracies of Support Vector Machine and Random Forest in Rice Mapping with Sentinel-1A, Landsat-8 and Sentinel-2A Datasets. Geocarto Int. 2020, 35, 1088–1108. [Google Scholar] [CrossRef]
  76. 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]
  77. Liu, Y.; Xiao, D.; Yang, W. An Algorithm for Early Rice Area Mapping from Satellite Remote Sensing Data in Southwestern Guangdong in China Based on Feature Optimization and Random Forest. Ecol. Inform. 2022, 72, 101853. [Google Scholar] [CrossRef]
  78. Majority Filter—Help|Documentation. Available online: https://desktop.arcgis.com/en/arcmap/10.7/tools/spatial-analyst-toolbox/majority-filter.htm (accessed on 8 May 2023).
  79. Zhang, C.; Zhang, H.; Du, J.; Zhang, L. Automated Paddy Rice Extent Extraction with Time Stacks of Sentinel Data: A Case Study in Jianghan Plain, Hubei, China. In Proceedings of the 2018 7th International Conference on Agro-geoinformatics (Agro-geoinformatics), Hangzhou, China, 6–9 August 2018; pp. 1–6. [Google Scholar]
  80. 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. [Google Scholar] [CrossRef]
  81. Ye, Q.; Yang, X.; Dai, S.; Chen, G.; Li, Y.; Zhang, C. Effects of Climate Change on Suitable Rice Cropping Areas, Cropping Systems and Crop Water Requirements in Southern China. Agric. Water Manag. 2015, 159, 35–44. [Google Scholar] [CrossRef]
  82. Zhao, H.; Fu, Y.H.; Wang, X.; Zhao, C.; Zeng, Z.; Piao, S. Timing of Rice Maturity in China Is Affected More by Transplanting Date than by Climate Change. Agric. For. Meteorol. 2016, 216, 215–220. [Google Scholar] [CrossRef]
  83. Ali, I.; Cao, S.; Naeimi, V.; Paulik, C.; Wagner, W. Methods to Remove the Border Noise From Sentinel-1 Synthetic Aperture Radar Data: Implications and Importance For Time-Series Analysis. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 777–786. [Google Scholar] [CrossRef]
  84. De Sá, N.C.; Baratchi, M.; Hauser, L.T.; Van Bodegom, P. Exploring the Impact of Noise on Hybrid Inversion of PROSAIL RTM on Sentinel-2 Data. Remote Sens. 2021, 13, 648. [Google Scholar] [CrossRef]
  85. Filipponi, F. Sentinel-1 GRD Preprocessing Workflow. Proceedings 2019, 18, 11. [Google Scholar] [CrossRef]
  86. Yan, S.; Xu, L.; Wu, R. Automatic Classification of Glaciers from Sentinel-2 Imagery Using A Novel Deep Learning Model. In Proceedings of the Proceedings of the 2019 3rd International Conference on Advances in Image Processing, Chengdu, China, 8 November 2019; pp. 155–159. [Google Scholar]
  87. Zhao, R.; Li, Y.; Chen, J.; Ma, M.; Fan, L.; Lu, W. Mapping a Paddy Rice Area in a Cloudy and Rainy Region Using Spatiotemporal Data Fusion and a Phenology-Based Algorithm. Remote Sens. 2021, 13, 4400. [Google Scholar] [CrossRef]
  88. Lei, T.C.; Wan, S.; Wu, Y.C.; Wang, H.-P.; Hsieh, C.-W. Multi-Temporal Data Fusion in MS and SAR Images Using the Dynamic Time Warping Method for Paddy Rice Classification. Agriculture 2022, 12, 77. [Google Scholar] [CrossRef]
  89. Pott, L.P.; Amado, T.J.C.; Schwalbert, R.A.; Corassa, G.M.; Ciampitti, I.A. Satellite-Based Data Fusion Crop Type Classification and Mapping in Rio Grande Do Sul, Brazil. ISPRS J. Photogramm. Remote Sens. 2021, 176, 196–210. [Google Scholar] [CrossRef]
  90. Stroppiana, D.; Villa, P.; Sona, G.; Ronchetti, G.; Candiani, G.; Pepe, M.; Busetto, L.; Migliazzi, M.; Boschetti, M. Early Season Weed Mapping in Rice Crops Using Multi-Spectral UAV Data. Int. J. Remote Sens. 2018, 39, 5432–5452. [Google Scholar] [CrossRef]
  91. Du, M.; Huang, J.; Wei, P.; Yang, L.; Chai, D.; Peng, D.; Sha, J.; Sun, W.; Huang, R. Dynamic Mapping of Paddy Rice Using Multi-Temporal Landsat Data Based on a Deep Semantic Segmentation Model. Agronomy 2022, 12, 1583. [Google Scholar] [CrossRef]
  92. Khan, H.R.; Gillani, Z.; Jamal, M.H.; Athar, A.; Chaudhry, M.T.; Chao, H.; He, Y.; Chen, M. Early Identification of Crop Type for Smallholder Farming Systems Using Deep Learning on Time-Series Sentinel-2 Imagery. Sensors 2023, 23, 1779. [Google Scholar] [CrossRef]
  93. Xu, J.; Yang, J.; Xiong, X.; Li, H.; Huang, J.; Ting, K.C.; Ying, Y.; Lin, T. Towards Interpreting Multi-Temporal Deep Learning Models in Crop Mapping. Remote Sens. Environ. 2021, 264, 112599. [Google Scholar] [CrossRef]
  94. Cámara, M.; López, F. Mathematical Morphology Applied to Raster Generalization of Urban City Block Maps. Cartographica 2000, 37, 33–48. [Google Scholar] [CrossRef]
  95. McDonald, T.; Chen, Y.R. Application of Morphological Image Processing in Agriculture. Trans. ASAE 1990, 33, 1346–1352. [Google Scholar] [CrossRef]
  96. You, N.; Dong, J.; Li, J.; Huang, J.; Jin, Z. Rapid Early-Season Maize Mapping without Crop Labels. Remote Sens. Environ. 2023, 290, 113496. [Google Scholar] [CrossRef]
  97. 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]
  98. Cui, J.; Zhang, X.; Wang, W.; Wang, L. Integration of Optical and SAR Remote Sensing Images for Crop-Type Mapping Based on a Novel Object-Oriented Feature Selection Method. Int. J. Agric. Biol. Eng. 2020, 13, 178–190. [Google Scholar] [CrossRef]
  99. Biswas, S. Importance of Chat GPT in Agriculture: According to Chat GPT. SSRN J. 2023. [Google Scholar] [CrossRef]
  100. Zhu, J.-J.; Jiang, J.; Yang, M.; Ren, Z.J. ChatGPT and Environmental Research. Environ. Sci. Technol. 2023. [Google Scholar] [CrossRef]
Figure 1. Study area (Monsoon region of the People’s Republic of China except for Hong Kong, Macao, and Taiwan), geographical boundaries, and zoning status. (Base map layer source: https://services.arcgisonline.com/arcgis/rest/services, accessed on 4 May 2023.)
Figure 1. Study area (Monsoon region of the People’s Republic of China except for Hong Kong, Macao, and Taiwan), geographical boundaries, and zoning status. (Base map layer source: https://services.arcgisonline.com/arcgis/rest/services, accessed on 4 May 2023.)
Remotesensing 15 04055 g001
Figure 2. Overview of the distribution of mid-season rice samples across China. (Base map layer source: https://services.arcgisonline.com/arcgis/rest/services, accessed on 4 May 2023.)
Figure 2. Overview of the distribution of mid-season rice samples across China. (Base map layer source: https://services.arcgisonline.com/arcgis/rest/services, accessed on 4 May 2023.)
Remotesensing 15 04055 g002
Figure 3. Mid-season rice sample collection based on visual interpretation of Google Earth images. (a) Partial sample points in Heilongjiang Province (128°24′E, 47°11′N). (b) Partial sample points in Liaoning Province (122°12′E, 40°55′N). (c) Partial sample points in Anhui Province (116°41′E, 33°20′N). (d) Partial sample points in Hunan Province (112°30′E, 28°50′N).
Figure 3. Mid-season rice sample collection based on visual interpretation of Google Earth images. (a) Partial sample points in Heilongjiang Province (128°24′E, 47°11′N). (b) Partial sample points in Liaoning Province (122°12′E, 40°55′N). (c) Partial sample points in Anhui Province (116°41′E, 33°20′N). (d) Partial sample points in Hunan Province (112°30′E, 28°50′N).
Remotesensing 15 04055 g003
Figure 4. Overall technology roadmap.
Figure 4. Overall technology roadmap.
Remotesensing 15 04055 g004
Figure 5. Four optical index profiles of a mid-season rice sample point (122°12′E, 40°55′N) in Liaoning Province in 2021. (The four curves in the figure are sixth-order polynomial fitted trend lines. Missing values mean that there is no valid pixel and the vegetation index is not available for calculation at that time.)
Figure 5. Four optical index profiles of a mid-season rice sample point (122°12′E, 40°55′N) in Liaoning Province in 2021. (The four curves in the figure are sixth-order polynomial fitted trend lines. Missing values mean that there is no valid pixel and the vegetation index is not available for calculation at that time.)
Remotesensing 15 04055 g005
Figure 6. Extraction of phenological characteristics of mid-season rice in southern provinces. (a) Time-series curves of VH polarization backscatter coefficients of various land classes of Sentinel I in a region of Anhui Province (116°55′E, 33°25′N) in 2021. (b) Growth patterns of rice at various phenological stages.
Figure 6. Extraction of phenological characteristics of mid-season rice in southern provinces. (a) Time-series curves of VH polarization backscatter coefficients of various land classes of Sentinel I in a region of Anhui Province (116°55′E, 33°25′N) in 2021. (b) Growth patterns of rice at various phenological stages.
Remotesensing 15 04055 g006
Figure 7. Schematic diagram of the process of majority filtering.
Figure 7. Schematic diagram of the process of majority filtering.
Remotesensing 15 04055 g007
Figure 8. Comparison of mid-season rice-mapping results in a certain district of Jiangxi (116°31′E, 29°8′N) before and after the process of majority filtering. (a) Before filtering. (b) After filtering.
Figure 8. Comparison of mid-season rice-mapping results in a certain district of Jiangxi (116°31′E, 29°8′N) before and after the process of majority filtering. (a) Before filtering. (b) After filtering.
Remotesensing 15 04055 g008
Figure 9. Overall distribution of mid-season rice cultivation areas of three typical provinces in northern China in 2021. (a) Heilongjiang Province. (b) Liaoning Province. (c) Jilin Province.
Figure 9. Overall distribution of mid-season rice cultivation areas of three typical provinces in northern China in 2021. (a) Heilongjiang Province. (b) Liaoning Province. (c) Jilin Province.
Remotesensing 15 04055 g009
Figure 10. Partial results of mapping the cultivation extent of mid-season rice in 2021 in three typical provinces in northern China against high-definition Google Earth satellite imagery. The blue blocks represent the extracted mid-season rice fields. (a) Heilongjiang Province. (b) Jilin Province. (c) Liaoning Province.
Figure 10. Partial results of mapping the cultivation extent of mid-season rice in 2021 in three typical provinces in northern China against high-definition Google Earth satellite imagery. The blue blocks represent the extracted mid-season rice fields. (a) Heilongjiang Province. (b) Jilin Province. (c) Liaoning Province.
Remotesensing 15 04055 g010
Figure 11. Overall distribution of mid-season rice cultivation areas of seven typical provinces in southern China in 2021. (a) Hunan Province. (b) Anhui Province. (c) Hubei Province. (d) Jiangsu Province. (e) Jiangxi Province. (f) Zhejiang Province. (g) Guangxi Zhuang Autonomous Region.
Figure 11. Overall distribution of mid-season rice cultivation areas of seven typical provinces in southern China in 2021. (a) Hunan Province. (b) Anhui Province. (c) Hubei Province. (d) Jiangsu Province. (e) Jiangxi Province. (f) Zhejiang Province. (g) Guangxi Zhuang Autonomous Region.
Remotesensing 15 04055 g011
Figure 12. Partial results of mapping the cultivation extent of mid-season rice in 2021 in seven typical provinces in southern China against high-definition Google Earth satellite imagery. The green blocks represent the extracted mid-season rice fields, and the blue blocks represent non-rice fields. (a) Hunan Province. (b) Hubei Province. (c) Jiangxi Province. (d) Jiangsu Province. (e) Anhui Province. (f) Zhejiang Province. (g) Guangxi Zhuang Autonomous Region.
Figure 12. Partial results of mapping the cultivation extent of mid-season rice in 2021 in seven typical provinces in southern China against high-definition Google Earth satellite imagery. The green blocks represent the extracted mid-season rice fields, and the blue blocks represent non-rice fields. (a) Hunan Province. (b) Hubei Province. (c) Jiangxi Province. (d) Jiangsu Province. (e) Anhui Province. (f) Zhejiang Province. (g) Guangxi Zhuang Autonomous Region.
Remotesensing 15 04055 g012
Table 1. Main differences between monsoon and non-monsoon areas in China.
Table 1. Main differences between monsoon and non-monsoon areas in China.
Area/CharacteristicMonsoon AreaNon-Monsoon Area
PrecipitationOver 400 mm per yearBelow 400 mm per year
Relationship between
Humidity and Temperature
SynchronousOut of sync
Terrain TypePlains, basins, and hillsPlateau, mountain, and basin
Vegetation TypeForests and grasslandsGrasslands and deserts
River TypeExoreic rivers
mainly recharged by rainwater
Internal streams
recharged by melting of snow and ice
Agricultural Production ModeFarming-dominatedLivestock-dominated
Table 2. Characteristics of the Sentinel series data obtained from the GEE.
Table 2. Characteristics of the Sentinel series data obtained from the GEE.
Satellite/PropertySentinel-1Sentinel-2
Product NameSentinel-1 SAR GRDHarmonized Sentinel-2 MSI, Level-1C
Time Coverage1 January–31 December 20211 January–31 December 2021
Data UnitdBNumerical value
Spatial Resolution10 m10–20 m
Temporal Resolution12 d5 d
Polarization/BandVH (IW Mode)B2 (Blue)
B4 (Red)
B8 (NIR)
B3 (Green)
B11, B12 (SWIR)
B5-B7, B8A (Red Edge)
Pre-processing StepsThermal noise removal,
Radiometric calibration,
Terrain correction
Preliminary cloud removal
Table 3. Optical indexes based on Sentinel-2 data for rice mapping in the northern region.
Table 3. Optical indexes based on Sentinel-2 data for rice mapping in the northern region.
Optical IndexCalculation Formula 1Phenological Stages
Bare Soil Index (BSI) B S I = S W I R + R e d N I R + B l u e S W I R + R e d + N I R + B l u e Bare Soil Stage
Land Surface Water Index (LSWI) L S W I = N I R S W I R N I R + S W I R Transplanting Stage
Normalized Difference Vegetation Index (NDVI) N D V I = N I R R e d N I R + R e d Growing Stage
Plant Senescence Reflectance Index (PSRI) P S R I = R e d B l u e R e d E d g e Mature Stage
Enhanced Vegetation Index (EVI) E V I = 2.5 × N I R R e d N I R + 6 × R e d 7.5 × B l u e + 1 Growing Stage
Green Chlorophyll Vegetation Index (GCVI) G C V I = N I R G r e e n 1 Transplanting Stage
1 Red, Blue, Green, NIR, SWIR, and Red Edge indicate the reflectance of red band, blue band, green band, near-infrared band, short-wavelength infrared band, and red edge band, respectively.
Table 4. Rice phenology in southern region of China in 2021 1.
Table 4. Rice phenology in southern region of China in 2021 1.
Region/StageJianghuai 2Jiangnan 3Huanan 4Xinan 5
Sowing Stage- 64.1–4.301.21–2.203.11–3.20
Seedling Stage4.11–5.315.1–5.312.21–3.203.21–3.31
Transplanting Stage6.1–6.106.1–6.104.1–4.204.1–5.20
Greening Stage6.11–6.206.11–6.20--
Tillering Stage6.21–7.316.21–7.314.21–5.315.21–7.10
Booting Stage8.1–8.10-6.1–6.107.11–7.31
Heading Stage8.11–8.208.1–8.106.11–6.208.1–8.10
Grouting Stage8.21–8.318.11–8.20-8.11–8.20
Milk-ripe Stage9.1–9.108.21–10.316.21–6.308.21–9.20
Mature Stage9.11–9.20-7.11–7.319.21–10.20
1 From the Crop Information Calendar on the website of the Ministry of Agriculture and Rural Affairs of the People’s Republic of China (http://zdscxx.moa.gov.cn:8080/nyb/pc/calendar.jsp, accessed on 4 May 2023). 2 Including Anhui, Hubei, and Jiangsu. 3 Including Zhejiang, Shanghai, Jiangxi, Hunan, Fujian, and Taiwan. 4 Including Guangdong, Guangxi, Hainan, Hong Kong, and Macau. 5 Including Chongqing, Sichuan, Guizhou, and Yunnan. 6 “-” indicates that this phenological period does not exist for rice in the region or is too short to be measured.
Table 5. Accuracy evaluation of rice-mapping results in three typical northern provinces.
Table 5. Accuracy evaluation of rice-mapping results in three typical northern provinces.
Province/ IndexHeilongjiangJilinLiaoning
Buffer Size N (m)101010
Confusion Matrix[62, 2],
[0, 43]
[55, 7],
[0, 38]
[43, 5],
[0, 29]
User’s Accuracy (UA)0.95560.84440.8529
Producer’s Accuracy (PA)0.96880.88710.8958
Overall Accuracy (OA)0.980.930.94
Kappa Coefficient (KC)0.960.870.87
Mapping Result Area (km2)37,345.738139.265009.72
Reference Area (km2)38,67083735206
Relative Error (%)−3.42−2.79−3.77
Table 6. Accuracy evaluation of rice-mapping results in seven typical southern provinces.
Table 6. Accuracy evaluation of rice-mapping results in seven typical southern provinces.
ProvinceHunanHubeiJiangxiJiangsuAnhuiZhejiangGuangxi
Buffer Size N (m)3545532
Confusion Matrix[68, 22],
[4, 67]
[94, 12],
[12, 72]
[62, 16],
[9, 51]
[93, 7],
[5, 77]
[100, 6],
[5, 38]
[80, 7],
[10, 44]
[104, 3],
[6, 22]
User’s Accuracy (UA)[0.94, 0.75][0.89, 0.86][0.87, 0.76][0.95, 0.92][0.95, 0.86][0.89, 0.86][0.95, 0.88]
Producer’s Accuracy (PA)[0.76, 0.94][0.89, 0.86][0.79, 0.85][0.93, 0.94][0.94, 0.88][0.89, 0.86][0.97, 0.79]
Overall Accuracy (OA)0.880.840.870.930.930.880.93
Kappa Coefficient (KC)0.780.740.830.870.820.740.79
F-score[0.88, 0.75][0.90, 0.90][0.90, 0.89][0.94, 0.93][0.95, 0.87][0.90, 0.84][0.96, 0.83]
Mapping Result Area (km2)13,298.6121,230.768963.9223,793.6026,755.345035.267285.36
Reference Area (km2)14761201929400221922512153148075
Relative Error (%)−9.915.14−4.647.226.51−5.25−9.78
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Huang, C.; You, S.; Liu, A.; Li, P.; Zhang, J.; Deng, J. High-Resolution National-Scale Mapping of Paddy Rice Based on Sentinel-1/2 Data. Remote Sens. 2023, 15, 4055. https://doi.org/10.3390/rs15164055

AMA Style

Huang C, You S, Liu A, Li P, Zhang J, Deng J. High-Resolution National-Scale Mapping of Paddy Rice Based on Sentinel-1/2 Data. Remote Sensing. 2023; 15(16):4055. https://doi.org/10.3390/rs15164055

Chicago/Turabian Style

Huang, Chenhao, Shucheng You, Aixia Liu, Penghan Li, Jianhua Zhang, and Jinsong Deng. 2023. "High-Resolution National-Scale Mapping of Paddy Rice Based on Sentinel-1/2 Data" Remote Sensing 15, no. 16: 4055. https://doi.org/10.3390/rs15164055

APA Style

Huang, C., You, S., Liu, A., Li, P., Zhang, J., & Deng, J. (2023). High-Resolution National-Scale Mapping of Paddy Rice Based on Sentinel-1/2 Data. Remote Sensing, 15(16), 4055. https://doi.org/10.3390/rs15164055

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