|
Size: 556
Comment:
|
Size: 3533
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 1: | Line 1: |
| In this tutorial we will compare the performance of the BEM, FEM forward computation in basic model (realistic head model 3 layer) and also within in a spherical head model with an analytical solution. | == Introduction == In this tutorial, we compare the head models computation using the methods available in brainstorm. At this time, the FEM implementation is available for the EEG and the MEG computation. We will integrate and test soon the FEM for sEEG and ECOG. == EEG within a spherical model == === The volume conductor or the head model: === In this part we used the spherical head model distributed by the duneuro team ([[https://gitlab.dune-project.org/duneuro/duneuro-tests/-/blob/master/src/test_sphere_tet.msh|sphere]]). The following figure shows on the left the surface model and on right the tetrahedral mesh. {{attachment:DuneuroModelMesh.JPG||height="300",width="600"}} The model has three layers, the brain (inner skull), the outer skull and the scalp. === The source space : === For the source space or dipole position, we will use a realistic cortex distributed with the ICBM default subject of brainstorm. {{attachment:DuneuroModelMeshAnCortex.JPG||height="300",width="600"}} === The sensor model === Regarding the electrode's positions, we use the same position defined in this [[https://gitlab.dune-project.org/duneuro/duneuro-tests/-/blob/master/src/test_sphere_electrodes.txt|file]]. The total number of electrodes is 200 regularly distributed on the outer layer. {{attachment:DuneuroModelElectrode.JPG||height="300",width="350"}} === The forward model === Now, we have all the componenet of the model, we can start the process to compute the head model. For the EEG, brainstorm offers three methods. We perform these three computations on this model. Cond: 0.330 0.004 0.330 1. 3-shell sphere : best fitting sphere based pn scalp, then OK, less 1. OpenMeeg BEM : use the conductivity 0.33, 0.004 and 0.33, and keep the default options. 1. DUNEuro FEM : use the default option with the same conductivity value as the previous method. === Quantitave comparaison of the methods === In this tutorial we use the display of the leadfield vector to compare visually the performances of the the three methods. Right click on the headmodel, and then select 'View leadfield vector'. You can select more than headmodel in order to overlay the vectors. In this figure we show the 3-Shell vs OpenMeeg and 3-Shell vs Duneuro and all of these method. the BEM, FEM forward computation in basic model (realistic head model 3 layer) and also within in a spherical head model with an analytical solution. |
| Line 5: | Line 66: |
| Compare the performance of the mesh resolution vs the lead field approximation v10, v1, v0.1 and v0.01 {{attachment:femMeshModels.JPG||height="400",width="700"}} v0.1 and V0.01 are almost the same .... difference appears on v10 {{attachment:sourceVsMeshResolution.JPG||height="400",width="700"}} in spherical head model same observation on the lf view {{attachment:lfview_v01andv001.JPG||height="400",width="700"}} difference appears on the v1 and worse for v10 |
|
| Line 6: | Line 83: |
| In this example we compute the head model for the ICBM head model using the openmeeg method and then we generate the fem mesh for the same subject, then we compute the the leadfield using the Duneuro method. | In this example we compute the head model for the ICBM head model using the openmeeg method and then we generate the fem mesh for the same subject, then we compute the the leadfield using the Duneuro method v = 0.1 default value. MNorm {{attachment:sourcelocalisation.JPG||height="400",width="700"}} |
| Line 8: | Line 85: |
| {{attachment:sourcelocalisation.JPG||height="400",width="700"}} | sloreta {{attachment:sourcelocalisationloreta.JPG||height="400",width="700"}} {{attachment:sourcelocalisationloreta2.JPG||height="400",width="700"}} |
Introduction
In this tutorial, we compare the head models computation using the methods available in brainstorm.
At this time, the FEM implementation is available for the EEG and the MEG computation. We will integrate and test soon the FEM for sEEG and ECOG.
EEG within a spherical model
The volume conductor or the head model:
In this part we used the spherical head model distributed by the duneuro team (sphere). The following figure shows on the left the surface model and on right the tetrahedral mesh. 
The model has three layers, the brain (inner skull), the outer skull and the scalp.
The source space :
For the source space or dipole position, we will use a realistic cortex distributed with the ICBM default subject of brainstorm.
The sensor model
Regarding the electrode's positions, we use the same position defined in this file. The total number of electrodes is 200 regularly distributed on the outer layer.
The forward model
Now, we have all the componenet of the model, we can start the process to compute the head model. For the EEG, brainstorm offers three methods. We perform these three computations on this model.
Cond: 0.330 0.004 0.330
- 3-shell sphere : best fitting sphere based pn scalp, then OK, less
OpenMeeg BEM : use the conductivity 0.33, 0.004 and 0.33, and keep the default options.
- DUNEuro FEM : use the default option with the same conductivity value as the previous method.
Quantitave comparaison of the methods
In this tutorial we use the display of the leadfield vector to compare visually the performances of the the three methods.
Right click on the headmodel, and then select 'View leadfield vector'. You can select more than headmodel in order to overlay the vectors.
In this figure we show the 3-Shell vs OpenMeeg and 3-Shell vs Duneuro and all of these method.
the BEM, FEM forward computation in basic model (realistic head model 3 layer) and also within in a spherical head model with an analytical solution.
Describe femVSbem here.
Compare the performance of the mesh resolution vs the lead field approximation
v10, v1, v0.1 and v0.01 
v0.1 and V0.01 are almost the same ....
difference appears on v10
in spherical head model
same observation on the lf view 
difference appears on the v1 and worse for v10
EEG forward within 3 layers head model :
In this example we compute the head model for the ICBM head model using the openmeeg method and then we generate the fem mesh for the same subject, then we compute the the leadfield using the Duneuro method v = 0.1 default value. MNorm 
sloreta 
