Wiener Deconvolution through Fast Fourier Transform.
deconv(x, y, index, dt, pt = 0.2, wiener = TRUE)z deconvolved/un-bioturbed series. length(z) = length(x)
Time series that needs to be deconvolved.
Green function/Impulse response function.
index in the impulse response function that corresponds to the deposition time.
time resolution for the input series.
Cosine-tapered window parameter: pt is the percent of the data series tapered (choose 0-1). When pt=1, this is equivalent to a Hann taper.
Apply Wiener filter? (T or F)
function 'deconv' is used by function 'unbioturb'. A cosine taper is applied to remove edge effects. The signal-to-noise ratio is chosen to be 0.05, and gamma is also chosen to be 0.05. x and y do not need to be of the same length. For additional information see Liu et al. (2021)
Liu, H., Meyers, S.R., and Marcott, S.A., 2021, Unmixing dee-sea paleoclimate records: A study on bioturbation effects through convolution and deconvolution, Earth and Planetary Science Letters, 564, 116883.