/
test_MNE_ROIs_Definition.py
35 lines (35 loc) · 1.83 KB
/
test_MNE_ROIs_Definition.py
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
30
31
32
33
34
35
import glob, os
from dirs_manage import reset_directory
from MNE_ROIs_Definition01 import apply_create_noise_covariance, apply_inverse
from MNE_ROIs_Definition01 import apply_rois, merge_rois, stan_rois, group_rois
evt_st, evt_rt = 'LLst', 'LLrt'
stmin, stmax = 0.0, 0.4
rtmin, rtmax = -0.4, 0.0
subjects_dir = os.environ['SUBJECTS_DIR']
method = 'dSPM'
#calculate noise cov from empty room file
#emp_list = glob.glob(subjects_dir+'/*/MEG/*empty-raw.fif')
#apply_create_noise_covariance(emp_list)
#inverse epochs into the source space
epo_st_list = glob.glob(subjects_dir+'/*/MEG/*evtW_%s_bc-epo.fif' %evt_st)
epo_rt_list = glob.glob(subjects_dir+'/*/MEG/*evtW_%s_bc-epo.fif' %evt_rt)
apply_inverse(epo_st_list[:], method=method, event=evt_st)
apply_inverse(epo_rt_list[:], method=method, event=evt_rt)
#make ROIs for special event
stc_st_list = glob.glob(subjects_dir+'/fsaverage/%s_ROIs/*/*,evtW_%s_bc-lh.stc' % (method, evt_st))
stc_rt_list = glob.glob(subjects_dir+'/fsaverage/%s_ROIs/*/*,evtW_%s_bc-lh.stc' % (method, evt_rt))
apply_rois(stc_st_list, event=evt_st, tmin=stmin, tmax=stmax)
apply_rois(stc_rt_list, event=evt_rt, tmin=rtmin, tmax=rtmax)
#merge kinds of ROIs together for each subject
labels_path = glob.glob(subjects_dir+'/fsaverage/%s_ROIs/*[0-9]' %method)
merge_rois(labels_path, group=False, evelist=['LLst','LLrt'])
#standardize the size of ROIs and interegrate all the subjects ROIs
stan_path = subjects_dir+'/fsaverage/%s_ROIs/standard/' %method
reset_directory(stan_path)
stan_rois(stc_st_list, stan_path, size=8.0)
#merge ROIs across subjects, and select the common ROIs
labels_path = subjects_dir+'/fsaverage/%s_ROIs/' %method
merge_rois(labels_path, group=True)
mer_path = subjects_dir+'/fsaverage/%s_ROIs/merged/' %method
com_path = subjects_dir+'/fsaverage/%s_ROIs/common/' %method
group_rois(am_sub=8, com_path=com_path, mer_path=mer_path)