deconv {astrochron} | R Documentation |
Wiener Deconvolution through Fast Fourier Transform
Description
Wiener Deconvolution through Fast Fourier Transform.
Usage
deconv(x, y, index, dt, pt = 0.2, wiener = TRUE)
Arguments
x |
Time series that needs to be deconvolved. |
y |
Green function/Impulse response function. |
index |
index in the impulse response function that corresponds to the deposition time. |
dt |
time resolution for the input series. |
pt |
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. |
wiener |
Apply Wiener filter? (T or F) |
Details
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)
Value
z deconvolved/un-bioturbed series. length(z) = length(x)
References
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.