mne.preprocessing.annotate_amplitude#
- mne.preprocessing.annotate_amplitude(raw, peak=None, flat=None, bad_percent=5, min_duration=0.005, picks=None, *, verbose=None)[source]#
Annotate raw data based on peak-to-peak amplitude.
Creates annotations
BAD_peak
orBAD_flat
for spans of data where consecutive samples exceed the threshold inpeak
or fall below the threshold inflat
for more thanmin_duration
. Channels where more thanbad_percent
of the total recording length should be annotated with eitherBAD_peak
orBAD_flat
are returned inbads
instead. Note that the annotations and the bads are not automatically added to theRaw
object; useset_annotations()
andinfo['bads']
to do so.- Parameters
- rawinstance of
Raw
The raw data.
- peak
float
|dict
|None
Annotate segments based on maximum peak-to-peak signal amplitude (PTP). Valid keys can be any channel type present in the object. The values are floats that set the maximum acceptable PTP. If the PTP is larger than this threshold, the segment will be annotated. If float, the minimum acceptable PTP is applied to all channels.
- flat
float
|dict
|None
Annotate segments based on minimum peak-to-peak signal amplitude (PTP). Valid keys can be any channel type present in the object. The values are floats that set the minimum acceptable PTP. If the PTP is smaller than this threshold, the segment will be annotated. If float, the minimum acceptable PTP is applied to all channels.
- bad_percent
float
The percentage of the time a channel can be above or below thresholds. Below this percentage,
Annotations
are created. Above this percentage, the channel involved is return inbads
. Note the returnedbads
are not automatically added toinfo['bads']
. Defaults to5
(5%).- min_duration
float
The minimum duration (sec) required by consecutives samples to be above
peak
or belowflat
thresholds to be considered. to consider as above or below threshold. For some systems, adjacent time samples with exactly the same value are not totally uncommon. Defaults to0.005
(5 ms).- picks
str
|list
|slice
|None
Channels to include. Slices and lists of integers will be interpreted as channel indices. In lists, channel type strings (e.g.,
['meg', 'eeg']
) will pick channels of those types, channel name strings (e.g.,['MEG0111', 'MEG2623']
will pick the given channels. Can also be the string values “all” to pick all channels, or “data” to pick data channels. None (default) will pick good data channels. Note that channels ininfo['bads']
will be included if their names or indices are explicitly provided.- verbosebool |
str
|int
|None
Control verbosity of the logging output. If
None
, use the default verbosity level. See the logging documentation andmne.verbose()
for details. Should only be passed as a keyword argument.
- rawinstance of
- Returns
- annotationsinstance of
Annotations
The annotated bad segments.
- bads
list
The channels detected as bad.
- annotationsinstance of
Notes
This function does not use a window to detect small peak-to-peak or large peak-to-peak amplitude changes as the
reject
andflat
argument fromEpochs
does. Instead, it looks at the difference between consecutive samples.When used to detect segments below
flat
, at leastmin_duration
seconds of consecutive samples must respectabs(a[i+1] - a[i]) ≤ flat
.When used to detect segments above
peak
, at leastmin_duration
seconds of consecutive samples must respectabs(a[i+1] - a[i]) ≥ peak
.
Thus, this function does not detect every temporal event with large peak-to-peak amplitude, but only the ones where the peak-to-peak amplitude is supra-threshold between consecutive samples. For instance, segments experiencing a DC shift will not be picked up. Only the edges from the DC shift will be annotated (and those only if the edge transitions are longer than
min_duration
).This function may perform faster if data is loaded in memory, as it loads data one channel type at a time (across all time points), which is typically not an efficient way to read raw data from disk.
New in version 1.0.