Rice Mapping Research Articles

Development of a new index for automated mapping of ratoon rice areas using time-series normalized difference vegetation index imagery

Over recent decades, the global demand for food has continued to grow, owing to population growth and the loss of arable land. Rice ratooning offers new opportunities for increasing rice production and has received renewed interest because of the minimal additional labor input required for its adoption. Regular, regional-scale monitoring of the spatial patterns of both traditional and ratoon rice cropping systems provides essential information for agricultural resource management and food security studies. However, the similar phenological characteristics of traditional double rice and ratoon rice cropping systems make it challenging to accurately classify these cropping practices based on satellite observations alone. In this study, we first proposed an improved phenology-based rice cropping area detection algorithm using moderate resolution imaging spectroradiometer (MODIS) normalized difference vegetation index (NDVI) imagery. A new index, ratoon rice index, was then developed to automatically delineate ratoon rice cropping areas with the aid of a base map of rice in Hubei Province, China. The accuracy assessment using ground truth data showed that our approach could map both traditional and ratoon rice cropping areas with high user accuracy (91.25% and 91.43%, respectively). The MODIS-retrieved rice cropping areas were validated using annual agricultural census data, and coefficient of determination ( R 2 ) values of 0.60 and 0.41 were recorded for traditional and ratoon rice cropping systems, respectively. The total area of ratoon rice was estimated to be 1 283.6 km 2 , 5.0% of the total rice cropping area, in Hubei Province in 2016. These demonstrated the feasibility of extracting the spatial patterns of both traditional and ratoon rice cropping systems solely from time-series NDVI and field survey data and strides made in facilitating the timely and routine monitoring of traditional and ratoon rice distribution at subnational level. Given sufficient historical satellite and phenology records, the proposed algorithm had the potential to enhance rice cropping area mapping efforts across a broad temporal scale ( e.g ., from the 1980s to the present).

Mapping paddy rice and rice phenology with Sentinel-1 SAR time series using a unified dynamic programming framework

AbstractMonitoring rice planting areas and their phenological phases is crucial for yield estimation and informed decision-making. This study proposed a unified method for mapping rice field and rice phenology with a dynamic time wrapping (DTW) distance-based classifier and its variant sub-DTW algorithm using Sentinel-1’s synthetic aperture radar (SAR) VH band. Field samplings were conducted for broad landcover types in one of the areas of interest (AOIs). We implemented a pixel-wisek-nearest neighbor classification model with DTW distance to identify paddy rice pixels. Standard rice phenological profiles of the SAR VH band were defined by ground monitoring of a sample rice field. Based on rice planting maps and the standard phenological profiles, rice phenological phases were estimated by pattern matching strategy with the sub-DTW algorithm. Experiments on six counties in Northeast China presented promising results. The overall producer and user accuracy reached 92.9 and 91.9% for rice mapping, respectively. The mean root mean square error (RMSE) for phenology estimation was 3.5 days. Rice planting and rice phenology maps were generated for the six AOIs. The phenological variances of the AOIs implied the effects of climate and rice cultivars on phenological development.

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

Rice is the staple crop for more than half the world’s population, but there is a lack of high-resolution maps outlining rice areas and their growth stages. Most remote sensing studies map the rice extent; however, in tropical regions, rice is grown throughout the year with variable planting dates and cropping frequency. Thus, mapping rice growth stages is more useful than mapping only the extent. This study addressed this challenge by developing a phenology-based method. The hypothesis was that the unsupervised classification (k-means clustering) of Sentinel-1 and 2 time-series data could identify rice fields and growth stages, because (1) the presence of flooding during transplanting can be identified by Sentinel-1 VH backscatter; and (2) changes in the canopy of rice fields during growth stages (vegetative, generative, and ripening phases) up to the point of harvesting can be identified by Normalized Difference Vegetation Index (NDVI) time series. Using the proposed method, this study mapped rice field extent and cropping calendars across Peninsular Malaysia (131,598 km2) on the Google Earth Engine (GEE) platform. The Sentinel-1 and 2 monthly time series data from January 2019 to December 2020 were classified using k-means clustering to identify areas with similar phenological patterns. This approach resulted in 10-meter resolution maps of rice field extent, intensity, and cropping calendars. Validation using very high-resolution street view images from Google Earth showed that the predicted map had an overall accuracy of 95.95%, with a kappa coefficient of 0.92. In addition, the predicted crop calendars agreed well with the local government’s granary data. The results show that the proposed phenology-based method is cost-effective and can accurately map rice fields and growth stages over large areas. The information will be helpful in measuring the achievement of self-sufficiency in rice production and estimates of methane emissions from rice cultivation.

