# DOsgram: Gabor transform In RSEIS: Seismic Time Series Analysis Tools

## 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.

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