plotevol: Plot Spectrogram

In RSEIS: Seismic Time Series Analysis Tools

Plot Spectrogram

Description

Plot Spectrogram

Usage

plotevol(DEVOL, log = 0, fl = 0, fh = 10, col = col, ylog = FALSE, ygrid
= FALSE, AXE = c(1, 2, 3, 4), CSCALE = FALSE, WUNITS = "Volts", STAMP =
NULL, STYLE = "fft")

plotevol2(DEVOL, log = 0, fl = 0, fh = 10, col = col, ylog = FALSE, ygrid
= FALSE, AXE = c(1, 2, 3, 4), CSCALE = FALSE, WUNITS = "Volts", STAMP =
NULL, STYLE = "fft", add=FALSE, IMAGE=TRUE, WIG=TRUE )

blankevol(DEVOL, log=0,  fl=0, fh=10 , col=col, ylog=FALSE, ygrid=FALSE,
AXE=c(1,2,3,4),
CSCALE=FALSE, WUNITS="Volts", STAMP=NULL, STYLE="fft", WIG=TRUE  )

Arguments

DEVOL	spectrogram structure
log	scale by logarithm
fl	low frequency
fh	high frequency
col	color palette
ylog	scale Y-axis by log
ygrid	logical, TRUE=add grid
AXE	sides to add axis
CSCALE	logical, TRUE=add color scale
WUNITS	character string for units
STAMP	character string for identification
STYLE	Plotting style. Default, "fft"=plot half the spectrum image , else plot whole spectrum
add	logical, add to existing plot, default=FALSE
IMAGE	logical, whether to plot the image or not
WIG	logical, whether to plot the wiggle or not

Details

Plot Spectrogram. Because the fft function returns positive and negative frequencies, ff STYLE="fft" then the image matrix is reduced IMAT = t(DSPEC[1:(numfreqs/2),]) otherwise IMAT = t(DSPEC).

plotevol2 is used to add secondary spectra to ones already plotted, or to manage graphical paramters, or create other plots that match the graphical presentation of the spectrogram (plots of frequency versus time, but not images)

Value

Graphical Side Effects

Author(s)

Jonathan M. Lees<jonathan.lees.edu>

See Also

evolfft

Examples

data(CE1)

Xamp <- CE1$y

DT <- CE1$dt

 tsecs <- DT*(length(Xamp)*.02)
    multi <- 2
scale.def <- 1
  TWOSEC <- tsecs*(1/DT)
    
  NS <- floor(multi*TWOSEC)
  NOV <- floor(multi*(TWOSEC-.2*TWOSEC))

    
    Nfft<-4096

    pal <- rainbow(100)

    fl  <- 0
   fh <- 1/(2*DT)

   flshow <-  .5
   fhshow <- 120

   DEV <- evolfft(Xamp,DT , Nfft=Nfft, Ns=NS , Nov=NOV,  fl=fl, fh=fh  )

  PE <-  plotevol(DEV, log=scale.def, fl=flshow, fh=fhshow,
                 col=pal, ygrid=FALSE, STAMP="HITHERE", STYLE="fft")

RSEIS documentation built on Aug. 19, 2023, 5:07 p.m.