Harmonized Landsat 8 and Sentinel-2 Time Series Data to Detect Irrigated Areas: An Application in Southern Italy

Lack of accurate and up-to-date data associated with irrigated areas and related irrigation amounts is hampering the full implementation and compliance of the Water Framework Directive (WFD). In this paper, we describe the framework that we developed and implemented within the DIANA project to map the actual extent of irrigated areas in the Campania region (Southern Italy) during the 2018 irrigation season. For this purpose, we considered 202 images from the Harmonized Landsat Sentinel-2 (HLS) products (57 images from Landsat 8 and 145 images from Sentinel-2). Such data were preprocessed in order to extract a multitemporal Normalized Difference Vegetation Index (NDVI) map, which was then smoothed through a gap-filling algorithm. We further integrated data coming from high-resolution (4 km) global satellite precipitation Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks (PERSIANN)-Cloud Classification System (CCS) products. We collected an extensive ground truth in the field represented by 2992 data points coming from three main thematic classes: bare soil and rainfed (class 0), herbaceous (class 1), and tree crop (class 2). This information was exploited to generate irrigated area maps by adopting a machine learning classification approach. We compared six different types of classifiers through a cross-validation approach and found that, in general, random forests, support vector machines, and boosted decision trees exhibited the best performances in terms of classification accuracy and robustness to different tested scenarios. We found an overall accuracy close to 90% in discriminating among the three thematic classes, which highlighted promising capabilities in the detection of irrigated areas from HLS products

Mapping Irrigated and Rainfed Wheat Areas Using Multi-Temporal Satellite Data

Irrigation is crucial to agriculture in arid and semi-arid areas and significantly contributes to crop development, food diversity and the sustainability of agro-ecosystems. For a specific crop, the separation of its irrigated and rainfed areas is difficult, because their phenology is similar and therefore less distinguishable, especially when there are phenology shifts due to various factors, such as elevation and latitude. In this study, we present a simple, but robust method to map irrigated and rainfed wheat areas in a semi-arid region of China. We used the Normalized Difference Vegetation Index (NDVI) at a 30 × 30 m spatial resolution derived from the Chinese HJ-1A/B (HuanJing(HJ) means environment in Chinese) satellite to create a time series spanning the whole growth period of wheat from September 2010 to July 2011. The maximum NDVI and time-integrated NDVI (TIN) that usually exhibit significant differences between irrigated and rainfed wheat were selected to establish a classification model using a support vector machine (SVM) algorithm. The overall accuracy of the Google-Earth testing samples was 96.0%, indicating that the classification results are accurate. The estimated irrigated-to-rainfed ratio was 4.4:5.6, close to the estimates provided by the agricultural sector in Shanxi Province. Our results illustrate that the SVM classification model can effectively avoid empirical thresholds in supervised classification and realistically capture the magnitude and spatial patterns of rainfed and irrigated wheat areas. The approach in this study can be applied to map irrigated/rainfed areas in other regions when field observational data are available

Characterizing and mapping cropping patterns in a complex agro-ecosystem: An iterative participatory mapping procedure using machine learning algorithms and MODIS vegetation indices

