Next Article in Journal
Blind Remote Sensing Image Deblurring Using Local Binary Pattern Prior
Previous Article in Journal
Deep-Learning-Based Multispectral Image Reconstruction from Single Natural Color RGB Image—Enhancing UAV-Based Phenotyping
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of a Phenology-Based Method for Identifying Sugarcane Plantation Areas in China Using High-Resolution Satellite Datasets

1
School of Atmospheric Sciences, Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Sun Yat-sen University, Zhuhai 519082, China
2
Key Laboratory of Environment Change and Resources Use in Beibu Gulf, Ministry of Education, Nanning Normal University, Nanning 530001, China
3
College of Geomatics & Municipal Engineering, Zhejiang University of Water Resources and Electric Power, Hangzhou 310018, China
*
Author to whom correspondence should be addressed.
Submission received: 19 January 2022 / Revised: 25 February 2022 / Accepted: 1 March 2022 / Published: 5 March 2022

Abstract

:
Sugarcane is an important sugar and biofuel crop with high socio-economic importance, and its planted area has increased rapidly in recent years. China is the world’s third or fourth sugarcane producer. However, to our knowledge, no study has investigated the mapping of sugarcane cultivation areas across entire China. In this study, we developed a phenology-based method to identify sugarcane plantations in China at 30-m spatial resolution from 2016–2020 using the time-series of Landsat and Sentinel-1/2 images derived from Google Earth Engine (GEE) platform. The method worked by comparing the phenological similarity in normalized difference vegetation index (NDVI) series between unknown pixels and sugarcane samples. The phenological similarity was assessed using the time-weighted dynamic time warping method (TWDTW), which has less sensitivity to training samples than machine learning methods and therefore can be easily applied to large areas with limited samples. More importantly, our method introduced multiple and moving time standard phenological curves of sugarcane to the TWDTW by fully considering the variable crop life-cycle of sugarcane, particularly its long harvest season spanning from December to March of the following year. Validations showed the method performed well in 2019, with overall accuracies of 93.47% and 92.74% for surface reflectance (SR) and top of atmosphere reflectance (TOA) data, respectively. The sugarcane maps agreed well with the agricultural statistical areas from 2016–2020. The mapping accuracies using TOA data were comparable to SR data in 2019–2020, but outperformed SR data in 2016–2018 when SR data had lower availability on GEE. The sugarcane maps produced in this study can be used to monitor growing conditions and production of sugarcane and, therefore, can benefit sugarcane management, sustainable sugarcane production, and national food security.

1. Introduction

Sugarcane (Saccharum officinarum) is a semi-perennial crop with high socio-economic importance planted across tropical and subtropical regions of the world [1]. Sugarcane can be used to produce sugar and provides over 70% of the world’s sugar production [2,3]. Sugarcane residues are rich in cellulose and hemicellulose, which can be used to produce bioethanol and other liquid transportation biofuels [4]. Sugarcane biofuel is an alternative energy source that can contribute to reducing greenhouse gases (GHG) emissions associated with fossil fuel combustion [5,6,7]. For example, sugarcane bioethanol could reduce fossil fuel GHG emissions by 85% in Brazil [8]. Sugarcane expanded rapidly in the past decades by replacing other crops (e.g., soybean, rice, and maize), pasture, and forest [9,10,11]. On the one hand, the replacement of other food crops with sugarcane may directly affect agriculture and food security [7,12]. On the other hand, sugarcane expansion may influence the local climate by altering surface albedo and evapotranspiration [13]. Timely and accurate estimates of the planted areas of sugarcane can be used for sugarcane growth condition monitoring and yield estimation and, therefore, is important for sugarcane management and sustainable sugarcane production.
China is the world’s fourth largest sugarcane producing country, accounting for 5.3% of the world sugarcane plantation after Brazil (37.6%), India (18.9%), and Thailand (6.9%) in 2019 [14]. During the past decades, over 90% of the sugar production in China came from sugarcane [15]. However, there are still no sugarcane cultivation maps covering all of China. In China, sugarcane generally begins to grow in March to April, matures from November to December, and can be harvested from December to March of the following year [16,17]. Sugarcane has a long harvest season, spanning three to four months; therefore, sugarcane maps can help detect the optimal harvest time and coordinate harvest operations in terms of laborers, machines, transportation, and sugar factories. In addition, sugarcane growth can be impacted by several biotic and abiotic stresses, such as diseases, pests, weeds, and natural disasters [15]. Especially, in south China with most sugarcane plantations, typhoon is a frequently occurred natural disaster during the sugarcane grand growth period. Typhoon and associated storms may lead to sugarcane lodging, causing serious plant damage and reducing sugarcane yield [18]. Therefore, there is an urgent need to develop an efficient method for sugarcane identification that can be easily applied to different regions and years in China.
Recently, several studies have used satellite based optical data and synthetic aperture radar (SAR) data to identify sugarcane plantations [1,19]. For example, the optical data include Landsat [20,21,22], Sentinel-2 [23], the Moderate Resolution Imaging Spectroradiometer (MODIS) [24], SPOT [25,26], Chinese HJ-1 CCD [27], Indian Remote Sensing Satellite (IRS-P6) [28]. The SAR data include ENVISAT Advanced SAR (ASAR) [16], TerraSAR-X [29], PALSAR [30], and Sentinel-1 [31]. Machine learning, especially random forest (RF), is a most current and popular method used in image classification and identification [32,33,34]. For example, Zhou et al. [27] identified sugarcane plantations in Suixi County in China using a decision tree and object-oriented classification technique based on optical images. Using sugarcane maps produced by the Canasat Project as ground reference samples [35], a space-time generalized random forest classifier was calibrated at local sites to identify sugarcane areas at the regional scale in São Paulo State in different years [36,37]. However, machine learning methods strongly rely on large volumes of ground training data, which are difficult to obtain at large scales, such as at country to continental scales [38].
Phenology-based methods can offer a potential alternative solution to identify crop plantations at large scales with low volumes of training data [39,40,41]. Compared with machine learning, phenology-based methods are developed by analyzing the life cycle of crops to construct rules for classification, which requires fewer training samples and are easier to be applied to other regions and years [42]. For example, Li et al. [29] extracted the polarization features of sugarcane from TerraSAR-X data in different tillering stages from early April to late May on the Leizhou Peninsula in China, and found tillering to be a suitable growing phase for sugarcane mapping. Sugarcane in Guangxi Province, China, was identified by using the differences in phenological thresholds between sugarcane and other landcover types, such as evergreen vegetation, double season crops, and triple season crops [40]. These studies were developed based on several phenological thresholds (e.g., green-up date, senescence date, and length of growing season) of crops in different growing stages, which may vary by region and year and need to be calibrated when extended to other regions or years. Therefore, traditional phenology-based methods are still insufficient to perform ideal mapping and may be limited by multiple thresholds when applied to large scales, such as the country to continental scales.
The time-weighted dynamic time warping method (TWDTW) was shown to be an efficient phenology-based method able to classify vegetation types at regional scales with limited samples [42,43,44]. The TWDTW method works by comparing the similarities in seasonal variations between the unknown pixels and the target classes. Compared with the dynamic time warping method (DTW), the TWDTW has superiority in time series classification by adding a temporal weight accounting for phenological changes of vegetation types besides shape matching in DTW [45]. The TWDTW method was shown to perform well in identifying crops (e.g., winter wheat, sugarcane, maize, and paddy rice), pasture, and forest with limited training samples at local to regional scales [23,42,44,46,47]. Moreover, the TWDTW method is less sensitive to training samples than machine learning methods, therefore can be easily applied to large scales and different years with limited samples [23].
In this study, we developed a phenology-based method by incorporating the TWDTW method into the sugarcane cultivation area mapping of China using Landsat and Sentinel-1/2 images. Additionally, we examined the mapping performance with both surface reflectance (SR) and top of atmosphere reflectance (TOA) data. Finally, sugarcane plantation area maps for China were produced at 30 m for the period 2016–2020 with both SR and TOA data.

