Related Experiment Video
Updated: Aug 29, 2025

Cereal Crop Ear Counting in Field Conditions Using Zenithal RGB Images
Published on: February 2, 2019
Optimizing Gaussian Process Regression for Image Time Series Gap-Filling and Crop Monitoring
Santiago Belda1, Luca Pipia1, Pablo Morcillo-Pallarés1
1Image Processing Laboratory (IPL), Parc Científic, University of Valencia, Paterna, 46980 Valencia, Spain.
This study introduces a faster method for processing satellite image time series using Gaussian process regression (GPR) for crop phenology monitoring. The new approach significantly reduces computation time with minimal impact on accuracy for essential agricultural insights.
Area of Science:
- Remote Sensing and Geospatial Analysis
- Agricultural Science and Crop Monitoring
- Artificial Intelligence and Machine Learning Applications
Background:
- Satellite image time series analysis is crucial for monitoring crop phenology, including start and end of season estimations.
- Gaussian Process Regression (GPR) is a powerful Bayesian algorithm for time series gap-filling, offering uncertainty estimates.
- Standard per-pixel GPR processing is computationally intensive, hindering large-scale application.
Purpose of the Study:
- To develop a computationally efficient GPR approach for satellite image time series gap-filling.
- To assess the accuracy and performance of the proposed method compared to standard per-pixel GPR.
- To evaluate the impact of the optimized GPR on crop phenology indicator retrieval.
Main Methods:
- Proposed a novel method involving cropland-based precalculations of GPR hyperparameters (θ) to replace per-pixel optimization.
- Utilized Sentinel-2 Leaf Area Index (LAI) time series data from an agricultural region in Spain for validation.
- Compared image reconstruction accuracy (RMSE) and processing time against standard per-pixel GPR and a scene-wide hyperparameter approach.
Main Results:
- The proposed method achieved comparable accuracy to per-pixel GPR (12% RMSE degradation) while accelerating processing time by approximately 90 times.
- Using scene-wide hyperparameters yielded similar accuracies and computational performance for crop areas.
- Calculated crop phenology indicators (start and end of season) showed analogous temporal patterns across all tested methods, with differences within five days.
Conclusions:
- The proposed GPR hyperparameter precalculation significantly enhances computational efficiency for satellite time series analysis without substantial accuracy loss.
- This optimized approach is suitable for large-scale crop monitoring applications, providing reliable phenological insights.
- The developed gap-filling and phenology retrieval techniques are integrated into the freely available DATimeS GUI toolbox.
Related Concept Videos
Linear Approximation in Time Domain
For a simple pendulum with a mass evenly distributed along its length and the center of mass located at half the pendulum's length,...
Calibration Curves: Linear Least Squares
For data that follow a straight line, the standard method for fitting is the linear...
Improving Translational Accuracy
Linear Approximation in Frequency Domain
In contrast, nonlinear systems do not inherently possess these properties. However, for small deviations around an operating point, a nonlinear system can often be approximated as linear....
Residuals and Least-Squares Property
If the observed data point lies above the line, the residual is positive, and the line underestimates the actual data value for y. If the observed data point lies below the line, the residual is negative, and the line overestimates the actual data value for y.
The process of fitting the best-fit...
Reducing Line Loss
With a step-up transformer at the source, the voltage is increased, thereby reducing the current in the transmission lines since power loss...

