Early monitoring of within-field yield variability and forecasting yield potential is critical for farmers and other key stakeholders such as policymakers. Remote sensing techniques are progressively being used in yield prediction studies due to easy access and affordability. Despite the increasing use of remote sensing techniques for yield prediction in agriculture, there is still a need for medium-resolution satellite imagery when predicting canola yield using a combination of crop and soil information. In this study, we investigated the utility of remotely sensed flowering information from PlanetScope (at 4 m) satellite imagery combined with derived soil and topography parameters to predict canola yield. Our yield prediction model was trained and validated using data from 21 fields managed under variable rate seed and fertilizer application, including cleaned harvester yield maps, soil, and topography maps. To quantify the flowering intensity of canola, 9 vegetation indices (VIs) were calculated using spectral bands from PlanetScope imagery acquired for the reproductive stages of canola. We created five random forest regression models using different subsets of covariates, including VIs, soil, and topography features, to predict canola yield within the season. Using a random forest regression algorithm, we recorded accuracies ranging from poor to best performing using coefficient of determination and root mean squared error (R2: 0.47 to 0.66, RMSE: 325 to 399 kg ha−1). The optimal subset of covariates identified electrical conductivity (EC), Normalized Difference Yellowness Index, and Canola Index as the key variables explaining within-spatial variability in canola yield. Our final model exhibited a validation R2 of 0.46 (RMSE = 730 kg ha−1), demonstrating the potential of medium-resolution satellite imagery during the flowering stage to detect and quantify sub-field spatial and temporal floral phenology changes when predicting canola yield.
