hilbert: Hilbert transform of stratigraphic series

In astrochron: A Computational Tool for Astrochronology

Hilbert transform of stratigraphic series

Description

Calculate instantaneous amplitude (envelope) and phase via Hilbert Transform of stratigraphic series

Usage

hilbert(dat,phase=F,padfac=2,demean=T,detrend=F,output=T,addmean=F,genplot=T,check=T,
        verbose=T)

Arguments

dat	Stratigraphic series to Hilbert Transform. First column should be location (e.g., depth), second column should be data value.
phase	Calculate instantaneous phase? (T or F)
padfac	Pad with zeros to (padfac*npts) points, where npts is the original number of data points.
demean	Remove mean from data series? (T or F)
detrend	Remove linear trend from data series? (T or F)
output	Return results as new data frame? (T or F)
addmean	Add mean value to instantaneous amplitude? (T or F)
genplot	Generate summary plots? (T or F)
check	Conduct compliance checks before processing? (T or F) In general this should be activated; the option is included for Monte Carlo simulation.
verbose	Verbose output? (T or F)

Examples

# generate example series with 3 precession terms and noise
ex <- cycles(noisevar=.0004,dt=5)
# bandpass precession terms using cosine-tapered window 
res_ex <- bandpass(ex,flow=0.038,fhigh=0.057,win=2,p=.4)
# hilbert transform
hil_ex <- hilbert(res_ex)

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