2. Materials and Methods

2.1. Study Area

Sugarcane in China is majorly planted in the southern portion of the country, from the north of the Qinling Mountains to the south of Hainan Province. As shown in Figure 1, we identified sugarcane over 95% of the cultivation areas in China, including Guangxi (64.01%), Yunnan (17.70%), Guangdong (12.20%), and Hainan provinces (1.32%) [48]. Specially, Guangxi is the dominant sugarcane producer in China, accounting for about 65% of the national sugarcane production. Of these provinces, the dominant production areas are in the central-south of Guangxi, southwest of Yunnan, west of Guangdong, and west of Hainan [15].

2.2. Datasets

2.2.1. Landsat Data and Sentinel-2 Data

To test the capability of the model in the identification of sugarcane by using atmospherically corrected optical data from the SR and TOA data, we used all the available Landsat and Sentinel-2 TOA data and SR data at the Google Earth Engine (GEE) platform to generate 16-day composite 30 m normalized difference vegetation index (NDVI) series from 2016 to 2020. We used the Landsat-7 Enhanced Thematic Mapper (ETM+) and Landsat-8 Operational Land Imager (OLI) Collection 1 SR and TOA data at 16-day intervals from 2016–2020. The quality bands pixel_qa and BQA were used to remove cloud and cloud shadow contaminations to obtain high-quality SR and TOA images, respectively. We also used the Sentinel-2A/B Multi-Spectral Instrument (MSI) Level-1C (TOA) and Level-2A (SR) data at 5-day intervals from 2016–2020. The quality band QA60 in Level-1C and Level-2A were used to remove cloud and cloud shadow contaminations. The Level-2A Sentinel-2 data on GEE were derived from the Sentinel Hub generated from Level-1C using the Sen2Cor processor. These data were produced by the European Space Agency (ESA) over Europe since March 2018 and extended to global since December 2018. The Level-2A data can also be derived from Level-1C using the Sentinel Toolbox provided by ESA. However, the toolbox is not directly available in GEE and the process needs huge computation and storage when applied to a large scales. Therefore, in China, the Level-2A (SR) data cannot be obtained directly from either GEE or ESA at large scales before December 2018.
The 16-day composite NDVI were produced by selecting the maximum value with cloud-free observation in each 16-day period. NDVI can be calculated from the reflectance of near-infrared band (ρNIR) and red band (ρRED):
NDVI = ρ NIR ρ RED ρ NIR + ρ RED
We further filled the cloud gaps in the 16-day NDVI using the linear interpolation method, and filtered the filled NDVI using the Whittaker smoother to eliminate noises [49,50]. The number of high-quality observations (cloud-free observations) for the 16-day composite SR and TOA data for each year from 2016 to 2020 are shown in Figure 2. For TOA data, the majority of the pixels (over 60%) had more than 17 high-quality observations for each year from 2018 to 2020, while the high-quality observations were fewer for each pixel in the years 2016 and 2017 for the absence of Sentinel-2B data (Sentinel-2B was launched on March 2018) (Figure 2a–e). For SR data, the number of high-quality observations was comparable to the TOA data from 2019 to 2020, while it greatly decreased in 2016–2018 with most of the pixels had fewer than 14 high-quality observations (Figure 2f–j). The main cause for this is that, at the time of this study, the Sentinel-2A/B SR data for China were not available on the GEE platform before December 2018.

2.2.2. Sentinel-1 Data

SAR data from Sentinel-1 were used to separate banana from sugarcane (see Section 2.3.3). The Sentinel-1 mission provides data from a dual-polarization C-band SAR instrument at 5.405 GHz at a spatial resolution of 10 m. The data were processed using the Sentinel-1 Toolbox to generate a calibrated, ortho-corrected product by ESA. The refined Lee filter method was used to remove the speckle noises on the SAR images [51]. Finally, we obtained the averaged value of the VH (vertical transmit/horizontal receive) polarization band during the growing season of sugarcane (from June to November) for each year in 2016–2020 from the GEE platform.

2.2.3. Field Data and Agricultural Statistical Data

Field data were used for training and validating the proposed sugarcane identification method. We collected the field data from two sources. First, field surveys were conducted in the four major sugarcane production provinces in the year 2019, and the vegetation types and their locations were marked using a hand-held GPS (G120, UniStrong, Beijing, China). Second, sample points were selected and labelled as sugarcane or non-sugarcane according to the color and textures on the high-resolution images from Google Earth. To further guarantee the accuracy of the labels, we cross-checked each sample using their NDVI series curve to verify the crop type (i.e., sugarcane or non-sugarcane). Finally, for the year 2019, we collected a total of 5764 samples, of which 1596 were for sugarcane and 4168 for non-sugarcane (Figure 1).
Agricultural statistical data for the planted areas of sugarcane was collected from the National Bureau of Statistics of China at the province, municipal, and county levels in 2016–2019 [48]. The statistical data for the year 2020 were not yet published at the time of this study and were replaced with the year 2019. The sugarcane planted area at the county level could not be found for Yunnan. For Hainan Province, the administrative regions of municipalities and counties are not separated, therefore we used all the statistical data at the municipalities and counties in the validation at both the municipal and county levels.

2.3. Methods

As a typical and commonly used vegetation index in landcover and land use classification, NDVI has been successfully applied to the TWDTW in crop identification in several studies [23,42]. In our study, we used NDVI time series to represent the phenological changes of vegetation in TWDTW. The workflow includes the following steps: (1) pre-process of the satellite data, including optical data (i.e., Landsat 7/8 and Sentinel-2) and SAR data (Sentinel-1) (Section 2.2.1; Section 2.2.2); (2) extract and generate the standard NDVI curves (i.e., the averaged NDVI temporal patterns) of randomly selected 100 sugarcane samples from the 16-day composite NDVI time series for SR and TOA data; (3) develop the mapping method to identify the sugarcane planted area in China based on the TWDTW; (4) validate the sugarcane maps at point and regional levels in different provinces and years.

