Multi-Temporal and Spectral Analysis of High-Resolution Hyperspectral Airborne Imagery for Precision Agriculture: Assessment of Wheat Grain Yield and Grain Protein Content

This study evaluates the potential of high resolution hyperspectral airborne imagery to capture within-field variability of durum wheat grain yield (GY) and grain protein content (GPC) in two commercial fields in the Yaqui Valley (northwestern Mexico). Through a weekly/biweekly airborne flight campaign, we acquired 10 mosaics with a micro-hyperspectral Vis-NIR imaging sensor ranging from 400–850 nanometres (nm). Just before harvest, 114 georeferenced grain samples were obtained manually. Using spectral exploratory analysis, we calculated narrow-band physiological spectral indices—normalized difference spectral index (NDSI) and ratio spectral index (RSI)—from every single hyperspectral mosaic using complete two by two combinations of wavelengths. We applied two methods for the multi-temporal hyperspectral exploratory analysis: (a) Temporal Principal Component Analysis (tPCA) on wavelengths across all images and (b) the integration of vegetation indices over time based on area under the curve (AUC) calculations. For GY, the best R2 (0.32) were found using both the spectral (NDSI—Ri, 750 to 840 nm and Rj, ±720–736 nm) and the multi-temporal AUC exploratory analysis (EVI and OSAVI through AUC) methods. For GPC, all exploratory analysis methods tested revealed (a) a low to very low coefficient of determination (R2 <= 0.21), (b) a relatively low overall prediction error (RMSE: 0.45–0.49%), compared to results from other literature studies, and (c) that the spectral exploratory analysis approach is slightly better than the multi-temporal approaches, with early season NDSI of 700 with 574 nm and late season NDSI of 707 with 523 nm as the best indicators. Using residual maps from the regression analyses of NDSIs and GPC, we visualized GPC within-field variability and showed that up to 75% of the field area could be mapped with relatively good predictability (residual class: -0.25 to 0.25%), therefore showing the potential of remote sensing imagery to capture the within-field variation of GPC under conventional agricultural practices