TDD (version 0.4)

DeconSeis: Deconvolve discrete instrument response from many traces

Description

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

Usage

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

Value

GH, with the instrument response removed from every trace.

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.

Examples

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)