The utility of unmanned aerial vehicles (UAV) imagery in retrieving phenotypic data to support plant breeding research has been a topic of increasing interest in recent years. The advantages of image-based phenotyping are related to the high spatial and temporal resolution of the retrieved data and the non-destructive and rapid method of data acquisition. This study trains parametric and nonparametric regression models to retrieve leaf area index (LAI), fraction of absorbed photosynthetically active radiation (fAPAR), fractional vegetation cover (fCover), leaf chlorophyll content (LCC), canopy chlorophyll content (CCC), and grain yield (GY) of winter durum wheat breeding experiment from four-bands UAV images. A ground dataset, collected during two field campaigns and complemented with data from a previous study, is used for model development. The dataset is split at random into two parts, one for training and one for testing the models. The tested parametric models use the vegetation index formula and parametric functions. The tested nonparametric models are partial least square regression (PLSR), random forest regression (RFR), support vector regression (SVR), kernel ridge regression (KRR), and Gaussian processes regression (GPR). The retrieved biophysical variables along with traditional phenotypic traits (plant height, yield, and tillering) are analysed for detection of genetic diversity, proximity, and similarity in the studied genotypes. Analysis of variance (ANOVA), Duncan’s multiple range test, correlation analysis, and principal component analysis (PCA) are performed with the phenotypic traits. The parametric and nonparametric models show close results for GY retrieval, with parametric models indicating slightly higher accuracy (R² = 0.49; RMSE = 0.58 kg/plot; rRMSE = 6.1%). However, the nonparametric model, GPR, computes per pixel uncertainty estimation, making it more appealing for operational use. Furthermore, our results demonstrate that grain filling was better than flowering phenological stage to predict GY. The nonparametric models show better results for biophysical variables retrieval, with GPR presenting the highest prediction performance. Nonetheless, robust models are found only for LAI (R² = 0.48; RMSE = 0.64; rRMSE = 13.5%) and LCC (R² = 0.49; RMSE = 31.57 mg m⁻²; rRMSE = 6.4%) and therefore these are the only remotely sensed phenotypic traits included in the statistical analysis for preliminary assessment of wheat productivity. The results from ANOVA and PCA illustrate that the retrieved remotely sensed phenotypic traits are a valuable addition to the traditional phenotypic traits for plant breeding studies. We believe that these preliminary results could speed up crop improvement programs; however, stronger interdisciplinary research is still needed, as well as uncertainty estimation of the remotely sensed phenotypic traits.

The monitoring of crops is of vital importance for food and environmental security in a global and European context. The main goal of this study was to assess the crop mapping performance provided by the 100 m spatial resolution of PROBA-V compared to coarser resolution data (e.g., PROBA-V at 300 m) for a 2250 km² test site in Bulgaria. The focus was on winter and summer crop mapping with three to five classes. For classification, single- and multi-date spectral data were used as well as NDVI time series. Our results demonstrate that crop identification using 100 m PROBA-V data performed significantly better in all experiments compared to the PROBA-V 300 m data. PROBA-V multispectral imagery, acquired in spring (March) was the most appropriate for winter crop identification, while satellite data acquired in summer (July) was superior for summer crop identification. The classification accuracy from PROBA-V 100 m compared to PROBA-V 300 m was improved by 5.8% to 14.8% depending on crop type. Stacked multi-date satellite images with three to four images gave overall classification accuracies of 74%–77% (PROBA-V 100 m data) and 66%–70% (PROBA-V 300 m data) with four classes (wheat, rapeseed, maize, and sunflower). This demonstrates that three to four image acquisitions, well distributed over the growing season, capture most of the spectral and temporal variability in our test site. Regarding the PROBA-V NDVI time series, useful results were only obtained if crops were grouped into two broader crop type classes (summer and winter crops). Mapping accuracies decreased significantly when mapping more classes. Again, a positive impact of the increased spatial resolution was noted. Together, the findings demonstrate the positive effect of the 100 m resolution PROBA-V data compared to the 300 m for crop mapping. This has important implications for future data provision and strengthens the arguments for a second generation of this mission originally designed solely as a “gap-filler mission”.

