Mapping of peanut crops in Queensland, Australia using time-series PROBA-V 100-m normalized difference vegetation index imagery

Mapping of peanut crops is essential in supporting peanut production, yield prediction, and commodity forecasting. While ground-based surveys can be used over small areas, the development of remote-sensing technologies could provide rapid and inexpensive crop area estimates with high accuracy over large regions. Some of these recent earth observation satellite systems, such as the Project for On-Board Autonomy Vegetation (PROBA-V), have the advantage of increased spatial and temporal resolution. With a study area located in the South Burnett region, Queensland, Australia, the primary aim of this study was to assess the ability of time-series PROBA-V 100-m normalized difference vegetation index (NDVI) for peanut crop mapping. Two datasets, i.e., PROBA-V NDVI time-series imagery and the corresponding phenological parameters generated from TIMESAT data analysis technique, were classified using maximum likelihood classification, spectral angle mapper, and minimum distance classification algorithms. The results show that among all methods used, the application of MLC in PROBA-V NDVI time series produced very good overall accuracy, i.e., 92.75%, with producer and user accuracy of each class ≥78.79  %  . For all algorithms tested, the mapping of peanut cropping areas produced satisfactory classification results, i.e., 75.95% to 100%. Our study confirmed that the use of finer resolution 100 m of PROBA-V imagery (i.e., relative to MODIS 250-m data) has contributed to the success of mapping peanut and other crops in the study area

Mapping of multitemporal rice (Oryza sativa L.) growth stages using remote sensing with multi-sensor and machine learning : a thesis dissertation presented in partial fulfilment of the requirements for the degree of Doctor of Philosophy in Earth Science at Massey University, Manawatū, New Zealand

Figure 2.1 is adapted and re-used under a Creative Commons Attribution 4.0 International (CC BY 4.0) license.Rice (Oryza Sativa) plays a pivotal role in food security for Asian countries, especially in Indonesia. Due to the increasing pressure of environmental changes, such as land use and climate, rice cultivation areas need to be monitored regularly and spatially to ensure sustainable rice production. Moreover, timely information of rice growth stages (RGS) can lead to more efficient of inputs distribution from water, seed, fertilizer, and pesticide. One of the efficient solutions for regularly mapping the rice crop is using Earth observation satellites. Moreover, the increasing availability of open access satellite images such as Landsat-8, Sentinel-1, and Sentinel-2 provides ample opportunities to map continuous and high-resolution rice growth stages with greater accuracy. The majority of the literature has focused on mapping rice area, cropping patterns and relied mainly on the phenology of vegetation. However, the mapping process of RGS was difficult to assess the accuracy, time-consuming, and depended on only one sensor. In this work, we discuss the use of machine learning algorithms (MLA) for mapping paddy RGS with multiple remote sensing data in near-real-time. The study area was Java Island, which is the primary rice producer in Indonesia. This study has investigated: (1) the mapping of RGS using Landsat-8 imagery and different MLAs, and their rigorous performance was evaluated by conducting a multitemporal analysis; (2) the temporal consistency of predicting RGS using Sentinel-2, MOD13Q1, and Sentinel-1 data; (3) evaluating the correlation of local statistics data and paddy RGS using Sentinel-2, PROBA-V, and Sentinel-1 with MLAs. The ground truth datasets were collected from multi-year web camera data (2014-2016) and three months of the field campaign in different regions of Java (2018). The study considered the RGS in the analysis to be vegetative, reproductive, ripening, bare land, and flooding, and MLAs such as support vector machines (SVMs), random forest (RF), and artificial neural network (ANN) were used. The temporal consistency matrix was used to compare the classification maps within three sensor datasets (Landsat-8 OLI, Sentinel-2, and Sentinel-2, MOD13Q1, Sentinel-1) and in four periods (5, 10, 15, 16 days). Moreover, the result of the RGS map was also compared with monthly data from local statistics within each sub-district using cross-correlation analysis. The result from the analysis shows that SVM with a radial base function outperformed the RF and ANN and proved to be a robust method for small-size datasets (< 1,000 points). Compared to Sentinel-2, Landsat-8 OLI gives less accuracy due to the lack of a red-edge band and larger pixel size (30 x 30 m). Integration of Sentinel-2, MOD13Q1, and Sentinel-1 improved the classification performance and increased the temporal availability of cloud-free maps. The integration of PROBA-V and Sentinel-1 improved the classification accuracy from the Landsat-8 result, consistent with the monthly rice planting area statistics at the sub-district level. The western area of Java has the highest accuracy and consistency since the cropping pattern only relied on rice cultivation. In contrast, less accuracy was noticed in the eastern area because of upland rice cultivation due to limited irrigation facilities and mixed cropping. In addition, the cultivation of shallots to the north of Nganjuk Regency interferes with the model predictions because the cultivation of shallots resembles the vegetative phase due to the water banks. One future research idea is the auto-detection of the cropping index in the complex landscape to be able to use it for mapping RGS on a global scale. Detection of the rice area and RGS using Google Earth Engine (GEE) can be an action plan to disseminate the information quickly on a planetary scale. Our results show that the multitemporal Sentinel-1 combined with RF can detect rice areas with high accuracy (>91%). Similarly, accurate RGS maps can be detected by integrating multiple remote sensing (Sentinel-2, Landsat-8 OLI, and MOD13Q1) data with acceptable accuracy (76.4%), with high temporal frequency and lower cloud interference (every 16 days). Overall, this study shows that remote sensing combined with the machine learning methodology can deliver information on RGS in a timely fashion, which is easy to scale up and consistent both in time and space and matches the local statistics. This thesis is also in line with the existing rice monitoring projects such as Crop Monitor, Crop Watch, AMIS, and Sen4Agri to support disseminating information over a large area. To sum up, the proposed workflow and detailed map provide a more accurate method and information in near real-time for stakeholders, such as governmental agencies against the existing mapping method. This method can be introduced to provide accurate information to rice farmers promptly with sufficient inputs such as irrigation, seeds, and fertilisers for ensuring national food security from the shifting planting time due to climate change

