Unmanned aerial vehicle (UAV) multispectral imagery has been applied in the remote sensing of wheat SPAD (Soil and Plant Analyzer Development) values. However, existing research has yet to consider the influence of different growth stages and UAV flight altitudes on the accuracy of SPAD estimation. This study aims to optimize UAV flight strategies and incorporate multiple feature selection techniques and machine learning algorithms to enhance the accuracy of the SPAD value estimation of different wheat varieties across growth stages. This study sets two flight altitudes (20 and 40 m). Multispectral images were collected for four winter wheat varieties during the green-up and jointing stages. Three feature selection methods (Pearson, recursive feature elimination (RFE), and correlation-based feature selection (CFS)) and four machine learning regression models (elastic net, random forest (RF), backpropagation neural network (BPNN), and extreme gradient boosting (XGBoost)) were combined to construct SPAD value estimation models for individual growth stages as well as across growth stages. The CFS-RF (40 m) model achieved satisfactory results (green-up stage: R² = 0.7270, RPD = 2.0672, RMSE = 1.1835, RRMSE = 0.0259; jointing stage: R² = 0.8092, RPD = 2.3698, RMSE = 2.3650, RRMSE = 0.0487). For cross-growth stage modeling, the optimal prediction results for SPAD values were achieved at a flight altitude of 40 m using the Pearson-XGBoost model (R² = 0.8069, RPD = 2.3135, RMSE = 2.0911, RRMSE = 0.0442). These demonstrate that the flight altitude of UAVs significantly impacts the estimation accuracy, and the flight altitude of 40 m (with a spatial resolution of 2.12 cm) achieves better SPAD value estimation than that of 20 m (with a spatial resolution of 1.06 cm). This study also showed that the optimal combination of feature selection methods and machine learning algorithms can more accurately estimate winter wheat SPAD values. In addition, this study includes multiple winter wheat varieties, enhancing the generalizability of the research results and facilitating future real-time and rapid monitoring of winter wheat growth.

Leaf Area Index (LAI) is a fundamental indicator of plant growth status in agronomy and environmental research. With the rapid development of drone technology, the estimation of crop LAI based on drone imagery and vegetation indices is becoming increasingly popular. However, there is still a lack of detailed research on the feasibility of using image texture to estimate LAI and the impact of combining texture indices with vegetation indices on LAI estimation accuracy. In this study, two key growth stages of winter wheat (i.e., the stages of green-up and jointing) were selected, and LAI was calculated using digital hemispherical photography. The feasibility of predicting winter wheat LAI was explored under three conditions: vegetation index, texture index, and a combination of vegetation index and texture index, at flight heights of 20 m and 40 m. Two feature selection methods (Lasso and recursive feature elimination) were combined with four machine learning regression models (multiple linear regression, random forest, support vector machine, and backpropagation neural network). The results showed that during the vegetative growth stage of winter wheat, the model combining texture information with vegetation indices performed better than the models using vegetation indices alone or texture information alone. Among them, the best prediction result based on vegetation index was RFECV-MLR at a flight height of 40 m (R² = 0.8943, RMSE = 0.4139, RRMSE = 0.1304, RPD = 3.0763); the best prediction result based on texture index was RFECV-RF at a flight height of 40 m (R² = 0.8894, RMSE = 0.4236, RRMSE = 0.1335, RPD = 3.0063); and the best prediction result combining texture and index was RFECV-RF at a flight height of 40 m (R² = 0.9210, RMSE = 0.3579, RRMSE = 0.1128, RPD = 3.5575). The results of this study demonstrate that combining vegetation indices and texture from multispectral drone imagery can improve the accuracy of LAI estimation during the vegetative growth stage of winter wheat. In addition, selecting a flight height of 40 m can improve efficiency in large-scale agricultural field monitoring, as this study showed that drone data at flight heights of 20 m and 40 m did not significantly affect model accuracy.

