psdPlot: Generate plots from a set of PSDs

In IRISSeismic: Classes and Methods for Seismic Data Analysis

Generate plots from a set of PSDs

Description

The psdPlot function is used to generate plots from the data in a psdList or psdDF dataframe.

Usage

psdPlot(PSDs, style='psd', evalresp=NULL, ylo=-200, yhi=-50, showNoiseModel=TRUE, 
        showMaxMin=TRUE, showMode=TRUE, showMean=FALSE, showMedian=FALSE, ...)

Arguments

PSDs	either a list as returned by psdList or a dataframe of PSD values obtained from the BSS
style	character identifier of plot type: 'psd' plots PSD lines, 'pdf' plots the pdfMatrix
evalresp	dataframe of freq, amp, phase information matching output of getEvalresp, optional
ylo	numeric setting lower limit of plot y-axis (default=-200)
yhi	numeric setting upper limit of plot y-axis (default=-50)
showNoiseModel	logical controlling plotting of noise model lines (default=TRUE)
showMaxMin	logical controlling plotting of PSD max and min lines (default=TRUE)
showMode	logical controlling plotting of PDF mode line (default=TRUE)
showMean	logical controlling plotting of PSD mean line (default=FALSE)
showMedian	logical controlling plotting of PSD median line (default=FALSE)
...	arguments to be passed to plotting methods

Details

The psdPlot function creates visualizations for sets of PSDs. Plots generated with style='pdf' mimic the plots presented in the McNamara paper.

Author(s)

Jonathan Callahan [e-mail]

References

Seismic Noise Analysis System Using Power Spectral Density Probability Density Functions (McNamara and Boaz 2005)

See Also

McNamaraPSD, psdList, psdStatistics

Examples

## Not run: 
# Create a new IrisClient
iris <- new("IrisClient", debug=TRUE)

# Get seismic data
starttime <- as.POSIXct("2011-05-05", tz="GMT") # 2011.125
endtime <- starttime + 1*24*3600
st <- getDataselect(iris,"IU","GRFO","--","BHE",starttime,endtime)

# Generate power spectral density for each hour long segment
psdList <- psdList(st)

# 'psd' line plot
psdPlot(psdList,style='psd',type='l',col=adjustcolor('black',0.3))

# McNamara 'pdf' plot
psdPlot(psdList,style='pdf')

## End(Not run)

IRISSeismic documentation built on Oct. 16, 2022, 1:09 a.m.