Improved Estimates of Arctic Land Surface Phenology Using Sentinel-2 Time Series

The high spatial resolution and revisit time of Sentinel-2A/B tandem satellites allow a potentially improved retrieval of land surface phenology (LSP). The biome and regional characteristics, however, greatly constrain the design of the LSP algorithms. In the Arctic, such biome-specific characteristics include prolonged periods of snow cover, persistent cloud cover, and shortness of the growing season. Here, we evaluate the feasibility of Sentinel-2 for deriving high-resolution LSP maps of the Arctic. We extracted the timing of the start and end of season (SoS and EoS, respectively) for the years 2019 and 2020 with a simple implementation of the threshold method in Google Earth Engine (GEE). We found a high level of similarity between Sentinel-2 and PhenoCam metrics; the best results were observed with Sentinel-2 enhanced vegetation index (EVI) (root mean squared error (RMSE) and mean error (ME) of 3.0 d and −0.3 d for the SoS, and 6.5 d and −3.8 d for the EoS, respectively), although other vegetation indices presented similar performances. The phenological maps of Sentinel-2 EVI compared well with the same maps extracted from the Moderate Resolution Imaging Spectroradiometer (MODIS) in homogeneous landscapes (RMSE and ME of 9.2 d and 2.9 d for the SoS, and 6.4 and −0.9 d for the EoS, respectively). Unreliable LSP estimates were filtered and a quality flag indicator was activated when the Sentinel-2 time series presented a long period (>40 d) of missing data; discontinuities were lower in spring and early summer (9.2%) than in late summer and autumn (39.4%). The Sentinel-2 high-resolution LSP maps and the GEE phenological extraction method will support vegetation monitoring and contribute to improving the representation of Artic vegetation phenology in land surface models.

Download Full-text

An updated national-scale mangrove forest map in China Using Sentinel-1 and Sentinel-2 Time-Series Data with Google Earth Engine

10.5194/egusphere-egu21-5310 ◽

2021 ◽

Author(s):

Luojia Hu ◽

Wei Yao ◽

Zhitong Yu ◽

Yan Huang

Keyword(s):

Time Series ◽

High Resolution ◽

Mangrove Forest ◽

Google Earth ◽

Series Data ◽

National Scale ◽

Mangrove Area ◽

Google Earth Engine ◽

Forest Maps ◽

Sentinel 2

<p>A high resolution mangrove map (e.g., 10-m), which can identify mangrove patches with small size (< 1 ha), is a central component to quantify ecosystem functions and help government take effective steps to protect mangroves, because the increasing small mangrove patches, due to artificial destruction and plantation of new mangrove trees, are vulnerable to climate change and sea level rise, and important for estimating mangrove habitat connectivity with adjacent coastal ecosystems as well as reducing the uncertainty of carbon storage estimation. However, latest national scale mangrove forest maps mainly derived from Landsat imagery with 30-m resolution are relatively coarse to accurately characterize the distribution of mangrove forests, especially those of small size (area < 1 ha). Sentinel imagery with 10-m resolution provide the opportunity for identifying these small mangrove patches and generating high-resolution mangrove forest maps. Here, we used spectral/backscatter-temporal variability metrics (quantiles) derived from Sentinel-1 SAR (Synthetic Aperture Radar) and sentinel-2 MSI (Multispectral Instrument) time-series imagery as input features for random forest to classify mangroves in China. We found that Sentinel-2 imagery is more effective than Sentinel-1 in mangrove extraction, and a combination of SAR and MSI imagery can get a better accuracy (F1-score of 0.94) than using them separately (F1-score of 0.88 using Sentinel-1 only and 0.895 using Sentinel-2 only). The 10-m mangrove map derived by combining SAR and MSI data identified 20,003 ha mangroves in China and the areas of small mangrove patches (< 1 ha) was 1741 ha, occupying 8.7% of the whole mangrove area. The largest area (819 ha) of small mangrove patches is located in Guangdong Province, and in Fujian the percentage of small mangrove patches in total mangrove area is the highest (11.4%). A comparison with existing 30-m mangrove products showed noticeable disagreement, indicating the necessity for generating mangrove extent product with 10-m resolution. This study demonstrates the significant potential of using Sentinel-1 and Sentinel-2 images to produce an accurate and high-resolution mangrove forest map with Google Earth Engine (GEE). The mangrove forest maps are expected to provide critical information to conservation managers, scientists, and other stakeholders in monitoring the dynamics of mangrove forest.</p>