Accurate and up-to-date spatial agricultural information is essential for applications including agro-environmental assessment, crop management, and appropriate targeting of agricultural technologies. There is growing research interest in spatial analysis of agricultural ecosystems applying satellite remote sensing technologies. However, usability of information generated from many of remotely sensed data is often constrained by accuracy problems. This is of particular concern in mapping complex agro-ecosystems in countries where small farm holdings are dominated by diverse crop types. This study is a contribution to the ongoing efforts towards overcoming accuracy challenges faced in remote sensing of agricultural ecosystems. We applied time-series analysis of vegetation indices (Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI)) derived from the Moderate Resolution Imaging Spectrometer (MODIS) sensor to detect seasonal patterns of irrigated and rainfed cropping patterns in five townships in the Central Dry Zone of Myanmar, which is an important agricultural region of the country has been poorly mapped with respect to cropping practices. To improve mapping accuracy and map legend completeness, we implemented a combination of (i) an iterative participatory approach to field data collection and classification, (ii) the identification of appropriate size and types of predictor variables (VIs), and (iii) evaluation of the suitability of three Machine Learning algorithms: Support Vector Machine (SVM), Random Forest (RF), and C5.0 algorithms under varying training sample sizes. Through these procedures, we were able to progressively improve accuracy and achieve maximum overall accuracy of 95% When a small sized training dataset was used, accuracy achieved by RF was significantly higher compared to SVM and C5.0 (P < 0.01), but as sample size increased, accuracy differences among the three machine learning algorithms diminished. Accuracy achieved by use of NDVI was consistently better than that of EVI (P < 0.01). The maximum overall accuracy was achieved using RF and 8-days NDVI composites for three years of remote sensing data. In conclusion, our findings highlight the important role of participatory classification, especially in areas where cropping systems are highly diverse and differ over space and time. We also show that the choice of classifiers and size of predictor variables are essential and complementary to the participatory mapping approach in achieving desired accuracy of cropping pattern mapping in areas where other sources of spatial information are scarce

Use and Improvement of Remote Sensing and Geospatial Technologies in Support of Crop Area and Yield Estimations in the West African Sahel

In arid and semi-arid West Africa, agricultural production and regional food security depend largely on small-scale subsistence farming and rainfed crops, both of which are vulnerable to climate variability and drought. Efforts made to improve crop monitoring and our ability to estimate crop production (areas planted and yield estimations by crop type) in the major agricultural zones of the region are critical paths for minimizing climate risks and to support food security planning. The main objective of this dissertation research was to contribute to these efforts using remote sensing technologies. In this regard, the first analysis documented the low reliability of existing land cover products for cropland area estimation (Chapter 2). Then two satellite remote sensing-based datasets were developed that 1) accurately map cropland areas in the five countries of Sahelian West Africa (Senegal, Mauritania, Mali, Burkina Faso and Niger; Chapter 3), and 2) focus on the country of Mali to identify the location and prevalence of the major subsistence crops (millet, sorghum, maize and non-irrigated rice; Chapter 4). The regional cropland area product is distributed as the West African Sahel Cropland area at 30 m (WASC30). The development of the new dataset involved high density training data (380,000 samples) developed by USGS in collaboration with CILSS for training about 200 locally optimized random forest (RF) classifiers using Landsat 8 surface reflectances and vegetation indices and the Google Earth Engine platform. WASC30 greatly improves earlier estimates through inclusion of cropland information for both rainfed and irrigated areas mapped with a class-specific accuracy of 79% across the West Africa Sahel. Used as a mask in crop monitoring systems, the new cropland area data could bring critical insights by reducing uncertainties in xv identification of croplands as crop growth condition metrics are extracted. WASC30 allowed us to derive detailed statistics on cultivated areas in the Sahel, at country and agroclimatic scales. Intensive agricultural zones were highlighted as well. The second dataset, mapping crop types for the country of Mali, is meant to separate signals of different crop types for improved crop yield estimation. The crop type map was used to derive detailed agricultural statistics (e.g. acreage by crop types, spatial distribution) at finer administrative scales than has previously been possible. The crop fraction information by crop type extracted from the map, gives additional details on farmers preferences by regions, and the natural adaptability of different crop types. The final analysis of this dissertation explores the use of ensemble machine learning techniques to predict maize yield in Mali (Chapter 5). Climate data (precipitation and temperature), and vegetation indices (Normalized Difference Vegetation Index, NDVI, the Enhanced Vegetation Index, EVI, and the Normalized Difference Water Index, NDWI) are used as predictors, while actual yields collected in 2017 by the Malian Ministry of Agriculture are the reference data. Random forest presented better predictive performance as compared to boosted regression trees (BRT). Results showed that climate variables have more predictive power for maize yield compared to vegetation indices. Among vegetation indices, the NDWI appeared to be the most influential predictor, maybe because of water requirement of maize and the sensitivity of this index to water in semi-arid regions. Tested with two different independent datasets, one constituted by 20% of the reference information, and another including observed yields for year 2018 (a one-year-left analysis), maize yield predictions were promising for year 2017 (RMSE = 362 kg/ha), but showed higher error for 2018 (RMSE = 707 kg/ha). That is, the fitted model may not capture accurately year to year variabilities in predicted maize yield. In this analysis, predictions were limited to field samples (~600 fields) across the country of Mali. It would be valuable in the future to predict maize yield for each pixel of the new developed crop type map. That will lead to a detailed spatial analysis of maize yield, allowing identification of low yielding regions for targeted interventions which could improve food security. Keywords: Agricultural identification of croplands as crop growth condition metrics are extracted. WASC30 allowed us to derive detailed statistics on cultivated areas in the Sahel, at country and agroclimatic scales. Intensive agricultural zones were highlighted as well. The second dataset, mapping crop types for the country of Mali, is meant to separate signals of different crop types for improved crop yield estimation. The crop type map was used to derive detailed agricultural statistics (e.g. acreage by crop types, spatial distribution) at finer administrative scales than has previously been possible. The crop fraction information by crop type extracted from the map, gives additional details on farmers preferences by regions, and the natural adaptability of different crop types. The final analysis of this dissertation explores the use of ensemble machine learning techniques to predict maize yield in Mali (Chapter 5). Climate data (precipitation and temperature), and vegetation indices (Normalized Difference Vegetation Index, NDVI, the Enhanced Vegetation Index, EVI, and the Normalized Difference Water Index, NDWI) are used as predictors, while actual yields collected in 2017 by the Malian Ministry of Agriculture are the reference data. Random forest presented better predictive performance as compared to boosted regression trees (BRT). Results showed that climate variables have more predictive power for maize yield compared to vegetation indices. Among vegetation indices, the NDWI appeared to be the most influential predictor, maybe because of water requirement of maize and the sensitivity of this index to water in semi-arid regions. Tested with two different independent datasets, one constituted by 20% of the reference information, and another including observed yields for year 2018 (a one-year-left analysis), maize yield predictions were promising for year 2017 (RMSE = 362 kg/ha), but showed higher error for 2018 (RMSE = 707 kg/ha). That is, the fitted model may not capture accurately year to year variabilities in predicted maize yield. In this analysis, predictions were limited to field samples (~600 fields) across the country of Mali. It would be valuable in the future to predict maize yield for each pixel of the new developed crop type map. That will lead to a detailed spatial analysis of maize yield, allowing identification of low yielding regions for targeted interventions which could improve food security