Regional water-saving potential calculation method for paddy rice based on remote sensing

To improve the calculation applicability and operability of regional water-saving potential (RWSP) for paddy rice, a calculation method based on remote sensing (RWSP-RS) was proposed. RWSP-RS consists of three sections: (a) paddy rice mapping by the decision tree algorithm, (b) rice evapotranspiration (ET) inversion under different irrigation modes by the surface energy balance algorithm for land (SEBAL), and (c) WSP based on ET (WSPE) and irrigation (WSPI) calculation by coupling water balance models for paddy fields. The RWSP-RS was applied in the Zhanghe Irrigation District in southern China in 2018 and 2019. The results showed that the three sections of RWSP-RS had high precision: paddy rice mapping errors ranged from 2% to 16%; WSPE of paddy rice errors were 26 mm and 5 mm for 2018 and 2019, respectively; and WSPI errors were 5 mm and 23 mm for 2018 and 2019, respectively. The WSPI of paddy rice in the whole region was 44.52 million m3 and 99.12 million m3 for 2018 and 2019, respectively. RWSP-RS has the characteristics of solid operability, good regional applicability, and time and labor savings, making it a recommended method for calculating the RWSP of paddy rice and contributing to regional water resource management.

Mapping Paddy Rice Distribution and Cropping Intensity in China from 2014 to 2019 with Landsat Images, Effective Flood Signals, and Google Earth Engine

Paddy rice cropping systems play a vital role in food security, water use, gas emission estimates, and grain yield prediction. Due to alterations in the labor structure and the high cost of paddy rice planting, the paddy rice cropping systems (single or double paddy rice) have drastically changed in China in recent years; many double-cropping paddy rice fields have been converted to single-cropping paddy rice or other crops, especially in southern China. Few maps detect single and double paddy rice and cropping intensity for paddy rice (CIPR) in China with a 30 m resolution. The Landsat-based and effective flooding signal-based phenology (EFSP) method, which distinguishes CIPR with the frequency of the effective flooding signal (EFe), was proposed and tested in China. The cloud/ice/shadow was excluded by bit arithmetic, generating a good observation map, and several non-paddy rice masks were established to improve the classification accuracy. Threshold values for single and double paddy rice were calculated through the mapped data and agricultural census data. Image processing (more than 684,000 scenes) and algorithm implementation were accomplished by a cloud computing approach with the Google Earth Engine (GEE) platform. The resultant maps of paddy rice from 2014 to 2019 were evaluated with data from statistical yearbooks and high-resolution images, with producer (user) accuracy and kappa coefficients ranging from 0.92 to 0.96 (0.76–0.87) and 0.67–0.80, respectively. Additionally, the determination coefficients for mapped and statistical data were higher than 0.88 from 2014 to 2019. Maps derived from EFSP illustrate that the single and double paddy rice systems are mainly concentrated in the Cfa (warm, fully humid, and hot summer, 49% vs. 56%) climate zone in China and show a slightly decreasing trend. The trend of double paddy rice is more pronounced than that of single paddy rice due to the high cost and shortages of rural household labor. However, single paddy rice fields expanded in Dwa (cold, dry winter, and hot summer, 11%) and Dwb (cold, dry winter, and warm summer, 9%) climate zones. The regional cropping intensity for paddy rice coincides with the paddy rice planting area but shows a significant decrease in south China, especially in Hunan Province, from 2014 to 2019. The results demonstrate that EFSP can effectively support the mapping of single and double paddy rice fields and CIPR in China, and the combinations of Landsat 7 and 8 provide enough good observations for EFSP to monitor paddy rice agriculture.

Large-Scale Rice Mapping Using Multi-Task Spatiotemporal Deep Learning and Sentinel-1 SAR Time Series

