## page was renamed from Tutorials/DuneuroMunsterData ''' [TUTORIAL UNDER CONSTRUCTION]''' https://neuroimage.usc.edu/brainstorm/Tutorials/tutoFEMadvancedFinal#preview = FEM tutorial: MEG/EEG Median nerve stimulation = ''Authors: Takfarinas Medani, Juan Garcia-Prieto, Wayne Mead '' This tutorial introduces the FEM modeling in the Brainstorm environment. Note that the operations used here are not detailed, the goal of this tutorial is not to introduce Brainstorm to new users. For in-depth explanations of the interface and theoretical foundations, please refer to the [[http://neuroimage.usc.edu/brainstorm/Tutorials#Get_started|introduction tutorials]]. In this tutorial, we describe the full FEM process as described in the SPIE paper <> == License (to check with Juan/John) == This tutorial dataset (MEG/EEG and MRI data) remains proprietary of xxxxyyyy. Its use and transfer outside the Brainstorm tutorial, e.g. for research purposes, is xxxx yyy. (allowed /prohibited /request) == Requirements == You have already followed all the introduction tutorials and you have a working copy of Brainstorm installed on your computer. If you want to reproduce this tutorial on your computer you need also to have these tools installed on your computer: [[https://simnibs.github.io/simnibs/build/html/index.html|SimNibs]]and [[http://brainsuite.org/|BrainSuite]]. Check [[https://neuroimage.usc.edu/brainstorm/Tutorials/FemMesh|this tutorial]] for further information about these tools. == Description of the experiment (todo: Juan) only one file? == The experiment consists of two stimulation protocols being conducted during a single scanning session in a MEG laboratory with an Elekta Triux (Megin, Finland) scanner. The subject is a right-handed 46 years old male. The two stimulation protocols consist of unilateral median nerve stimulation and an eyes-closed resting-state recording. Median nerve stimulation: * The stimulation signal was a square-wave pulse with 2Hz frequency and duration of 0.2ms. * An ISI (inter-stimulus interval) of 500ms with a variation of ±20ms in order to be able to average out time-locked noise to the stimulation, while remaining unnoticeable by the subject. * The stimulation was performed on both hands/wrists, independently, with a Digitimer DS7A stimulator. An electrode was placed on each current value were tuned to match the motor threshold of the subject on the stimulated hand, with a result of 10~12mA approximately. * Initially, the left wrist was stimulated for approximately 2 minutes (this corresponds with the file containing L1 in its name). After a 2-minute rest, while sitting in the MSR and with his head in the helmet, the stimulation was repeated (L2 file). After this stimulation, the subject was asked to have a 10 minutes rest, during which he was allowed to sit calmly with his head out of the helmet, although remaining hooked to the scanner at all times. Finally, two subsequent runs of right wrist stimulation (R1 and R2 files) with an intermediate 2-minute rest were performed. * Recordings were performed with a 1kHz sampling rate. Continuous HPI was disabled during these recordings. And high-pass filters were set to DC for MEG channels and 0.03Hz for EEG channels. * All files underwent a MaxFilter (version 2.3.13) tsss post-processing. Resting-State protocol: * The subject was recorded for approximately 80 minutes. Due to a 2GB maximum-size limitation for FIFF files, this translates into the recording being saved in files 'epi1' to 'epi4'. * The subject had his eyes closed. * Recordings were performed with a 1kHz sampling rate. Continuos HPI was enabled. High-pass filters were set to DC for MEG channels and 0.03Hz for EEG channels. * All files underwent a MaxFilter (version 2.3.13) tsss post-processing, with motion compensation active. == Download and installation (todo Francois upload data) == * Go to the [[http://neuroimage.usc.edu/bst/download.php|Download]] page of this website, and download the file: '''sample_tutoFEMadvance.zip''' * Unzip it in a folder that is not in any of the Brainstorm folders (program or database folder) * Start Brainstorm (Matlab scripts or stand-alone version) * Select the menu File > Create new protocol. Name it "'''tutoFEMadvanced'''" and select the options: * "'''No, use individual anatomy'''", * "'''No, use one channel file per condition'''". == Import the MRIs == * Switch to the "anatomical data" view, the left button in the toolbar above the database explorer. * Right-click on the tutoFEMadvanced folder > New subject > '''Subject01''' * Leave the default options you set for the protocol * Right-click on the subject node > Import anatomy MRI: Navigate to the selected folder and then select the T1 MRI, with the name "T1.nii" * Set the 6 required fiducial points (indicated in MRI coordinates). Fiducials point not always have to be located where we recommend in our [[https://neuroimage.usc.edu/brainstorm/Tutorials/ImportAnatomy#Fiducial_points|introductory]] and our [[https://neuroimage.usc.edu/brainstorm/CoordinateSystems?highlight=(auricular)|coordinate systems]] tutorials. Every subject is different and sometimes things get in the way, and it is not possible to follow canonical recommendations. That is fine! As long as there is a correct correspondence between fiducial points in the fif file and in the MRI, everything works fine. * For the following subject you can use these values :<
> . . - NAS: [121, 214, 50] . - LPA: [45, 111, 31] . - RPA: [198, 109, 33] . - AC: [120, 126, 75] . - PC: [120, 101, 73] . - IH: [120, 110, 132] . {{attachment:editMRIfiducials.jpg||height="150",width="600"}} * Click on the link "''' Click here to compute MNI transformation'''" and keep the default values. * Load the T2, by following the same process as previous, then select the T2 MRI with the name "T2.nii". {{attachment:importT2.jpg||height="150",width="600"}} Brainstorm will ask you to register the T2 to the T1, select the SPM option (require to have SPM installed), then reslice the T2, select "Yes". This process will take few minutes, just be patient. == Head model construction == === FEM head model === The first step requires the generation of the FEM head model, where the MRIs are segmented into the main issues and then tesselated into hexahedral or tetrahedral elements. The available methods within Brainstorm are listed in [[https://neuroimage.usc.edu/brainstorm/Tutorials/FemMesh|this page]]. Note: In the following, the SimNibs method is used. If you don't have this software installed you can import the SimNibs results (TODO: add the simNibsOutput to the data) * Select both the T1 and the T2 MRIs, using the mouse and the "Shift" key * Right-click on and then select the "Generate FEM mesh" * Select the SimNIbs method from the panel and keep the default options. In order to use the call SimNibs you need to have it installed on your computer, please follow the instruction as explained in [[https://neuroimage.usc.edu/brainstorm/Tutorials/FemMesh#SimNIBS|here]]. Keep all the options to their default values. Depending on your computer performance, this process can take 2 to 4 hours, so be patient. {{attachment:generateFEMmeshSimNibs2.png||height="150",width="600"}} At the end of this computation, Brainstorm will populate the windows with the following nodes * tissues: the segmentation output of the MRIs * FEM 72358V (simnibs, 5 layers): the FEM mesh of the MRIs with 72353 vertices, generated with the method SimNibs and the model has 5 layers (wm, gm, CSF, skull, and scalp). * Head mask: the outer surface of the scalp (or the head) * cortex_fem: the outer surface mesh of the grey matter * cortex_280753V: the pial surface generated by spm/cat within the SimNibs process * cortex_15002V : low version of the cortex_280753 At the end of the process, make sure that the file "cortex_15002V" is selected (downsampled pial surface, which will be used for the source estimation). If it is not, double-click on it to select it as the default cortex surface. <> === Remesh the head model === In some cases, from our experiences, the FEM mesh generated with SimNibs may have some issues with the DUNEuro FEM solver, this issue is mainly related to the air cavities in the head model as well as the hole. In order to avoid these possible issues in the next steps, we will correct the mesh using the FEM mesh tools integrated within Brainstorm. * Richt click on the FEM head model generated by simNibs > "Extract Surfaces". This option will extract the outer surfaces of the five tissues available in the FEM head model. * Hold the "Shif" keyboard and then select all the extracted surfaces (white, gray, CSF, skull, and scalp), then right-click and select "Generate FEM mesh" again, select "MergeMesh", and set the mesh resolution as shown in the following figure. Note that the more the mesh resolution is higher the better is the accuracy, however, it requires higher memory and computation time. * {{attachment:remeshHeadIso2mesh.jpg||height="150",width="600"}} The following figure shows the two model, left is the initial model obtained with SimNibs, right is the second model obtained from the Iso2Mesh re-mesh {{attachment:femMeshSimNibsVSiso2mesh2.jpg||height="150",width="600"}} You can also use this process to generate FEM model with high mesh densities. === FEM tensors === The FEM has the ability to incorporate anisotropic conductivity, es[acially for the white matter. Brainstorm offers the known methods to estimate the tensors from the DWI data. For a more detailed example please [[https://neuroimage.usc.edu/brainstorm/FemTensors|refere to this page.]] There are two main phases to compute the tensors, the first is the computation of the DTI from the DWI. The second is the estimation of the conductivity tensors from the DTI on each of the FEM mesh elements. __'''Step one: ''' __ * Right-click on the '''Subject01 '''> "Convert DWI to DTI" * Select the DWI data from the folder, "dwi.nii", the associated files "dwi.bvec" and "dwi.bval" must be in the same folder (same name but different extensions) {{attachment:compute&DisplayTensors.png|compute&DisplayTensors.png|height="150",width="600"}} BrainStorm will call internally the Brainsuite Software to compute the DTI. This process can take up to 20 minutes. At the end of this process, a node will appear in the Brainstorm database explorer under the name "DTI-EIG", this is a volume data that contains the 12 values of the eigenvalues and eigenvectors at each voxel. ''__Explanation of the options:__'' Brainstorm recognizes the tissue listed on the FEM had and assigns the default isotropic conductivities, as shown on the panel, users can change and use their own values. When the DWI data are computed, the conductivity tensors can be estimated on the white matter tissues using the Effective Medium Approach (EMA), Brainstorm offers two option, the EMA with a fixed factor k=0.736, or the EMA with the volume constraint (EMA+VC), please refer to this tutorial and the cited publication for further information. __'''Step two:''' __ Once the DTIs are computed, * Richt click on the FEM mesh and then '''"Compute FEM tensors"''' * From the panel select the options as shown in the following figure, then hit "Ok" {{attachment:computeDTI.png||height="150",width="600"}} This process can take up to 5min, depending on the resolution of the FEM mesh. === Visualisation of the FEM volume conductor === Brainstorms include the possibilities to display the FEM head models and the tensors, users can also overlay the display with the MRI as well as with the different surfaces. {{attachment:overlayModalities.png||height="150",width="600"}} The FEM tensors can be also displayed either on the mesh or on the MRI, to do so, right-click on the FEM mesh, then "Display FEM tensors", you can choose the displaying mode, the tensors can be displayed either as arrows (line) on the main eigenvector or as ellipsoids, on each FEM element (tetrahedron). {{attachment:dispTensorMenu.png||height="150",width="600"}} {{attachment:tensorsOnBrain.jpg||height="150",width="600"}} {{attachment:tensorsOnMRI.jpg||height="150",width="600"}} The size of the displayed tensors can be changed from the keyboard with the "Up" or "Down " arrows keys. You can also switch the display mode (from lines to ellipsoids to line or inversely) by using the shortcut "Shift + Space". ---- <> === On the hard disc === The "DTI-EIG" is saved structure as the Brainstorm MRI format, but it contains 12 values, from 1 to 9: components of the three eigenvectors, and from 10 to 12, the values of their norm to the eigenvalue. The FEM mesh, as shown in the figure, contains the following fields . Comment: 'FEM 722358V (simnibs, 5 layers)' . Vertices: number of vertices x3, the coordinates of the nodes Elements: number of element x4, list of the connectivity between nodes . Tissue: number of element x1, indexes or labels of the tissues, in Brainstorm the tissues are numbered from the inner to the outer layer (wm=1, gm=2, csf=3, skull=4, scalp=5) . TissueLabels: the labels of the tissues, {'white' 'gray' 'csf' 'skull' 'scalp'} . Tensors: number of elements ×12 if tensors are computed, otherwise it's empty. the 12 values are the eigenvalues and eigenvectors interpolated on each element of the mesh. . History: history of the processes applied on this file {{attachment:femDataOnHardDisc.png||height="150",width="600"}} === BEM head model === For the purpose of comparison between the FEM and the BEM, we will generate also the BEM surfaces for this subject and we will follow the same step as explained in [[https://neuroimage.usc.edu/brainstorm/Tutorials/TutBem?highlight=(bem)#Brainstorm|this page]]. The obtained surfaces will be used later for the BEM source computation. Richt-click on the subject01, then "MRI Segmentation" and then "Generate BEM surfaces", then keep the default options, Brainstorm, ('''1922 vertices) ''' for each layer and 4mm for the skull. Three new nodes appear in the database explorer, bem_head_1922V, bem_outerskull_1922V, and bem_innerskull_1922V. These surfaces will be used for the BEM forward computation. Now we have two head geometry available for this subject, BEM and FEM {{attachment:fem&bemModels.jpg|fem&bemModels.jpg|height="150",width="600"}} == Access the recordings == === Link the recordings === * Switch to the "functional data" view, the middle button in the toolbar above the database explorer. * You notice a folder with "eeg_postions" that contains an EEG cap, this folder is generated by the simnibs process. we will not use it in this tutorial. * Right-click on the subject folder > '''Review raw file''': * Select the file type: "MEG/EEG: Electa-Neuromag(*.fif)" * Navigate to the downloaded folder and select the file "subject01_median_nerve_L2_raw_tsss_anonymized.fif" * On the FIF event file, select "Event Channel", this will read the events markers * Apply the refine registration to the head points, click on yes. * A figure is opened to show the current registration MRI/MEG. It is already quite good for this dataset. The new file "Link to raw file" lets you access directly the contents of the MEG/EEG recordings. The "Neuromag channels(374)" contains the name of the channels and the position of the corresponding sensors (MEG/EEG) {{attachment:reviewMEEGsensors.jpg||height="150",width="600"}} === The channel file === * The recordings contain signals coming from different types of electrodes: * 306 MEG channels (102 magnetometers and 204 gradiometers ) * 64 EEG channels * 2 EOG channels: EO1 and EO2 * 1 Stim channel that contain the stimulation signals === Refine the MRI registration === * Right-click on the channel file > MRI registration > Edit... (EEG) * The white points are the electrodes, the green points are the additional digitized head points. To display the label of the electrodes, click on the [LABEL] button in the toolbar. To see what the other buttons in the toolbar are doing and how to use them, leave your mouse over them for a few seconds and read the description. * Now try to manipulate the position of the EEG+MEG sensors using '''rotations '''and '''translations '''only (no "resize" or individual electrodes adjustments). The objective is to have all the points close to the surface and the three forehead points inside the little peaks on the surface (due to markers in the MRI). * The rotation+translation is going to be applied both to the EEG and the MEG sensors. After you are done with this solid registration part, you can click on the button '''"Project electrodes on scalp surface"''', it will help for the source modeling later. The green points (digitized) stay in place, the white points (electrodes) are now projected on the skin of the subject. * If you feel like you didn't do this correctly, close the figure and cancel the modifications, then try again. It takes a few trials to get used to this rotation/translation interface. * Click on '''[OK]''' when done. * Answer '''YES''' to save the modifications. * Answer '''YES''' again to apply the solid transformation (rotation+translation) to the MEG sensors. * Now all the electrodes are projected on the surface. * (refine head with digitized points) {{attachment:MEEGsensors.jpg||height="150",width="600"}} === Pre-processing === * Drag and drop the "Link to raw file" into the Process1 list. * Select the process "'''Frequency > Power spectrum density'''", configure it as in the figure * After the computation, double-click on the new PSD file to display it. {{attachment:psd.jpg||height="500",width="600"}} * The lines on the top represent the EEG electrodes, the lines at the bottom the MEG sensors(MAG and GRAD). If you want to get clearer plots, you can calculate separately the spectrum for the two types of sensors separately, by running twice the process "Power spectrum density", once with sensor types = "MEG" and once with "EEG", instead of running in on both at the same time as we did. * Observations (below 250Hz): * Peak around 10Hz: Alpha waves from the subject's brain * Peaks at 60Hz, 120Hz, 180Hz, 240Hz on EEG + MEG: Power lines (60Hz+harmonics) * Most of the signal looks clean and no bad channel is observed * If we review quickly the EEG and EOG signals, we notice that they are quite clean and no artifacts are observed. The experiment is well designed and the subject was not moving the eyes (blinks and slow movements), maybe because there was a fixation cross for this experiment. However, we will apply at least a high-pass filter to make the signals easier to process (we are not interested in very low frequencies in this experiment). === Frequency filters === * In Process1, select the "Link to raw file". * Select process '''Pre-process > Band-pass filter''': Frequency='''[20, 250]Hz''', Sensors='''MEG,EEG''' * Add processs '''Pre-process > Notch filter''': Frequencies='''[60 120 180]Hz''', Sensors='''MEG,EEG''' {{attachment:freqFilters.jpg||height="500",width="600"}} * In Process1, select the filtered file "Raw | Band | Notch". * Run the process "'''Frequency > Power spectrum density'''", with the same options as before. And compare with the previous figure. You may find that the power line frequency and its harmonics are removed and the data is filtered between 20 and 250Hz. {{attachment:spectreClean.jpg||height="500",width="600"}} === EEG: Average reference === * Right-click on the filtered file "Raw|band|notch" > EEG > Display time series. * In the Record tab, menu '''Artifacts > Re-reference EEG''' > "AVERAGE". {{attachment:avgRef.jpg||height="500",width="600"}} At the end, the window "Select active projectors" is open to show the new re-referencing projector. Just close this window. To get it back, use the menu Artifacts > Select active projectors. To keep the interface clean and eeasy to follow, you may need also to remove the previous and keep only the last folder that inlude all the post-processing. === Epoching === In this experiment, the electric stimulation is sent with a frequency of 2Hz, meaning that the inter-stimulus interval is 500ms. We are going to import epochs of 300ms around the stimulation events (-100 to 200ms). * Right-click on the filtered file "Raw|band|notch" > '''Import in database''': * Check "'''Use events'''" and select "'''2'''", this is the name that brainstorm assigns for the stimulus onset, it may have different name on your computer. Set epoch time: '''[-100, 200]''' '''ms ''' and Apply SSP/ICA in order to apply the average reference. {{attachment:importEpochs.jpg||height="500",width="600"}} In the end, you are asked whether you want to ignore one epoch that is shorter than the others. This happens because the acquisition of the MEG signals was stopped less than 300ms after the last stimulus trigger was sent. Therefore, the last epoch cannot have the full [-100,200]ms time definition. This shorter epoch would prevent us from averaging all the trials easily. As we already have enough repetitions in this experiment (240), we can just ignore it. The total number of epochs is then 239. === Averaging === * Drag and drop all the Trigger01 trials to the Process1 tab. * Run the process '''Average > Average files''': By trial group (folder average) {{attachment:avrgEpochs.jpg||height="500",width="600"}} Review the average for the '''MEG''' and '''EEG''' as a topography plot -Right-click on the averaged signal -EEG (then MEG) > 2D disc {{attachment:megEegTopography.jpg||height="700",width="800"}} AT the selected time, an ERP is visible with a nice dipolar pattern on the sensors, both for the EEG and MEG. == Source estimation == We are going to use the realistic FEM head model, the one previously generated from the MRI as explained in sections above. Go to the "Anatomy" view, and make sure that the FEM head model is highlighted in green color (it should be the case if you have only one model). You may also highlight the cortex to use for the computation (select the cortex_15002V). '''IMPORTANT:''' You can compute the forward model both for EEG and MEG simultaneously, however, using the high mesh resolution model we recommend to compute separately the head model for each modality(EEG and then MEG). The time required for EEG is aound one hour for ~70 channles, for the MEG with 306 sensors it can take up to 4 hours or more (with the integrations points). The EEG/MEG FEM computation depends on the computation of the FEM transfer matrix, and the computation of this matrix is related to the resolution of the FEM head mesh (number of vertices) and the number of sensors. In most of the case, the number of EEG sensors is lower than the number of MEG sensors. Furthermore, internally the MEG sensors modeling uses the integrations points, which increase the number of computation points (~ multiplied by 4 for the magnetometers and by 8 for the gradiometers). Therefore the MEG requires more time than the EEG. To reduce the MEG computation time, there are some tips : 1. Use only the inners tissues (wm, gm, and CSF) ==> reduce the number of vertices 1. Do not use the integration points ==> reduce the number of virtual sensors 1. These parameters can be tuned from the DUNEuro Options panel (see the advanced panel) In this tutorial, we will separate the modalities. === EEG Head model with DTI tensors === * Go back to the "Functional data" view, right-click on the channel file > '''Compute head model''' * Select the DUNEuro FEM, change the comment to "DUNEuro FEM EEG DTI tensors", '''select only the EEG, ''' as shown in the panel * On the DUNEuro Options panel, click on the "Show details" button for more option * Set the DUNEuro option as shown in the following figure, and then click on "Ok" {{attachment:eegDuneuroPanels.jpg||height="300",width="600"}} The FEM EEG head modeling will start and depending on the performance and the workload of your computer, it can last between 30min to 2 hours. Once the computation is done, a new node with the name "DUNEuro FEM EEG DTI tensors" appear on the Brainstorm database navigator, which can be used for the source localization process. === MEG Head model with DTI tensors === * Go back to the "Functional data" view, right-click on the channel file > '''Compute head model''' * Select the DUNEuro FEM, change the comment to "DUNEuro FEM MEG DTI tensors", '''select only the MEG, ''' as shown in the panel * On the DUNEuro Options panel, click on the "Show details" button for more option * For this case (MEG), in order to reduce the computation, we will check only the inner layers (white, gray, and CSF). Not that this is true only for the MEG [ref]. * Set the DUNEuro option as shown in the following figure, and then click on "Ok" {{attachment:megDuneuroPanelsDTI.jpg||height="300",width="600"}} The FEM MEG head modeling will start and depending on the performance and the workload of your computer, it can last between 1hour to 4 hours (even more if you select all the tissues). Once the computation is done, a new node with the name "DUNEuro FEM MEG DTI tensors" appears on the Brainstorm database navigator, which can be used for the source localization process. === EEG Head model with isotropic conductivity === In the case where the DTIs are not available, or if the users want to use the isotropic conductivity, the following section explains how to do it. * First you need to make sure that the tensors are not computed * If you have the tensor computed, as in this tutorial, we need to clear these values. Go to the "Anatomy" view, right-click on the FEM head model, and then "Clear tensor" * Go back to the "Functional data" view, right-click on the channel file > '''Compute head model,''' and set the value as in the figures. * Compare to the previous panels, in this case, the DUNEuro options allow you to change the isotropic conductivity for each tissue. In this example keep the * Click on "Ok", then the FEM computation will start with the new value of the conductivities. {{attachment:eegDuneuroPanelsIso.jpg||height="300",width="600"}} === MEG Head model with isotropic conductivity === Same step as before, first you need to make sure that the tensors are not computed, * right-click on the channel file > '''Compute head model,''' and set the value as in the figures. * Compare to the previous panels, in this case, the DUNEuro options allow you to change the isotropic conductivity for each tissue. In this example keep only the inner tissues (wm, gm and csf) * Click on "Ok", then the FEM computation will start with the new value of the conductivities. {{attachment:megDuneuroPanelsIso.jpg||height="300",width="600"}} In this tutorial, for each modelity we have computed two forward models (isotropic and anisotrpic). This is a good option to compare if there is any difference between the isotropic model and the anisotropic model. === Noise covariance matrix === * We will use the baseline of the single epochs to calculate the noise covariance matrix. * Right-click on the Trigger01 epochs group > Noise covariance > '''Compute from recordings'''. Enter the same baseline interval we used for removing the DC offset: '''[-100, -10] ms''' {{attachment:noiseCovariance.jpg||height="300",width="600"}} === Inverse model === * Right-click on the head model > '''Compute sources'''. * Select "Dipole modeling option" and '''EEG'''. Repeat the same operation for '''MEG''' and make sure to select the appropriate head model at each step. {{attachment:computeDipolesInverseAll.jpg||height="300",width="600"}} It is better to study separately the two modalities because the method for combining MEG and EEG are not working well yet within Brainstorm. This operation creates a shared inversion kernel and one source link for each block of recordings in the folder. {{attachment:sharedKernelDipolesEegMeg.jpg||height="300",width="600"}} Drag and drop the Dipoles kernels in the Process table and select the process > Run> Dipole Scanning> {{attachment:dipoleScanningProcess.jpg||height="300",width="600"}} If you are not familiar with those concepts, please refer to the [[https://neuroimage.usc.edu/brainstorm/Tutorials/SourceEstimation|Source estimation tutorial]]. Once the dipoles are computed, we display both EEG and MEG separately in the following figures. We merge the MEG ISO and the MEG DTI dipoles as well as the EEG dipoles. {{attachment:dipolesMergeEEG&MEG.jpg|dipolesMergeEEG&MEG.jpg|height="300",width="600"}} === P20/N20 source localization === The following figures (left EEG, right MEG), show the localization of the dipoles on the MRI at the 21ms. The ISO model is colored in greeen wherease the DTI model is colored in red. {{attachment:eegMegDipolesSagittal.jpg||height="300",width="600"}} {{attachment:eegMegDipolesAxial.jpg||height="300",width="600"}} {{attachment:eegMegDipolesCoronal.jpg||height="300",width="600"}} In this experiment, the anisotropy does not show a significant effect on the source localization, whereas it shows a slich difference in the orientation. When we compare between the EEG and the MEG, there is a difference of 20mm between the localization between the MEG and the EEG as well as a difference in the orientation in the coronal view. If we compare the difference between the EEG and MEG dipoles, the following figures (left anisotropy, right isotropic), show the localization of the dipoles on the MRI at the 21ms. The EEG dipole is colored in green wherease the MEG dipole is colored in red. {{attachment:eegMegDipolesSagittal2.jpg||height="300",width="600"}} {{attachment:eegMegDipolesAxial2.jpg||height="300",width="600"}} {{attachment:eegMegDipolesCoronal2.jpg||height="300",width="600"}} Aniso: Eeg [-7.5 -35.7 89.8] vs Meg [-0.2 -33.6 88.6] ==> '''distance = 7.7mm ''' Iso: Eeg [-7.5 -35.7 89.8] vs Meg [-0.2 -33.6 88.6] ==> '''distance = 7.7mm ''' In both modalities, the dipoles are located exactely in the same positions, and there is 7.7mm between the two dipoles. Further more, we notice a difference on the orientation on the coronal and sagital views === Comparaison with BEM for EEG and with OS for the MEG === In order to compare qualitatively the results of the FEM models, we will use the previous methods already implemented in Brainstorm. For the EEG: compute the forward using the BEM For the MEG: compute the forward using overlapping spheres, For all the methods, we will perform the dipole scanning as done previously but with the new forward models, the results are shown in the following figures: In the following figures, left is the dipoles computed from the EEG (BEM in red and FEM in green), right are the dipoles computed from the MEG (OS in red and FEM in green) {{attachment:eegMegDipolesSagittal3.jpg||height="300",width="600"}} {{attachment:eegMegDipolesAxial3.jpg||height="300",width="600"}} {{attachment:eegMegDipolesCoronal3.jpg||height="300",width="600"}} EEG: fem [-7.5 -35.7 89.8] ; bem [-10.7 -35.0 89.5] ==> '''3.3mm''' MEG: fem [-0.2 -33.6 88.6] ; os [-10.7 -39.3 83.2] ==> '''13.1 mm''' <> === Resect the neck and the bottom part of the head model === The FEM computation can be time and memory-consuming, on way to reduce this issue is to remove the part of the FEM mesh that is not required. In some cases where the field of the MRI is large, you can have the mesh of the neck and even the shoulders. In Brainstorm, it is possible to remove these parts. You can this by following these steps : right-click on the FEM-mesh => Resect Neck ==> keep the default value ==> Ok ---- {{attachment:resectNeck.JPG||height="150",width="600"}} ---- Once the process is finished, a new FEM head model appears in the database with the same name and a prefix "resect" which is the resected head. The following figure shows the model before and after re-section. {{attachment:FemMeshAllandResect.JPG||height="150",width="600"}} Note that the initial model has 743828 nodes and 4079587 and the resected model has 613955 nodes and 3400957 elements. == Scripting == This section is under devellopement == DUNEuro advanced options == This section is under development A set of advanced options are made available and can be easily changed. Here a short explanation is given for each option. * '''FEM''' layers & conductivities: same explanation as in the previous section. Moreover, in the case where the conductivity tensor are previously computed, Brainstorm detects these tensors and load them. In this case, the users can't change the conductivities scalar values, since they are not used. The following is displayed panel is displayed. * '''FEM''' solver type: * CG or Continuous Galerkin: This is the standard Lagrangian method. * DC or the Discontinuous Galerkin: In this version, only the Fitted '''FEM''' approaches are integrated, that require the '''FEM''' mesh of the head model. Unfitted methods are also available within DUNEuro, and will be integrated soon in Brainstorm. * '''FEM''' source model: The list of the available source are * Venant * Subtraction * Partial Integration For more information about these methods, users can check this thesis (Vorwek thesis). * Source space * Shrink source space: the location of dipoles are moved inward by the specified value in this field(in mm). * Force source space: this is required in the case where the dipoles are not within the GM matter. * Outputs options * Save transfer matrix: All these parameters are stored and passed to the DUNEuro as a text file. This file is the main interface that passes the parameters from Brainstorm to DUNEuro. More details about the integration can be found in these links: * GitHub repository for the [[https://github.com/brainstorm-tools/bst-duneuro|Brainstorm-DUNEuro]] compilation and integration * GitHub repository for the [[https://github.com/tmedani/duneuro_interface|matlab-duneuro interface]] * Brainstorm-DUNEuro integration discussions: * [[https://github.com/brainstorm-tools/brainstorm3/issues/185|Brainstorm-simbio/DUNEuro implementation]]/head model generation * [[https://github.com/brainstorm-tools/brainstorm3/issues/242|Integrate the DUNEuro FEM computation]] * [[https://github.com/brainstorm-tools/brainstorm3/pull/282|Genration of the FEM tensor]] * [[https://github.com/brainstorm-tools/brainstorm3/issues/302|Generation of the FEM source space]] == Troubleshooting and solution(TODO) == Dipoles outside of the grid for the MEG: remove this dipole or replace it with its spatial neighbors (not the neighbor in the list) ==> not easy to handle write small tuto?