Content area
This study aims to address the challenge of monitoring Plant Height (PH), SPAD, Leaf Area Index (LAI), and Above-Ground Biomass (AGB) in Gerbera under greenhouse cultivation conditions. We initially gathered multi-spectral images and corresponding ground truth data of these parameters at various growth stages using a low-altitude UAV. From the collected images, we derived five Vegetation Indices (VIs): NDVI, GNDVI, LCI, NDRE, and OSAVI, and extracted their textural features as fusion features. An adaptive ensemble model, OBM-RFEcv, was then developed by integrating six base models (Linear Regression, Decision Tree Regressor, Random Forest Regressor, XGBoost Regressor, and Support Vector Regressor) with Recursive Feature Elimination (RFE) to predict the key growth indicators. The results indicate that the OBM-RFEcv model outperforms the other models when using the fusion of the five VIs, particularly in the test dataset, where it achieved the highest accuracy for PH (NDVI), SPAD (GNDVI), LAI (GNDVI), and AGB (NDRE) with R2 values of 0.92, 0.90, 0.89, and 0.93, respectively. The root mean square error (RMSE) values were 0.04, 0.07, 0.08, and 0.07, respectively, showing improvements over the best individual model by 0.01, 0.03, 0.01, and 0.09 in R2, and reductions in RMSE by 0.01, 0.07, 0.08, and 0.03, respectively. These findings confirm that the OBM-RFEcv model, based on multi-spectral image fusion, effectively monitors key growth indicators in Gerbera, providing a non-invasive and precise method for greenhouse crop monitoring.
Introduction
Gerbera, commonly known as the African daisy, is one of the world’s top five cut flowers [1]. In China, Gerbera has gained popularity as a cut flower due to its high yield, extensive cultivation, and significant economic returns [1,2]. To ensure robust growth and high-quality flower production, Gerbera requires stringent environmental controls [3]. Key indicators of crop health, such as Plant Height (PH), SPAD, Leaf Area Index (LAI), and Above-Ground Biomass (AGB), are essential for assessing the developmental status of Gerbera [4–6]. Accurate and efficient monitoring of these indicators in greenhouse settings is vital for effective field management [4].
Traditional methods for measuring PH, SPAD, LAI, and AGB include field measurements or laboratory analyses [7]. These methods, particularly those for LAI and AGB, often require destructive sampling at different growth stages, which is both time-consuming and labor-intensive [8,9]. Therefore, there is a pressing need for developing a method that can provide rapid, non-destructive, and accurate monitoring of these critical parameters during the growth cycle of greenhouse crops.
In the context of crop growth monitoring, the most common traditional approach involves the use of tools like DSSAT and WOFOST to construct crop growth models. For example, Mulla et al. [10] calibrated the DSSAT CERES-Maize and CROPGRO-PFM models with high nitrogen application rates to simulate intermediate wheat grass and predict yield, biomass, soil water balance, and soil nitrogen balance under medium and low nitrogen application rates. Wang et al. [11] used the WOFOST model to evaluate the impact of varying model parameters on goji berry yield under different climatic conditions. Tang et al. [12] developed a WOFOST-Chili model to quantify dry matter production as a function of fertilizer management, climate, and soil conditions. Amiri [13] employed the WOFOST model to simulate maize growth, LAI, biomass, grain yield, yield gap, and soil moisture under both irrigation and rainfall conditions, evaluating different growth periods and the entire growing season. Raghavendra et al. [14] used the DSSAT model to estimate potential soybean yields in different regions of India, validating the model under various agricultural ecological conditions and deriving the correlation between soybean yield and different agricultural ecological conditions. These studies have improved crop cultivation practices and and contributed to monitoring of the crop growth process. However, traditional crop growth models rely on a large number of environmental parameter inputs, including meteorological conditions, soil properties, etc. The collection of these data often requires complex equipment and professional technology, leading to high data acquisition costs and low efficiency. Moreover, most of the aforementioned studies have focused on field crops, with relatively few studies on horticultural crops or crops in special environments, which limits the application scope of the models. By using spectral sensors mounted on Unmanned Aerial Vehicles (UAVs), it is possible to efficiently obtain real-time image data of field crops, not only improving the efficiency of data collection but also expanding the monitoring range, which contributes to more accurately assessing crop growth status. Optimizing model parameters through machine learning algorithms can improve the prediction accuracy of the models while reducing the heavy reliance on specific regional and crop data.
In recent years, satellite and Unmanned Aerial Vehicle (UAV)-based multi-spectral remote sensing technologies have gained popularity for monitoring crop growth processes. These technologies primarily use satellite remote sensing platforms to monitor large areas of crop growth and estimate key crop indicators. For instance, Sun et al. [15] used time-series Sentinel-1A satellite data to extract the Dual-Polarization Radar Vegetation index (DPRVI) and integrated it with empirical regression models to invert sugarcane plant height, effectively capturing the changes in plant height during the early growth stage. Ma et al. [16] assimilated the LAI derived from Sentinel-2 satellite data into the SAFY (Simple Algorithm for Yield Estimation crop growth model for wheat yield estimation. Integrating remote sensing data into crop growth models aids in estimating LAI, biomass, and yield. Gupta et al. [17] used LISS-IV satellite data to spatially model SPAD (Soil-Plant Analysis Development) values, enabling effective monitoring of chlorophyll content distribution in agricultural areas. Tang et al. [18] built and validated a customized model for grassland biomass in the Yellow River source region using AGB ground truth data, environmental variables, and vegetation indices, analyzing the spatiotemporal distribution and dynamic trends of grassland biomass over 20 years. Yang et al. [19] used multi-source, multi-angle image data from Landsat-5 and small satellite multispectral sensors to construct a LAI inversion model combining radiative transfer models, mean clustering, and artificial neural networks, improving the average accuracy of LAI inversion by approximately 30% compared to single-angle data. However, in current research, satellite remote sensing images mostly have lower resolution, which limits the ability to accurately monitor small plots or orchards. Moreover, they predominantly focus on specific crop growth indicators, lacking comprehensive monitoring methods that can fully reflect the crop growth status and its environmental impact factors. Combining UAVs with multispectral sensors for low-altitude flights can acquire higher-resolution data, enabling more precise monitoring of crop growth conditions within small plots or specific agricultural areas. Utilizing machine learning or deep learning algorithms to analyze integrated features allows for a deeper assessment of crop growth status while also enhancing the diversity and accuracy of monitoring indicators.
Concurrently, UAV-based multi-spectral remote sensing technology, with its advantages of fast, real-time, and high-resolution imaging, has become a crucial tool in agricultural remote sensing. Zu et al. [20] used multi-spectral images obtained by UAVs to build machine learning and Convolutional Neural Network (CNN) models for estimating winter wheat LAI, demonstrating that the CNN model had the highest estimation accuracy. Ji et al. [21] constructed 26 multi-spectral vegetation indices combined with ground-measured SPAD values of wheat, using four machine learning algorithms to develop SPAD value estimation models for different stages of wheat. Xie et al. [22] used UAV multi-spectral data, combined with vegetation indices and texture features, to predict the SPAD value of lychee fruits during the growth and autumn sprouting periods. Prakriti Sharma et al. [23] used multi-spectral data before and after oat heading obtained by UAVs, combined with four machine learning algorithms, to construct oat biomass estimation models, exploring the potential for estimating AGB in oat breeding nurseries.
Machine learning methods can enhance the prediction accuracy of crop growth traits through data fusion based on UAV multi-spectral images [24]. Yu et al. [25] developed a simple and effective LAI estimation method, using Kalman Filtering (KF) to integrate the estimation results of multiple vegetation index models for precise LAI monitoring during different growth stages of rice. Zhou et al. [26] used machine learning algorithms to develop Canopy Cover Density (CCD), Leaf Chlorophyll Density (LCD), and LAI estimation models by integrating spectral features, texture features, and wavelet coefficient features obtained from multi-spectral images. Ma et al. [27] combined vegetation indices (VIs), color indices (CIs), and texture features to improve the accuracy of assessing crop disease severity from UAV images. Wen-jing [28] formed a multi-source feature dataset using UAV spectral features, thermal features, structural features, and chlorophyll content data, building chlorophyll content estimation models for corn at the jointing, tassel, and large tassel stages using four machine learning algorithms, proving that multi-source feature fusion and stacked ensemble learning can effectively improve the accuracy of chlorophyll content estimation. Naveed et al. [29] used a Predictive Coding/Biased Competition-Divisive Input Modulation (PC/BC-DIM) neural network to determine multi-spectral fusion saliency maps for predicting prominent crops and detecting weeds, achieving higher average accuracy compared to unfused data.
In current research, although machine learning methods have been used to fuse data from UAV multispectral images, thereby improving the prediction accuracy of crop growth traits, most models and methods lack an adaptive mechanism. When faced with changes in environmental conditions, crop types, or growth stages, they cannot automatically adapt to new situations. Despite the fact that data fusion can improve the prediction accuracy of crop growth traits, these methods often involve complex algorithms and require substantial computational resources, increasing the difficulty and cost of practical applications. By simplifying data fusion strategies, reducing the complexity of data processing, and eliminating unnecessary complex computational steps, while incorporating adaptive mechanisms into machine learning models, it is possible to enable these models to automatically adjust their structure based on real-time acquired data. This allows them to adapt to different environmental conditions and crop types.
The aforementioned studies have utilized various methods, including crop growth models, satellite remote sensing, UAV multispectral remote sensing technology, and multispectral image feature fusion, to achieve the inversion, estimation, and prediction of key characteristic values during the crop growth process. Currently, researchers mainly focus on field crops, with multispectral image data primarily collected via high-altitude UAVs. There is very limited research reported on economic crops grown in greenhouses. Moreover, in terms of feature fusion, it often involves complex calculations or simply stacking different features, and there is also limited research on ensemble methods for regression models. Therefore, the main objective of this study is to propose a feature fusion method based on extracting texture features from vegetation index images using multispectral images captured at low altitudes within greenhouse environments, combined with ground-measured data such as PH, SPAD, LAI, and AGB. By employing an adaptive ensemble model for model and feature variable selection, we aim to improve the monitoring accuracy of growth indicators for Gerbera at different growth stages and construct a monitoring model for the growth process of Gerbera. The specific objectives are: (i) to evaluate the performance of the adaptive ensemble model with model and feature variable selection in multi-indicator monitoring during the growth process of Gerbera; (ii) to assess the application potential of the proposed feature fusion method that extracts texture features from different vegetation index images in monitoring key growth indicators of Gerbera. This provides a reference for using multispectral UAVs for low-altitude monitoring of economic crops grown in greenhouses, further promoting the application of multispectral technology in the growth models of economic crops.
Paper structure
The structure of the other sections of this paper is as follows: it is divided into four sections—Materials and methods, Results, Discussion, and Conclusion. Materials and methods mainly introduce the experimental site and planting conditions for the study subject Gerbera, followed by the methods for collecting growth indicators and multispectral image data of Gerbera and the data processing procedures. It concludes with the construction process of the adaptive ensemble model OBM-RFEcv and the evaluation metrics used for the models. In the Results section, it showcases the application effects of the OBM-RFEcv model in monitoring PH, SPAD, LAI, and AGB of Gerbera, comparing the performance between base models and the OBM-RFEcv model. Through charts and data analysis, it reveals the relationships among various parameters and their impacts on Gerbera’s growth status, while also incorporating SHAP visualizations to analyze the interpretability of the OBM-RFEcv model. In the Discussion section, based on the previous results, it delves into the significance of feature fusion in monitoring Gerbera growth indicators, comparing and analyzing these findings with existing literature. It also discusses in-depth the reasons behind the superior performance of the OBM-RFEcv model, its limitations, and potential directions for future research. In the final Conclusion section, it summarizes the main conclusions based on the research findings, emphasizing the practical value of the OBM-RFEcv model, and offers insights into potential future development trends.
Materials and methods
Experimental site and growing conditions
The experimental area is located in the greenhouse of zone 2–22 at the New Agricultural Science Comprehensive Practical Teaching Base of Yunnan Agricultural University (longitude 102°45’49’‘E to 102°45’50’‘E, latitude 25°13’61’‘N to 25°13’63’‘N, see Fig 1). The dimensions of the greenhouse are approximately 35 meters in length, 10 meters in width, and 5 meters in internal height. Four varieties of Gerbera, namely “Yunnanhong”, “Lixiang”, “Taiyang”, and “Hexie”, were cultivated. The layout of the planting areas is illustrated in Fig 1. The planting period was from December 20, 2023, to April 10, 2024. Each variety was planted in an area of approximately 40 square meters, with the entire experimental plot covering about 160 square meters. The planting density was 6–9 plants per square meter, and the average plant spacing was 40 centimeters. The research area’s average growing environment included a carbon dioxide concentration of 313 ppm, soil moisture content of 16.18%, soil temperature of 21.87°C, air temperature of 22°C, air humidity of 70.99%, and solar radiation of 3.57 kilolux (kLx). To ensure healthy growth of Gerbera, organic fertilizer was applied once every 20 days, watering was conducted twice a week, and weeding as well as removal of yellow leaves were performed once a week during the planting period of Gerbera.
[Figure omitted. See PDF.]
Data collection
Data collection was conducted during the growth, early flowering, and full bloom stages of the Gerbera. The data were collected from January 20, 2024, to April 10, 2024, with a total of seven collections: three during the growth stage, two during the early flowering stage, and two during the full bloom stage. The ground data collection for Gerbera growth indicators primarily included PH, LAI, SPAD, and AGB.
First, the experimental area was divided into smaller sections, with each variety being divided into six relatively uniform small sections, resulting in a total of 24 small sections. In each small section, four evenly growing plants were randomly selected as samples for each data collection. After measuring the growth indicators of the four sample plants, the average value was taken as one set of sample data. The methods for collecting and measuring the various growth indicators are as follows:
PH measurement.
The PH of Gerbera was measured as the height of the highest leaf in its natural state relative to the ground. A height ruler was used, held vertically to the ground, for the measurement.
SPAD measurement.
For the collection of SPAD data, the SPAD-502 Plus instrument was used. Functional leaves with a width greater than 3 cm or a length greater than 5 cm were selected from each sample plant, and their SPAD values were measured. The average value of these measurements was taken as the SPAD value for that sample plant.
LAI measurement.
The specific leaf weight method, which uses the ratio of unit leaf area to dry leaf weight to obtain the leaf area index, was employed. This method is simple and relatively accurate, often used for comparative measurements of crops, although the measurement process is complex. To calculate LAI, four sample plants from each small section (A) (area approximately 6 m²) were collected and measured in the laboratory. Before measurement, the total number of plants (N) in the sampling area was recorded. All expanded green leaves were removed from the samples, and 2 cm long segments were cut from parts of the leaves where the width was consistent. The total width was measured using a ruler to calculate the area (S). The segments were then dried and weighed (M1), and the remaining green leaves were also dried and weighed (M2) [30]. The formula for calculating LAI is:
(1)
The average LAI value of the four sample plants was calculated to represent the LAI of the small section.
AGB measurement.
For AGB measurement, four sample plants were collected from each small section and placed in a vacuum drying oven at a constant temperature of 100°C until they reached a constant weight (approximately 6 hours) [31]. The samples were then weighed on an electronic scale with a precision of 0.01 g. The average value of the four sample plants was taken as the AGB value for the small section.
Multi-spectral image data collection.
Multispectral images of the Gerbera at different growth stages across the entire experimental area were collected using a DJI Mavic 3M (DJI, Shenzhen, China). This device integrates four sensors with different central wavelengths (green: 560 ± 16 nm, red: 650 ± 16 nm, red edge: 730 ± 16 nm, near-infrared: 840 ± 26 nm), with a resolution of 2594*1944 pixels. The frequency of multispectral image collection was synchronized with the ground data collection, and flights were conducted between 9:00 and 15:00 on the same day, before the ground data collection. To reduce light interference within the greenhouse, the shading curtains were opened before the flight to ensure uniform lighting across all areas. During the acquisition of multispectral images, the UAV was manually flown. Based on the height of the greenhouse, the safe flying distance for the UAV, and the efficiency of image collection, the flight altitude was maintained at 3 meters with a spatial resolution of 1.57 cm, ensuring an image overlap greater than 80%. This was done to ensure the acquisition of complete high-resolution multispectral images.
Data processing
Multi-spectral image processing.
After the multispectral image data collection, images from each of the four bands collected during each flight were selected according to the number of collections. Using DJI Terra 4.0.1 (which includes image correction and aerial triangulation optimization), 2D reconstructions of visible light images and Vegetation Index (VIs) images were generated for each collection. A total of five VIs images were generated for each collection, as shown in Table 1: Normalized Difference Vegetation Index (NDVI) images, Green Normalized Difference Vegetation Index (GNDVI) images, Chlorophyll Index (LCI) images, Normalized Difference Red Edge (NDRE) images, and Optimized Soil-Adjusted Vegetation Index (OSAVI) images. After generating the VIs images, they were cropped based on the predefined experimental area, resulting in 24 images corresponding to the 24 small sections of the experimental area. Over the entire growing cycle, a total of 840 VIs images (NDVI, GNDVI, LCI, NDRE, OSAVI) were obtained, which were used for the next stage of feature extraction.
[Figure omitted. See PDF.]
Fusion feature extraction.
The fusion features in this study are based on a deep integration of VIs images with texture features. The Gray-Level Co-occurrence Matrix (GLCM) method was used to calculate Contrast, Homogeneity, Energy, Correlation, and Entropy [37] for the five VIs images. These features help represent the complexity of the leaf surface structure, health status, and boundaries between plants. In this experiment, calculations were performed using the skimage package in Python. Firstly, the pixel values of the images were linearly scaled to a range between 0 and 255. Then, the GLCM was calculated for adjacent pixels at angles of horizontal, 45 degrees to the upper right, vertical, and 135 degrees to the upper left, with the image grayscale levels set to 256. Finally, the five features of the GLCM were averaged. A total of 840 fusion feature samples of texture features were extracted from five VIs during three growth stages of Gerbera.
Data preprocessing.
To convert the 840 fusion features and growth indicators to the same scale, we used the Z-Score method to eliminate unit and scale differences before modeling. This method also helps detect and handle potential outliers in the dataset. Z-Score is a statistical measure that describes how many standard deviations an element is from the mean. It is calculated by subtracting the mean from the value and dividing by the standard deviation [38].The formula for Z-Score is as follows:
(2)
Where is the data point, is the mean of all data points, is the standard deviation, and is the Z-Score value [39].
In this study, Z-Score was calculated for each column of data to detect outliers (using a default threshold of 3; if the absolute value of a data point’s Z-Score exceeds 3, it is considered an outlier). These outliers were then replaced with the median of that column to mitigate their impact on the analysis. After handling the outliers, Z-Score normalization was performed on the processed data. This converts the data into a form where the mean is 0 and the standard deviation is 1, facilitating subsequent data analysis or model training. An example of the processed data after Z-Score normalization is shown in Table 2. We obtained 840 data points, including five fusion features and growth indicators, for subsequent analysis and modeling.
[Figure omitted. See PDF.]
Model evaluation metrics.
In this study, the Coefficient of Determination (R2) and Root Mean Square Error (RMSE) were chosen to evaluate the accuracy of model training and validation. The R2 represents the proportion of the variance in the dependent variable that is predictable from the independent variable. The formula for R2 is as follows:
(3)
where is the observed value, is the predicted value, is the mean of the observed values, and is the number of observations [40]. The value of R2 ranges from 0 to 1, where 0 indicates that the model does not explain any of the variability in the data, and 1 indicates that the model perfectly explains all the variability in the data.
The RMSE is a measure of the differences between the predicted and actual values. It provides a standard size of prediction error and has the same units as the original data, making it easy to understand and interpret. The formula for RMSE is as follows:
(4)
where is the observed value, is the predicted value, and is the number of observations [40]. In model evaluation, a smaller RMSE value indicates better model performance.
Construction of the OBM-RFEcv adaptive ensemble model.
In this study, we extracted five texture features (Contrast, Homogeneity, Energy, Correlation, and Entropy) from the five VIs images (NDVI, GNDVI, LCI, NDRE, OSAVI) using the GLCM method. These features were used to build the models. Six machine learning regression algorithms were selected as base models and comparison models: Linear Regression (LR) [41], Decision Tree Regression (DTR) [42], Random Forest Regression (RFR) [43], Gradient Boosting Regression (GBR) [44], Extreme Gradient Boosting Regression (XGBR) [45], and Support Vector Regression (SVR) [46]. First, the fusion features were used as model inputs to construct six base models: LR, DTR, RFR, GBR, XGBR, and SVR. During the model training and evaluation phase, 5-fold cross-validation was incorporated to ensure model stability. After calculating the R2 and RMSE for these six base models, the Python ‘max’ function was used to automatically select the base model with the highest R2 value as the Optimal Base Model (OBM). This OBM also serves as the base model for the ensemble model.
To recursively remove the least important features during model training, we used the cross-validated version of Recursive Feature Elimination (RFEcv) for feature selection. RFEcv repeatedly builds models and retains the best features until a predetermined number of features is reached. The core idea of RFE is to judge the importance of features based on model weights and gradually eliminate less important features to optimize model performance [47]. Compared to RFE, RFEcv can automatically determine the optimal number of features through cross-validation. The RFEcv algorithm process is shown in Fig 2. Since RFEcv is model-based, it can be integrated into various types of models [48].
[Figure omitted. See PDF.]
As shown in Fig 2, in this experiment, the Random Forest Regressor algorithm was used as the estimator for RFEcv. In each iteration, one feature was removed, and 5-fold cross-validation was used to determine the optimal number of features. The negative mean squared error (NMSE) of each feature was calculated as the evaluation metric for automatic feature selection.
After selecting the OBM, we integrated the RFEcv algorithm into the OBM to form an adaptive ensemble model, OBM-RFEcv, for monitoring the growth indicators of Gerbera. The technical route of this study is shown in Fig 3.
[Figure omitted. See PDF.]
Results
This study developed monitoring models for the growth stages (growth, early flowering, and full bloom) of Gerbera, specifically for PH, SPAD, LAI, and AGB, based on the fusion features extracted from low-altitude UAV multispectral images, including VIs and texture features. We used the fusion feature variables extracted from five VIs images across three growth stages of Gerbera to build 140 monitoring models (4*5*7) for four growth indicators using seven models. We compared and analyzed the fusion features of five different VIs images for PH, SPAD, LAI, and AGB in Gerbera using six base models and the OBM-RFEcv adaptive ensemble model. The performance of the base models and the adaptive ensemble model was compared, and the model performance was tested on a test set using R2 and RMSE as evaluation metrics. Regarding the dataset division, 80% of the samples were used for training and 20% for testing in this experiment. For model parameters, n_estimators = 100 and random_state = 42 were set, with all other parameters being default values. The cross-validation was performed using 5-fold cross-validation.
Monitoring results of Gerbera PH indicators with different VIs image fusion features
As shown in Table 3, the combination of the OBM-RFEcv model and the NDVI fusion feature performed best in the PH indicator monitoring results, with R2 and RMSE values of 0.92 and 0.04, respectively. Compared to the optimal base models RFR and GBR, the R2 value increased by 0.08, and the RMSE decreased by 0.01. Although the base models did not perform optimally with NDVI compared to other VIs, the accuracy improved significantly after incorporating RFEcv, with R2 values 0.2 to 0.5 higher and the highest RMSE reduction of 0.2. For OSAVI, although the base models already showed high accuracy, the inclusion of RFEcv reduced the accuracy, indicating that the fusion features based on OSAVI do not require feature selection, but the accuracy is still lower than the combination of OBM-RFEcv and NDVI.
[Figure omitted. See PDF.]
From Fig 4, it can also be seen that the OBM-RFEcv model, compared to the base models, not only showed a significant improvement in R2 and RMSE but also made the evaluation results of the fusion features of the five different VIs more balanced. The most significant improvement was observed with the NDVI feature, indicating that the combination of the OBM-RFEcv model and the NDVI fusion feature is the most effective in capturing and monitoring plant height. Based on the PH monitoring model with higher R2 and lower RMSE, managers can better understand the crop’s growth rate and developmental stages. Accurate PH data can also assist in planning more reasonable planting layouts and densities.
[Figure omitted. See PDF.]
Fig 5 illustrates the impact of two texture features, Energy and Contrast, on the model’s prediction results after selection by the OBM-RFEcv model. For the Energy feature, most points are located on the right side, with more red points than blue ones, while the left side is mostly blue points. As the feature value increases (from blue to red), its positive contribution to the model’s prediction gradually strengthens. For the Contrast feature, the points are mainly concentrated on the left side, with colors ranging from red to purple, indicating a significant negative contribution to the model’s prediction. As the feature value increases (from red to purple), its negative contribution to the model’s prediction gradually decreases.
[Figure omitted. See PDF.]
Monitoring results of Gerbera SPAD indicators with different VIs image fusion features
According to the SPAD monitoring results in Table 4, the OBM-RFEcv model also showed higher accuracy than the base models, especially when combined with the GNDVI, where R2 and RMSE were 0.90 and 0.07, respectively. Compared to the optimal base model GBR, the R2 value increased by 0.03, and the RMSE decreased by 0.01. Compared to other VIs, the R2 value increased by 0.02 to 0.17, and the RMSE decreased by up to 0.03. Similarly, in the GNDVI, the average performance of the base models was better than in other VIs, indicating that the GNDVI feature is more sensitive to SPAD in Gerbera. The lowest performance was observed with the LCI, with R2 and RMSE values of 0.73 and 0.10, respectively, indicating a lower correlation between SPAD and the LCI index.
[Figure omitted. See PDF.]
Based on Fig 6, it is evident that the OBM-RFEcv model achieved a significant increase in R2 compared to the base models for all five VIs image fusion features, and the RMSE values were the lowest. The performance of all fusion features was relatively balanced, with the most significant improvement observed in GNDVI. Higher R2 and lower RMSE allow managers to gain a deeper understanding of the crop’s health status and nutrient absorption efficiency. By monitoring changes in SPAD values, it is possible to quickly detect whether crops suffer from nutrient deficiencies such as nitrogen deficiency. Accurate SPAD values can help determine the optimal timing for fertilization, ensuring that necessary nutrients are provided precisely when the crops need them most, thereby promoting smooth transitions through all growth stages.
[Figure omitted. See PDF.]
Fig 7 illustrates the impact of five selected texture features—Contrast, Homogeneity, Energy, Entropy, and Correlation—on the model’s prediction results after selection by the OBM-RFEcv model. Except for the Entropy feature, the points for the other four features are concentrated on the right side, indicating that these four features generally have a positive contribution to the model’s predictions. Conversely, the Entropy feature shows an opposite trend. As the value of the Contrast feature increases (from blue to red), its negative contribution to the model’s prediction gradually strengthens. In Contrast, for the other four features (Homogeneity, Energy, Correlation), as their values increase (from blue to red), their positive contributions to the model’s prediction gradually strengthen. Notably, the Homogeneity feature exhibits a higher positive contribution compared to the other three features (Energy, Correlation, and Contrast).
[Figure omitted. See PDF.]
Monitoring results of Gerbera LAI indicators with different VIs image fusion features
As shown in Table 5, the OBM-RFEcv model combined with the GNDVI performed better and had higher accuracy than other combinations for LAI monitoring. The R2 and RMSE values were 0.89 and 0.08, respectively, which were 0.12 higher in R2 and 0.02 lower in RMSE compared to the optimal base model XGBR. In the OBM-RFEcv model, the R2 value for the GNDVI feature was 0.01 to 0.14 higher than for other VI combinations, and the RMSE values also showed good performance, decreasing by 0.01 to 0.04. The worst performance was observed with the OSAVI, indicating that it is not sensitive enough to the LAI of Gerbera, and there were significant differences in the performance of different models.
[Figure omitted. See PDF.]
From Fig 8, it can also be seen that the GNDVI showed significantly better R2 and RMSE values than other indices, and the OBM-RFEcv model’s performance was more accurate than the base models. This indicates that the OBM-RFEcv model has good generalization ability for different VIs features, and combining it with the GNDVI is more suitable for predicting the LAI of Gerbera. In LAI monitoring, higher R2 and lower RMSE allow managers to more accurately understand changes in crop canopy structure and photosynthesis efficiency. Changes in LAI can serve as early warning signals, indicating potential pest infestations or nutrient imbalances.
[Figure omitted. See PDF.]
Fig 9 illustrates the impact of three selected texture features—Entropy, Contrast, and Correlation—on the model’s prediction results after selection by the OBM-RFEcv model. The points for these three features are mostly concentrated on the left side, indicating that most of their feature values have a slight negative contribution to the model’s predictions.As the feature values increase (from blue to red), the positive contributions of the Entropy and Correlation features to the model’s predictions gradually strengthen. Among them, the Entropy feature has a greater positive contribution compared to the Correlation feature. Conversely, the Contrast feature shows an opposite trend, where its negative contribution gradually strengthens as its feature values increase.
[Figure omitted. See PDF.]
Monitoring results of Gerbera AGB indicators with different VIs image fusion features
From Table 6, the OBM-RFEcv model combined with the NDRE performed best in AGB monitoring, with R2 and RMSE values of 0.93 and 0.07, respectively. The R2 value was the highest among all growth indicators. Compared to the optimal base model RFR, the R2 value of the OBM-RFEcv model increased by 0.11, and the RMSE decreased by 0.04, showing a significant improvement. Although the base models did not perform optimally with the NDRE index, the accuracy improved significantly after using the OBM-RFEcv model, indicating that not all features are sensitive to AGB in the NDRE, and feature selection is necessary.
[Figure omitted. See PDF.]
In Fig 10, the OBM-RFEcv model performed best among all VIs image features, and the results of all features based on the OBM-RFEcv model were more balanced, indicating that the OBM-RFEcv model can significantly improve the prediction accuracy of various vegetation features, with stronger overall performance, making it very suitable for monitoring the AGB of Gerbera. By accurately assessing the biomass accumulation of crops, managers can adjust fertilization plans according to the actual needs of the crops, promoting healthy growth while reducing the use of chemical fertilizers and minimizing environmental impact. Utilizing improved AGB monitoring results, farmers can predict the final crop yield in advance.
[Figure omitted. See PDF.]
Fig 11 illustrates the impact of three selected texture features—Entropy, Contrast, and Correlation—on the model’s prediction results after selection by the OBM-RFEcv model. The points for these three features are mostly concentrated on the left side, indicating that most of their feature values have a negative contribution to the model’s predictions.
[Figure omitted. See PDF.]
As the feature values increase (from blue to red), the positive contributions of the Entropy and Correlation features to the model’s predictions gradually strengthen. Similar to the LAI monitoring model, the Entropy feature has a greater positive contribution compared to the Correlation feature. Conversely, the Contrast feature shows an opposite trend, where its negative contribution to the model’s predictions gradually strengthens as its feature values increase.
Discussion
Impact of fusion features based on different VIs images on estimating Gerbera growth indicators
By utilizing high-resolution multispectral images obtained from UAVs, we extracted five texture features—Contrast, Homogeneity, Energy, Correlation, and Entropy—from five VIs images using the GLCM technique. These features, when combined with traditional VIs, provided a richer information source for estimating PH, SPAD, LAI, and AGB in Gerbera. Wei et al. [49] monitored chlorophyll content in winter wheat canopies under different nitrogen treatments using UAV multispectral images, selecting 16 spectral vegetation indices to establish SPAD estimation models for different growth stages, with an R2 of 0.73 in the optimal stepwise regression model. Guomin et al. [50] used five VIs as input parameters and employed the random forest regression algorithm to establish the relationship between canopy vegetation indices and LAI for field maize under different irrigation conditions, with an R2 of 0.87. Guo et al. [51] used images acquired from a multispectral imaging spectrometer mounted on a UAV and extracted plant height from Digital Surface Model (DSM) information to establish an improved AGB estimation model using back propagation neural networks, with an R2 of 0.88. Despite the innovations in spectral data processing and the relatively high model accuracy achieved in the aforementioned studies, most of them rely solely on vegetation indices or directly use raw spectral data for crop growth status assessment during feature extraction. Additionally, these studies do not employ model ensemble techniques in model construction, which may result in weaker generalization ability and poorer applicability of the models. In contrast, this study employs a feature fusion strategy that extracts texture features from VI images and constructs an adaptive ensemble model. This approach can more accurately capture subtle changes during the plant growth process and enhances the model’s generalization ability.
Additionally, by using various regression models to select the fusion features extracted from different VIs images, we further ensured that the selected features have strong correlations with Gerbera growth indicators. This approach not only effectively captures the complex relationships between different VIs and target variables, improving model prediction accuracy, but also reduces unnecessary computational burden. For economic crops like Gerbera, rapid and accurate monitoring of growth status is crucial for optimizing field management decisions.
Constructing the adaptive ensemble model OBM-RFEcv to improve overall model performance
To determine the effectiveness of the adaptive ensemble model OBM-RFEcv in monitoring key growth indicators of Gerbera, we compared the performance of multiple single regression models, including LR, DTR, RFR, GBR,XGBR and SVR. The experimental results showed that while some single models performed well under specific conditions, their overall predictive capabilities were limited.
In contrast, the proposed OBM-RFEcv adaptive ensemble model, which combines the feature selection method RFEcv with the optimal model among the various regression models, improved the estimation accuracy of Gerbera PH, SPAD, LAI, and AGB across different VIs, showing better overall performance. Especially in estimating above-ground biomass, the OBM-RFEcv + NDRE combination achieved R2 = 0.93 and RMSE = 0.07, significantly outperforming the best-performing model without feature selection. This indicates that the RFEcv method, by recursively removing the least important features and retraining the model, is capable of identifying the most predictive subset of features. This approach not only reduces model complexity but also avoids overfitting, making the model more efficient and accurate when handling complex nonlinear relationships. The OBM-RFEcv enhances this by adaptively selecting the optimal base models for ensemble learning, improving the fitness of different data across various models, thereby enhancing overall predictive capability. This strategy is particularly suitable for agricultural data characterized by high heterogeneity and complexity, aiding in capturing subtle changes under different growth stages and environmental conditions.Future research could further expand the number of base models to better match different datasets, improving the universality and practicality of the models.
Whether for field crops or economic crops, research related to the adaptive ensemble model strategy is still relatively scarce. Although the adaptive ensemble model proposed in this study has shown good performance and has been tested on multiple varieties of Gerbera, the number of varieties covered by the current study remains limited compared to the rich diversity of Gerbera species. Significant genetic differences, morphological characteristics, and varying responses to external environments among different varieties may affect the prediction accuracy of the model. In addition to increasing the number of Gerbera varieties, it is also necessary to validate the model on other types of economic crops to test its generalizability across a broader range of crop types.
Even though the model has been tested on various Gerbera varieties, these tests were mostly conducted within specific geographic regions and under relatively consistent planting conditions. Different climatic zones, soil conditions, and agricultural management practices can have significant impacts on crop growth. Therefore, testing the model’s performance in more diverse environments is essential.
In addition to utilizing data acquired via UAVs, future work could consider integrating information from other sources. Combining with micro-environmental data such as soil moisture and temperature provided by ground sensor networks could offer richer input information for crop growth models. For instance, real-time monitoring of soil moisture can help adjust irrigation schedules, while temperature data can assist in predicting the risk of pest and disease outbreaks. Integrating data from meteorological stations, such as precipitation, wind speed, and light intensity, could further refine crop growth models, enabling them to better reflect the impact of external environments on crops. This comprehensive approach would lead to the construction of more thorough and precise models for monitoring the crop growth process.
Conclusions
This study explored the impact of fusion features based on different VIs images and texture features on the monitoring performance of Gerbera growth indicators. The OBM-RFEcv adaptive ensemble model was used to analyze and evaluate the monitoring results of Gerbera PH, SPAD, LAI, and AGB, and comparisons were made with base models.
The results showed that the OBM-RFEcv model exhibited high accuracy in estimating Gerbera PH, SPAD, LAI, and AGB, yielding good prediction results for fusion features of all VIs images. Compared to the base models, the OBM-RFEcv model achieved an approximately 2% improvement in R2 and RMSE, with more balanced performance. Specifically, the OBM-RFEcv model performed best in monitoring PH with NDVI, achieving R2 = 0.92 and RMSE = 0.04; in monitoring SPAD, GNDVI performed best, with R2 = 0.90 and RMSE = 0.07; in monitoring LAI, GNDVI performed best, with R2 = 0.89 and RMSE = 0.08; and finally, in monitoring AGB, NDRE performed best, with R2 = 0.93 and RMSE = 0.07. In summary, the feature fusion method based on extracting texture features from different VIs images has high application potential and practicality for monitoring growth indicators during the growth stages of Gerbera. The OBM-RFEcv model demonstrated more balanced and superior performance than the base models across different VIs image fusion features. The advantage of the OBM-RFEcv model lies in its use of the RFEcv feature selection method to optimize input variables, reducing the influence of irrelevant features and thus improving prediction accuracy. Furthermore, the OBM-RFEcv model more effectively captured the complex relationships between different VIs and growth indicators, significantly enhancing the model’s generalization ability, particularly under NDVI, GNDVI, and NDRE. These findings provide important references for future research and applications in monitoring the growth indicators of economic crops in greenhouse environments using low-altitude multispectral UAVs.
Supporting information
S1 File. The original data used in this study.
https://doi.org/10.1371/journal.pone.0322851.s001
(XLSX)
References
1. 1. Brandler D, Divensi LJ, Tonin RJ, Pilla TP, Rezendes I, Milanesi PM. Evaluation of biological control of fusarium wilt in gerbera with Trichoderma asperellum. OH. 2017;23(3):234.
* View Article
* Google Scholar
2. 2. Lim M-J, Han J-E, Murthy HN, Song H-Y, Lee S-Y, Park S-Y. A Comparison of semi-solid, liquid, and temporary immersion bioreactor systems for effective plant regeneration of Gerbera jamesonii “Shy Pink”. Horticulturae. 2024;10(8):836.
* View Article
* Google Scholar
3. 3. JNeog PP, Deka P. Screening of Gerbera cultivars against root-knot nematode meloidogyne incognita. J Exp Agric Int. 2024;46(9):899–905.
* View Article
* Google Scholar
4. 4. Su X, Nian Y, Shaghaleh H, Hamad A, Yue H, Zhu Y, et al. Combining features selection strategy and features fusion strategy for SPAD estimation of winter wheat based on UAV multispectral imagery. Front Plant Sci. 2024;15:1404238. pmid:38799101
* View Article
* PubMed/NCBI
* Google Scholar
5. 5. Song E, Shao G, Zhu X, Zhang W, Dai Y, Lu J. Estimation of plant height and biomass of rice using unmanned aerial vehicle. Agronomy. 2024;14(1):145.
* View Article
* Google Scholar
6. 6. Zhang X, Zhang K, Sun Y, Zhao Y, Zhuang H, Ban W, et al. Combining spectral and texture features of UAS-based multispectral images for maize leaf area index estimation. Remote Sensing. 2022;14(2):331.
* View Article
* Google Scholar
7. 7. Guo Y, Chen S, Li X, Cunha M, Jayavelu S, Cammarano D, et al. Machine learning-based approaches for predicting SPAD values of maize using multi-spectral images. Remote Sensing. 2022;14(6):1337.
* View Article
* Google Scholar
8. 8. Bajocco S, Ginaldi F, Savian F, Morelli D, Scaglione M, Fanchini D, et al. On the use of NDVI to estimate LAI in field crops: implementing a conversion equation library. Remote Sensing. 2022;14(15):3554.
* View Article
* Google Scholar
9. 9. Liu Y, Feng H, Yue J, Jin X, Fan Y, Chen R, et al. Improved potato AGB estimates based on UAV RGB and hyperspectral images. Computers and Electronics in Agriculture. 2023;214:108260.
* View Article
* Google Scholar
10. 10. Mulla DJ, Tahir M, Jungers JM. Comparative simulation of crop productivity, soil moisture and nitrate-N leaching losses for intermediate wheatgrass and maize in Minnesota using the DSSAT model. Front Sustain Food Syst. 2023;7.
* View Article
* Google Scholar
11. 11. Wang Z, Shi Y, Hou C, Zhang P. Sensitivity analysis of simulated Lycium barbarum L. yield in the WOFOST model under different climate conditions. Ecological Modelling. 2024;488:110602.
* View Article
* Google Scholar
12. 12. Tang R, Supit I, Hutjes R, Zhang F, Wang X, Chen X, et al. Modelling growth of chili pepper (Capsicum annuum L.) with the WOFOST model. Agricultural Systems. 2023;209:103688.
* View Article
* Google Scholar
13. 13. Ebrahim A, Suat I, Hadis Y. Performance of WOFOST model for simulating maize growth, leaf area index, biomass, grain yield, yield gap, and soil water under irrigation and rainfed conditions. J Irrigation and Drainage Engineering. 2022; 148(2):05021005.
* View Article
* Google Scholar
14. 14. Nargund R, Bhatia VS, Sinha NK, Mohanty M, Jayaraman S, Dang YP, et al. Assessing soybean yield potential and yield gap in different agroecological regions of india using the DSSAT model. Agronomy. 2024;14(9):1929.
* View Article
* Google Scholar
15. 15. Sheng S, Lilu L, Zhongwen H, Xu Y. Inversion method of sugarcane plant height based on Sentinel-1A dual-polarization time series data. Transactions of the Chinese society for agricultural machinery. 2021; 53(186-194):1000–298.
* View Article
* Google Scholar
16. 16. Ma C, Liu M, Ding F, Li C, Cui Y, Chen W, et al. Wheat growth monitoring and yield estimation based on remote sensing data assimilation into the SAFY crop growth model. Sci Rep. 2022;12(1):5473. pmid:35361910
* View Article
* PubMed/NCBI
* Google Scholar
17. 17. Gupta DK, Prasad R, Kumar P, Mishra VN, Vishwakarma AK, Singh RS, et al. Spatial modeling of SPAD values for different type of crops using LISS-IV satellite imagery. 2015 International Conference on Microwave, Optical and Communication Engineering (ICMOCE). IEEE, 2015:212–215.
18. 18. Tang R, Zhao Y, Lin H. Spatio-temporal variation characteristics of aboveground biomass in the headwater of the yellow river based on machine learning. Remote Sensing. 2021;13(17):3404.
* View Article
* Google Scholar
19. 19. Yang G, Zhao C, Liu Q, Huang W, Wang J. Inversion of a radiative transfer model for estimating forest lai from multisource and multiangular optical remote sensing data. IEEE Trans Geosci Remote Sensing. 2011;49(3):988–1000.
* View Article
* Google Scholar
20. 20. Zu J, Yang H, Wang J, Cai W, Yang Y. Inversion of winter wheat leaf area index from UAV multispectral images: classical vs. deep learning approaches. Front Plant Sci. 2024;15:1367828. pmid:38550285
* View Article
* PubMed/NCBI
* Google Scholar
21. 21. Ji J, Li N, Cui H, Li Y, Zhao X, Zhang H, et al. Study on Monitoring SPAD Values for Multispatial Spatial Vertical Scales of Summer Maize Based on UAV Multispectral Remote Sensing. Agriculture. 2023; 13(5):1004.
* View Article
* Google Scholar
22. 22. Xie J, Wang J, Chen Y, Gao P, Yin H, Chen S, et al. Estimating the SPAD of litchi in the growth period and autumn shoot period based on UAV multi-spectrum. Remote Sensing. 2023;15(24):5767.
* View Article
* Google Scholar
23. 23. Sharma P, Leigh L, Chang J, Maimaitijiang M, Caffé M. Above-ground biomass estimation in oats using UAV remote sensing and machine learning. Sensors (Basel). 2022;22(2):601. pmid:35062559
* View Article
* PubMed/NCBI
* Google Scholar
24. 24. Yang Z, Yu Z, Wang X, Yan W, Sun S, Feng M, et al. Estimation of millet aboveground biomass utilizing multi-source UAV image feature fusion. Agronomy. 2024;14(4):701.
* View Article
* Google Scholar
25. 25. Yu M, He J, Li W, Zheng H, Wang X, Yao X, et al. Estimation of rice leaf area index utilizing a kalman filter fusion methodology based on multi-spectral data obtained from Unmanned Aerial Vehicles (UAVs). Remote Sensing. 2024;16(12):2073.
* View Article
* Google Scholar
26. 26. Zhou L, Nie C, Su T, Xu X, Song Y, Yin D, et al. Evaluating the canopy chlorophyll density of maize at the whole growth stage based on multi-scale UAV image feature fusion and machine learning methods. Agriculture. 2023;13(4):895.
* View Article
* Google Scholar
27. 27. Ma R, Zhang N, Zhang X, Bai T, Yuan X, Bao H, et al. Cotton Verticillium wilt monitoring based on UAV multispectral-visible multi-source feature fusion. Computers and Electronics in Agriculture. 2024;217:108628.
* View Article
* Google Scholar
28. 28. Wen-jing Z, Zhan-kang F, Shi-yuan D, Zhang P-P, Wen J, Ai-chen W, et al. Multi-feature fusion detection of wheat lodging information based on UAV multispectral images. Spectroscopy and Spectral Analysis. 2024; 44:197-206.
* View Article
* Google Scholar
29. 29. Naveed A, Muhammad W, Irshad MJ, Aslam MJ, Manzoor SM, Kausar T, et al. Saliency-based semantic weeds detection and classification using UAV multispectral imaging. IEEE Access. 2023;11:11991–2003.
* View Article
* Google Scholar
30. 30. Daughtry CST, Hollinger SE. Costs of measuring leaf area index of corn1. Agronomy Journal. 1984;76(5):836–41.
* View Article
* Google Scholar
31. 31. Zheng H, Cheng T, Zhou M, Li D, Yao X, Tian Y, et al. Improved estimation of rice aboveground biomass combining textural and spectral analysis of UAV imagery. Precision Agric. 2018;20(3):611–29.
* View Article
* Google Scholar
32. 32. JMao D, Wang Z, Luo L, Ren C. Integrating AVHRR and MODIS data to monitor NDVI changes and their relationships with climatic parameters in Northeast China. Int J Appl Earth Observation and Geoinformation. 2012;18:528–36.
* View Article
* Google Scholar
33. 33. Gitelson AA, Kaufman YJ, Merzlyak MN. Use of a green channel in remote sensing of global vegetation from EOS-MODIS. Remote Sensing of Environment. 1996;58(3):289–98.
* View Article
* Google Scholar
34. 34. Hydrology UA, MD RSLB, Hydrology UA, MD RSLB, ID UoIGLfEDM, OR UACPCRCP. Remote sensing leaf chlorophyll content using a visible band index. Agronomy Journal. 2011; 103(4):1090–99.
* View Article
* Google Scholar
35. 35. Li W, Wang J, Zhang Y, Yin Q, Wang W, Zhou G, et al. Combining texture, color, and vegetation index from unmanned aerial vehicle multispectral images to estimate winter wheat leaf area index during the vegetative growth stage. Remote Sensing. 2023; 15(24).
* View Article
* Google Scholar
36. 36. Rondeaux G, Steven M, Baret F. Optimization of soil-adjusted vegetation indices. Remote Sensing of Environment. 1996;55(2):95–107.
* View Article
* Google Scholar
37. 37. JMutlag WK, Ali SK, Aydam ZM, Taher BH. Feature extraction methods: a review. J Phys: Conf Ser. 2020;1591(1):012028.
* View Article
* Google Scholar
38. 38. Curtis AE, Smith TA, Ziganshin BA, Elefteriades JA. The mystery of the Z-score. Aorta (Stamford). 2016;4(4):124–30. pmid:28097194
* View Article
* PubMed/NCBI
* Google Scholar
39. 39. JWu H, Hayes MJ, Weiss A, Hu Q. An evaluation of the standardized precipitation index, the China‐Z Index and the statistical Z‐Score. Intl J Climatology. 2001;21(6):745–58.
* View Article
* Google Scholar
40. 40. Kumar C, Mubvumba P, Huang Y, Dhillon J, Reddy K. Multi-stage corn yield prediction using high-resolution UAV multispectral data and machine learning models. Agronomy. 2023;13(5):1277.
* View Article
* Google Scholar
41. 41. Forkuor G, Hounkpatin OKL, Welp G, Thiel M. High resolution mapping of soil properties using remote sensing variables in south-western burkina faso: a comparison of machine learning and multiple linear regression models. PLoS One. 2017;12(1):e0170478. pmid:28114334
* View Article
* PubMed/NCBI
* Google Scholar
42. 42. XU M, WATANACHATURAPORN P, VARSHNEY P, ARORA M. Decision tree regression for soft classification of remote sensing data. Remote Sensing of Environment. 2005;97(3):322–36.
* View Article
* Google Scholar
43. 43. Biau G, Scornet E. A random forest guided tour. TEST. 2016;25(2):197–227.
* View Article
* Google Scholar
44. 44. eNatekin A, eKnoll A. Gradient Boosting Machines, A Tutorial. Frontiers in Neurorobotics. 2013; 7:21.
* View Article
* Google Scholar
45. 45. Chen T, Guestrin C. XGBoost: A scalable tree boosting system. Proceedings of the 22nd acm sigkdd international conference on knowledge discovery and data mining. 2016:785–794.
46. 46. Ben-Hur A, Ong CS, Sonnenburg S, Schölkopf B, Rätsch G. Support vector machines and kernels for computational biology. PLoS Comput Biol. 2008;4(10):e1000173. pmid:18974822
* View Article
* PubMed/NCBI
* Google Scholar
47. 47. Khaing KT. Enhanced features ranking and selection using recursive feature elimination(rfe) and k-nearest neighbor algorithms in support vector machine for intrusion detection system. International J Network and Mobile Technologies. 2010; 1(1):8.
* View Article
* Google Scholar
48. 48. Yan K, Zhang D. Feature selection and analysis on correlated gas sensor data with recursive feature elimination. Sensors and Actuators B: Chemical. 2015;212:353–63.
* View Article
* Google Scholar
49. 49. Wei Q, Zhang B-Z, Wei Z, Han X, Duan C-F. Estimation of canopy chlorophyll content in winter wheat by UAV multispectral remote sensing. CABI Digital Library. 2020:365–372.
* View Article
* Google Scholar
50. 50. Guomin S, Yajie W, Wenting H. Estimation method of leaf area index for summer maize using UAV-based multispectral remote sensing. Smart Agriculture. 2020; 2(3):118.
* View Article
* Google Scholar
51. 51. Guo Y, He J, Zhang H, Shi Z, Wei P, Jing Y, et al. Improvement of winter wheat aboveground biomass estimation using digital surface model information extracted from unmanned-aerial-vehicle-based multispectral images. Agriculture. 2024;14(3):378.
* View Article
* Google Scholar
Citation: Wang X, Shen Y, Tian P, Wu M, Li Z, Zhao J, et al. (2025) OBM-RFEcv: An adaptive ensemble model for monitoring key growth indicators of Gerbera using multi-spectral image fusion features. PLoS One 20(5): e0322851. https://doi.org/10.1371/journal.pone.0322851
About the Authors:
Xinrui Wang
Roles: Conceptualization, Formal analysis, Investigation, Methodology, Validation, Writing – original draft
Affiliation: College of Big Data, Yunnan Agricultural University, Kunming, Yunnan, China
Yingming Shen
Roles: Writing – original draft, Writing – review & editing
Affiliation: International Cooperation Office, Yunnan Provincial Academy of Science and Technology, Kunming, Yunnan, China
Peng Tian
Roles: Validation, Writing – original draft
Affiliation: College of Big Data, Yunnan Agricultural University, Kunming, Yunnan, China
Mengyao Wu
Roles: Validation, Writing – original draft
Affiliation: College of Big Data, Yunnan Agricultural University, Kunming, Yunnan, China
Zhaowen Li
Roles: Data curation, Writing – original draft
Affiliation: College of Big Data, Yunnan Agricultural University, Kunming, Yunnan, China
Jiawei Zhao
Roles: Data curation
Affiliation: College of Big Data, Yunnan Agricultural University, Kunming, Yunnan, China
Jihong Sun
Roles: Conceptualization, Formal analysis, Investigation, Methodology, Validation, Writing – original draft
E-mail: [email protected]; [email protected]
Affiliations: College of Agronomy and Biotechnology, Yunnan Agricultural University, Kunming, Yunnan, China, Yunnan Key Laboratory of Computer Technologies Application, Kunming University of Science and Technology, Kunming, Yunnan, China
ORICD: https://orcid.org/0009-0000-2058-2675
Ye Qian
Roles: Conceptualization, Formal analysis, Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing
E-mail: [email protected]; [email protected]
Affiliations: College of Big Data, Yunnan Agricultural University, Kunming, Yunnan, China, The Key Laboratory for Crop Production and Smart Agriculture of Yunnan Province, Yunnan Agricultural University, Kunming, Yunnan, China, College of Agronomy and Biotechnology, Yunnan Agricultural University, Kunming, Yunnan, China
ORICD: https://orcid.org/0009-0006-6620-0971
[/RAW_REF_TEXT]
[/RAW_REF_TEXT]
1. Brandler D, Divensi LJ, Tonin RJ, Pilla TP, Rezendes I, Milanesi PM. Evaluation of biological control of fusarium wilt in gerbera with Trichoderma asperellum. OH. 2017;23(3):234.
2. Lim M-J, Han J-E, Murthy HN, Song H-Y, Lee S-Y, Park S-Y. A Comparison of semi-solid, liquid, and temporary immersion bioreactor systems for effective plant regeneration of Gerbera jamesonii “Shy Pink”. Horticulturae. 2024;10(8):836.
3. JNeog PP, Deka P. Screening of Gerbera cultivars against root-knot nematode meloidogyne incognita. J Exp Agric Int. 2024;46(9):899–905.
4. Su X, Nian Y, Shaghaleh H, Hamad A, Yue H, Zhu Y, et al. Combining features selection strategy and features fusion strategy for SPAD estimation of winter wheat based on UAV multispectral imagery. Front Plant Sci. 2024;15:1404238. pmid:38799101
5. Song E, Shao G, Zhu X, Zhang W, Dai Y, Lu J. Estimation of plant height and biomass of rice using unmanned aerial vehicle. Agronomy. 2024;14(1):145.
6. Zhang X, Zhang K, Sun Y, Zhao Y, Zhuang H, Ban W, et al. Combining spectral and texture features of UAS-based multispectral images for maize leaf area index estimation. Remote Sensing. 2022;14(2):331.
7. Guo Y, Chen S, Li X, Cunha M, Jayavelu S, Cammarano D, et al. Machine learning-based approaches for predicting SPAD values of maize using multi-spectral images. Remote Sensing. 2022;14(6):1337.
8. Bajocco S, Ginaldi F, Savian F, Morelli D, Scaglione M, Fanchini D, et al. On the use of NDVI to estimate LAI in field crops: implementing a conversion equation library. Remote Sensing. 2022;14(15):3554.
9. Liu Y, Feng H, Yue J, Jin X, Fan Y, Chen R, et al. Improved potato AGB estimates based on UAV RGB and hyperspectral images. Computers and Electronics in Agriculture. 2023;214:108260.
10. Mulla DJ, Tahir M, Jungers JM. Comparative simulation of crop productivity, soil moisture and nitrate-N leaching losses for intermediate wheatgrass and maize in Minnesota using the DSSAT model. Front Sustain Food Syst. 2023;7.
11. Wang Z, Shi Y, Hou C, Zhang P. Sensitivity analysis of simulated Lycium barbarum L. yield in the WOFOST model under different climate conditions. Ecological Modelling. 2024;488:110602.
12. Tang R, Supit I, Hutjes R, Zhang F, Wang X, Chen X, et al. Modelling growth of chili pepper (Capsicum annuum L.) with the WOFOST model. Agricultural Systems. 2023;209:103688.
13. Ebrahim A, Suat I, Hadis Y. Performance of WOFOST model for simulating maize growth, leaf area index, biomass, grain yield, yield gap, and soil water under irrigation and rainfed conditions. J Irrigation and Drainage Engineering. 2022; 148(2):05021005.
14. Nargund R, Bhatia VS, Sinha NK, Mohanty M, Jayaraman S, Dang YP, et al. Assessing soybean yield potential and yield gap in different agroecological regions of india using the DSSAT model. Agronomy. 2024;14(9):1929.
15. Sheng S, Lilu L, Zhongwen H, Xu Y. Inversion method of sugarcane plant height based on Sentinel-1A dual-polarization time series data. Transactions of the Chinese society for agricultural machinery. 2021; 53(186-194):1000–298.
16. Ma C, Liu M, Ding F, Li C, Cui Y, Chen W, et al. Wheat growth monitoring and yield estimation based on remote sensing data assimilation into the SAFY crop growth model. Sci Rep. 2022;12(1):5473. pmid:35361910
17. Gupta DK, Prasad R, Kumar P, Mishra VN, Vishwakarma AK, Singh RS, et al. Spatial modeling of SPAD values for different type of crops using LISS-IV satellite imagery. 2015 International Conference on Microwave, Optical and Communication Engineering (ICMOCE). IEEE, 2015:212–215.
18. Tang R, Zhao Y, Lin H. Spatio-temporal variation characteristics of aboveground biomass in the headwater of the yellow river based on machine learning. Remote Sensing. 2021;13(17):3404.
19. Yang G, Zhao C, Liu Q, Huang W, Wang J. Inversion of a radiative transfer model for estimating forest lai from multisource and multiangular optical remote sensing data. IEEE Trans Geosci Remote Sensing. 2011;49(3):988–1000.
20. Zu J, Yang H, Wang J, Cai W, Yang Y. Inversion of winter wheat leaf area index from UAV multispectral images: classical vs. deep learning approaches. Front Plant Sci. 2024;15:1367828. pmid:38550285
21. Ji J, Li N, Cui H, Li Y, Zhao X, Zhang H, et al. Study on Monitoring SPAD Values for Multispatial Spatial Vertical Scales of Summer Maize Based on UAV Multispectral Remote Sensing. Agriculture. 2023; 13(5):1004.
22. Xie J, Wang J, Chen Y, Gao P, Yin H, Chen S, et al. Estimating the SPAD of litchi in the growth period and autumn shoot period based on UAV multi-spectrum. Remote Sensing. 2023;15(24):5767.
23. Sharma P, Leigh L, Chang J, Maimaitijiang M, Caffé M. Above-ground biomass estimation in oats using UAV remote sensing and machine learning. Sensors (Basel). 2022;22(2):601. pmid:35062559
24. Yang Z, Yu Z, Wang X, Yan W, Sun S, Feng M, et al. Estimation of millet aboveground biomass utilizing multi-source UAV image feature fusion. Agronomy. 2024;14(4):701.
25. Yu M, He J, Li W, Zheng H, Wang X, Yao X, et al. Estimation of rice leaf area index utilizing a kalman filter fusion methodology based on multi-spectral data obtained from Unmanned Aerial Vehicles (UAVs). Remote Sensing. 2024;16(12):2073.
26. Zhou L, Nie C, Su T, Xu X, Song Y, Yin D, et al. Evaluating the canopy chlorophyll density of maize at the whole growth stage based on multi-scale UAV image feature fusion and machine learning methods. Agriculture. 2023;13(4):895.
27. Ma R, Zhang N, Zhang X, Bai T, Yuan X, Bao H, et al. Cotton Verticillium wilt monitoring based on UAV multispectral-visible multi-source feature fusion. Computers and Electronics in Agriculture. 2024;217:108628.
28. Wen-jing Z, Zhan-kang F, Shi-yuan D, Zhang P-P, Wen J, Ai-chen W, et al. Multi-feature fusion detection of wheat lodging information based on UAV multispectral images. Spectroscopy and Spectral Analysis. 2024; 44:197-206.
29. Naveed A, Muhammad W, Irshad MJ, Aslam MJ, Manzoor SM, Kausar T, et al. Saliency-based semantic weeds detection and classification using UAV multispectral imaging. IEEE Access. 2023;11:11991–2003.
30. Daughtry CST, Hollinger SE. Costs of measuring leaf area index of corn1. Agronomy Journal. 1984;76(5):836–41.
31. Zheng H, Cheng T, Zhou M, Li D, Yao X, Tian Y, et al. Improved estimation of rice aboveground biomass combining textural and spectral analysis of UAV imagery. Precision Agric. 2018;20(3):611–29.
32. JMao D, Wang Z, Luo L, Ren C. Integrating AVHRR and MODIS data to monitor NDVI changes and their relationships with climatic parameters in Northeast China. Int J Appl Earth Observation and Geoinformation. 2012;18:528–36.
33. Gitelson AA, Kaufman YJ, Merzlyak MN. Use of a green channel in remote sensing of global vegetation from EOS-MODIS. Remote Sensing of Environment. 1996;58(3):289–98.
34. Hydrology UA, MD RSLB, Hydrology UA, MD RSLB, ID UoIGLfEDM, OR UACPCRCP. Remote sensing leaf chlorophyll content using a visible band index. Agronomy Journal. 2011; 103(4):1090–99.
35. Li W, Wang J, Zhang Y, Yin Q, Wang W, Zhou G, et al. Combining texture, color, and vegetation index from unmanned aerial vehicle multispectral images to estimate winter wheat leaf area index during the vegetative growth stage. Remote Sensing. 2023; 15(24).
36. Rondeaux G, Steven M, Baret F. Optimization of soil-adjusted vegetation indices. Remote Sensing of Environment. 1996;55(2):95–107.
37. JMutlag WK, Ali SK, Aydam ZM, Taher BH. Feature extraction methods: a review. J Phys: Conf Ser. 2020;1591(1):012028.
38. Curtis AE, Smith TA, Ziganshin BA, Elefteriades JA. The mystery of the Z-score. Aorta (Stamford). 2016;4(4):124–30. pmid:28097194
39. JWu H, Hayes MJ, Weiss A, Hu Q. An evaluation of the standardized precipitation index, the China‐Z Index and the statistical Z‐Score. Intl J Climatology. 2001;21(6):745–58.
40. Kumar C, Mubvumba P, Huang Y, Dhillon J, Reddy K. Multi-stage corn yield prediction using high-resolution UAV multispectral data and machine learning models. Agronomy. 2023;13(5):1277.
41. Forkuor G, Hounkpatin OKL, Welp G, Thiel M. High resolution mapping of soil properties using remote sensing variables in south-western burkina faso: a comparison of machine learning and multiple linear regression models. PLoS One. 2017;12(1):e0170478. pmid:28114334
42. XU M, WATANACHATURAPORN P, VARSHNEY P, ARORA M. Decision tree regression for soft classification of remote sensing data. Remote Sensing of Environment. 2005;97(3):322–36.
43. Biau G, Scornet E. A random forest guided tour. TEST. 2016;25(2):197–227.
44. eNatekin A, eKnoll A. Gradient Boosting Machines, A Tutorial. Frontiers in Neurorobotics. 2013; 7:21.
45. Chen T, Guestrin C. XGBoost: A scalable tree boosting system. Proceedings of the 22nd acm sigkdd international conference on knowledge discovery and data mining. 2016:785–794.
46. Ben-Hur A, Ong CS, Sonnenburg S, Schölkopf B, Rätsch G. Support vector machines and kernels for computational biology. PLoS Comput Biol. 2008;4(10):e1000173. pmid:18974822
47. Khaing KT. Enhanced features ranking and selection using recursive feature elimination(rfe) and k-nearest neighbor algorithms in support vector machine for intrusion detection system. International J Network and Mobile Technologies. 2010; 1(1):8.
48. Yan K, Zhang D. Feature selection and analysis on correlated gas sensor data with recursive feature elimination. Sensors and Actuators B: Chemical. 2015;212:353–63.
49. Wei Q, Zhang B-Z, Wei Z, Han X, Duan C-F. Estimation of canopy chlorophyll content in winter wheat by UAV multispectral remote sensing. CABI Digital Library. 2020:365–372.
50. Guomin S, Yajie W, Wenting H. Estimation method of leaf area index for summer maize using UAV-based multispectral remote sensing. Smart Agriculture. 2020; 2(3):118.
51. Guo Y, He J, Zhang H, Shi Z, Wei P, Jing Y, et al. Improvement of winter wheat aboveground biomass estimation using digital surface model information extracted from unmanned-aerial-vehicle-based multispectral images. Agriculture. 2024;14(3):378.
© 2025 Wang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License: http://creativecommons.org/licenses/by/4.0/ (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.