This function plots dynamically a sliding spectrogram along a time wave.
dynspectro(wave, f, channel = 1, slidframe = 10,
wl = 512, wn = "hanning", zp = 0, ovlp = 75,
fftw = FALSE, dB = TRUE, plot = TRUE,
title = TRUE, osc = FALSE,
tlab = "Time (s)", flab = "Frequency (kHz)", alab = "Amplitude",
from = NULL, to = NULL,
collevels = NULL, palette = spectro.colors,
envt = NULL, msmooth = NULL, ksmooth = NULL,
coltitle = "black", colbg = "white", colline = "black",
colaxis = "black", collab = "black", cexlab = 1,
fontlab = 1, colwave = "black",
coly0 = "lightgrey", colcursor = "red", bty = "l")
an R object.
sampling frequency of wave
(in Hz). Does not need to
be specified if embedded in wave
.
channel of the R object, by default left channel (1).
size of the sliding frame (in percent of the wave duration).
if at
is not null, length of the window for the analysis
(even number of points, by defaults = 512).
window name, see ftwindow
(by default "hanning"
).
zero-padding (even number of points), see Details
.
overlap between two successive windows (in % ).
if TRUE
calls the function FFT
of the
library fftw
. See Notes of the spectro
.
a logical, if TRUE
then uses dB values
logical, if TRUE
plots in an ew graphics device the successive
spectrograms sliding along the time wave (by default TRUE
).
logical, if TRUE
adds a title with the time position of the current
spectrogram along the time wave.
logical, if TRUE
plots an oscillogram beneath
the sliding spectrogram with a cursor showing the position of the
current spectrum (by default FALSE
).
title of the time axis.
title of the frequency axis.
title of the amplitude axis.
start mark where to compute the sliding spectrogram (in s).
end mark where to compute the sliding spectrogram (in s).
a set of levels which are used to partition the amplitude range of the spectrogram.
a color palette function to be used to assign colors in the plot.
the type of envelope to be plooted:
either "abs" for absolute amplitude envelope or "hil" for Hilbert amplitude envelope.
See env
.
when env
is not NULL
,
a vector of length 2 to smooth the amplitude envelope with a
mean sliding window. The first component is the window length
(in number of points). The second component is the overlap between
successive windows (in %). See env
.
if title
is TRUE
, colour of the title.
background colour.
colour of axes line.
colour of the axes.
colour of axes title.
character size for axes title.
font for axes title.
colour of the oscillogram or of the envelope (only when osc
is TRUE
).
colour of the y=0 line (only when osc
is TRUE
).
colour of oscillogram cursor (only when osc
is TRUE
).
the type of box to be drawn around the oscillogram (only
when osc
is TRUE
).
This function returns a list of three items:
a numeric vector corresponding to the time axis.
a numeric vector corresponding to the frequency axis.
a numeric matrix corresponding to the amplitude values.
Each column is a Fourier transform of length wl/2
.
Use the slider panel to move along the time wave.
The function requires the package rpanel that is based on the
package tcltk.
The function is mainly written for inspecting long sounds.
The function is based on image
for fast display when
spectro
is based on filled.contour
.
Displaying the amplitude envelope with the argument envt
can
slow down significantly the display.
# NOT RUN {
data(sheep)
require(rpanel)
dynspectro(sheep, ovlp=95, osc=TRUE)
# }
Run the code above in your browser using DataLab