“…Within more heterogeneous environments such as the Frasassi gorge, the boxplot function performs better for certain indices when compared to the tsclean function. In literature, different techniques have been used to carry out the smoothing and correction phase of time-series satellite data, such as the curve fitting (Pickers and Manning, 2015), the Fourier decomposition (Mingwei et al, 2008), the asymmetric Gaussian function (Jonsson and Eklundh, 2002), the double logistic functions (Atkinson et al, 2012, Eklundh andJönsson, 2015), the Whittaker smoother (Shao et al, 2016, Kandasamy et al, 2013, the Savitzky-Golay filter (Huang et al, 2021), the high order spline with roughness damping (Hermance et al, 2007), the spatio-temporal tensor completion method (Chu et al, 2021) and other spatio-temporal combination methods such as the adaptive spatio-temporal weighted method (Li et al, 2017) and hybrid Generalised Additive Model (GAM)geostatistical space-time model (Poggio et al, 2012) wich are even useful to fill temporal gaps. GAM utilization for regression model fitting widely demonstrated in literature (Hua et al, 2021).…”