deconv: Wiener Deconvolution through Fast Fourier Transform

In astrochron: A Computational Tool for Astrochronology

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.

astrochron documentation built on Sept. 30, 2024, 9:14 a.m.