PSTLTcurve: Short Term/Long Term curve

In RSEIS: Seismic Time Series Analysis Tools

Short Term/Long Term curve

Description

ST/LT ratio curve for sutomated picking routines

Usage

PSTLTcurve(y, dt = 0.008, fwlen = 125, bwlen = 125,
 perc = 0.05, stretch = 1000, MED = 255, PLOT = FALSE)

Arguments

y	signal
dt	deltaT (s)
fwlen	forward window
bwlen	backward window
perc	percent cut-off
stretch	stretch curve
MED	Median smoothing parameter
PLOT	logical, TRUE=PLOT

Value

list(flag=1, ind=ix, eye=eye, mix=mix, SNR=SNR, s2=s2, rat=therat)

flag	flag on success
ind	index of pick estimate 1
eye	index of pick estimate 2
mix	index of pick estimate 3
SNR	Signal/Noise ratio
s2	sum squared
rat	ratio curve

Author(s)

Jonathan M. Lees<jonathan.lees.edu>

Examples

require(stats)

data(CE1)
plot(CE1$x, CE1$y, type='l')

z <- CE1$y[ CE1$x>5.352622 & CE1$x<5.589836]
x <- CE1$x[ CE1$x>5.352622 & CE1$x<5.589836]

G <- PSTLTcurve(z, dt = CE1$dt, fwlen = 10,
bwlen = 10, perc = 0.05,
stretch = 10, MED = 11, PLOT = FALSE)

### get time from beginning of trace
tpick <- x[G$ind]
abline(v=x[G$ind], col='red', lty=2)

RSEIS documentation built on Sept. 13, 2024, 1:09 a.m.