DeconSeis: Deconvolve discrete instrument response from many traces In TDD: Time-Domain Deconvolution of Seismometer Response

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.

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) ```

TDD documentation built on May 29, 2017, 5:53 p.m.