Deconvolve discrete instrument response from many traces

Share:

Description

Deconvolves instrument responses from a seismogram structure from the RSEIS package.

Usage

1
2
DeconSeis(GH, inst, L, fl = 0.1, fh = NaN, bitweight = NULL, dec =
rep(1, length(GH$JSTR)))

Arguments

GH

Seismogram structure

inst

Vector of indices of instrument responses within L to deconvolve from each trace in GH

L

List in which each element is a discrete instrument response (from MakeDPZ, for example)

fl

Low corner of filter (NaN for no high-pass filtering) (Hz)

fh

High corner of filter (NaN for no low-pass filtering) (Hz)

bitweight

Vector of optional counts-to-volts factors for data in counts (volts/counts)–NULL if data are already in volts

dec

Oversampling/decimation factor (optional); vector equal to the number of traces considered

Details

Discrete instrument responses are specific to a given sampling rate. If the response you give has a different sample rate (given by DPZ$dt) from the trace x, you will get incorrect results. DPZ$dt * dec should be equal to the sample intervals of the traces.

Value

GH, with the instrument response removed from every trace.

Author(s)

Jake Anderson

Examples

1
2
3
4
5
6
7
data(COLOC)
swig(COLOC)
L = GetDPZ(c(4, 14), c(0.01, 0.01)) # get responses for 3T and 40T-1s
inst = c(1,1,1,2,2,2) # deconvolve 3T response from channel 1-3,
                      # 40T-1 response from channel 4-6
D = DeconSeis(COLOC, inst, L)
swig(D)