2.3.1. Time-Weighted Dynamic Time Warping (TWDTW) Method

The DTW was originally proposed for speech recognition [52]. The DTW compares the similarity between two sequences, such as the unknown sequence (Y) and target sequence (X), through warping the unknown sequence (Y) to search for their optimal path and obtain their minimum distance, namely the similarity (dissimilarity). The TWDTW is an improved version of the DTW method by adding a temporal weight accounting for phenological changes [43]. The TWDTW method includes three main steps in land cover and land use classification: (1) generation of the standard phenological curves (temporal pattern) of the field samples of crops to be classified using NDVI time series; (2) calculation of TWDTW distances between unknown pixels and the standard NDVI curves; (3) identification or classification of the unknown pixels according to the TWDTW distances, in which pixels with lower distances indicate higher similarities and, therefore, higher probabilities of being the target class.

2.3.2. Incorporating the TWDTW Method into Sugarcane Mapping in China

In China, newly planted sugarcane and ratoon sugarcane have similar growing stages and NDVI variations. Most of the sugarcane in China begins to grow from March to April and can be harvested from December to March of the following year, with about a one year life cycle [16,17]. The sugarcane crop cycle in China generally can be divided into four phenological stages, including (1) seeding or germination phase (including ratoon cane), from March to early June with low NDVI; (2) tillering phase, from May to August, with NDVI increasing quickly; (3) grand growth phase, from July to the end of November, with NDVI reaching its peak; (4) ripening and harvesting phase, from November to March of the following year, with NDVI sharply decreasing. The temporal patterns of NDVI, namely standard NDVI curves, for sugarcane in China for both SR and TOA data are shown in Figure 3. The standard NDVI curves were produced as following. (1) Randomly selecting 100 sugarcane samples across China from the field samples. (2) Extracting their NDVI curves from the 16-day composited NDVI time series in 2019/2020 crop year for both SR and TOA data. (3) Adjusting the timeline of each NDVI curves manually by moving forward or back 16 days (or 32 days) to ensure all the 100 NDVI curves generally in the same growing period. Because the majority of these NDVI curves are consistent with each other, we only need to adjust a very small number of them manually to let their germination phase, harvesting phase, and peak of growing season are relatively in the same 16-day period with the majorities. (4) Calculating their averaged NDVI values to obtain the standard NDVI curves for both SR and TOA data (Figure 3).
The TWDTW method was applied to sugarcane planted area identification in China with the following processes for both the SR and TOA data. (1) Generation of multiple standard curves from the standard NDVI curve. Because the germination and harvest times may vary in different areas and fields, we tested the result in each province using multiple and moving time standard curves. Namely, the standard curves of sugarcane in Figure 3 were moved forward and back by 16 days and 32 days, and therefore five standard curves were derived from the original standard curves for both SR and TOA data (Figure 4). (2) Calculation of the minimum TWDTW distances. For each unknown pixel, all the five NDVI distance (NDVIdis) and their minimum distance were calculated using the TWDTW method based on the five standard NDVI curves in Figure 4 for each year in 2016–2020. (3) Calculation of the NDVI difference (NDVIdiff) between sugarcane grand growth period and harvest season. The NDVI of sugarcane is high in its grand growth period and decreases quickly during harvest season, which can be a good feature to help further separate sugarcane from other vegetation types that have no discernible NDVI decrease during the harvest season of sugarcane (e.g., evergreen forest). Therefore, for each pixel, we calculated the “difference” between the maximum NDVI (NDVImax) in sugarcane grand growth period and the minimum NDVI (NDVImin) in sugarcane harvest season. The NDVImax was calculated as the averaged value of the three maximum 16-day NDVI from July to December (namely from the 161-day to the 321-day, shadowed with orange in Figure 4); and the NDVImin was calculated as the two minimum 16-day NDVI form December to April in the next year (namely from the 337-day to the 113-day in the next year, shadowed with grey in Figure 4). From 100 randomly selected sugarcane samples, we found most of the sugarcane samples had higher NDVIdiff than 0.36. Therefore, pixels with NDVIdiff less than 0.36 were labelled as non-sugarcane, and pixels with NDVIdiff greater than 0.36 were selected for further identification. (4) Identification and classification to obtain sugarcane maps. Pixels with lower dissimilarities (NDVIdis) have a higher probability of being sugarcane. Agricultural statistical areas of sugarcane at the province level were used to determine the threshold of NDVIdis (T-NDVIdis). Pixels with NDVIdis less than T-NDVIdis were recognized as sugarcane, with the condition that the identified area of sugarcane should be equal to the statistical area of sugarcane in the studied province. For each province, except Guangxi, municipalities with small areas of sugarcane plantation (less than 1000 ha or 1% plantation areas of the studied province) were identified separately to obtain higher accuracies.

2.3.3. Separating Banana from Sugarcane

Banana is a multi-year growing fruit exhibiting NDVI curves similar to those of sugarcane. Both banana and sugarcane have low NDVI at the start of the growing season (from March to May) and the end of the growing season (from December to March of the following year) because of harvest (Figure 5a). Where both crops are grown, some banana pixels may be misclassified as sugarcane. While compared with sugarcane, banana has larger leaves and branches with a rougher canopy surface and therefore exhibits larger backscatter coefficients (e.g., in VH polarization) during the growing season. Therefore, we calculated the mean value of VH from Sentinel-1 images from June to November for field samples of sugarcane and banana. The histogram shows that most of the sugarcane has an averaged VH less than −13.3 and most of the banana has an averaged VH greater than −13.3 (Figure 5b). Therefore, we set the threshold of the averaged VH from June to November as −13.3 to separate banana from sugarcane. Pixels with averaged VH greater than −13.3 were identified as non-sugarcane.

2.3.4. Accuracy Assessment

The model accuracy for sugarcane mapping was validated at both pixel and regional levels. At the pixel level, we assessed the sugarcane map by comparing it with the field samples derived from field surveys and high-resolution images from Google Earth using the indices of producer’s accuracy (PA), user’s accuracy (UA), and overall accuracy (OA). At the regional level, the identified sugarcane areas were calculated and compared with the agricultural statistical area at both municipal and county levels. The coefficient of determination (R2), RMSE (root mean squared error), and RMAE (relative mean absolute error) between the estimated areas and agricultural statistical areas of sugarcane were adopted to evaluate the mapping accuracy. The RMSE can be calculated as:
RMSE = 1 n i = 1 n ( S i S ^ i ) 2
The RMAE is the ratio of MAE (mean absolute error) to the averaged value of the agricultural statistical areas for all the administrative regions:
RMAE = MAE i = 1 n S i / n
where S ^ i and S i are the estimated area and agricultural statistical area of sugarcane for the i th administrative region, respectively; n is the number of administrative regions included in the validation. The MAE can be calculated as:
MAE = 1 n i = 1 n | S i S ^ i |
The MAE represents the accumulated absolute errors, and the RMSE reflects the deviation. Compared with MAE, RMSE is more sensitive to outliers.