Download Full-text

Fusing Geostationary Satellite Observations with Harmonized Landsat-8 and Sentinel-2 Time Series for Monitoring Field-Scale Land Surface Phenology

Remote Sensing ◽

10.3390/rs13214465 ◽

2021 ◽

Vol 13 (21) ◽

pp. 4465

Author(s):

Yu Shen ◽

Xiaoyang Zhang ◽

Weile Wang ◽

Ramakrishna Nemani ◽

Yongchang Ye ◽

...

Keyword(s):

Time Series ◽

Land Surface ◽

Cloud Cover ◽

Landsat 8 ◽

Field Scale ◽

Spatiotemporal Resolution ◽

Essential Information ◽

Land Surface Phenology ◽

The Earth ◽

Sentinel 2

Accurate and timely land surface phenology (LSP) provides essential information for investigating the responses of terrestrial ecosystems to climate changes and quantifying carbon and surface energy cycles on the Earth. LSP has been widely investigated using daily Visible Infrared Imaging Radiometer Suite (VIIRS) or Moderate Resolution Imaging Spectroradiometer (MODIS) observations, but the resultant phenometrics are frequently influenced by surface heterogeneity and persistent cloud contamination in the time series observations. Recently, LSP has been derived from Landsat-8 and Sentinel-2 time series providing detailed spatial pattern, but the results are of high uncertainties because of poor temporal resolution. With the availability of data from Advanced Baseline Imager (ABI) onboard a new generation of geostationary satellites that observe the earth every 10–15 min, daily cloud-free time series could be obtained with high opportunities. Therefore, this study investigates the generation of synthetic high spatiotemporal resolution time series by fusing the harmonized Landsat-8 and Sentinel-2 (HLS) time series with the temporal shape of ABI data for monitoring field-scale (30 m) LSP. The algorithm is verified by detecting the timings of greenup and senescence onsets around north Wisconsin/Michigan states, United States, where cloud cover is frequent during spring rainy season. The LSP detections from HLS-ABI are compared with those from HLS or ABI alone and are further evaluated using PhenoCam observations. The result indicates that (1) ABI could provide ~3 times more high-quality observations than HLS around spring greenup onset; (2) the greenup and senescence onsets derived from ABI and HLS-ABI are spatially consistent and statistically comparable with a median difference less than 1 and 10-days, respectively; (3) greenup and senescence onsets derived from HLS data show sharp boundaries around the orbit-overlapped areas and shifts of ~13 days delay and ~15 days ahead, respectively, relative to HLS-ABI detections; and (4) HLS-ABI greenup and senescence onsets align closely to PhenoCam observations with an absolute average difference of less than 2 days and 5 days, respectively, which are much better than phenology detections from ABI or HLS alone. The result suggests that the proposed approach could be implemented the monitor of 30 m LSP over regions with persistent cloud cover.

Download Full-text

Advancing the Mapping of Mangrove Forests at National-Scale Using Sentinel-1 and Sentinel-2 Time-Series Data with Google Earth Engine: A Case Study in China

Remote Sensing ◽

10.3390/rs12193120 ◽

2020 ◽

Vol 12 (19) ◽

pp. 3120

Author(s):

Luojia Hu ◽

Nan Xu ◽

Jian Liang ◽

Zhichao Li ◽

Luzhen Chen ◽

...

Keyword(s):

Time Series ◽

High Resolution ◽

Mangrove Forest ◽

Time Series Data ◽

Landsat Imagery ◽

Google Earth ◽

Series Data ◽

National Scale ◽

Google Earth Engine ◽

Sentinel 2

