|
Size: 2076
Comment:
|
Size: 9033
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 1: | Line 1: |
| = Realistic head model: FEM with Duneuro = ''Authors: Takfarinas MEDANI,'' |
= Realistic head model: FEM with DUNEuro = '''[TUTORIAL UNDER DEVELOPMENT: NOT READY FOR PUBLIC USE]''' |
| Line 4: | Line 4: |
| This forward model uses a finite element element method (FEM) and was developed by the '''...(l''''''ink)'''. It uses the volume mesh of the realistic head model. The goal of this forward solution is mostly for '''EEG users''', to provide more accurate results than the spherical models and more complex geometry than the BEM. It is not necessary for MEG users, as the "overlapping spheres" method gives similar results but much faster. This method is illustrated using the tutorial [[http://neuroimage.usc.edu/brainstorm/Tutorials/Epilepsy|EEG and epilepsy]] or similar data ('''todo'''). | ''Authors: Takfarinas, Juan, François ...Sophie, Maria, Christian, Carsten, John, Richard ? '' |
| Line 6: | Line 6: |
| refe to these point : http://duneuro.org/ | This tutorial explains how to use brainstorm-dueneuro to compute the head model using the finite element method, implemented in the Duneuro software. |
| Line 8: | Line 8: |
| and to the documentation from here : | This forward model uses a finite element element method (FEM) implemented within the [[http://duneuro.org/|Duneuro]] software. More documentation could be found on this [[https://gitlab.dune-project.org/duneuro/duneuro/wikis/home|gitlab page]]. |
| Line 10: | Line 10: |
| https://gitlab.dune-project.org/duneuro/duneuro/wikis/home | The FEM methods use the realistic volume mesh of the head model genertaed from the segmentation of the magnetic reasonance images (MRI). The goal of the forward solution is, mostly for '''EEG users (and for sEEG/ECOG)''', to provide more accurate results than the spherical models and more realistic geometry than the BEM. |
| Line 12: | Line 12: |
| Here we will describe the FEM and run some examples | In this tutorial we will describe the software requirement to use the DUNEuro computation from brainstrom and examples on the FEM computation through the graphical interface. |
| Line 14: | Line 14: |
| == DUNEuro == [[http://duneuro.org/|DUNEuro]] is an open-source C++ software library for solving partial differential equations (PDE) in neurosciences using mesh bases methods. It is based on the [[https://www.dune-project.org/|DUNE library]] and its main features that include solving the electroencephalography (EEG) and magnetoencephalography (MEG) forward problem and providing simulations for brain stimulation. |
|
| Line 15: | Line 17: |
== Duneuro == |
{{attachment:duneuroFromDune.png|height="100",width="350"}} <<BR>><<BR>> |
| Line 19: | Line 20: |
| Duneuro is build on top of DUNE Library, at this time (sept 2019) these libraries work only on Linux systems. However we are able to run generate binaries for windows and can be easily used from Matlab for all platforms. | === FEM Computation & DUNEuro === The FEM computation is done with DUNEuro. From Brainstrom, the DUNEuro binaries files are distributed throgh the [[https://github.com/brainstorm-tools/bst-duneuro|bst-duneuro toolbox]]. These binaries are adapted and tested for Windows, MAC and Ubunto 64bit systems. Therefore you don't need to install any extra software package to run the FEM computation'''.''' |
| Line 21: | Line 23: |
| However if you want to use the source code of Duneuro you can visite: http://duneuro.org/ | '''<<TAG(Advanced)>>''' |
| Line 23: | Line 25: |
| We also put here : | DUNEuro is build on top of DUNE Library, the source code works only on Linux operating systems. However, we are able to generate the binaries for the main platforms (windows 64, Linux and Mac), therefor it can be easily used from Matlab without the need to install and compile the tidiuos Dune and Duneuro libraries on your computer. |
| Line 25: | Line 27: |
| The main steps you need in order to compile for windows are listed here : | We have also develloped an alternative way to install and generates the binaries for the main platforms. For the advanced user users, the main steps you need in order to compile DUNEuro for are explained [[https://github.com/brainstorm-tools/bst-duneuro|here]]. |
| Line 27: | Line 29: |
| https://github.com/svdecomposer/brainstorm-duneuro | === FEM Mesh generation & processing === In order to used the FEM to compute the electromagnetic field, the volume mesh of the head is required. |
| Line 29: | Line 32: |
| We distribute the binary files for Windows/Ubunto 64bit systems, you don't need to installan extra Microsoft software package to run '''Duneuro''' | Brain |
| Line 31: | Line 34: |
| Microsoft Visual C++ 2010 Redistributable Package (x64). |
The FEM mesh visualisation and some of the mesh processing integrated with brainstorm requires the installation of the iso2mesh toolbox. If you want to work on offline, you can download it and add it to your matlab path. If it's not installed, Brainstrom will download and istall this toolbox when it needed. |
| Line 38: | Line 39: |
| === SIMNIBS === https://simnibs.github.io/simnibs/build/html/index.html |
|
| Line 39: | Line 42: |
| === SimNIBS Installation === Please follow the instructions on this webapge: |
|
| Line 40: | Line 45: |
| https://simnibs.github.io/simnibs/build/html/installation/simnibs_installer.html | |
| Line 41: | Line 47: |
| === Volume mesh generation from Brainstorm === To generate the volume mesh head model we use the "headreco" : https://simnibs.github.io/simnibs/build/html/documentation/command_line/headreco.html This function is part of the SimNibs software: https://simnibs.github.io/simnibs/build/html/index.html === headreco === The headreco function is fully integrated to brainstorm. With this option, brainstorm can reconstructs a tetrahedral head mesh from T1- and T2-weighted structural MR images. It runs also with only a T1w image, but it will achieve more reliable skull segmentations when a T2w image is supplied. === Surface mesh generation from volume mesh === |
|
| Line 42: | Line 61: |
| In this part you can generate your FEM mesh from surfaces that you can get fron the segmentation software (brainSuite, FreeSurfer ....). | |
| Line 43: | Line 63: |
| === Volume generation from T1 MRI data === | This process will |
| Line 45: | Line 65: |
| - merge the surfaces, | |
| Line 46: | Line 67: |
| - check the self intersecting | |
| Line 47: | Line 69: |
| - fixe the size of the mesh | |
| Line 48: | Line 71: |
| - generate the volum mesh | |
| Line 49: | Line 73: |
| == Forward model == | - visual checking ... |
| Line 51: | Line 75: |
| === Head model === === Source model === === Advanced paramaters === == Additional documentation == == Reported Errors & alternative solution == |
- TODO : may be we can add some mesh auqlity measures ?? |
| Line 75: | Line 79: |
| === Volume generation from T1/T2 MRI data === You can also generate your own FEM head model and then load it to brainstorm. However the automatic head model generation from from imaging techniques are not accurate and most of the time visual checking are needed and manual correction are required. ==> this depends lagely on the quality of the T1/T2 MRI image(https://simnibs.github.io/simnibs/build/html/tutorial/head_meshing.html). This step is based on the "roast" toolbox (link to roast : https://github.com/andypotatohy/roast ) that we adapted for the MEEG forward computation. If you want to generate your own FEM head model from an MRI, you will need to download these file (link), then run the bst process as explained here. * f there is a MRI file with the string "T2" in the subject anatomy folder, it will use it * Otherwise, if you select explicitly two MRI files with CTRL+Click, it will use the first one as the T1 and the second one as the T2 (this needs to be documented in the tutorial) === FEM Head model generation with SimNibs === This method used the SimNibs software. So to call this process, you need to download and install the SimNibs software, the process of the installation is explained in the SimNibs webpage : https://simnibs.github.io/simnibs/build/html/installation/simnibs_installer.html. When you have installed SimNibs, Brainstorm can call the main function used for the mesh generation frm the main graphical interface. Depemding on your computer performances, this process will take between 2 to 5 hours. We highly recommend to close all other running process and application on our computer in order to speed this process. - Explain here the main steps with screenshots : 1- Create new subject within the current protocole 2- Load the T1 of the subject to the brainstorm database. 3- Associate a T2 mri to the subject if it's available (this is better for csf/skull/scalp segmentation) 4- Right click on the subject, select the "Generate FEM mesh" . Select "SIMNIBS", and choose "Tetrahedral element" and keep the other options to the default value. 5- When this process is finished, a new node will appear in the data base, which hase he name "FEM xxxV, (simNibs, yLayers)". This is the FEM mesh model generated from the T1. === FEM Head model template === - Load the FEM volumic mesh (template created from ICBM T1 MRI using SimNibs) - Load the surface mesh (template created also from ICBM using ICBM ) and then generates the volume mesh (either tetra or hexa) by calling the tetgen process cia iso2mesh toolbox (if hexa are desired, the tetra mesh will be converted to hexa ... ) https://github.com/brainstorm-tools/brainstorm3/issues/185#issuecomment-576749612 === Head model based on the level set approach === TODO and Validate if users want to use freesurfer and simnibs, the priority is given to freesurfer, then load the fs files, then call simnibs headmodel generation from t1/t2. == FEM computation and interface to DUNEuro == === Head model === Number of layers, conductivity value, isotropy/anisotropy/ mesh resolution/ === Electrode model === Check the position of the electrodes and align to the head model (projection if needed) === Source model === Similarly to the spherical nad BEM head model, the source position are defined on the cortex surface vertices. We can either define a contraned or not constrained orientation. However, for the FEM model, more paramters could be tuned for the source model. Choice of the source model : PI, Venant, Subtraction, Whitney Panel of the options choice that the user can select. (other wise we will set to default ) '''<<TAG(Advanced)>>''' === Advanced paramaters === - Solver parameters - Electrodes projection - maybe explain here the relevant option of the mini file ?? == Additional documentation == refer to : http://duneuro.org/ https://www.dune-project.org/ https://simnibs.github.io/simnibs/build/html/index.html == Reported Errors & alternative solution == '''<<TAG(Advanced)>>''' simnibs pblm : https://simnibs.github.io/simnibs/build/html/installation/throubleshooting.html |
|
| Line 76: | Line 164: |
| == License == == Reference == |
|
| Line 77: | Line 167: |
| == License == | == Additional documentation == https://github.com/brainstorm-tools/bst-duneuro/issues/1 |
| Line 79: | Line 170: |
| https://github.com/brainstorm-tools/brainstorm3/issues/242 | |
| Line 80: | Line 172: |
== Reference == |
https://github.com/brainstorm-tools/brainstorm3/issues/185 |
Realistic head model: FEM with DUNEuro
[TUTORIAL UNDER DEVELOPMENT: NOT READY FOR PUBLIC USE]
Authors: Takfarinas, Juan, François ...Sophie, Maria, Christian, Carsten, John, Richard ?
This tutorial explains how to use brainstorm-dueneuro to compute the head model using the finite element method, implemented in the Duneuro software.
This forward model uses a finite element element method (FEM) implemented within the Duneuro software. More documentation could be found on this gitlab page.
The FEM methods use the realistic volume mesh of the head model genertaed from the segmentation of the magnetic reasonance images (MRI). The goal of the forward solution is, mostly for EEG users (and for sEEG/ECOG), to provide more accurate results than the spherical models and more realistic geometry than the BEM.
In this tutorial we will describe the software requirement to use the DUNEuro computation from brainstrom and examples on the FEM computation through the graphical interface.
DUNEuro
DUNEuro is an open-source C++ software library for solving partial differential equations (PDE) in neurosciences using mesh bases methods. It is based on the DUNE library and its main features that include solving the electroencephalography (EEG) and magnetoencephalography (MEG) forward problem and providing simulations for brain stimulation.
Requirements
FEM Computation & DUNEuro
The FEM computation is done with DUNEuro. From Brainstrom, the DUNEuro binaries files are distributed throgh the bst-duneuro toolbox. These binaries are adapted and tested for Windows, MAC and Ubunto 64bit systems. Therefore you don't need to install any extra software package to run the FEM computation.
DUNEuro is build on top of DUNE Library, the source code works only on Linux operating systems. However, we are able to generate the binaries for the main platforms (windows 64, Linux and Mac), therefor it can be easily used from Matlab without the need to install and compile the tidiuos Dune and Duneuro libraries on your computer.
We have also develloped an alternative way to install and generates the binaries for the main platforms. For the advanced user users, the main steps you need in order to compile DUNEuro for are explained here.
FEM Mesh generation & processing
In order to used the FEM to compute the electromagnetic field, the volume mesh of the head is required.
Brain
The FEM mesh visualisation and some of the mesh processing integrated with brainstorm requires the installation of the iso2mesh toolbox. If you want to work on offline, you can download it and add it to your matlab path. If it's not installed, Brainstrom will download and istall this toolbox when it needed.
FEM surfaces / Volume generation
SIMNIBS
https://simnibs.github.io/simnibs/build/html/index.html
SimNIBS Installation
Please follow the instructions on this webapge:
https://simnibs.github.io/simnibs/build/html/installation/simnibs_installer.html
Volume mesh generation from Brainstorm
To generate the volume mesh head model we use the "headreco" :
https://simnibs.github.io/simnibs/build/html/documentation/command_line/headreco.html
This function is part of the SimNibs software:
https://simnibs.github.io/simnibs/build/html/index.html
headreco
The headreco function is fully integrated to brainstorm. With this option, brainstorm can reconstructs a tetrahedral head mesh from T1- and T2-weighted structural MR images. It runs also with only a T1w image, but it will achieve more reliable skull segmentations when a T2w image is supplied.
Surface mesh generation from volume mesh
Volume generation from surface files
In this part you can generate your FEM mesh from surfaces that you can get fron the segmentation software (brainSuite, FreeSurfer ....).
This process will
- merge the surfaces,
- check the self intersecting
- fixe the size of the mesh
- generate the volum mesh
- visual checking ...
- TODO : may be we can add some mesh auqlity measures ??
Volume generation from T1/T2 MRI data
You can also generate your own FEM head model and then load it to brainstorm. However the automatic head model generation from from imaging techniques are not accurate and most of the time visual checking are needed and manual correction are required.
==> this depends lagely on the quality of the T1/T2 MRI image(https://simnibs.github.io/simnibs/build/html/tutorial/head_meshing.html).
This step is based on the "roast" toolbox (link to roast : https://github.com/andypotatohy/roast
) that we adapted for the MEEG forward computation. If you want to generate your own FEM head model from an MRI, you will need to download these file (link), then run the bst process as explained here.
- f there is a MRI file with the string "T2" in the subject anatomy folder, it will use it
- Otherwise, if you select explicitly two MRI files with CTRL+Click, it will use the first one as the T1 and the second one as the T2 (this needs to be documented in the tutorial)
FEM Head model generation with SimNibs
This method used the SimNibs software. So to call this process, you need to download and install the SimNibs software, the process of the installation is explained in the SimNibs webpage : https://simnibs.github.io/simnibs/build/html/installation/simnibs_installer.html.
When you have installed SimNibs, Brainstorm can call the main function used for the mesh generation frm the main graphical interface. Depemding on your computer performances, this process will take between 2 to 5 hours. We highly recommend to close all other running process and application on our computer in order to speed this process.
- Explain here the main steps with screenshots :
1- Create new subject within the current protocole
2- Load the T1 of the subject to the brainstorm database.
3- Associate a T2 mri to the subject if it's available (this is better for csf/skull/scalp segmentation)
4- Right click on the subject, select the "Generate FEM mesh"
- Select "SIMNIBS", and choose "Tetrahedral element" and keep the other options to the default value.
5- When this process is finished, a new node will appear in the data base, which hase he name "FEM xxxV, (simNibs, yLayers)". This is the FEM mesh model generated from the T1.
FEM Head model template
- Load the FEM volumic mesh (template created from ICBM T1 MRI using SimNibs)
- Load the surface mesh (template created also from ICBM using ICBM ) and then generates the volume mesh (either tetra or hexa) by calling the tetgen process cia iso2mesh toolbox (if hexa are desired, the tetra mesh will be converted to hexa ... )
https://github.com/brainstorm-tools/brainstorm3/issues/185#issuecomment-576749612
Head model based on the level set approach
TODO and Validate
if users want to use freesurfer and simnibs, the priority is given to freesurfer, then load the fs files, then call simnibs headmodel generation from t1/t2.
FEM computation and interface to DUNEuro
Head model
Number of layers, conductivity value, isotropy/anisotropy/ mesh resolution/
Electrode model
Check the position of the electrodes and align to the head model (projection if needed)
Source model
Similarly to the spherical nad BEM head model, the source position are defined on the cortex surface vertices. We can either define a contraned or not constrained orientation.
However, for the FEM model, more paramters could be tuned for the source model.
Choice of the source model : PI, Venant, Subtraction, Whitney
Panel of the options choice that the user can select. (other wise we will set to default )
Advanced paramaters
- Solver parameters
- Electrodes projection
- maybe explain here the relevant option of the mini file ??
Additional documentation
refer to :
https://simnibs.github.io/simnibs/build/html/index.html
Reported Errors & alternative solution
simnibs pblm :
https://simnibs.github.io/simnibs/build/html/installation/throubleshooting.html
The MEEG forward problem with the FEM
License
Reference
Additional documentation
https://github.com/brainstorm-tools/bst-duneuro/issues/1