3. Results

Annual sugarcane maps in the four major production provinces (Guangxi, Yunnan, Guangdong, and Hainan) in China from 2016–2020 were produced using the TWDTW method. We show year 2019 with SR data as an example in Figure 6. Based on the validations at pixel level using field samples from field surveys and Google Earth in 2019, the producer’s, user’s, and overall accuracies were 86.76%, 88.30%, and 93.47% with SR data; and 85.63%, 86.73%, and 92.74% with TOA data for sugarcane in China, respectively. The performance in different provinces decreased with the area of sugarcane, that is province with larger sugarcane area or proportion generally showed higher performance (Table 1). For example, Guangxi, the province with the largest sugarcane plantation (accounting for about 64.01% of the sugarcane in China), showed high producer’s, user’s, and overall accuracies of 87.66%, 89.25%, and 94.46% with SR data; and 86.31%, 87.89%, and 93.81% with TOA data, respectively (Table 1). Yunnan (17.70%) and Guangdong (12.20%) showed moderate performance, with an overall accuracy of 89.32% derived from SR data and 88.64% derived from TOA data for Yunnan, and 90.43% derived from SR data and 88.78% derived from TOA data for Guangdong (Table 1). Hainan, accounting for only 1.32% of the sugarcane in China, exhibited low performance (with an overall accuracy of 87.33% derived from SR data and 86.33% derived from TOA data) compared with other provinces. All the investigated provinces exhibited overall accuracies (OA) over 86% with both SR and TOA data (Table 1).
In addition, our method accurately estimated the planted area of sugarcane compared with agricultural statistical data at the municipal and county levels. The R2 were 0.97 and 0.92 at the municipal and county levels, respectively, and the respective RMAEs were 23.7% and 38.2% with SR data (Figure 7a,c). The model performance by using TOA data was comparable to that with SR data, with R2 of 0.95 and 0.91 at the municipal and county levels, respectively, and the respective RMAEs were 26.6% and 40.0% (Figure 7b,d).
Correlations between the estimated and agricultural statistical planted areas of sugarcane exhibited differently in different provinces at the municipal and county levels. At the municipal level, the R2 ranged from 0.80 (Yunnan) to 0.99 (Guangdong) and the RMAE ranged from 19.4% (Guangxi) to 42.7% (Hainan) with SR data. And the R2 ranged from 0.73 (Yunnan) to 0.98 (Guangdong) and the RMAE ranged from 21.4% (Guangxi) to 43.2% (Hainan) with TOA data (Figure 8a). At the county level, the R2 ranged from 0.82 (Hainan) to 0.96 (Guangdong) and the RMAE ranged from 35.7% (Guangxi) to 51.1% (Guangdong) with SR data. And the R2 ranged from 0.83 (Hainan) to 0.94 (Guangdong) and the RMAE ranged from 37.1% (Guangxi) to 55.0% (Guangdong) with TOA data (Figure 8b).
The performance of the method also varied in different years (Figure 9 and Figure 10). For the SR data, the R2 generally increased with the year from 2016–2020, with lower R2 in 2016–2018 (the lowest R2 was 0.83 in 2016) and higher R2 in 2019 and 2020 (the highest R2 was 0.97 in 2019) for the entire China at the municipal level (Figure 9a). The RMAE showed the opposite trend in value to R2, with higher RMAE in 2016–2018 (the highest RMAE was 45.65% in 2016) and lower RMAE in 2019 and 2020 (the lowest RMAE was 23.66% in 2019) for the entire China at the municipal level (Figure 9e). Similar to SR data, for the TOA data, the R2 generally increased with the year from 2016–2020, with the lowest R2 (0.85) in 2016 and the highest R2 (0.95) in 2019, and the RMAE generally decreased with the year from 2016–2020, with the highest RMAE (42.87%) in 2016 and the lowest RMAE (26.57%) in 2019 (Figure 9b,f). Additionally, at the county level, the identification using TOA data generally showed good performance from 2016–2020. However, the accuracy using SR data were high in 2019 and 2020 but a bit lower from 2016–2018 (Figure 10).

4. Discussion

Sugarcane is distributed in tropical and subtropical areas of the world, with different phenological characteristics in different fields, regions, and countries. For example, sugarcane in Brazil has a long harvest season, spanning from April to December in the south-central Brazil and spanning from September to April of the following year in northeast Brazil [35]. In China, on the other hand, sugarcane enters the ripening stage in the November and can be harvested thereafter until March of the following year [31]. The long harvest season and variable phenological characteristics of sugarcane make it difficult to map sugarcane plantation areas at large scales. With the rapid expansion of sugarcane in recent years, studies appeared on sugarcane identification at local or regional scale. The first high resolution (30 m) sugarcane map in Brazil was generated for São Paulo for the 2003/04 crop year and updated annually using visual/manual image interpretation [35,53].
China is the world’s fourth largest sugarcane-producing country, accounting for about 5.64% of the world’s sugarcane production in 2019 [14]. In the recent ten years (2010–2019), the planted area of sugarcane in China decreased by about 14.35%, while its production increased by about 3.21% as a result of increasing yield (20.51%) [48]. In the recent four years, the planted area of sugarcane in China remained relatively stable (decreasing by about 0.78%), while its production and yield significantly increased by 5.98% and 6.18%, respectively [48]. Therefore, mapping sugarcane plantations in China is helpful to detect the mechanism of sugarcane plantation and yield variations in China. Some studies attempted to investigate sugarcane plantations in China at the local or regional scales. Tan et al. [17] extracted the sugarcane planted areas in Guangxi Province using maximum likelihood method and MODIS data. Jiang et al. [31] mapped the sugarcane in Zhanjiang City in Guangdong using RF and XGBoost classification methods based on Sentinel-1A and Sentinel-2 data. Wang et al. [40] identified the sugarcane in Guangxi Province using the thresholds of NDVI and phenology periods in different growing stages. However, to our knowledge, there are no moderate to high resolution sugarcane distribution maps all over the entire China. This is probably because machine learning methods strongly require extensive training samples, which are difficult to obtain [20,36,37]. Traditional phenology-based methods need fewer samples for training, but need local training samples to calibrate several thresholds when applied to other regions or years [40].
We proposed a phenology-based sugarcane identification method using the TWDTW method and our results showed a good performance in identifying planted areas of sugarcane. We also compared the performance of the method by using SR and TOA optical data from 2016–2020. Classifications based on high resolution optical data from the GEE platform often use TOA data instead of SR data mainly because the Sentinel-2A/B SR data before 2019 are not available on the GEE platform [40,54]. Therefore, it is necessary to examine whether it is suitable to use TOA data instead of SR data in sugarcane mapping on a large scale in our study. Study indicated that, in many applications involving classification and change detection, the atmospheric correction was unnecessary as long as the training data and the data to be classified were in the same relative scale [55]. Our study exhibited that the mapping accuracy of sugarcane was comparable by using SR and TOA data in 2019–2020, when there were sufficient data with high-quality observations (Figure 2 and Figure 11). However, the mapping accuracy of sugarcane using SR data was lower than TOA data in 2016–2018, with lower R2 and slope and higher RMSE and RMAE, when there were not enough SR images (Figure 2 and Figure 11). Our study demonstrated that, TOA data can be used to identify sugarcane plantation areas based on the proposed method when there are not enough available SR data.
Overall, the method developed in our study can identify the sugarcane planted areas automatically and effectively with limited training samples for China, but there remain potential uncertainties in our proposed method. For example, the sugarcane maps produced in our study contain some isolated pixels. As shown in Figure 12, sugarcane patches with only one pixel account for 5.52% (Guangxi), 8.85% (Hainan), 13.37% (Guangdong), 14.74% (Yunnan) of the sugarcane planted area with SR data, and account for 6.12% (Guangxi), 9.23% (Hainan), 11.41% (Guangdong), 11.68% (Yunnan) of the sugarcane planted area with TOA data in the year 2019. On the one hand, south China is located in the tropical and subtropical areas charactered with frequent rainy and cloudy weather, therefore speckle salt-and-pepper effects may exist in the sugarcane maps due to poor quality optical images [16]. In future study, applying an object-based image segmentation technique to our proposed method instead of the pixel-based method may reduce the salt-and-pepper effects in the sugarcane maps [23]. Furthermore, image fusion with high revisit MODIS data or cloud-free SAR data to reconstruct high quality NDVI time series may help to reduce the salt-and-pepper effects and improve mapping performance [56,57,58]. On the other hand, the landscape in south China is also highly heterogeneous and fragmented, with complex sugarcane plantation patterns. Some of the sugarcane is planted by small farm holders and is scattered in small patches. For example, the most sugarcane planted municipality in Guangdong is Zhanjiang, accounting for about 81.98% (139,100 ha) of the sugarcane plantations in Guangdong, while the remaining 18.02% (30,580 ha) of sugarcane plantations are dispersed in other municipalities [48]. The mapping accuracy at 30 m spatial resolution could be limited by the spectral mixing with other types of vegetation in the heterogeneous or scattered planted areas [27]. Additionally, directly comparing with other prevailing methods, such as machine learning, can help to deepen understanding of the method in our study and therefore reduce uncertainties. But the machine learning methods rely on plenty of samples which are difficult to apply to such a large region in current study.