Timely and accurate cropland information at large spatial scales can improve crop management and support the government in decision making. Mapping the spatial extent and distribution of crops on a large spatial scale is challenging work due to the spatial variability. A multi-task spatiotemporal deep learning model, named LSTM-MTL, was developed in this study for large-scale rice mapping by utilizing time-series Sentinel-1 SAR data. The model showed a reasonable rice classification accuracy in the major rice production areas of the U.S. (OA = 98.3%, F1 score = 0.804), even when it only utilized SAR data. The model learned region-specific and common features simultaneously, and yielded a significant improved performance compared with RF and AtBiLSTM in both global and local training scenarios. We found that the LSTM-MTL model achieved a regional F1 score up to 10% higher than both global and local baseline models. The results demonstrated that the consideration of spatial variability via LSTM-MTL approach yielded an improved crop classification performance at a large spatial scale. We analyzed the input-output relationship through gradient backpropagation and found that low VH value in the early period and high VH value in the latter period were critical for rice classification. The results of in-season analysis showed that the model was able to yield a high accuracy (F1 score = 0.746) two months before rice maturity. The integration between multi-task learning and multi-temporal deep learning approach provides a promising approach for crop mapping at large spatial scales.

Selection of Candidate Genes Conferring Blast Resistance and Heat Tolerance in Rice through Integration of Meta-QTLs and RNA-Seq.

Due to global warming, high temperature is a significant environmental stress for rice production. Rice (Oryza sativa L.), one of the most crucial cereal crops, is also seriously devastated by Magnaporthe oryzae. Therefore, it is essential to breed new rice cultivars with blast and heat tolerance. Although progress had been made in QTL mapping and RNA-seq analysis in rice in response to blast and heat stresses, there are few reports on simultaneously mining blast-resistant and heat-tolerant genes. In this study, we separately conducted meta-analysis of 839 blast-resistant and 308 heat-tolerant QTLs in rice. Consequently, 7054 genes were identified in 67 blast-resistant meta-QTLs with an average interval of 1.00 Mb. Likewise, 6425 genes were obtained in 40 heat-tolerant meta-QTLs with an average interval of 1.49 Mb. Additionally, using differentially expressed genes (DEGs) in the previous research and GO enrichment analysis, 55 DEGs were co-located on the common regions of 16 blast-resistant and 14 heat-tolerant meta-QTLs. Among, OsChib3H-c, OsJAMyb, Pi-k, OsWAK1, OsMT2b, OsTPS3, OsHI-LOX, OsACLA-2 and OsGS2 were the significant candidate genes to be further investigated. These results could provide the gene resources for rice breeding with excellent resistance to these 2 stresses, and help to understand how plants response to the combination stresses of blast fungus and high temperature.

Spatio-Temporal Estimation of Rice Height Using Time Series Sentinel-1 Images

Rice height, as the fundamental biophysical attribute, is a controlling factor in crop phenology estimation and yield estimation. The aim of this study was to use time series Sentinel-1A images to estimate the spatio-temporal distribution of rice height. In this study, a particle filter (PF) was applied for the real-time estimation of rice height compared with a simplified water cloud model (SWCM) on the basis of rice mapping and transplanting date. It was found that the VH backscatter (σvho) can potentially be applied to accurately estimate rice height compared with VV backscatter (σvvo), the σvho/σvv0 ratio, and the Radar Vegetation Index (RVI, 4* σvho/(σvho+σvvo)). The results show that the rice height estimation by PF generated a better result with a root-mean-square error (RMSE) equal to 7.36 cm and a determination factor (R2) of 0.95 compared with SWCM (RMSE = 12.59 cm and R2 = 0.86). Moreover, rice height in the south and east of the study area was higher than in the north and west. The reason for this is that the south and east are near to the South China Sea, and there are higher temperatures and earlier transplanting. Altogether, our results demonstrate the potential of PF and σvho to study the spatio-temporal distribution of crop height estimation. As a result, the PF method can contribute greatly to improvements in crop monitoring.

Whole-Genome Sequencing of 117 Chromosome Segment Substitution Lines for Genetic Analyses of Complex Traits in Rice