Comparative geospatial approach for agricultural crops identification in interfluvial plain - A case study of Sahiwal district, Pakistan

Agricultural crop cover identification is a major issue and time-consuming effort to verify the crop type through surveys of the individual field or using prehistoric methods. To establish the scenario of crop identification, the stage of crop provides diverse spatial information about the variety of crops due to its spectral changes. The main aim of this study was to the identify the crop types and their behavior using remote sensing and geographical information system-based approach. Moreover, two main methods were applied to the Sentinel-2 satellite data in which one is random forest based supervised classification and another was Normalize Difference Vegetation Index (NDVI) density estimation method through the google earth engine to procure the data in time-efficient way. This study also established the comparison between classified and vegetation index based seasonal compositional datasets for wheat, cotton, maize, and fodder crops. Study discussed the best fit technique for crops identification in the light of observed methods. Furthermore, the vegetation index ranges by the zonal statistics of the field samples were established according to crop precision. Results showed that -22.94, -43.72, 20.61, and 32.49 % dissimilarities existed in wheat, fodder, cotton, and maize results respectively, after comparison of both techniques. Although, the accuracy assessment was performed on the classified dataset for validation of results by confusion matrix accuracy assessment process using field sample data. Moreover, the vegetation index was used to evaluate crop land surface temperature to estimate the crop growth stage valuation that revealed noticeably enthralling outcomes. The results determined that the classified accuracies of wheat, cotton, maize and fodder were 84, 80, 81 and 71 % respectively. This study also revealed that the random forest classifier has used more features and information potentially during the classifier trainings but vegetation index just implies the limited number of features such as crop growing status