Copernicus Global Land Cover Layers—Collection 2

In May 2019, Collection 2 of the Copernicus Global Land Cover layers was released. Next to a global discrete land cover map at 100 m resolution, a set of cover fraction layers is provided depicting the percentual cover of the main land cover types in a pixel. This additional continuous classification scheme represents areas of heterogeneous land cover better than the standard discrete classification scheme. Overall, 20 layers are provided which allow customization of land cover maps to specific user needs or applications (e.g., forest monitoring, crop monitoring, biodiversity and conservation, climate modeling, etc.). However, Collection 2 was not just a global up-scaling, but also includes major improvements in the map quality, reaching around 80% or more overall accuracy. The processing system went into operational status allowing annual updates on a global scale with an additional implemented training and validation data collection system. In this paper, we provide an overview of the major changes in the production of the land cover maps, that have led to this increased accuracy, including aligning with the Sentinel 2 satellite system in the grid and coordinate system, improving the metric extraction, adding better auxiliary data, improving the biome delineations, as well as enhancing the expert rules. An independent validation exercise confirmed the improved classification results. In addition to the methodological improvements, this paper also provides an overview of where the different resources can be found, including access channels to the product layer as well as the detailed peer-review product documentation

Potential and Limitations of Open Satellite Data for Flood Mapping

Satellite remote sensing is a powerful tool to map flooded areas. In recent years, the availability of free satellite data significantly increased in terms of type and frequency, allowing the production of flood maps at low cost around the world. In this work, we propose a semi-automatic method for flood mapping, based only on free satellite images and open-source software. The proposed methods are suitable to be applied by the community involved in flood hazard management, not necessarily experts in remote sensing processing. As case studies, we selected three flood events that recently occurred in Spain and Italy. Multispectral satellite data acquired by MODIS, Proba-V, Landsat, and Sentinel-2 and synthetic aperture radar (SAR) data collected by Sentinel-1 were used to detect flooded areas using different methodologies (e.g., Modified Normalized Difference Water Index, SAR backscattering variation, and supervised classification). Then, we improved and manually refined the automatic mapping using free ancillary data such as the digital elevation model-based water depth model and available ground truth data. We calculated flood detection performance (flood ratio) for the different datasets by comparing with flood maps made by official river authorities. The results show that it is necessary to consider different factors when selecting the best satellite data. Among these factors, the time of the satellite pass with respect to the flood peak is the most important. With co-flood multispectral images, more than 90% of the flooded area was detected in the 2015 Ebro flood (Spain) case study. With post-flood multispectral data, the flood ratio showed values under 50% a few weeks after the 2016 flood in Po and Tanaro plains (Italy), but it remained useful to map the inundated pattern. The SAR could detect flooding only at the co-flood stage, and the flood ratio showed values below 5% only a few days after the 2016 Po River inundation. Another result of the research was the creation of geomorphology-based inundation maps that matched up to 95% with official flood maps

Revisiting the coupling between NDVI trends and cropland changes in the Sahel drylands:a case study in western Niger

