Hi,
I'm interested in creating an MRI volume head model with the cerebellum, and scripting it.
I found this very helpful tutorial on how to create a new surface: http://neuroimage.usc.edu/brainstorm/Tutorials/DeepAtlas#Source_estimation
Using the GUI:
I successfully created a cortex+cerebellum surface by following the tutorial, computed the head model for an MRI volume, and computed sources. This gives a nice source localization that includes the entire cortex as well as the cerebellum. This is perfect and exactly what I want to do, but I want to script it.
Using Scripts:
I'm having trouble merging the cortex with the cerebellum. I can import anatomy, compute the head model, and compute sources but for the cortex only. I don't know what to do in order to merge the cortex ans cerebellum. See code below.
Any help is appreciated!
Thanks,
Elizabeth
I'm importing anatomy:
% Process: Import anatomy folder
sFiles = bst_process('CallProcess', 'process_import_anatomy', sFiles, [], ...
'subjectname', SubjectName.Name, ...
'mrifile', {FsDir, 'FreeSurfer'}, ...
'nvertices', 15000, ...;
'nas', NAS, ...
'lpa', LPA, ...
'rpa', RPA, ...
'ac', AC, ...
'pc', PC, ...
'ih', IH;
This is where I think the cerebellum should be added, but I don't know what code to use
Then I compute the head model:
% Process: Compute head model
sFiles = bst_process('CallProcess', 'process_headmodel', sFiles, [], ...
'Comment', '', ...
'sourcespace', 2, ... % 1 Cortex surface; 2 = MRI volume
'volumegrid', [], ...
'meg', 3, ... % Overlapping spheres
'eeg', 1, ... %
'ecog', 1, ... %
'seeg', 1, ... %
'openmeeg', struct(...
'BemFiles', {{}}, ...
'BemNames', {{'Scalp', 'Skull', 'Brain'}}, ...
'BemCond', [1, 0.0125, 1], ...
'BemSelect', [1, 1, 1], ...
'isAdjoint', 0, ...
'isAdaptative', 1, ...
'isSplit', 0, ...
'SplitLength', 4000));