Learn R Programming

astrochron (version 0.4.3)

hilbert: Hilbert Transform of Stratigraphic Series

Description

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

Usage

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

Arguments

dat
Stratigraphic series to Hilbert Transform. First column should be location (e.g., depth), second column should be data value.
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)
verbose
Verbose output? (T or F)

Examples

Run this code
# 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)

Run the code above in your browser using DataLab