R: Velocity Seismogram

VELOCITY.SEISN {RSEIS} | R Documentation

## Velocity Seismogram

### Description

Removes seismic instrument response and corrects for sensitivity of seismoc instrument, returning units of m/s rather than volts.

### Usage

    VELOCITY.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"))

### 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

### Details

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

### Value

Same as input list with new traces representing velocity versus volts

### Author(s)

Jonathan M. Lees<jonathan.lees@unc.edu>

### See Also

DISPLACE.SEISN, deconinst

### Examples

    ## Not run:
    Kal  <-  PreSet.Instr()
    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

    VH  <-     VELOCITY.SEISN(TH, sel = 1:length(TH$JSTR), inst = 1,
    Kal = Kal, FILT = list(ON = FALSE, fl = 1/30, fh = 7,
    type = "HP", proto = "BU"))

    ## End(Not run)

[Package RSEIS version 3.7-4 Index]