A high resolution mangrove map (e.g., 10-m), including mangrove patches with small size, is urgently needed for mangrove protection and ecosystem function estimation, because more small mangrove patches have disappeared with influence of human disturbance and sea-level rise. However, recent national-scale mangrove forest maps are mainly derived from 30-m Landsat imagery, and their spatial resolution is relatively coarse to accurately characterize the extent of mangroves, especially those with small size. Now, Sentinel imagery with 10-m resolution provides an opportunity for generating high-resolution mangrove maps containing these small mangrove patches. Here, we used spectral/backscatter-temporal variability metrics (quantiles) derived from Sentinel-1 SAR (Synthetic Aperture Radar) and/or Sentinel-2 MSI (Multispectral Instrument) time-series imagery as input features of random forest to classify mangroves in China. We found that Sentinel-2 (F1-Score of 0.895) is more effective than Sentinel-1 (F1-score of 0.88) in mangrove extraction, and a combination of SAR and MSI imagery can get the best accuracy (F1-score of 0.94). The 10-m mangrove map was derived by combining SAR and MSI data, which identified 20003 ha mangroves in China, and the area of small mangrove patches (<1 ha) is 1741 ha, occupying 8.7% of the whole mangrove area. At the province level, Guangdong has the largest area (819 ha) of small mangrove patches, and in Fujian, the percentage of small mangrove patches is the highest (11.4%). A comparison with existing 30-m mangrove products showed noticeable disagreement, indicating the necessity for generating mangrove extent product with 10-m resolution. This study demonstrates the significant potential of using Sentinel-1 and Sentinel-2 images to produce an accurate and high-resolution mangrove forest map with Google Earth Engine (GEE). The mangrove forest map is expected to provide critical information to conservation managers, scientists, and other stakeholders in monitoring the dynamics of the mangrove forest.

Download Full-text

High-Resolution Urban Land Mapping in China from Sentinel 1A/2 Imagery Based on Google Earth Engine

Remote Sensing ◽

10.3390/rs11070752 ◽

2019 ◽

Vol 11 (7) ◽

pp. 752 ◽

Cited By ~ 21

Author(s):

Zhongchang Sun ◽

Ru Xu ◽

Wenjie Du ◽

Lei Wang ◽

Dengsheng Lu

Keyword(s):

High Resolution ◽

Vegetation Index ◽

Urban Land ◽

Google Earth ◽

Economic Research ◽

Optical Data ◽

Large Area ◽

Google Earth Engine ◽

Sar Data ◽

Sentinel 2

Accurate and timely urban land mapping is fundamental to supporting large area environmental and socio-economic research. Most of the available large-area urban land products are limited to a spatial resolution of 30 m. The fusion of optical and synthetic aperture radar (SAR) data for large-area high-resolution urban land mapping has not yet been widely explored. In this study, we propose a fast and effective urban land extraction method using ascending/descending orbits of Sentinel-1A SAR data and Sentinel-2 MSI (MultiSpectral Instrument, Level 1C) optical data acquired from 1 January 2015 to 30 June 2016. Potential urban land (PUL) was identified first through logical operations on yearly mean and standard deviation composites from a time series of ascending/descending orbits of SAR data. A Yearly Normalized Difference Vegetation Index (NDVI) maximum and modified Normalized Difference Water Index (MNDWI) mean composite were generated from Sentinel-2 imagery. The slope image derived from SRTM DEM data was used to mask mountain pixels and reduce the false positives in SAR data over these regions. We applied a region-specific threshold on PUL to extract the target urban land (TUL) and a global threshold on the MNDWI mean, and slope image to extract water bodies and high-slope regions. A majority filter with a three by three window was applied on previously extracted results and the main processing was carried out on the Google Earth Engine (GEE) platform. China was chosen as the testing region to validate the accuracy and robustness of our proposed method through 224,000 validation points randomly selected from high-resolution Google Earth imagery. Additionally, a total of 735 blocks with a size of 900 × 900 m were randomly selected and used to compare our product’s accuracy with the global human settlement layer (GHSL, 2014), GlobeLand30 (2010), and Liu (2015) products. Our method demonstrated the effectiveness of using a fusion of optical and SAR data for large area urban land extraction especially in areas where optical data fail to distinguish urban land from spectrally similar objects. Results show that the average overall, producer’s and user’s accuracies are 88.03%, 94.50% and 82.22%, respectively.

Download Full-text

Influence of Varying Solar Zenith Angles on Land Surface Phenology Derived from Vegetation Indices: A Case Study in the Harvard Forest

Remote Sensing ◽

10.3390/rs13204126 ◽

2021 ◽

Vol 13 (20) ◽

pp. 4126

Author(s):

Yang Li ◽

Ziti Jiao ◽

Kaiguang Zhao ◽

Yadong Dong ◽

Yuyu Zhou ◽

...

Keyword(s):

