-
Notifications
You must be signed in to change notification settings - Fork 1.4k
/
Copy pathmne_surf2bem.py
53 lines (38 loc) · 1.18 KB
/
mne_surf2bem.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
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
r"""Convert surface to BEM FIF file.
Examples
--------
.. code-block:: console
$ mne surf2bem --surf ${SUBJECTS_DIR}/${SUBJECT}/surf/lh.seghead \
--fif ${SUBJECTS_DIR}/${SUBJECT}/bem/${SUBJECT}-head.fif \
--id=4
"""
# Authors: The MNE-Python contributors.
# License: BSD-3-Clause
# Copyright the MNE-Python contributors.
import sys
import mne
def run():
"""Run command."""
from mne.commands.utils import get_optparser
parser = get_optparser(__file__)
parser.add_option(
"-s", "--surf", dest="surf", help="Surface in Freesurfer format", metavar="FILE"
)
parser.add_option(
"-f", "--fif", dest="fif", help="FIF file produced", metavar="FILE"
)
parser.add_option(
"-i",
"--id",
dest="id",
default=4,
help=("Surface Id (e.g. 4 for head surface)"),
)
options, args = parser.parse_args()
if options.surf is None:
parser.print_help()
sys.exit(1)
print(f"Converting {options.surf} to BEM FIF file.")
surf = mne.bem._surfaces_to_bem([options.surf], [int(options.id)], sigmas=[1])
mne.write_bem_surfaces(options.fif, surf)
mne.utils.run_command_if_main()