Rice is one of the most important food crops in Asia. Genetic analyses of complex traits and molecular breeding studies in rice greatly rely on the construction of various genetic populations. Chromosome segment substitution lines (CSSLs) serve as a powerful genetic population for quantitative trait locus (QTL) mapping in rice. Moreover, CSSLs containing target genomic regions can be used as improved varieties in rice breeding. In this study, we developed a set of CSSLs consisting of 117 lines derived from the recipient ‘Huanghuazhan’ (HHZ) and the donor ‘Basmati Surkb 89–15’ (BAS). The 117 lines were extensively genotyped by whole-genome resequencing, and a high-density genotype map was constructed for the CSSL population. The 117 CSSLs covered 99.78% of the BAS genome. Each line contained a single segment, and the average segment length was 6.02 Mb. Using the CSSL population, we investigated three agronomic traits in Shanghai and Hangzhou, China, and a total of 25 QTLs were detected in both environments. Among those QTLs, we found that RFT1 was the causal gene for heading date variance between HHZ and BAS. RFT1 from BAS was found to contain a loss-of-function allele based on yeast two-hybrid assay, and its causal variation was a P to S change in the 94th amino acid of the RFT1 protein. The combination of high-throughput genotyping and marker-assisted selection (MAS) is a highly efficient way to construct CSSLs in rice, and extensively genotyped CSSLs will be a powerful tool for the genetic mapping of agronomic traits and molecular breeding for target QTLs/genes.

Rice Mapping in Training Sample Shortage Regions Using a Deep Semantic Segmentation Model Trained on Pseudo-Labels

A deep semantic segmentation model-based method can achieve state-of-the-art accuracy and high computational efficiency in large-scale crop mapping. However, the model cannot be widely used in actual large-scale crop mapping applications, mainly because the annotation of ground truth data for deep semantic segmentation model training is time-consuming. At the operational level, it is extremely difficult to obtain a large amount of ground reference data by photointerpretation for the model training. Consequently, in order to solve this problem, this study introduces a workflow that aims to extract rice distribution information in training sample shortage regions, using a deep semantic segmentation model (i.e., U-Net) trained on pseudo-labels. Based on the time series Sentinel-1 images, Cropland Data Layer (CDL) and U-Net model, the optimal multi-temporal datasets for rice mapping were summarized, using the global search method. Then, based on the optimal multi-temporal datasets, the proposed workflow (a combination of K-Means and random forest) was directly used to extract the rice-distribution information of Jiangsu (i.e., the K–RF pseudo-labels). For comparison, the optimal well-trained U-Net model acquired from Arkansas (i.e., the transfer model) was also transferred to Jiangsu to extract local rice-distribution information (i.e., the TF pseudo-labels). Finally, the pseudo-labels with high confidences generated from the two methods were further used to retrain the U-Net models, which were suitable for rice mapping in Jiangsu. For different rice planting pattern regions of Jiangsu, the final results showed that, compared with the U-Net model trained on the TF pseudo-labels, the rice area extraction errors of pseudo-labels could be further reduced by using the U-Net model trained on the K–RF pseudo-labels. In addition, compared with the existing rule-based rice mapping methods, he U-Net model trained on the K–RF pseudo-labels could robustly extract the spatial distribution information of rice. Generally, this study could provide new options for applying a deep semantic segmentation model to training sample shortage regions.

From Cascade to Parallel: Automated Phenology-Guided Supervised Paddy Rice Mapping with Landsat and Sentinel Imagery in Google Earth Engine

From Cascade to Parallel: Automated Phenology-Guided Supervised Paddy Rice Mapping with Landsat and Sentinel Imagery in Google Earth Engine

Semantic Segmentation Based on Temporal Features: Learning of Temporal–Spatial Information From Time-Series SAR Images for Paddy Rice Mapping

Synthetic aperture radar (SAR) can be used to obtain remote sensing images of different growth stages of crops under all weather conditions. Such time-series SAR images can provide an abundance of temporal and spatial features for use in large-scale crop mapping and analysis. In this study, we propose a temporal feature-based segmentation (TFBS) model for accurate crop mapping using time-series SAR images. This model first extracts deep-seated temporal features and then learns the spatial context of the extracted temporal features for crop mapping. The results indicate that the TFBS model significantly outperforms traditional long short-term memory (LSTM), U-network, and convolutional LSTM models in crop mapping based on time-series SAR images. TFBS demonstrates better generalizability than other models in the study area, which makes it more transferable, and the results show that data augmentation can significantly improve this generalizability. The visualization of the temporal features extracted by the TFBS shows that there is a high degree of intraclass homogeneity among rice fields and interclass heterogeneity between rice fields and other features. TFBS also achieved the highest accuracy of the four deep learning models for multicrop classification in the study area. This study presents a feasible way of producing high-accuracy large-scale crop maps based on the proposed model.

