Hello. I would like to perform MEG source reconstruction of Cerebellum sources located on the (outer, simple-shaped) surface of the Cerebellum and I donāt understand how to do it.
There is an example on the MNE website allowing to get volume sources for Cerebellum (but not surface once). I want to reduce the number of sources and restrict myself to the surface ones. But for the surface sources one has to do something different because the setup_source_space and setup_volume_source_space have different arguments (in particular one cannot just provide a volume label).
If I understand well, the way to go it to provide a surface that contains Cerebellum to setup_source_space and then select Cerebellum-related source using an appropriate atlas. The problem is that my knowledge of freesurfer and its interaction with MNE is not solid enough so I donāt know
what is the right surface to use (is there a ready one or I have to recreate a new one?) and
what is the right atlas to use (if I understand well, many atlases donāt contain Cerebellum).
Could you help me with these two issues?
Iām using latest stable version MNE 1.2 on Ubuntu 20 (though perhaps it is not important for this question).
Sorry for confusion, I did not mean surface of the cerebellum in the sense of caring about all of its super complex geometry. I just want the surface of the outer convex hull of the cerebellum, especially its part that is closest to the MEG sensors (I realize I have not mentioned that I am working with MEGā¦).
But if I were to follow your advise on increasing volumetric sampling size, do you know how could I do it? Just looking at documentation for setup_volume_source_space function I donāt see which parameter should control it.
The pos variable in setup_volume_source_space sets the grid sampling size, just increase it to whatever you need (the default 5mm is a typical good choice though).
If you want the outer surface of the cerebellum - you can probably do this in 3Dslicer software from the labelled voxels. I just donāt know how you would import those surfaces back into MNE. I would try the volume based methods first, since it will be easier to get to an end result and you wonāt be gaining any benefit from the surface, since you cant use the surface orientation in your model.
Hello again. I have tried to follow your suggestion but unfortunately pos argument of setup_volume_source_space seems to have no effect. The number of sources I get is always the same, whether it is 4, 40 or 400. Do you know what am I doing wrong?
conductivity = (0.3,) # for single layer
model = mne.make_bem_model(subject=subject, ico=4,
conductivity=conductivity,
subjects_dir=subjects_dir)
bem = mne.make_bem_solution(model)
volume_label = 'Left-Cerebellum-Cortex'
lh_cereb = mne.setup_volume_source_space(
subject, volume_label=volume_label, bem=bem,
subjects_dir=subjects_dir, pos=CB_spacing)
volume_label = 'Right-Cerebellum-Cortex'
rh_cereb = mne.setup_volume_source_space(
subject, volume_label=volume_label, bem=bem,
subjects_dir=subjects_dir, pos=CB_spacing)
In the first example its a 4mm sampling and I get 690 sources, in the second its a 10mm sampling and I get 46 sources. So it does appear to be downsampling the source space using the pos flag. I tested this on mne version 1.2.1, but it should work on earlier versions as well. Not sure why you are seeing different results.
I have found a problem that was not related to MNE (I have copied in the message only a part of the code inside an āifā clause that did not get actually executed).
After solving it I have gotten results similar to yours.