computing forward solutions separately for MEG and EEG

✉️ Mailing List Archive (read-only)
1

Hi and thanks Alex for the quick answer,

I already tried "EEG 064" but got the same message.

I followed your advice and inspect the first fiff file of my list: this is
what I got

In [2]: import mne

In [3]: info =
mne.io.read_info('/neurospin/meg/meg_tmp/MTT_MEG_Baptiste/MEG/sd130343/raw_sss/run3_DG_trans_sss.fif')

In [4]: info
Out[4]:
<Info | 22 non-empty fields
    acq_pars : unicode | 20983 items
    bads : list | 0 items
    ch_names : list | MEG0113, MEG0112, MEG0111, MEG0122, MEG0123,
MEG0121, ...
    chs : list | 393 items
    comps : list | 0 items
    description : unicode | 32 items
    dev_head_t : dict | 3 items
    dig : list | 69 items
    experimenter : unicode | 6 items
    file_id : dict | 4 items
    highpass : float | 0.0299999993294
    line_freq : float | 50.0
    lowpass : float | 330.0
    meas_date : numpy.ndarray | 2013-10-25 10:50:49
    meas_id : dict | 4 items
    nchan : int | 393
    orig_blocks : dict | 2 items
    proj_id : numpy.ndarray | 1 items
    proj_name : unicode | 7 items
    projs : list | 0 items
    sfreq : float | 1000.0
    subject_info : dict | 6 items
    acq_stim : NoneType
    ctf_head_t : NoneType
    dev_ctf_t : NoneType

In [11]: info["ch_names"]

[...]

'EEG059',

'EEG060',
'EOG061',
'EOG062',
'ECG063',
'EEG064',

[...]

So it seems that EEG064 name is not the problem here.

and finally...
I got it by looking at the rest of info dict: "projs : list | 0 items"
I was using not fully preprocessed data...
Then, when using fully preprocessed data:

arr2=($listRunPerSubj)

Dir=/neurospin/meg/meg_tmp/MTT_MEG_Baptiste/MEG/

for i in {0..87}; do
    echo \{arr1\[i\]\} &nbsp;&nbsp;&nbsp;&nbsp;echo {arr2[i]}

mne_process_raw --raw Dir{arr1[i]}/raw_sss/${arr2[i]}_trans_sss.fif
--lowpass 40 --save
Dir{arr1[i]}/mne_python/preproc_${arr2[i]}_trans_sss.fif --projon --proj
Dir{arr1[i]}/raw_sss/ecg_eogv/PCA_MEG_ECG_GRADS.fif --proj
Dir{arr1[i]}/raw_sss/ecg_eogv/PCA_MEG_ECG_MAGS.fif --proj
Dir{arr1[i]}/raw_sss/ecg_eogv/PCA_MEG_EOG61_GRADS.fif --proj
/Dir{arr1[i]}/raw_sss/ecg_eogv/PCA_MEG_EOG61_MAGS.fif --proj
Dir{arr1[i]}/raw_sss/ecg_eogv/PCA_EEG_EOG61.fif --proj
Dir{arr1[i]}/raw_sss/ecg_eogv/PCA_EEG_ECG.fif

done

mne_check_eeg_locations --file

Dir{arr1[i]}/mne_python/preproc_run3_DG_trans_sss_raw.fif --dig
/neurospin/meg/meg_tmp/MTT_MEG_Baptiste/MEG/wl130316/raw_sss/run3_DG.fif
--fix

I got the correct fix:

mne_check_eeg_locations version 1.2 compiled at Jan 5 2015 04:22:49

input file :
/neurospin/meg/meg_tmp/MTT_MEG_Baptiste/MEG/sd130343/mne_python/preproc_run3_DG_trans_sss_raw.fif
digitizer data file :
/neurospin/meg/meg_tmp/MTT_MEG_Baptiste/MEG/wl130316/raw_sss/run3_DG.fif

Checking...
EEG064 needs to be fixed (r = 0.0 mm).
1 channels need fixing.

Fixing...
Loaded 60 EEG locations + ref from
/neurospin/meg/meg_tmp/MTT_MEG_Baptiste/MEG/wl130316/raw_sss/run3_DG.fif
60 EEG locations applied

Sorry then for not having noticed that before and thanks for the advice
I will use this quick info checking in the future

Best,

Baptiste