def test_source_space():
"Test SourceSpace dimension"
data_dir = mne.datasets.sample.data_path()
subjects_dir = os.path.join(data_dir, 'subjects')
annot_path = os.path.join(subjects_dir, '%s', 'label', '%s.%s.annot')
for subject in ['fsaverage', 'sample']:
mne_src = datasets._mne_source_space(subject, 'ico-4', subjects_dir)
vertices = [mne_src[0]['vertno'], mne_src[1]['vertno']]
ss = SourceSpace(vertices, subject, 'ico-4', subjects_dir)
# labels
for hemi_vertices, hemi in zip(ss.vertices, ('lh', 'rh')):
labels, _, names = read_annot(annot_path %
(subject, hemi, 'aparc'))
start = 0 if hemi == 'lh' else len(ss.lh_vertices)
hemi_tag = '-' + hemi
for i, v in enumerate(hemi_vertices, start):
label = labels[v]
if label == -1:
assert ss.parc[i] == 'unknown' + hemi_tag
else:
assert ss.parc[i] == names[label].decode() + hemi_tag
# connectivity
conn = ss.connectivity()
mne_conn = mne.spatial_src_connectivity(mne_src)
assert_array_equal(conn, _matrix_graph(mne_conn))
# sub-space connectivity
sssub = ss[ss._array_index('superiortemporal-rh')]
ss2 = SourceSpace(vertices, subject, 'ico-4', subjects_dir, 'aparc')
ss2sub = ss2[ss2._array_index('superiortemporal-rh')]
assert_array_equal(sssub.connectivity(), ss2sub.connectivity())
def test_source_space():
"Test SourceSpace dimension"
data_dir = mne.datasets.sample.data_path()
subjects_dir = os.path.join(data_dir, 'subjects')
annot_path = os.path.join(subjects_dir, '%s', 'label', '%s.%s.annot')
for subject in ['fsaverage', 'sample']:
mne_src = datasets._mne_source_space(subject, 'ico-4', subjects_dir)
vertno = [mne_src[0]['vertno'], mne_src[1]['vertno']]
ss = SourceSpace(vertno, subject, 'ico-4', subjects_dir)
# labels
for hemi_vertices, hemi in izip(ss.vertno, ('lh', 'rh')):
labels, _, names = read_annot(annot_path % (subject, hemi, 'aparc'))
start = 0 if hemi == 'lh' else len(ss.lh_vertno)
hemi_tag = '-' + hemi
for i, v in enumerate(hemi_vertices, start):
label = labels[v]
if label == -1:
eq_(ss.parc[i], 'unknown' + hemi_tag)
else:
eq_(ss.parc[i], names[label] + hemi_tag)
# connectivity
conn = ss.connectivity()
mne_conn = mne.spatial_src_connectivity(mne_src)
assert_array_equal(conn, _matrix_graph(mne_conn))
# sub-space connectivity
sssub = ss[ss.dimindex('superiortemporal-rh')]
ss2 = SourceSpace(vertno, subject, 'ico-4', subjects_dir, 'aparc')
ss2sub = ss2[ss2.dimindex('superiortemporal-rh')]
assert_array_equal(sssub.connectivity(), ss2sub.connectivity())
def test_source_space():
"Test SourceSpace dimension"
for subject in ['fsaverage', 'sample']:
mne_src = datasets._mne_source_space(subject, 'ico-4', subjects_dir)
vertno = [mne_src[0]['vertno'], mne_src[1]['vertno']]
ss = SourceSpace(vertno, subject, 'ico-4', subjects_dir, 'aparc')
# connectivity
conn = ss.connectivity()
mne_conn = mne.spatial_src_connectivity(mne_src)
assert_array_equal(conn, _matrix_graph(mne_conn))
# sub-space connectivity
sssub = ss[ss.dimindex('superiortemporal-rh')]
ss2 = SourceSpace(vertno, subject, 'ico-4', subjects_dir, 'aparc')
ss2sub = ss2[ss2.dimindex('superiortemporal-rh')]
assert_array_equal(sssub.connectivity(), ss2sub.connectivity())
def test_source_space():
"Test SourceSpace dimension"
for subject in ['fsaverage', 'sample']:
mne_src = datasets._mne_source_space(subject, 'ico-4', subjects_dir)
vertno = [mne_src[0]['vertno'], mne_src[1]['vertno']]
ss = SourceSpace(vertno, subject, 'ico-4', subjects_dir, 'aparc')
# connectivity
conn = ss.connectivity()
mne_conn = mne.spatial_src_connectivity(mne_src)
assert_array_equal(conn, connectivity_from_coo(mne_conn))
# sub-space connectivity
sssub = ss[ss.dimindex('superiortemporal-rh')]
ss2 = SourceSpace(vertno, subject, 'ico-4', subjects_dir, 'aparc')
ss2sub = ss2[ss2.dimindex('superiortemporal-rh')]
assert_array_equal(sssub.connectivity(), ss2sub.connectivity())
def test_source_space():
"Test SourceSpace dimension"
for subject in ['fsaverage', 'sample']:
path = os.path.join(subjects_dir, subject, 'bem', subject + '-ico-4-src.fif')
mne_src = mne.read_source_spaces(path)
vertno = [mne_src[0]['vertno'], mne_src[1]['vertno']]
ss = SourceSpace(vertno, subject, 'ico-4', subjects_dir, 'aparc')
# connectivity
conn = ss.connectivity()
mne_conn = mne.spatial_src_connectivity(mne_src)
assert_array_equal(conn, connectivity_from_coo(mne_conn))
# sub-space connectivity
sssub = ss[ss.dimindex('superiortemporal-rh')]
ss2 = SourceSpace(vertno, subject, 'ico-4', subjects_dir, 'aparc')
ss2sub = ss2[ss2.dimindex('superiortemporal-rh')]
assert_array_equal(sssub.connectivity(), ss2sub.connectivity())
