deconv {astrochron}

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.