5. Conclusions

Large area sugarcane identification is difficult because of the geographically variable phenology characteristics of sugarcane. In our study, we developed a phenology-based method to identify sugarcane plantation areas in China using Landsat-7/8, Sentinel-1, and Sentinel-2 images. The proposed method can identify sugarcane cultivation areas effectively and accurately with limited training sample data. Sugarcane maps were produced for China at 30 m spatial resolution from 2016 to 2020 with both SR and TOA data. Validations against field samples and agricultural statistical data showed the maps were in high accuracy. The sugarcane maps produced in our study are applicable to monitor the growing conditions and yields of sugarcane, and contribute to sustainable sugarcane production. Of course, there are uncertainties in our sugarcane maps. In the future, integrating other techniques to the proposed method, such as the object-based image segmentation technique or image fusion with high revisit MODIS data or cloud-free SAR data, may help to improve the model performance and obtain sugarcane maps with higher accuracy and longer time period.

Author Contributions

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

Funding

This research was funded by the Guangdong Major Project of Basic and Applied Basic Research, grant number 2020B0301030004; and the National Natural Science Foundation of China, grant number 31860144.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The sugarcane maps presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to thank the editors and the reviewers for their constructive comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Abdel-Rahman, E.M.; Ahmed, F.B. The application of remote sensing techniques to sugarcane (Saccharum spp. Hybrid) production: A review of the literature. Int. J. Remote Sens. 2008, 29, 3753–3767. [Google Scholar] [CrossRef]
  2. Lakshmanan, P.; Geijskes, R.J.; Aitken, K.S.; Grof, C.L.P.; Bonnett, G.D.; Smith, G.R. Sugarcane biotechnology: The challenges and opportunities. Vitr. Cell. Dev. Biol.-Plant 2005, 41, 345–363. [Google Scholar] [CrossRef]
  3. Brar, N.S.; Dhillon, B.S.; Saini, K.; Sharma, P. Agronomy of sugarbeet cultivation—A review. Agric. Rev. 2015, 36, 184–197. [Google Scholar] [CrossRef]
  4. Sindhu, R.; Gnansounou, E.; Binod, P.; Pandey, A. Bioconversion of sugarcane crop residue for value added products—An overview. Renew. Energy 2016, 98, 203–215. [Google Scholar] [CrossRef]
  5. Bordonal, R.d.O.; de Figueiredo, E.B.; Aguiar, D.A.; Adami, M.; Theodor Rudorff, B.F.; La Scala, N. Greenhouse gas mitigation potential from green harvested sugarcane scenarios in sao paulo state, Brazil. Biomass Bioenergy 2013, 59, 195–207. [Google Scholar] [CrossRef]
  6. Bordonal, R.d.O.; Lal, R.; Aguiar, D.A.; de Figueiredo, E.B.; Perillo, L.I.; Adami, M.; Theodor Rudorff, B.F.; La Scala, N. Greenhouse gas balance from cultivation and direct land use change of recently established sugarcane (Saccharum officinarum) plantation in South-Central Brazil. Renew. Sustain. Energy Rev. 2015, 52, 547–556. [Google Scholar] [CrossRef]
  7. Jaiswal, D.; De Souza, A.P.; Larsen, S.; LeBauer, D.S.; Miguez, F.E.; Sparovek, G.; Bollero, G.; Buckeridge, M.S.; Long, S.P. Brazilian sugarcane ethanol as an expandable green alternative to crude oil use. Nat. Clim. Change 2017, 7, 788–792. [Google Scholar] [CrossRef]
  8. Borjesson, P. Good or bad bioethanol from a greenhouse gas perspective—What determines this? Appl. Energy 2009, 86, 589–594. [Google Scholar] [CrossRef]
  9. Adami, M.; Theodor Rudorff, B.F.; Freitas, R.M.; Aguiar, D.A.; Sugawara, L.M.; Mello, M.P. Remote sensing time series to evaluate direct land use change of recent expanded sugarcane crop in brazil. Sustainability 2012, 4, 574–585. [Google Scholar] [CrossRef] [Green Version]
  10. Defante, L.R.; Vilpoux, O.F.; Sauer, L. Rapid expansion of sugarcane crop for biofuels and influence on food production in the first producing region of brazil. Food Policy 2018, 79, 121–131. [Google Scholar] [CrossRef]
  11. Ferreira, M.P.; Alves, D.S.; Shimabukuro, Y.E. Forest dynamics and land-use transitions in the brazilian atlantic forest: The case of sugarcane expansion. Reg. Environ. Chang. 2015, 15, 365–377. [Google Scholar] [CrossRef]
  12. Mello, F.F.C.; Cerri, C.E.P.; Davies, C.A.; Holbrook, N.M.; Paustian, K.; Maia, S.M.F.; Galdos, M.V.; Bernoux, M.; Cerri, C.C. Payback time for soil carbon and sugar-cane ethanol. Nat. Clim. Chang. 2014, 4, 605–609. [Google Scholar] [CrossRef]
  13. Loarie, S.R.; Lobell, D.B.; Asner, G.P.; Mu, Q.; Field, C.B. Direct impacts on local climate of sugar-cane expansion in brazil. Nat. Clim. Change 2011, 1, 105–109. [Google Scholar] [CrossRef]
  14. FAOSTAT, Food and Agriculture Organization of the United Nations (FAO). FAO Statistical Databases; FAO: Rome, Italy, 2020. [Google Scholar]
  15. Zhang, M.; Govindaraju, M. Sugarcane Production in China. In Sugarcane: Technology and Research; de Oliveira, A., Ed.; IntechOpen: London, UK, 2018; pp. 49–66. [Google Scholar]
  16. Lin, H.; Chen, J.; Pei, Z.; Zhang, S.; Hu, X. Monitoring sugarcane growth using envisat asar data. IEEE Trans. Geosci. Remote Sens. 2009, 47, 2572–2580. [Google Scholar] [CrossRef]
  17. Tan, Z.; Wu, L.; Ding, M.; Yang, X.; Ou, Z.; He, Y.; Kuang, Z. Study on the extraction of sugarcane planting areas from eos/modis data. Meteorol. Mon. 2007, 33, 76–81. [Google Scholar]
  18. Li, H.; Han, Y.; Chen, J. Capability of multidate radarsat-2 data to identify sugarcane lodging. J. Appl. Remote Sens. 2019, 13, 044514. [Google Scholar] [CrossRef]
  19. Som-ard, J.; Atzberger, C.; Izquierdo-Verdiguier, E.; Vuolo, F.; Immitzer, M. Remote sensing applications in sugarcane cultivation: A review. Remote Sens. 2021, 13, 4040. [Google Scholar] [CrossRef]
  20. Cechim Junior, C.; Johann, J.A.; Antunes, J.F.G. Mapping of sugarcane crop area in the parana state using landsat/tm/oli and irs/liss-3 images. Rev. Bras. Eng. Agric. E Ambient. 2017, 21, 427–432. [Google Scholar] [CrossRef] [Green Version]
  21. Johnson, B.A.; Scheyvens, H.; Shivakoti, B.R. An ensemble pansharpening approach for finer-scale mapping of sugarcane with landsat 8 imagery. Int. J. Appl. Earth Obs. Geoinf. 2014, 33, 218–225. [Google Scholar] [CrossRef]
  22. Mulyono, S.; Nadirah. Identifying sugarcane plantation using landsat-8 images with support vector machines. In Proceedings of the 2nd International Conference of Indonesian Society for Remote Sensing, Yogyakarta, Indonesia, 17–19 October 2016; Volume 47. [Google Scholar]
  23. Belgiu, M.; Csillik, O. Sentinel-2 cropland mapping using pixel-based and object-based time-weighted dynamic time warping analysis. Remote Sens. Environ. 2018, 204, 509–523. [Google Scholar] [CrossRef]
  24. Xavier, A.C.; Rudorff, B.F.T.; Berka, L.M.S.; Moreira, M.A. Multi-temporal analysis of modis data to classify sugarcane crop. Int. J. Remote Sens. 2006, 27, 755–768. [Google Scholar] [CrossRef]
  25. El Hajj, M.; Begue, A.; Guillaume, S.; Martine, J.-F. Integrating spot-5 time series, crop growth modeling and expert knowledge for monitoring agricultural practices—The case of sugarcane harvest on reunion island. Remote Sens. Environ. 2009, 113, 2052–2061. [Google Scholar] [CrossRef]
  26. Yang, C.; Everitt, J.H.; Murden, D. Evaluating high resolution spot 5 satellite imagery for crop identification. Comput. Electron. Agric. 2011, 75, 347–354. [Google Scholar] [CrossRef]
  27. Zhou, Z.; Huang, J.; Wang, J.; Zhang, K.; Kuang, Z.; Zhong, S.; Song, X. Object-oriented classification of sugarcane using time-series middle-resolution remote sensing data based on adaboost. PLoS ONE 2015, 10, e0142069. [Google Scholar] [CrossRef] [Green Version]
  28. Verma, A.K.; Garg, P.K.; Prasad, K.S.H. Sugarcane crop identification from liss iv data using isodata, mlc, and indices based decision tree approach. Arab. J. Geosci. 2017, 10, 16. [Google Scholar] [CrossRef]
  29. Li, H.; Chen, J.; Liang, S.; Li, Q. Sugarcane mapping in tillering period by quad-polarization terrasar-x data. IEEE Geosci. Remote Sens. Lett. 2015, 12, 993–997. [Google Scholar]
  30. Baghdadi, N.; Boyer, N.; Todoroff, P.; El Hajj, M.; Begue, A. Potential of sar sensors terrasar-x, asar/envisat and palsar/alos for monitoring sugarcane crops on reunion island. Remote Sens. Environ. 2009, 113, 1724–1738. [Google Scholar] [CrossRef]
  31. Jiang, H.; Li, D.; Jing, W.; Xu, J.; Huang, J.; Yang, J.; Chen, S. Early season mapping of sugarcane by applying machine learning algorithms to sentinel-1a/2 time series data: A case study in zhanjiang city, China. Remote Sens. 2019, 11, 861. [Google Scholar] [CrossRef] [Green Version]
  32. Virnodkar, S.S.; Pachghare, V.K.; Patil, V.C.; Jha, S.K. Application of Machine Learning on Remote Sensing Data for Sugarcane Crop Classification: A Review. In ICT Analysis and Applications; Springer: Berlin/Heidelberg, Germany, 2020; pp. 539–555. [Google Scholar]
  33. Boryan, C.; Yang, Z.; Mueller, R.; Craig, M. Monitoring us agriculture: The us department of agriculture, national agricultural statistics service, cropland data layer program. Geocarto Int. 2011, 26, 341–358. [Google Scholar] [CrossRef]
  34. 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]
  35. Theodor Rudorff, B.F.; de Aguiar, D.A.; da Silva, W.F.; Sugawara, L.M.; Adami, M.; Moreira, M.A. Studies on the rapid expansion of sugarcane for ethanol production in sao paulo state (Brazil) using landsat data. Remote Sens. 2010, 2, 1057–1076. [Google Scholar] [CrossRef] [Green Version]
  36. Dos Santos Luciano, A.C.; Araujo Picoli, M.C.; Rocha, J.V.; Junqueira Franco, H.C.; Sanches, G.M.; Lima Verde Leal, M.R.; le Maire, G. Generalized space-time classifiers for monitoring sugarcane areas in brazil. Remote Sens. Environ. 2018, 215, 438–451. [Google Scholar] [CrossRef]
  37. Dos Santos Luciano, A.C.; Araujo Picoli, M.C.; Rocha, J.V.; Duft, D.G.; Camargo Lamparelli, R.A.; Lima Verde Leal, M.R.; Le Maire, G. A generalized space-time obia classification scheme to map sugarcane areas at regional scale, using landsat images time-series and the random forest algorithm. Int. J. Appl. Earth Obs. Geoinf. 2019, 80, 127–136. [Google Scholar] [CrossRef]
  38. Dong, J.; Liu, W.; Han, W.; Xiang, K.; Lei, T.; Yuan, W. A phenology-based method for identifying the planting fraction of winter wheat using moderate-resolution satellite data. Int. J. Remote Sens. 2020, 41, 6892–6913. [Google Scholar] [CrossRef]
  39. Massey, R.; Sankey, T.T.; Congalton, R.G.; Yadav, K.; Thenkabail, P.S.; Ozdogan, M.; Meador, A.J.S. Modis phenology-derived, multi-year distribution of conterminous us crop types. Remote Sens. Environ. 2017, 198, 490–503. [Google Scholar] [CrossRef]
  40. Wang, J.; Xiao, X.; Liu, L.; Wu, X.; Qin, Y.; Steiner, J.L.; Dong, J. Mapping sugarcane plantation dynamics in guangxi, china, by time series sentinel-1, sentinel-2 and landsat images. Remote Sens. Environ. 2020, 247, 111951. [Google Scholar] [CrossRef]
  41. Wardlow, B.D.; Egbert, S.L.; Kastens, J.H. Analysis of time-series modis 250 m vegetation index data for crop classification in the US central great plains. Remote Sens. Environ. 2007, 108, 290–310. [Google Scholar] [CrossRef] [Green Version]
  42. Dong, J.; Fu, Y.; Wang, J.; Tian, H.; Fu, S.; Niu, Z.; Han, W.; Zheng, Y.; Huang, J.; Yuan, W. Early-season mapping of winter wheat in china based on landsat and sentinel images. Earth Syst. Sci. Data 2020, 12, 3081–3095. [Google Scholar] [CrossRef]
  43. Maus, V.; Camara, G.; Cartaxo, R.; Sanchez, A.; Ramos, F.M.; de Queiroz, G.R. A time-weighted dynamic time warping method for land-use and land-cover mapping. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 3729–3739. [Google Scholar] [CrossRef]
  44. Pan, B.; Zheng, Y.; Shen, R.; Ye, T.; Zhao, W.; Dong, J.; Ma, H.; Yuan, W. High resolution distribution dataset of double-season paddy rice in china. Remote Sens. 2021, 13, 4609. [Google Scholar] [CrossRef]
  45. Petitjean, F.; Inglada, J.; Gancarski, P. Satellite image time series analysis under time warping. IEEE Trans. Geosci. Remote Sens. 2012, 50, 3081–3095. [Google Scholar] [CrossRef]
  46. Guan, X.; Huang, C.; Liu, G.; Meng, X.; Liu, Q. Mapping rice cropping systems in vietnam using an ndvi-based time-series similarity measurement based on dtw distance. Remote Sens. 2016, 8, 19. [Google Scholar] [CrossRef] [Green Version]
  47. Li, M.; Bijker, W. Vegetable classification in indonesia using dynamic time warping of sentinel-1a dual polarization sar time series. Int. J. Appl. Earth Obs. Geoinf. 2019, 78, 268–280. [Google Scholar] [CrossRef]
  48. National Bureau of Statistics of China. National Statistical Yearbook 2017, 2018, 2019, 2020; China Statistics Press: Beijing, China, 2017–2020.
  49. Eilers, P.H.C. A perfect smoother. Anal. Chem. 2003, 75, 3631–3636. [Google Scholar] [CrossRef]
  50. Kong, D.; Zhang, Y.; Gu, X.; Wang, D. A robust method for reconstructing global modis evi time series on the google earth engine. ISPRS J. Photogramm. Remote Sens. 2019, 155, 13–24. [Google Scholar] [CrossRef]
  51. Abramov, S.; Rubel, O.; Lukin, V.; Kozhemiakin, R.; Kussul, N.; Shelestov, A.; Lavreniuk, M. Speckle reducing for sentinel-1 sar data. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Fort Worth, TX, USA, 23–28 July 2017; pp. 2353–2356. [Google Scholar]
  52. Sakoe, H.; Chiba, S. Dynamic programming algorithm optimization for spoken word recognition. IEEE Trans. Acoust. Speech Signal Process. 1978, 26, 43–49. [Google Scholar] [CrossRef] [Green Version]
  53. Theodor Rudorff, B.F.; Berka, L.M.S.; Moreira, M.A.; Duarte, V.; Xavier, A.C.; Rosa, V.G.C.; Shimabukuro, Y.E. Imagens de satélite no mapeamento e estimativa de área de cana-de-açúcar em São Paulo: Ano-safra 2003/04. Agric. São Paulo 2005, 52, 21–39. [Google Scholar]
  54. Jin, Z.; Azzari, G.; You, C.; Di Tommaso, S.; Aston, S.; Burke, M.; Lobell, D.B. Smallholder maize area and yield mapping at national scales with google earth engine. Remote Sens. Environ. 2019, 228, 115–128. [Google Scholar] [CrossRef]
  55. Song, C.; Woodcock, C.E.; Seto, K.C.; Lenney, M.P.; Macomber, S.A. Classification and change detection using landsat tm data: When and how to correct atmospheric effects? Remote Sens. Environ. 2001, 75, 230–244. [Google Scholar] [CrossRef]
  56. Gao, F.; Anderson, M.C.; Zhang, X.; Yang, Z.; Alfieri, J.G.; Kustas, W.P.; Mueller, R.; Johnson, D.M.; Prueger, J.H. Toward mapping crop progress at field scales through fusion of landsat and modis imagery. Remote Sens. Environ. 2017, 188, 9–25. [Google Scholar] [CrossRef] [Green Version]
  57. Hazaymeh, K.; Hassan, Q.K. Spatiotemporal image-fusion model for enhancing the temporal resolution of landsat-8 surface reflectance images using modis images. J. Appl. Remote Sens. 2015, 9, 96095. [Google Scholar] [CrossRef] [Green Version]
  58. Roy, D.P.; Ju, J.; Lewis, P.; Schaaf, C.; Gao, F.; Hansen, M.; Lindquist, E. Multi-temporal modis-landsat data fusion for relative radiometric normalization, gap filling, and prediction of landsat data. Remote Sens. Environ. 2008, 112, 3112–3130. [Google Scholar] [CrossRef]