The impact of human activities via land use/cover changes on NDVI trends is critical for an improved understanding of satellite-observed changes in vegetation productivity in drylands. The dominance of positive NDVI trends in the Sahel, the so-called re-greening, is sometimes interpreted as a combined effect of an increase in rainfall and cropland expansion or agricultural intensification. Yet, the impact of changes in land use has yet to be thoroughly tested and supported by empirical evidence. At present, no studies have considered the importance of the different seasonal NDVI signals of cropped and fallowed fields when interpreting NDVI trends, as both field types are commonly merged into a single ‘cropland’ class. We make use of the distinctly different phenology of cropped and fallowed fields and use seasonal NDVI curves to separate these two field types. A fuzzy classifier is applied to quantify cropped and fallowed areas in a case study region in the southern Sahel (Fakara, Niger) on a yearly basis between 2000 and 2014. We find that fallowed fields have a consistently higher NDVI than unmanured cropped fields and by using two seasonal NDVI metrics (the amplitude and the decreasing rate) derived from the MODIS time series, a clear separation between classes of fields is achieved (r = 0.77). The fuzzy classifier can compute the percentage of a pixel (250 m) under active cultivation, thereby alleviating the problem of small field sizes in the region. We find a predominant decrease in NDVI over the period of analysis associated with an increased area of cropped fields at the expense of fallowed fields. Our findings couple cropping abandonment (more frequent fallow years) with positive NDVI trends and an increase in the percentage of the cropped area (fallow period shortening) with negative trends. These findings profoundly impact our understanding of greening and browning trends in agrarian Sahelian drylands and in other drylands of developing countries characterized by limited use of fertilizers

Estimating the crop leaf area index using hyperspectral remote sensing

AbstractThe leaf area index (LAI) is an important vegetation parameter, which is used widely in many applications. Remote sensing techniques are known to be effective but inexpensive methods for estimating the LAI of crop canopies. During the last two decades, hyperspectral remote sensing has been employed increasingly for crop LAI estimation, which requires unique technical procedures compared with conventional multispectral data, such as denoising and dimension reduction. Thus, we provide a comprehensive and intensive overview of crop LAI estimation based on hyperspectral remote sensing techniques. First, we compare hyperspectral data and multispectral data by highlighting their potential and limitations in LAI estimation. Second, we categorize the approaches used for crop LAI estimation based on hyperspectral data into three types: approaches based on statistical models, physical models (i.e., canopy reflectance models), and hybrid inversions. We summarize and evaluate the theoretical basis and different methods employed by these approaches (e.g., the characteristic parameters of LAI, regression methods for constructing statistical predictive models, commonly applied physical models, and inversion strategies for physical models). Thus, numerous models and inversion strategies are organized in a clear conceptual framework. Moreover, we highlight the technical difficulties that may hinder crop LAI estimation, such as the “curse of dimensionality” and the ill-posed problem. Finally, we discuss the prospects for future research based on the previous studies described in this review

Soil map density and a nation's wealth and income

Little effort has been made to link soil mapping and soil data density to a nation’s welfare. Soil map density in 31 European countries and 44 low and middle income countries is linked to Gross Domestic Product (GDP) per capita and the number of soil scientists per country

Using Low Resolution Satellite Imagery for Yield Prediction and Yield Anomaly Detection

Low resolution satellite imagery has been extensively used for crop monitoring and yield forecasting for over 30 years and plays an important role in a growing number of operational systems. The combination of their high temporal frequency with their extended geographical coverage generally associated with low costs per area unit makes these images a convenient choice at both national and regional scales. Several qualitative and quantitative approaches can be clearly distinguished, going from the use of low resolution satellite imagery as the main predictor of final crop yield to complex crop growth models where remote sensing-derived indicators play different roles, depending on the nature of the model and on the availability of data measured on the ground. Vegetation performance anomaly detection with low resolution images continues to be a fundamental component of early warning and drought monitoring systems at the regional scale. For applications at more detailed scales, the limitations created by the mixed nature of low resolution pixels are being progressively reduced by the higher resolution offered by new sensors, while the continuity of existing systems remains crucial for ensuring the availability of long time series as needed by the majority of the yield prediction methods used today.JRC.H.4-Monitoring Agricultural Resource

Merging the Minnaert-k parameter with spectral unmixing to map forest heterogeneity with CHRIS/PROBA data

The Compact High Resolution Imaging Spectrometer (CHRIS) mounted onboard the Project for Onboard Autonomy (PROBA) spacecraft is capable of sampling reflected radiation at five viewing angles over the visible and near-infrared regions of the solar spectrum with high spatial resolution. We combined the spectral domain with the angular domain of CHRIS data in order to map the surface heterogeneity of an Alpine coniferous forest during winter. In the spectral domain, linear spectral unmixing of the nadir image resulted in a canopy cover map. In the angular domain, pixelwise inversion of the Rahman-Pinty-Verstraete (RPV) model at a single wavelength at the red edge (722 nm) yielded a map of the Minnaert-k parameter that provided information on surface heterogeneity at a subpixel scale. However, the interpretation of the Minnaert-k parameter is not always straightforward because fully vegetated targets typically produce the same type of reflectance anisotropy as non-vegetated targets. Merging both maps resulted in a forest cover heterogeneity map, which contains more detailed information on canopy heterogeneity at the CHRIS subpixel scale than is possible to realize from a single-source optical data set