DOsgram: Gabor transform

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/DOsgram.R

Description

Gabor transform with simple spectrum

Usage

1
2
DOsgram(Xamp, DT = 0.008, multi = 1, scale.def = 0,
 TWIN = 2, TSKIP = 0.2, PCTTAP = 0.05, PLOT=TRUE)

Arguments

Xamp

signal

DT

sample rate interval (s)

multi

Multiples of time window estimate

scale.def

scaling flag for plotting (0=raw, 1=log, 2=sqrt)

TWIN

time for window

TSKIP

time for skip

PCTTAP

percent of taper to apply to individual windows

PLOT

logical, TRUE=plot to device

Details

This is a non-interactive version of SPECT.drive.

Value

list

sig

input signal

dt

deltat

numfreqs

Number of frequencies output

wpars

input parameters list(Nfft=numfreqs, Ns=Ns, Nov=Nov, fl=fl, fh=fh)

DSPEC

spectrum image

freqs

output frequencies (y axis)

tims

output times (x-axis)

Author(s)

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

References

Lees, J. M. and Park, J., 1995: Multiple-taper spectral analysis: A stand-alone C-subroutine, Computers \& Geology, 21(2), 199-236.

See Also

evolMTM, evolfft, evolAR, plotevol

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
data(KH)
###   swig(KH)

Xamp  <-  KH$JSTR[[1]]
Xamp  <-  Xamp[57914:72989]

Nfft <- 1024   ###  fft length
Ns <- 512      ###  number of samples in a window
Nov <- 480    ###  number of samples of overlap per window
fl <- 0        ###  low frequency to return
fh <- 12     ###  high frequency to return


EV  <-  DOsgram(Xamp, DT = 0.008, multi = 1, scale.def = 0,
TWIN = 2, TSKIP = 0.2, PCTTAP = 0.05)

RSEIS documentation built on Aug. 2, 2021, 5:07 p.m.