Paddy Rice Mapping Using a Dual-Path Spatio-Temporal Network Based on Annual Time-Series Sentinel-2 Images

Paddy rice is one of the main foods of the global population. To guarantee paddy rice acreage is essential to ensure food security. Currently, techniques for large-area paddy field mapping rely mainly on complex rule-based machine learning algorithms. But it is difficult for them to achieve an optimal balance between discriminability and robustness. In this article, we proposed a novel deep learning-based approach for large-scale paddy rice mapping, termed dual-path interactive network (DPIN). An annual time-series Sentinel-2 remote sensing images are used as data source. Taking several areas of interest over the middle and lower Yangtze River plain of China as experimental fields, our model achieves an F1-score of 91.22% on the test dataset, which is 1.09% higher than the existing state-of-the-art predictive model, and its inference speed is 1.18 times faster than it. DPIN-Lite is a lightweight variant of DPIN, and while keeping a competitive mapping accuracy, its inference speed is 1.91 times faster than the compared method (with the best score except for DPIN and DPIN-Lite).

Mapping of Rice Field to Support Food Safety at Banyumas Regency

The rice area in Banyumas Regency needs attention in an effort to support food safety and avoid land conversion. The objective of this study was to update the information through mapping rice land area at Banyumas Regency based on the type of irrigation, altitude and soil type. Geographic Information System (GIS) ArcView 3.2 -based software was applied, then analyzed by using the overlay method and direct checking/surveying in the field to process of analyzing spatial data of paddy fields in Banyumas Regency. Mapping was carried out on three variables, namely irrigation type, altitude and soil type. The results showed that the irrigated rice field type is more dominant than the rainfed type, the altitude is lowland with an altitude of 25-250 m above sea level and the dominant soil types are andosols and latosols with an area of> 12,000 ha.

A Novel Weighted Amplitude Modulation (WAM) System for Ambiguity Suppression of Spaceborne Hybrid Quad-Pol SAR

Quadrature-polarimetric synthetic aperture radar (quad-pol SAR) has extensive applications, including climate zones classification, extraction of surface roughness, soil moisture mapping, moving target indication, and rice mapping. Hybrid quad-pol SAR ameliorates the range ambiguity performance of conventional quad-pol SAR. However, the azimuth ambiguity of its cross-polarized (cross-pol) echo signals is serious, limiting the swath width of SAR. Therefore, this paper proposes a spaceborne weighted amplitude modulation (WAM) full-polarimetric (full-pol) SAR system, and it can suppress the azimuth ambiguity of hybrid quad-pol SAR. The performance boost of the azimuth ambiguity by the two imaging modes of the proposed SAR system is detailed and evaluated with the L-band system parameters. Moreover, the chirp scaling algorithm (CSA) is adopted to execute scene simulations for the two imaging modes. The results indicate that the proposed SAR system can effectively suppress the azimuth ambiguity of hybrid quad-pol SAR and verify the theoretical analysis.

NESEA-Rice10: high-resolution annual paddy rice maps for Northeast and Southeast Asia from 2017 to 2019

Abstract. An accurate paddy rice map is crucial for ensuring food security, particularly for Southeast and Northeast Asia. MODIS satellite data are useful for mapping paddy rice at continental scales but have a mixed-pixel problem caused by the coarse spatial resolution. To reduce the mixed pixels, we designed a rule-based method for mapping paddy rice by integrating time series Sentinel-1 and MODIS data. We demonstrated the method by generating annual paddy rice maps for Southeast and Northeast Asia in 2017–2019 (NESEA-Rice10). We compared the resultant paddy rice maps with available agricultural statistics at subnational levels and existing rice maps for some countries. The results demonstrated that the linear coefficient of determination (R2) between our paddy rice maps and agricultural statistics ranged from 0.80 to 0.97. The paddy rice planting areas in 2017 were spatially consistent with the existing maps in Vietnam (R2=0.93) and Northeast China (R2=0.99). The spatial distribution of the 2017–2019 composite paddy rice map was consistent with that of the rice map from the International Rice Research Institute. The paddy rice planting area may have been underestimated in the region in which the flooding signal was not strong. The dataset is useful for water resource management, rice growth, and yield monitoring. The full product is publicly available at https://doi.org/10.5281/zenodo.5645344 (Han et al., 2021a). Small examples can be found from the following DOI: https://doi.org/10.17632/cnc3tkbwcm.1 (Han et al., 2021b).

