Add blink noise to raw data.
Raw
The raw instance to modify.
None
| str
| dict
| tuple
| array
Name of the position estimates file. Should be in the format of
the files produced by MaxFilter. If dict, keys should
be the time points and entries should be 4x4 dev_head_t
matrices. If None, the original head position (from
info['dev_head_t']
) will be used. If tuple, should have the
same format as data returned by head_pos_to_trans_rot_t
.
If array, should be of the form returned by
mne.chpi.read_head_pos()
.
str
Either ‘hann’, ‘cos2’ (default), ‘linear’, or ‘zero’, the type of forward-solution interpolation to use between forward solutions at different head positions.
int
| None
The number of jobs to run in parallel. If -1
, it is set
to the number of CPU cores. Requires the joblib
package.
None
(default) is a marker for ‘unset’ that will be interpreted
as n_jobs=1
(sequential execution) unless the call is performed under
a joblib.parallel_backend()
context manager that sets another
value for n_jobs
.
None
| int
| instance of RandomState
A seed for the NumPy random number generator (RNG). If None
(default),
the seed will be obtained from the operating system
(see RandomState
for details), meaning it will most
likely produce different output every time this function or method is run.
To achieve reproducible results, pass a value here to explicitly initialize
the RNG with a defined state.
The random generator state used for blink, ECG, and sensor noise
randomization.
str
| int
| None
Control verbosity of the logging output. If None
, use the default
verbosity level. See the logging documentation and
mne.verbose()
for details. Should only be passed as a keyword
argument.
Raw
The instance, modified in place.
See also
Notes
The blink artifacts are generated by:
Random activation times are drawn from an inhomogeneous poisson process whose blink rate oscillates between 4.5 blinks/minute and 17 blinks/minute based on the low (reading) and high (resting) blink rates from [1].
The activation kernel is a 250 ms Hanning window.
Two activated dipoles are located in the z=0 plane (in head coordinates) at ±30 degrees away from the y axis (nasion).
Activations affect MEG and EEG channels.
The scale-factor of the activation function was chosen based on visual inspection to yield amplitudes generally consistent with those seen in experimental data. Noisy versions of the activation will be stored in the first EOG channel in the raw instance, if it exists.
References
mne.simulation.add_eog
#Simulate raw data using subject anatomy