DISPLACE.SEISN: Displacement seismogram

In RSEIS: Seismic Time Series Analysis Tools

Displacement seismogram

Description

Removes seismic instrument response and integrates to displacement.

Usage

DISPLACE.SEISN(TH, sel = 1:length(TH$JSTR), inst = 1,
Kal = Kal,waterlevel = 1e-08,  FILT = list(ON = FALSE,
 fl = 1/30, fh = 7, type = "HP", proto = "BU",RM=FALSE, zp=TRUE))

Arguments

TH	list structure of seismic traces
sel	select which tracesin list to deconvolve
inst	index to instrument in Kal list for calibration and instrument response
Kal	list of instrument responses
waterlevel	waterlevel for low frequency division
FILT	filter output, after instrumentation, see butfilt

Details

Instrument responses are lists of poles and zeros for each instrument defined.

Value

Same as input list with new traces representing displacement versus velocity

Author(s)

Jonathan M. Lees<[e-mail]>

See Also

VELOCITY.SEISN, deconinst, butfilt

Examples

## Not run: 
Majorevent   <-  250+8/24+2250/(24*3600)

at1  <-  Majorevent - 1/24
at2  <-  Majorevent + 1/24


TH  <-   Mine.seis(at1, at2, DB99, NULL, NULL , kind = 1,
 Iendian=1, BIGLONG=FALSE, CHOP=TRUE, verbose=FALSE)
inst  <-  rep(0, length(TH$STNS))
inst[TH$COMPS %in% c("V", "N", "E")]  <-  1
inst[TH$STNS=="krm2" &  TH$COMPS %in% c("V", "N", "E")]  <-  2
inst[TH$STNS=="krm0" &  TH$COMPS %in% c("V", "N", "E")]  <-  2


DH  <-   DISPLACE.SEISN(TH, sel = which(inst>0) , inst = inst,
Kal = Kal, FILT = list(ON = FALSE, fl = 1/200, fh = 7,
type = "BP", proto = "BU"))

  SOUT  <-   swig(DH, PADDLAB=c("CENTER", "fspread", "HALF", "PREV"),
filters=thefilts  )



## End(Not run)

RSEIS documentation built on Aug. 19, 2023, 5:07 p.m.