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”.

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.

The aim of this study is to predict and map winter wheat yield in the Parvomay municipality, situated in the Upper Thracian Lowland of Bulgaria, utilizing satellite data from Sentinel-2. The main crops grown in the research area are winter wheat, rapeseed, sunflower, and maize. To distinguish winter wheat fields accurately, we evaluated classification methods such as Support Vector Machines (SVM) and Random Forest (RF). These methods were applied to satellite multispectral data acquired by the Sentinel-2 satellites during the growing season of 2020–2021. In accordance with their development cycles, temporal image composites were developed to identify suitable moments when each crop is most accurately distinguished from others. Ground truth data obtained from the integrated administration and control system (IACS) were used for training the classifiers and assessing the accuracy of the final maps. Winter wheat fields were masked using the crop mask created from the best-performing classification algorithm. Yields were predicted with regression models calibrated with in situ data collected in the Parvomay study area. Both SVM and RF algorithms performed well in classifying winter wheat fields, with SVM slightly outperforming RF. The produced crop maps enable the application of crop-specific yield models on a regional scale. The best predictor of yield was the green NDVI index (GNDVI) from the April monthly composite image.

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.

Barley is an ecologically adaptable crop widely used in agriculture and well suited for organic farming. Satellite imagery from Sentinel-2 can support crop monitoring and yield prediction, optimising production processes. This study compares two types of Sentinel-2 data—standard (S2) data with 10 m and 20 m resolution and Sentinel-2 Deep Resolution 3 (S2DR3) data with 1 m resolution—to assess their (i) relationship with yield in organically grown barley and (ii) utility for winter crop mapping. Vegetation indices were generated and analysed across different phenological phases to determine the most suitable predictors of yield. The results indicate that using 10 × 10 m data, the BBCH-41 phase is optimal for yield prediction, with the Green Chlorophyll Vegetation Index (GCVI; r = 0.80) showing the strongest correlation with yield. In contrast, S2DR3 data with a 1 × 1 m resolution demonstrated that Transformed the Chlorophyll Absorption in Reflectance Index (TCARI), TO, and Normalised Difference Red Edge Index (NDRE1) were consistently reliable across all phenological stages, except for BBCH-51, which showed weak correlations. These findings highlight the potential of remote sensing in organic barley farming and emphasise the importance of selecting appropriate data resolutions and vegetation indices for accurate yield prediction. With the use of three-date spectral band stacks, the Random Forest (RF) and Support Vector Classification (SVC) methods were used to differentiate between wheat, barley, and rapeseed. A five-fold cross-validation approach was applied, training data were stratified with 200 points per crop, and classification accuracy was assessed using the User’s and Producer’s accuracy metrics through pixel-by-pixel comparison with a reference raster. The results for S2 and S2DR3 were very similar to each other, confirming the significant potential of S2DR3 for high-resolution crop mapping.