Time Series ◽

Land Surface ◽

Vegetation Index ◽

Vegetation Indices ◽

Mean Difference ◽

Land Surface Phenology ◽

Harvard Forest ◽

Ndvi Time Series ◽

Bidirectional Reflectance Distribution

Vegetation indices are widely used to derive land surface phenology (LSP). However, due to inconsistent illumination geometries, reflectance varies with solar zenith angles (SZA), which in turn affects the vegetation indices, and thus the derived LSP. To examine the SZA effect on LSP, the MODIS bidirectional reflectance distribution function (BRDF) product and a BRDF model were employed to derive LSPs under several constant SZAs (i.e., 0°, 15°, 30°, 45°, and 60°) in the Harvard Forest, Massachusetts, USA. The LSPs derived under varying SZAs from the MODIS nadir BRDF-adjusted reflectance (NBAR) and MODIS vegetation index products were used as baselines. The results show that with increasing SZA, NDVI increases but EVI decreases. The magnitude of SZA-induced NDVI/EVI changes suggests that EVI is more sensitive to varying SZAs than NDVI. NDVI and EVI are comparable in deriving the start of season (SOS), but EVI is more accurate when deriving the end of season (EOS). Specifically, NDVI/EVI-derived SOSs are relatively close to those derived from ground measurements, with an absolute mean difference of 8.01 days for NDVI-derived SOSs and 9.07 days for EVI-derived SOSs over ten years. However, a considerable lag exists for EOSs derived from vegetation indices, especially from the NDVI time series, with an absolute mean difference of 14.67 days relative to that derived from ground measurements. The SOSs derived from NDVI time series are generally earlier, while those from EVI time series are delayed. In contrast, the EOSs derived from NDVI time series are delayed; those derived from the simulated EVI time series under a fixed illumination geometry are also delayed, but those derived from the products with varying illumination geometries (i.e., MODIS NBAR product and MODIS vegetation index product) are advanced. LSPs derived from varying illumination geometries could lead to a difference spanning from a few days to a month in this case study, which highlights the importance of normalizing the illumination geometry when deriving LSP from NDVI/EVI time series.

Download Full-text

Incorporating the Plant Phenological Trajectory into Mangrove Species Mapping with Dense Time Series Sentinel-2 Imagery and the Google Earth Engine Platform

Remote Sensing ◽

10.3390/rs11212479 ◽

2019 ◽

Vol 11 (21) ◽

pp. 2479 ◽

Cited By ~ 5

Author(s):

Huiying Li ◽

Mingming Jia ◽

Rong Zhang ◽

Yongxing Ren ◽

Xin Wen

Keyword(s):

Time Series ◽

Vegetation Index ◽

Google Earth ◽

Late Winter ◽

Mangrove Forests ◽

Mangrove Species ◽

Ndvi Time Series ◽

Google Earth Engine ◽

Species Mapping ◽

Sentinel 2

Information on mangrove species composition and distribution is key to studying functions of mangrove ecosystems and securing sustainable mangrove conservation. Even though remote sensing technology is developing rapidly currently, mapping mangrove forests at the species level based on freely accessible images is still a great challenge. This study built a Sentinel-2 normalized difference vegetation index (NDVI) time series (from 2017-01-01 to 2018-12-31) to represent phenological trajectories of mangrove species and then demonstrated the feasibility of phenology-based mangrove species classification using the random forest algorithm in the Google Earth Engine platform. It was found that (i) in Zhangjiang estuary, the phenological trajectories (NDVI time series) of different mangrove species have great differences; (ii) the overall accuracy and Kappa confidence of the classification map is 84% and 0.84, respectively; and (iii) Months in late winter and early spring play critical roles in mangrove species mapping. This is the first study to use phonological signatures in discriminating mangrove species. The methodology presented can be used as a practical guideline for the mapping of mangrove or other vegetation species in other regions. However, future work should pay attention to various phenological trajectories of mangrove species in different locations.

Download Full-text

Monitoring Broadscale Vegetational Diversity and Change across North American Landscapes Using Land Surface Phenology

Forests ◽

10.3390/f11060606 ◽

2020 ◽

Vol 11 (6) ◽

pp. 606

Author(s):

Bjorn-Gustaf J. Brooks ◽

Danny C. Lee ◽

Lars Y. Pomara ◽

