Intersecting MNE BEM meshes

Dear All,

I would like to compute nonintersecting decimated (~<=5000 triangles)
MNE BEM meshes, however using MNE on several patients with the default
mne_watershed seems to lead to intersecting surfaces. Is there a
workaround this ?
Thank you,

Octavian.

Hi Octavian,

There are several things to consider about this

Do you have any EEG data?

If NO

Just use a single layer BEM (This is supported by the literature).

If you do have EEG data,

First I would say you need to try to get better MR data FreeSurfer has
pages on their wiki about this. Also, I highly recommend using the
multi-echo flash sequences as these dramatically improve the BEM (I
would argue they are necessary, but I imagine many others would
disagree).

If it is too late for that, you can try to generate them by utilizing
the -less and -more options described on the FreeSurfer wiki for
mri_watershed (the mne_watershed_bem script can help you find a good
place to start). vary the names of the output files and use the best
ones you can get. This is slow and difficult, but could work.

D

hi,

when you do:

mne_setup_forward_model --surf --ico 4

as in :

https://github.com/mne-tools/mne-scripts/blob/master/sample-data/run_meg_tutorial.sh#L25

you properly decimate the bem meshes.

the -bem.fif file then contains the triangulations you want with
an acceptable resolution for BEM computation.

hope this helps,

Alex