IDENTIFICATION OF AGRICULTURAL LAND USE IN CALIFORNIA THROUGH REMOTE SENSING

Ground truthing actual crop types in an area can be expensive and time-consuming. The California Department of Water Resources attempts to ground truth land use in each county in California every five years. However, this is limited by budgetary constraints and often results in infrequent (more than every ten years) surveying of many counties. An accurate accounting of crops growing in a region is important for a variety of purposes including farm production estimates, groundwater and surface water modeling, evapotranspiration estimation, water planning, research applications, etc. Agricultural land use is continually changing due to development and environmental factors. Currently, USDA NASS provides georeferenced land use maps of regions throughout the U.S. While these are beneficial, the accuracy is not very high for California due to the wide variety of crops grown throughout the state. California has an increasingly complex agricultural system which includes multi-crops changing on an annual and even semiannual basis, long growing seasons, and complex and flexible irrigation schedules. Remotely sensed data from available satellites are used to more accurately classify crop types within the Madera and Merced Counties of California’s Central Valley. An initial classification approach utilizing a simplified decision tree for a data subset of the area considered is presented. In order to accommodate the larger dataset at hand, a computer based approach is applied using the Nearest Neighbor classification algorithm in the computer program eCognition. Iterative analyses were performed to consider a range of scenarios with varying spectral inputs. The results show the methods presented can be beneficial in discriminating 24 of the major crop types from multi-temporal spectral data

Landcover and crop type classification with intra-annual times series of sentinel-2 and machine learning at central Portugal

Dissertation submitted in partial fulfilment of the requirements for the degree of Master of Science in Geospatial TechnologiesLand cover and crop type mapping have benefited from a daily revisiting period of sensors such as MODIS, SPOT-VGT, NOAA-AVHRR that contains long time-series archive. However, they have low accuracy in an Area of Interest (ROI) due to their coarse spatial resolution (i.e., pixel size > 250m). The Copernicus Sentinel-2 mission from the European Spatial Agency (ESA) provides free data access for Sentinel 2-A(S2a) and B (S2b). This satellite constellation guarantees a high temporal (5-day revisit cycle) and high spatial resolution (10m), allowing frequent updates on land cover products through supervised classification. Nevertheless, this requires training samples that are traditionally collected manually via fieldwork or image interpretation. This thesis aims to implement an automatic workflow to classify land cover and crop types at 10m resolution in central Portugal using existing databases, intra-annual time series of S2a and S2b, and Random Forest, a supervised machine learning algorithm. The agricultural classes such as temporary and permanent crops as well as agricultural grasslands were extracted from the Portuguese Land Parcel Identification System (LPIS) of the Instituto de Financiamento da Agricultura e Pescas (IFAP); land cover classes like urban, forest and water were trained from the Carta de Ocupação do Solo (COS) that is the national Land Use and Land Cover (LULC) map of Portugal; and lastly, the burned areas are identified from the corresponding national map of the Instituto da Conservação da Natureza e das Florestas (ICNF). Also, a set of preprocessing steps were defined based on the implementation of ancillary data allowing to avoid the inclusion of mislabeled pixels to the classifier. Mislabeling of pixels can occur due to errors in digitalization, generalization, and differences in the Minimum Mapping Unit (MMU) between datasets. An inner buffer was applied to all datasets to reduce border overlap among classes; the mask from the ICNF was applied to remove burned areas, and NDVI rule based on Landsat 8 allowed to erase recent clear-cuts in the forest. Also, the Copernicus High-Resolution Layers (HRL) datasets from 2015 (latest available), namely Dominant Leaf Type (DLT) and Tree Cover Density (TCD) are used to distinguish between forest with more than 60% coverage (coniferous and broadleaf) such as Holm Oak and Stone Pine and between 10 and 60% (coniferous) for instance Open Maritime Pine. Next, temporal gap-filled monthly composites were created for the agricultural period in Portugal, ranging from October 2017 till September 2018. The composites provided data free of missing values in opposition to single date acquisition images. Finally, a pixel-based approach classification was carried out in the “Tejo and Sado” region of Portugal using Random Forest (RF). The resulting map achieves a 76% overall accuracy for 31 classes (17 land cover and 14 crop types). The RF algorithm captured the most relevant features for the classification from the cloud-free composites, mainly during the spring and summer and in the bands on the Red Edge, NIR and SWIR. Overall, the classification was more successful on the irrigated temporary crops whereas the grasslands presented the most complexity to classify as they were confused with other rainfed crops and burned areas