William W. Hargrove

Keyword(s):

Time Series ◽

Factor Analysis ◽

Coordinate Transformation ◽

Land Surface ◽

Vegetation Index ◽

Normalized Difference Vegetation Index ◽

Ecological Gradients ◽

Land Surface Phenology ◽

Polar Coordinate ◽

Vegetational Diversity

We describe a polar coordinate transformation of vegetation index profiles which permits a broad-scale comparison of location-specific phenological variability influenced by climate, topography, land use, and other factors. We apply statistical data reduction techniques to identify fundamental dimensions of phenological variability and to classify phenological types with intuitive ecological interpretation. Remote sensing-based land surface phenology can reveal ecologically meaningful vegetational diversity and dynamics across broad landscapes. Land surface phenology is inherently complex at regional to continental scales, varying with latitude, elevation, and multiple biophysical factors. Quantifying phenological change across ecological gradients at these scales is a potentially powerful way to monitor ecological development, disturbance, and diversity. Polar coordinate transformation was applied to Moderate Resolution Imaging Spectroradiometer (MODIS) normalized difference vegetation index (NDVI) time series spanning 2000-2018 across North America. In a first step, 46 NDVI values per year were reduced to 11 intuitive annual metrics, such as the midpoint of the growing season and degree of seasonality, measured relative to location-specific annual phenological cycles. Second, factor analysis further reduced these metrics to fundamental phenology dimensions corresponding to annual timing, productivity, and seasonality. The factor analysis explained over 95% of the variability in the metrics and represented a more than ten-fold reduction in data volume from the original time series. In a final step, phenological classes (‘phenoclasses’) based on the statistical clustering of the factor data, were computed to describe the phenological state of each pixel during each year, which facilitated the tracking of year-to-year dynamics. Collectively the phenology metrics, factors, and phenoclasses provide a system for characterizing land surface phenology and for monitoring phenological change that is indicative of ecological gradients, development, disturbance, and other aspects of landscape-scale diversity and dynamics.

Download Full-text

Time-Series of Cloud-Free Sentinel-2 NDVI Data Used in Mapping the Onset of Growth of Central Spitsbergen, Svalbard

Remote Sensing ◽

10.3390/rs13153031 ◽

2021 ◽

Vol 13 (15) ◽

pp. 3031

Author(s):

Stein Rune Karlsen ◽

Laura Stendardi ◽

Hans Tømmervik ◽

Lennart Nilsen ◽

Ingar Arntzen ◽

...

Keyword(s):

Climate Change ◽

Time Series ◽

Cloud Cover ◽

The Arctic ◽

Vegetation Types ◽

Threshold Method ◽

High Increase ◽

Temperature Changes ◽

Sentinel 2 ◽

Phenological Phases

The Arctic is a region that is expected to experience a high increase in temperature. Changes in the timing of phenological phases, such as the onset of growth (as observed by remote sensing), is a sensitive bio-indicator of climate change. In this paper, the study area was the central part of Spitsbergen, Svalbard, located between 77.28°N and 78.44°N. The goals of this study were: (1) to prepare, analyze and present a cloud-free time-series of daily Sentinel-2 NDVI datasets for the 2016 to 2019 seasons, and (2) to demonstrate the use of the dataset in mapping the onset of growth. Due to a short and intense period with greening-up and frequent cloud cover, all the cloud-free Sentinel-2 data were used. The onset of growth was then mapped by a NDVI threshold method, which showed significant correlation (r2 = 0.47, n = 38, p < 0.0001) with ground-based phenocam observation of the onset of growth in seven vegetation types. However, large bias was found between the Sentinel-2 NDVI-based mapped onset of growth and the phenocam-based onset of growth in a moss tundra, which indicates that the data in these vegetation types must be interpreted with care. In 2018, the onset of growth was about 10 days earlier compared to 2017.

Download Full-text

Assessing Spatiotemporal Variations of Landsat Land Surface Temperature and Multispectral Indices in the Arctic Mackenzie Delta Region between 1985 and 2018

Remote Sensing ◽

10.3390/rs11192329 ◽

2019 ◽

Vol 11 (19) ◽

pp. 2329 ◽

Cited By ~ 2

Author(s):

Leon Nill ◽

Tobias Ullmann ◽

Christof Kneisel ◽

