Hi Brainstorm team,
I am attempting to estimate the deep sources of our sleep EEG data, but am running into an issue I cannot solve and am hoping you can help me. This is new territory for me, so I apologize in advance if there is a basic answer to my issue that I am missing.
For starters, the data: I have T1 anatomical scans for each individual that have been recon-all'ed in FreeSurfer, and the EEG are preprocessed in our homegrown MATLAB software and imported into Brainstorm as single-timepoint scalp maps for each individual (two identical time samples per map actually, as BST did not like importing a single sample FieldTrip-format dataset). The cortical and subcortical surfaces from FreeSurfer look appropriate after import into BST, to my eye at least, as do the EEG scalp maps and electrode locations. I generated realistic mixed (cortical/subcortical) forward models using OpenMEEG as described in the DBA and other tutorials (thank you by the way, for having such fully featured tutorials), then estimated sources using MNE/sLORETA with the identity matrix for noise covariance.
The issue I am having is that for some individuals, the lead field vector maps seem to appear empty when I quality control the OpenMEEG output with the "View lead field vectors" menu option. Further, for these individuals the source estimates have extremely small magnitude, to the point of labelling the color axis with "invalid scale", and seemingly noise-driven widespread diffuse activation compared to other participants. This is not the case for all individuals though, as some data points have lead field vector maps that seem sensible to me- meaning that there are small dipoles throughout the cortex and subcortical structures, with larger vector magnitudes at cortical patches directly under the electrode site. The individuals with sensible seeming vector maps also have more sensible source estimates, with activations centered in cortical patches that I would a priori theoretically expect to see. Additionally, for the problem participants, using OpenMEEG to generate a forward model based only on their cortex (i.e., not the mixed surface with subcortical structures added) results in a more sensible seeming lead field vector map, which leads me to suspect something is going wrong (or more accurately, that I am doing something wrong) with the mixed model/deep source estimation.
Can you all offer any guidance on where to look for what might be going wrong for me? Or, is this the nature of the deep source modeling, that it does not work well for some anatomies? The hippocampal, cerebellar, and striatal contributions are of primary importance to our research question, so if there is any way to make this work I would like to find it. Moreover, performing a separate cortex-only analysis of these data based on the template ICBM152 brain (for troubleshooting) reveals reasonable significant source-domain group effects in our comparisons, so I do believe there must be room for improvement in my OpenMEEG/DBA workflow.
In case it is useful, here is a more detailed step-by-step of the procedures I have taken so far per subject in BST:
-
Imported FreeSurfer recon-all folder and preprocessed FieldTrip-format EEG dataset
-
Created a vertex-downsampled atlas from the imported subcortical structures containing hippocampus, cerebellum, brainstem, caudate, putamen, pallidum, thalamus, amygdala, and accumbens, merged with cortex into a mixed surface
-
Created a source modeling atlas for the mixed surface, and set all structures to deep brain modeling
-
Generated Scalp, Inner Skull, and Outer Skull BEM surfaces with 1922 vertices each
-
Used "Force inside skull" menu option to correct cerebellum and brainstem clipping outside skull (this was in response to warnings that vertices were outside the skull for most participants)
-
Edited the MRI registration of electrode locations using the "Project electrodes on surface" function (we do not have individually measured electrode locations e.g. Polhemus, only standard cap layout positions)
-
Computed the head model using OpenMEEG BEM with adaptive integration and all other default settings
-
Set the noise covariance to be the identity matrix
-
Estimated sources with MNE/sLORETA
If additional information or any screenshots would be helpful for diagnosis please let me know, and I will gladly provide them. Thank you in advance for any advice you can offer, and thanks again for your work on Brainstorm!
- Ahren