我正在做一些工作,包括分析波形并隔离超出幅度阈值的特定事件。那部分工作正常。 但是,完成分析之后,我希望能够绘制出结果波形,并在其上绘制振幅阈值(以及事件,但这是次要的)。
我已经使用grapics :: plot()做到了这一点,但是这花费了很长时间,并且使用tuneR包和tuneR :: plot()可以更快地吐出波形。但是,我不知道如何在其上绘制阈值。
这是我得到的最近的东西:
library(tuneR)
Wave <- readWave(filename = "Wave.wav", units = "seconds", from = 0, to = 1800)
Right_Threshold <- -6*sd(Wave@left)
Left_Threshold <- -6*sd(Wave@right)
tuneR::plot(Wave, panel.first={
lines(x = c(0,1800), y = c(Left_Threshold, Left_Threshold), col = 2)
})
哪个生成此内容:waveform plot with amplitude threshold overplotted
我希望能够在两个面板上绘制阈值,理想情况下,该解决方案还允许我使用graphics :: points()稍后绘制阈值穿越事件的记录。