1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
|
"""
=========================================
Reading/Writing a noise covariance matrix
=========================================
"""
# Author: Alexandre Gramfort <alexandre.gramfort@telecom-paristech.fr>
#
# License: BSD (3-clause)
print(__doc__)
import mne
from mne.datasets import sample
data_path = sample.data_path()
fname = data_path + '/MEG/sample/sample_audvis-cov.fif'
cov = mne.Covariance(fname)
print(cov)
###############################################################################
# Show covariance
# Note: if you have the measurement info you can use mne.viz.plot_cov
import matplotlib.pyplot as plt
plt.matshow(cov.data)
plt.title('Noise covariance matrix (%d channels)' % cov.data.shape[0])
plt.show()
|
