plot_ppsd {eseis}

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))