Land cover is one of the key terrestrial variables used for monitoring and as input for modelling in support of achieving the United Nations Strategical Development Goals. Global and Continental Land Cover Products (GCLCs) aim to provide the required harmonized information background across areas; thus, they are not being limited by national or other administrative nomenclature boundaries and their production approaches. Moreover, their increased spatial resolution, and consequently their local relevance, is of high importance for users at a local scale. During the last decade, several GCLCs were developed, including the Global Historical Land-Cover Change Land-Use Conversions (GLC), the Globeland-30 (GLOB), Corine-2012 (CLC) and GMES/ Copernicus Initial Operation High Resolution Layers (GIOS). Accuracy assessment is of high importance for product credibility towards incorporation into decision chains and implementation procedures, especially at local scales. The present study builds on the collaboration of scientists participating in the Global Observations of Forest Cover—Global Observations of Land Cover Dynamics (GOFC-GOLD), South Central and Eastern European Regional Information Network (SCERIN). The main objective is to quantitatively evaluate the accuracy of commonly used GCLCs at selected representative study areas in the SCERIN geographic area, which is characterized by extreme diversity of landscapes and environmental conditions, heavily affected by anthropogenic impacts with similar major socio-economic drivers. The employed validation strategy for evaluating and comparing the different products is detailed, representative results for the selected areas from nine SCERIN countries are presented, the specific regional differences are identified and their underlying causes are discussed. In general, the four GCLCs products achieved relatively high overall accuracy rates: 74–98% for GLC (mean: 93.8%), 79–92% for GLOB (mean: 90.6%), 74–91% for CLC (mean: 89%) and 72–98% for GIOS (mean: 91.6%), for all selected areas. In most cases, the CLC product has the lower scores, while the GLC has the highest, closely followed by GIOS and GLOB. The study revealed overall high credibility and validity of the GCLCs products at local scale, a result, which shows expected benefit even for local/regional applications. Identified class dependent specificities in different landscape types can guide the local users for their reasonable usage in local studies. Valuable information is generated for advancing the goals of the international GOFC-GOLD program and aligns well with the agenda of the NASA Land-Cover/Land-Use Change Program to improve the quality and consistency of space-derived higher-level products.

This paper presents the results of a sub-pixel classification of crop types in Bulgaria from PROBA-V 100 m normalized difference vegetation index (NDVI) time series. Two sub-pixel classification methods, artificial neural network (ANN) and support vector regression (SVR) were used where the output was a set of area fraction images (AFIs) at 100 m resolution with pixels containing estimated area fractions of each class. High-resolution maps of two test sites derived from Sentinel-2 classifications were used to obtain training data for the sub-pixel classifications. The estimated area fractions have a good correspondence with the true area fractions when aggregated to regions of 10 × 10 km², especially when the SVR method was used. For the five dominant classes in the test sites the R² obtained after the aggregation was 86% (winter cereals), 81% (sunflower), 92% (broad-leaved forest), 89% (maize), and 67% (grasslands) when the SVR method was used.

This study tested the potential of parametric and nonparametric regression modeling utilizing multispectral data from two different unoccupied aerial vehicles (UAVs) as a tool for the prediction of and indirect selection of grain yield (GY) in barley breeding experiments. The coefficient of determination (R²) of the nonparametric models for GY prediction ranged between 0.33 and 0.61 depending on the UAV and flight date, where the highest value was achieved with the DJI Phantom 4 Multispectral (P4M) image from 26 May (milk ripening). The parametric models performed worse than the nonparametric ones for GY prediction. Independent of the retrieval method and UAV, GY retrieval was more accurate in milk ripening than dough ripening. The leaf area index (LAI), fraction of absorbed photosynthetically active radiation (fAPAR), fraction vegetation cover (fCover), and leaf chlorophyll content (LCC) were modeled at milk ripening using nonparametric models with the P4M images. A significant effect of the genotype was found for the estimated biophysical variables, which was referred to as remotely sensed phenotypic traits (RSPTs). Measured GY heritability was lower, with a few exceptions, compared to the RSPTs, indicating that GY was more environmentally influenced than the RSPTs. The moderate to strong genetic correlation of the RSPTs to GY in the present study indicated their potential utility as an indirect selection approach to identify high-yield genotypes of winter barley.

In the present work we applied a recently developed procedure for multidimensional data clustering to processing of spectral satellite images. The core of our approach lays in projection of multidimensional image to a two dimensional one. The main aim is to discover points with similar characteristics. This was done by clustering of the resulting image. The processing technique exploits equilibrium states of a kind of recurrent neural network - Echo state network (ESN) - that are obtained after intrinsic plasticity (IP) tuning of the ESN using multidimensional data as inputs. The proposed in our previous work automated procedure for multidimensional data clustering is further refined and tested on the satellite image data. The obtained number and position of clusters of a multi-spectral image of a mountain region in Bulgaria is compared with the classification of the region landscape given by the Ministry of Regional Development and Public Works.