Fine Classification of Rice Paddy Based on RHSI-DT Method Using Multi-Temporal Compact Polarimetric SAR Data

In recent years, the compact polarimetric (CP) synthetic aperture radar (SAR) has become a hotspot of SAR Earth observation. Meanwhile, CP SAR provides both relatively rich polarization information and large swath-width for rice mapping. Fine classification of rice paddy plays an important role in growth monitoring, pest prevention and yield estimation of rice. In this study, the multi-temporal CP SAR data were firstly simulated by fully polarimetric RADARSAT-2 data, and 22 CP parameters from each of the six temporal CP SAR data were extracted. Then we built a rice height-sensitive index (RHSI). Furthermore, a decision tree (DT) method was established by using the optimal CP parameters based on RHSI. Finally, the classification results of rice paddy based on DT and support vector machine (SVM) methods were compared. Results showed that the RHSI-DT method could obtain better results, with an overall accuracy of 97.94% and a kappa coefficient of 0.973, which was 2% higher and 0.03 larger than those of the SVM method. Besides, we found that the surface scattering of m-χ decomposition (m-χ_s (0627)) and ΔShannon entropy intensity Hi (Hi (1015)-Hi (0627)) were highly effective parameters to distinguish paddies of transplanting hybrid rice (T-H) and direct-sown japonica rice (D-J).

Registration of two rice mapping populations using weedy rice ecotypes as a novel germplasm resource

AbstractTwo mapping populations were developed from crosses of the Asian indica rice (Oryza sativa L.) cultivar ‘Dee Geo Woo Gen’ (DGWG; PI 699210 Parent, PI 699212 Parent) and two weedy rice ecotypes, an early‐flowering straw hull (SH) biotype AR‐2000‐1135‐01 (PI 699209 Parent) collected in Arkansas and a late‐flowering black hull (BHA) biotype MS‐1996‐9 (PI 699211 Parent) collected in Mississippi. The weed and crop‐based rice recombinant inbred line (RIL) mapping populations have been used to identify genomic regions associated with weedy traits as well as resistance to sheath blight and rice blast diseases. The mapping population consists of 185 (DGWG/SH; Reg. no. MP‐9, NSL 541035 MAP) and 234 (BHA/DGWG; Reg. no. MP‐10, NSL 541036 MAP) F8 RILs, of which 175 (DGWG/SH) and 224 (BHA/DGWG) were used to construct two linkage maps using single nucleotide polymorphic markers to identify weedy traits, sheath blight, and blast resistance loci. These mapping populations and related datasets represent a valuable resource for basic rice evolutionary genomic research and applied marker‐assisted breeding efforts in disease resistance.

Characterization of Genes Associated with Salt Tolerance Using Transcriptome Analysis and Quantitative Trait Loci Mapping in Rice

Characterization of Genes Associated with Salt Tolerance Using Transcriptome Analysis and Quantitative Trait Loci Mapping in Rice

High Resolution Distribution Dataset of Double-Season Paddy Rice in China

Although China is the largest producer of rice, accounting for about 25% of global production, there are no high-resolution maps of paddy rice covering the entire country. Using time-weighted dynamic time warping (TWDTW), this study developed a pixel- and phenology-based method to identify planting areas of double-season paddy rice in China, by comparing temporal variations of synthetic aperture radar (SAR) signals of unknown pixels to those of known double-season paddy rice fields. We conducted a comprehensive evaluation of the method’s performance at pixel and regional scales. Based on 145,210 field surveyed samples from 2018 to 2020, the producer’s and user’s accuracy are 88.49% and 87.02%, respectively. Compared to county-level statistical data from 2016 to 2019, the relative mean absolute errors are 34.11%. This study produced distribution maps of double-season rice at 10 m spatial resolution from 2016 to 2020 over nine provinces in South China, which account for more than 99% of the planting areas of double-season paddy rice of China. The maps are expected to contribute to timely monitoring and evaluating rice growth and yield.