Displacement seismogram
Removes seismic instrument response and integrates to displacement.
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))
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 |
Instrument responses are lists of poles and zeros for each instrument defined.
Same as input list with new traces representing displacement versus velocity
Jonathan M. Lees<jonathan.lees@unc.edu>
VELOCITY.SEISN, deconinst, butfilt
## 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)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.