plot_ppsd: Plot a probabilistic power spectral density estimate (PPSD)

In coffeemuggler/eseis: Environmental Seismology Toolbox

Plot a probabilistic power spectral density estimate (PPSD)

Description

The function uses the output of signal_spectrogram() to plot a probabilistic power spectral density estimate.

Usage

plot_ppsd(data, res = c(500, 500), n, ...)

Arguments

data	List object, spectrogram to be plotted. Must be output of signal_spectrogram() or of equivalent structure.
res	Integer vector of length two, factors of image resolution in pixels, i.e. in time and frequency dimension. Default is c(100, 100).
n	Integer vector of length two, factors by which the image will be smoothend by a running average. n sets the filter window size, in x and y direction, respectively. By default, the window sizes are set to one percent of the input data set dimension.
...	Additional arguments passed to the plot function.

Value

Graphic output of a spectrogram.

Author(s)

Michael Dietze

See Also

signal_spectrogram

Examples

## load example data set
data(rockfall)

## deconvolve data set
r <- signal_deconvolve(data = rockfall_eseis)

## calculate PSD
p <- signal_spectrogram(data = r)

## plot PPSD
plot_ppsd(data = p$PSD)

## plot PPSD with lower resolution, more smoothing and other colour
ppsd_color <- colorRampPalette(c("white", "black", "red"))

plot_ppsd(data = p$PSD, 
          res = c(200, 200), 
          n = c(15, 20), 
          col = ppsd_color(200))

coffeemuggler/eseis documentation built on Aug. 19, 2023, 9:57 p.m.