Apply a 1D generic function to a time series
and specific methods for common tasks,
such as missing values removal and smoothing.
`sits_apply()` returns a sits tibble with the same samples
and new bands computed by `fun`, `fun_index` functions.
These functions must be defined inline; they are called by `sits_apply`
for each band, whose vector values is passed as the function argument.
The `fun` function may either return a vector or a list of vectors.
In the first case, the vector will be the new values
of the corresponding band.
In the second case, the returned list must have names,
and each element vector will generate a new band which name composed
by concatenating original band name and the corresponding list element name.
If a suffix is provided in `bands_suffix`, all resulting band
names will end with provided suffix separated by a ".".