In China’s second-largest wheat-producing region, the mid-lower Yangtze River area, cold stress impacts winter wheat production during the pre-heading growth stage. Previous research focused on specific growth stages, lacking a comprehensive approach. This study utilizes Unmanned Aerial Vehicle (UAV) multispectral imagery to monitor Soil-Plant Analysis Development (SPAD) values throughout the pre-heading stage, assessing crop stress resilience. Vegetation Indices (VIs) and Texture Indices (TIs) are extracted from UAV imagery. Recursive Feature Elimination (RFE) is applied to VIs, TIs, and fused variables (VIs + TIs), and six machine learning algorithms are employed for SPAD value estimation. The fused VIs and TIs model, based on Long Short-Term Memory (LSTM), achieves the highest accuracy (R² = 0.8576, RMSE = 2.9352, RRMSE = 0.0644, RPD = 2.6677), demonstrating robust generalization across wheat varieties and nitrogen management practices. This research aids in mitigating winter wheat frost risks and increasing yields.

Monitoring winter wheat Soil–Plant Analysis Development (SPAD) values using Unmanned Aerial Vehicles (UAVs) is an effective and non-destructive method. However, predicting SPAD values during the booting stage is less accurate than other growth stages. Existing research on UAV-based SPAD value prediction has mainly focused on low-altitude flights of 10–30 m, neglecting the potential benefits of higher-altitude flights. The study evaluates predictions of winter wheat SPAD values during the booting stage using Vegetation Indices (VIs) from UAV images at five different altitudes (i.e., 20, 40, 60, 80, 100, and 120 m, respectively, using a DJI P4-Multispectral UAV as an example, with a resolution from 1.06 to 6.35 cm/pixel). Additionally, we compare the predictive performance using various predictor variables (VIs, Texture Indices (TIs), Discrete Wavelet Transform (DWT)) individually and in combination. Four machine learning algorithms (Ridge, Random Forest, Support Vector Regression, and Back Propagation Neural Network) are employed. The results demonstrate a comparable prediction performance between using UAV images at 120 m (with a resolution of 6.35 cm/pixel) and using the images at 20 m (with a resolution of 1.06 cm/pixel). This finding significantly improves the efficiency of UAV monitoring since flying UAVs at higher altitudes results in greater coverage, thus reducing the time needed for scouting when using the same heading overlap and side overlap rates. The overall trend in prediction accuracy is as follows: VIs + TIs + DWT > VIs + TIs > VIs + DWT > TIs + DWT > TIs > VIs > DWT. The VIs + TIs + DWT set obtains frequency information (DWT), compensating for the limitations of the VIs + TIs set. This study enhances the effectiveness of using UAVs in agricultural research and practices.

Accurately estimating leaves’ relative chlorophyll contents (widely represented by Soil and Plant Analysis Development (SPAD) values) across growth stages is crucial for assessing crop health, particularly in regions characterized by varying sowing dates. Unlike previous studies focusing on high-resolution UAV imagery or specific growth stages, this research incorporates satellite-derived texture indices (TIs) into a SPAD value estimation model applicable across multiple growth stages (from tillering to grain-filling). Field experiments were conducted in Jiangsu Province, China, where winter wheat sowing dates varied significantly from field to field. Sentinel-2 imagery was employed to extract vegetation indices (VIs) and TIs. Following a two-step variable selection method, Random Forest (RF)-LassoCV, five machine learning algorithms were applied to develop estimation models. The newly developed model (SVR-RBF_VIs+TIs) exhibited robust estimation performance (R² = 0.8131, RMSE = 3.2333, RRMSE = 0.0710, and RPD = 2.3424) when validated against independent SPAD value datasets collected from fields with varying sowing dates. Moreover, this optimal model also exhibited a notable level of transferability at another location with different sowing times, wheat varieties, and soil types from the modeling area. In addition, this research revealed that despite the lower resolution of satellite imagery compared to UAV imagery, the incorporation of TIs significantly improved estimation accuracies compared to the sole use of VIs typical in previous studies.