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

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

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)

Example output

Loading required package: signal

Attaching package: 'signal'

The following objects are masked from 'package:stats':

    filter, poly

Loading required package: RSEIS
Loading required package: pracma

Attaching package: 'pracma'

The following objects are masked from 'package:RSEIS':

    detrend, hypot, logspace, peaks, trapz

The following objects are masked from 'package:signal':

    conv, ifft, interp1, pchip, polyval, roots


Attaching package: 'TDD'

The following object is masked from 'package:RSEIS':

    ReadInstr

$but
[1] "None"

$zloc
[1] 0

$pix
[1] 0

$but
[1] "None"

$zloc
[1] 0

$pix
[1] 0

TDD documentation built on May 2, 2019, 4:51 a.m.