Jennifer Sobiech-Wolf ◽

Roland Baumhauer

Keyword(s):

Land Surface ◽

Large Scale ◽

Vegetation Index ◽

Google Earth ◽

Influential Factors ◽

The Arctic ◽

Moisture Regime ◽

Mackenzie Delta ◽

Delta Region

Air temperatures in the Arctic have increased substantially over the last decades, which has extensively altered the properties of the land surface. Capturing the state and dynamics of Land Surface Temperatures (LSTs) at high spatial detail is of high interest as LST is dependent on a variety of surficial properties and characterizes the land–atmosphere exchange of energy. Accordingly, this study analyses the influence of different physical surface properties on the long-term mean of the summer LST in the Arctic Mackenzie Delta Region (MDR) using Landsat 30 m-resolution imagery between 1985 and 2018 by taking advantage of the cloud computing capabilities of the Google Earth Engine. Multispectral indices, including the Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI) and Tasseled Cap greenness (TCG), brightness (TCB), and wetness (TCW) as well as topographic features derived from the TanDEM-X digital elevation model are used in correlation and multiple linear regression analyses to reveal their influence on the LST. Furthermore, surface alteration trends of the LST, NDVI, and NDWI are revealed using the Theil-Sen (T-S) regression method. The results indicate that the mean summer LST appears to be mostly influenced by the topographic exposition as well as the prevalent moisture regime where higher evapotranspiration rates increase the latent heat flux and cause a cooling of the surface, as the variance is best explained by the TCW and northness of the terrain. However, fairly diverse model outcomes for different regions of the MDR (R2 from 0.31 to 0.74 and RMSE from 0.51 °C to 1.73 °C) highlight the heterogeneity of the landscape in terms of influential factors and suggests accounting for a broad spectrum of different factors when modeling mean LSTs. The T-S analysis revealed large-scale wetting and greening trends with a mean decadal increase of the NDVI/NDWI of approximately +0.03 between 1985 and 2018, which was mostly accompanied by a cooling of the land surface given the inverse relationship between mean LSTs and vegetation and moisture conditions. Disturbance through wildfires intensifies the surface alterations locally and lead to significantly cooler LSTs in the long-term compared to the undisturbed surroundings.

Download Full-text

Development of a New Phenology Algorithm for Fine Mapping of Cropping Intensity in Complex Planting Areas Using Sentinel-2 and Google Earth Engine

ISPRS International Journal of Geo-Information ◽

10.3390/ijgi10090587 ◽

2021 ◽

Vol 10 (9) ◽

pp. 587

Author(s):

Yan Guo ◽

Haoming Xia ◽

Li Pan ◽

Xiaoyang Zhao ◽

Rumeng Li ◽

...

Keyword(s):

Land Surface ◽

Vegetation Index ◽

Normalized Difference Vegetation Index ◽

Growth Cycle ◽

Google Earth ◽

Henan Province ◽

Cropping Intensity ◽

Water Index ◽

Google Earth Engine ◽

Sentinel 2

Cropping intensity is a key indicator for evaluating grain production and intensive use of cropland. Timely and accurately monitoring of cropping intensity is of great significance for ensuring national food security and improving the level of national land management. In this study, we used all Sentinel-2 images on the Google Earth Engine cloud platform, and constructed an improved peak point detection method to extract the cropping intensity of a heterogeneous planting area combined with crop phenology. The crop growth cycle profiles were extracted from the multi-temporal normalized difference vegetation index (NDVI) and land surface water index (LSWI) datasets. Results show that by 2020, the area of single cropping, double cropping, and triple cropping in the Henan Province are 52,236.9 km2, 74,334.1 km2, and 1927.1 km2, respectively; the corresponding producer accuracies are 86.12%, 93.72%, and 91.41%, respectively; the corresponding user accuracies are 88.99%, 92.29%, and 71.26%, respectively. The overall accuracy is 90.95%, and the Kappa coefficient is 0.81. Using the sown area in the statistical yearbook data of cities in the Henan Province to verify the extraction results of this paper, the R2 is 0.9717, and the root mean square error is 1715.9 km2. This study shows that using all the Sentinel-2 data, the phenology algorithm, and cloud computing technology has great potential in producing a high spatio-temporal resolution dataset for crop remote sensing monitoring and agricultural policymaking in complex planting areas.

Download Full-text