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)
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(RSEIS)
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/RSEIS/plotevol.Rd_%03d_medium.png", width=480, height=480)
> ### Name: plotevol
> ### Title: Plot Spectrogram
> ### Aliases: plotevol blankevol plotevol2
> ### Keywords: hplot
>
> ### ** 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 )
evolfft kcol= 117 krow= 2048 Ns 29 Nov 23
>
> PE <- plotevol(DEV, log=scale.def, fl=flshow, fh=fhshow,
+ col=pal, ygrid=FALSE, STAMP="HITHERE", STYLE="fft")
>
>
>
>
>
>
>
>
>
> dev.off()
null device
1
>