peaks_find | R Documentation |
Finds local maxima in sequential data.
peaks_find(object, ...)
## S4 method for signature 'GammaSpectrum'
peaks_find(object, method = c("MAD"), SNR = 2, span = NULL, ...)
object |
A GammaSpectrum object. |
... |
Extra parameters to be passed to internal methods. |
method |
A |
SNR |
An |
span |
An |
A local maximum has to be the highest one in the given window and has to be
higher than SNR \times noise
to be recognized as peak.
The following methods are available for noise estimation:
MAD
Median Absolute Deviation.
A PeakPosition object.
N. Frerebeau
Other signal processing:
baseline
,
peaks_search()
,
signal_integrate()
,
signal_slice()
,
signal_split()
,
signal_stabilize()
,
smooth()
## Import a LaBr spectrum
LaBr_file <- system.file("extdata/LaBr.TKA", package = "gamma")
LaBr_spc <- read(LaBr_file)
## Find peaks by channel
(LaBr_pks <- peaks_find(LaBr_spc)) # Ugly
plot(LaBr_spc, LaBr_pks)
## Search peaks by channel
(LaBr_pks <- peaks_search(LaBr_spc, index = c(86L, 207L, 496L), span = 7))
plot(LaBr_spc, LaBr_pks, split = TRUE)
## Import a BEGe spectrum
BEGe_file <- system.file("extdata/BEGe.CNF", package = "gamma")
BEGe_spc <- read(BEGe_file)
## Search peaks by energy
(BEGe_pks <- peaks_search(BEGe_spc, index = c(47, 63, 911, 1460)))
plot(BEGe_spc, BEGe_pks, split = TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.