DOsgram {RSEIS}R Documentation

Gabor transform

Description

Gabor transform with simple spectrum

Usage

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




data(KH)
###   PICK.GEN(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)





[Package RSEIS version 2.3-3 Index]