RSEIS (version 3.5-2)

PSTLTcurve: 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

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)