Learn R Programming

RSEIS (version 2.8-1)

plotevol: 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

Value

  • Graphical Side Effects

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)

See Also

evolfft

Examples

Run this code
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")

Run the code above in your browser using DataLab