Crop type identification and spatial mapping using Sentinel-2 satellite data with focus on field-level information

Accurate monitoring of croplands helps in making decisions (for insurance claims, crop management and contingency plans) at the macro-level, especially in drylands where variability in cropping is very high owing to erratic weather conditions. Dryland cereals and grain legumes are key to ensuring the food and nutritional security of a large number of vulnerable populations living in the drylands. Reliable information on area cultivated to such crops forms part of the national accounting of food production and supply in many Asian countries, many of which are employing remote sensing tools to improve the accuracy of assessments of cultivated areas. This paper assesses the capabilities and limitations of mapping cultivated areas in the Rabi (winter) season and corresponding cropping patterns in three districts characterized by small-plot agriculture. The study used Sentinel-2 Normalized Difference Vegetation Index (NDVI) 15-day time-series at 10m resolution by employing a Spectral Matching Technique (SMT) approach. The use of SMT is based on the well-studied relationship between temporal NDVI signatures and crop phenology. The rabi season in India, dominated by non-rainy days, is best suited for the application of this method, as persistent cloud cover will hamper the availability of images necessary to generate clearly differentiating temporal signatures. Our study showed that the temporal signatures of wheat, chickpea and mustard are easily distinguishable, enabling an overall accuracy of 84%, with wheat and mustard achieving 86% and 94% accuracies, respectively. The most significant misclassifications were in irrigated areas for mustard and wheat, in small-plot mustard fields covered by trees and in fragmented chickpea areas. A comparison of district-wise national crop statistics and those obtained from this study revealed a correlation of 96%

Mapping Small-Scale Irrigation Areas Using Expert Decision Rules and the Random Forest Classifier in Northern Ethiopia

The mapping of small-scale irrigation areas is essential for food security and water resource management studies. The identification of small-scale irrigation areas is a challenge, but it can be overcome using expert knowledge and satellite-derived high-spatial-resolution multispectral information in conjunction with monthly normalized difference vegetation index (NDVI) time series, and additional terrain information. This paper presents a novel approach to characterize small-scale irrigation schemes that combine expert knowledge, multi-temporal NDVI time series, multispectral high-resolution satellite images, and the random forest classifier in the Zamra catchment, North Ethiopia. A fundamental element of the approach is mapping small-scale irrigation areas using expert decision rules to incorporate the available water resources. We apply expert decision rules to monthly NDVI composites from September 2020 to August 2021 along with the digital elevation model (DEM) data on the slope, drainage order, and distance maps to derive the sample set. The samples werebased on the thresholds obtained by expert knowledge from field surveys. These data, along with the four spectral bands of a cloud-free Planet satellite image composite, 12 NDVI monthly composites, slope, drainage order, and distance map were used as input into a random forest classifier which was trained to classify pixels as either irrigated or non-irrigated. The results show that the analysis allows the mapping of small-scale irrigation areas with high accuracy. The classification accuracy for identifying irrigated areas showed a user accuracy ranging from 81% to 87%, along with a producer accuracy ranging from 64% to 79%. Furthermore, the classification accuracy and the kappa coefficient for the classified irrigation schemes were 80% and 0.70, respectively. As a result, these findings highlight a substantial level of agreement between the classification results and the reference data. The use of different expert knowledge-based decision rules, as a method, can be applied to extract small-scale and larger irrigation areas with similar agro-ecological characteristics.<br/