Figure 1. The study area included four provinces in the south of China. The solid dots represent the sample data used for calibration and validation derived from field surveys and Google Earth images for the year 2019.
Figure 1. The study area included four provinces in the south of China. The solid dots represent the sample data used for calibration and validation derived from field surveys and Google Earth images for the year 2019.
Remotesensing 14 01274 g001
Figure 2. Number of high-quality observations for the 16-day composite top of atmosphere reflectance (TOA) and surface reflectance (SR) data from 2016 to 2020.
Figure 2. Number of high-quality observations for the 16-day composite top of atmosphere reflectance (TOA) and surface reflectance (SR) data from 2016 to 2020.
Remotesensing 14 01274 g002
Figure 3. Temporal patterns of NDVI (standard NDVI curves) and growing stages for sugarcane in China for both SR and TOA data.
Figure 3. Temporal patterns of NDVI (standard NDVI curves) and growing stages for sugarcane in China for both SR and TOA data.
Remotesensing 14 01274 g003
Figure 4. Multiple standard NDVI curves for sugarcane in China for both SR and TOA data. The orange areas are the periods used to calculate maximum NDVI (NDVImax) in sugarcane grand growth period and the grey areas are the periods used to calculate the minimum NDVI (NDVImin) in sugarcane harvest season.
Figure 4. Multiple standard NDVI curves for sugarcane in China for both SR and TOA data. The orange areas are the periods used to calculate maximum NDVI (NDVImax) in sugarcane grand growth period and the grey areas are the periods used to calculate the minimum NDVI (NDVImin) in sugarcane harvest season.
Remotesensing 14 01274 g004
Figure 5. (a) Temporal patterns of the 16-day composite NDVI for sugarcane and banana (taking SR data as an example); (b) and distributions of averaged VH (vertical transmit/horizontal receive) polarization backscatter coefficient from June to November for sugarcane and banana.
Figure 5. (a) Temporal patterns of the 16-day composite NDVI for sugarcane and banana (taking SR data as an example); (b) and distributions of averaged VH (vertical transmit/horizontal receive) polarization backscatter coefficient from June to November for sugarcane and banana.
Remotesensing 14 01274 g005
Figure 6. Sugarcane identification map derived from SR data in China in 2019. (a) Sugarcane map in China; (bg) zoomed-in areas of the typical areas in subfigure (a).
Figure 6. Sugarcane identification map derived from SR data in China in 2019. (a) Sugarcane map in China; (bg) zoomed-in areas of the typical areas in subfigure (a).
Remotesensing 14 01274 g006
Figure 7. Correlations between the estimated and agricultural statistical planted areas of sugarcane at the (a,b) municipal and (c,d) county level in China in 2019.
Figure 7. Correlations between the estimated and agricultural statistical planted areas of sugarcane at the (a,b) municipal and (c,d) county level in China in 2019.
Remotesensing 14 01274 g007
Figure 8. Correlations between the estimated and agricultural statistical planted areas of sugarcane at the (a) municipal and (b) county level in China for each province in 2019.
Figure 8. Correlations between the estimated and agricultural statistical planted areas of sugarcane at the (a) municipal and (b) county level in China for each province in 2019.
Remotesensing 14 01274 g008
Figure 9. Correlations between the estimated and agricultural statistical planted areas of sugarcane (R2, slope, and RMAE) at the municipal level in China for each year from 2016 to 2020.
Figure 9. Correlations between the estimated and agricultural statistical planted areas of sugarcane (R2, slope, and RMAE) at the municipal level in China for each year from 2016 to 2020.
Remotesensing 14 01274 g009
Figure 10. Correlations between the estimated and agricultural statistical planted areas of sugarcane (R2, slope, and RMAE) at the county level in China for each year from 2016 to 2020.
Figure 10. Correlations between the estimated and agricultural statistical planted areas of sugarcane (R2, slope, and RMAE) at the county level in China for each year from 2016 to 2020.
Remotesensing 14 01274 g010
Figure 11. Comparisons of the regional performance (R2, slope, RMSE, and RMAE) of the method in identification sugarcane planted areas with SR and TOA data at the municipal level for China from 2016 to 2020.
Figure 11. Comparisons of the regional performance (R2, slope, RMSE, and RMAE) of the method in identification sugarcane planted areas with SR and TOA data at the municipal level for China from 2016 to 2020.
Remotesensing 14 01274 g011
Figure 12. Statistics of patches with different pixel numbers (from 1 to 100) on the sugarcane maps for both SR and TOA data in 2019.
Figure 12. Statistics of patches with different pixel numbers (from 1 to 100) on the sugarcane maps for both SR and TOA data in 2019.
Remotesensing 14 01274 g012
Table 1. Confusion matrix of the sugarcane identification maps in China in 2019 derived from SR and TOA data.
Table 1. Confusion matrix of the sugarcane identification maps in China in 2019 derived from SR and TOA data.
ProvinceReferenceMapProducer’s AccuracyUser’s AccuracyOverall Accuracy
SugarcaneNon-Sugarcane
SRGuangXi (GX)Sugarcane98013887.66%89.25%94.46%
Non-sugarcane118338596.63%96.08%
Yunnan (YN)Sugarcane1582685.87%88.27%89.32%
Non-sugarcane2123591.80%90.04%
Guangdong (GD)Sugarcane871585.29%86.14%90.43%
Non-sugarcane1418793.03%92.57%
Hainan (HN)Sugarcane731979.35%79.35%87.33%
Non-sugarcane1918990.87%90.87%
TOAGuangXi (GX)Sugarcane96515386.31%87.89%93.81%
Non-sugarcane133337096.20%95.66%
Yunnan (YN)Sugarcane1562884.78%87.64%88.64%
Non-sugarcane2223491.41%89.31%
Guangdong (GD)Sugarcane851783.33%83.33%88.78%
Non-sugarcane1718491.54%91.54%
Hainan (HN)Sugarcane751781.52%75.76%86.33%
Non-sugarcane2418488.46%91.54%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zheng, Y.; Li, Z.; Pan, B.; Lin, S.; Dong, J.; Li, X.; Yuan, W. Development of a Phenology-Based Method for Identifying Sugarcane Plantation Areas in China Using High-Resolution Satellite Datasets. Remote Sens. 2022, 14, 1274. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14051274

AMA Style

Zheng Y, Li Z, Pan B, Lin S, Dong J, Li X, Yuan W. Development of a Phenology-Based Method for Identifying Sugarcane Plantation Areas in China Using High-Resolution Satellite Datasets. Remote Sensing. 2022; 14(5):1274. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14051274

Chicago/Turabian Style

Zheng, Yi, Zhuoting Li, Baihong Pan, Shangrong Lin, Jie Dong, Xiangqian Li, and Wenping Yuan. 2022. "Development of a Phenology-Based Method for Identifying Sugarcane Plantation Areas in China Using High-Resolution Satellite Datasets" Remote Sensing 14, no. 5: 1274. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14051274

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