hilbert {astrochron}

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)