|
Size: 1290
Comment: more specific
|
Size: 10269
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 1: | Line 1: |
| ## page was renamed from brainsuiteBDP | |
| Line 2: | Line 3: |
| = Estimation of brain's tissues anisortopy using Brainsuite Diffusion Pipline = | = FEM tensors estimation with BrainSuite = ''Authors: [[https://neuroimage.usc.edu/brainstorm/AboutUs/tmedani|Takfarinas Medani]], Francois Tadel, Anand Joshi and Richard Leahy'' '''[TUTORIAL UNDER WRITING: NOT READY FOR PUBLIC USE]''' |
| Line 4: | Line 9: |
'''[TUTORIAL UNDER REVISION/CORRECTION: NOT READY FOR PUBLIC USE]''' ''Authors: [[https://neuroimage.usc.edu/brainstorm/AboutUs/tmedani|Takfarinas Medani]], Francois Tadel, Anand Joshi'' |
|
| Line 11: | Line 12: |
| --(model using the '''finite element method''' ('''FEM'''). The FEM methods use the realistic volume mesh of the head generated from the segmentation of the MRI. The FEM models provides more accurate results than the spherical forward models, and more realistic geometry and tissue propriety than the BEM methods.)-- | refer to this page |
| Line 13: | Line 14: |
| --(The scope of this page is limited to a '''basic example''' (head model with 3 layers), more advanced options for head model generation and forward model options are discussed in the tutorial about FEM mesh generation. We assume that you have already followed the introduction tutorials (or at least the head modeling tutorial), we will not discuss the general principles of forward modeling here.)-- | [[https://neuroimage.usc.edu/brainstorm/Tutorials/SegBrainSuite?highlight=(Anand)|https://neuroimage.usc.edu/brainstorm/Tutorials/SegBrainSuite?highlight=%28anand%29]] The realistic tensors are estimated from the Diffusion-Weighted Images (DWI). For this purpose, Brainstorm calls internally the BrainSuite Diffusion Pipeline to compute the diffusion tensors on each brain voxel. Afterward, the Effective Medium Approach is applied to convert the diffusion tensors to the conductivity tensors. The following section shows the users how to do it from the graphical interface. Only the NIfTI is supported. All the diffusion data, including the DWI file and direction and the value of the gradient files, respectively the *.nii, the *.bval and the *.bvec are required. Ideally, these files should have the same name and saved in the same folder. <<TableOfContents(2,2)>> == Requirement == * You have already followed all the introduction tutorials * You have a working copy of Brainstorm installed on your computer == Brainsuite Installation == 1. Download the latest version of BrainSuite from http://www.brainsuite.org/download. 1. Install it on your computer by following the instructions in [[http://brainsuite.bmap.ucla.edu/quickstart/installation/|BrainSuite's quick start installation guide]]. 1. Note that you will be using BrainSuite Diffusion Pipeline(BDP), so you need to install a compatible [[http://www.mathworks.com/products/compiler/mcr|MATLAB Compiler Runtime]](last version). 1. Start BrainSuite to check if the installation (It's not required to open BrainSuite to run this tutorial). 1. The BrainSuite installation folder should be informed in the Brainstorm preferences ---- /!\ '''Edit conflict - other version:''' ---- ---- /!\ '''Edit conflict - other version:''' ---- {{attachment:brainsuiteInstalation.JPG||height="250",width="650"}} ---- /!\ '''Edit conflict - your version:''' ---- {{attachment:brainsuiteInstalation.JPG||height="250",width="650"}} ---- /!\ '''End of edit conflict''' ---- ---- /!\ '''Edit conflict - your version:''' ---- {{attachment:brainsuiteInstalation.JPG||height="250",width="650"}} ---- /!\ '''End of edit conflict''' ---- == Dataset == In this tutorial, we use the Brainsuite dataset example available on the Brainsuite tutorial webpage http://brainsuite.org/tutorials/. The T1w of the subject can be download from this link [[http://brainsuite.org/WebTutorialData/BrainSuiteTutorialSVReg_Sept16.zip|BrainSuiteTutorialSVReg_Sept16.zip]]. The DWI from this link [[http://brainsuite.org/WebTutorialData/DWI_Feb15.zip|DWI_Feb15.zip]] The first link contains the T1 MRI, with the name '2523412.nii.gz'. as in this figure, {{attachment:vierMriFolder.JPG||height="130",width="580"}} The second file is the DWI and should contain at least three files {{attachment:vierDWIFolder.JPG||height="150",width="580"}} Where the *.bval is a text file that contains the value of the gradient, and the *.bvec is also a text file that contains the orientation of the gradient. The nii.gz file is the NifTi file of the DWI where the images are stored. == Realistic condctivity tensors == First, you need to create a new subject in your protocol, let call it the 'BrainSuiteSubject'. Then import the T1 MRI of the subject and set the fiducials points as explained in the previous tutorial. {{attachment:mri.JPG||height="100",width="250"}} === Diffusion tensor generation (DTI) from DWI === Right-click on the subject and then select the item "Convert DWI to DTI". ---- /!\ '''Edit conflict - other version:''' ---- ---- /!\ '''Edit conflict - other version:''' ---- Then follow the popup windows by selecting the DWI, bval, and bval. If these files are in the same folder, Brainstorm will detect them automatically, otherwise, the user will be asked to browse the files one by one (as it's the case in this tutorial). ---- /!\ '''Edit conflict - your version:''' ---- Then follow the popup windows by selecting the DWI, bval, and bval. If these files are in the same folder, Brainstorm will detect them automatically, otherwise, the user will be asked to browse the files one by one (as it's the case in this tutorial). ---- /!\ '''End of edit conflict''' ---- ---- /!\ '''Edit conflict - your version:''' ---- Then follow the popup windows by selecting the DWI, bval, and bval. If these files are in the same folder, Brainstorm will detect them automatically, otherwise, the user will be asked to browse the files one by one (as it's the case in this tutorial). ---- /!\ '''End of edit conflict''' ---- In this step, Brainstorm calls the Brainsuite internally, and the diffusion tensors are computed. At the end of this process, a new node will appear in the Brainstorm database with the name 'DTI-EIT'. This name refers to, DTI: diffusion tensors images, and EIG for eigenvalue, since the eigenvalues and eigenvectors are computed at voxel and stored in Brainstorm database. {{attachment:mriAndDTI-EIG.JPG||height="100",width="250"}} === Conductivity tensor generation from DTI === The Effective Medium Approach is applied to convert the diffusion tensors to the conductivity tensors. This step requires the FEM mesh of the head model. You can generate the FEM head model from the MRI as explained on this page (add link). ---- /!\ '''Edit conflict - other version:''' ---- ---- /!\ '''Edit conflict - other version:''' ---- For the following, we used the SimNibs FEM mesh generation. The following figure shows the FEM mesh obtained from using the T1 MRI. Note that this mesh is obtained only from the T1, the use of the T2 is highly recommended if it's available. ---- /!\ '''Edit conflict - your version:''' ---- For the following, we used the SimNibs FEM mesh generation. The following figure shows the FEM mesh obtained from using the T1 MRI. Note that this mesh is obtained only from the T1, the use of the T2 is highly recommended if it's available. ---- /!\ '''End of edit conflict''' ---- ---- /!\ '''Edit conflict - your version:''' ---- For the following, we used the SimNibs FEM mesh generation. The following figure shows the FEM mesh obtained from using the T1 MRI. Note that this mesh is obtained only from the T1, the use of the T2 is highly recommended if it's available. ---- /!\ '''End of edit conflict''' ---- The FEM head model to use for tensors should be selected and highlighted with the green color (double click on the FEM mesh node to select it) When this is done, then right-click on the subject > Convert DWI to DTI, ---- /!\ '''Edit conflict - other version:''' ---- ---- /!\ '''Edit conflict - other version:''' ---- Brainstorm will load the available tissue in the FEM head model and the following windows appear. Select the WM anisotropy and kee all the other tissues as isotropic. The process of conversion from DWI to Conductivity tensors uses the EMA, furthermore, brainstorm proposes the option to use the adaptative EMA with the volume constraint option [ref]. In this example, we select the EMA with the VC. The process will take around 10 min, and then the FEM tensors are computed and stored in the FEM structure. [explain how it is organized and how to use it outside brainstorm ] ---- /!\ '''Edit conflict - your version:''' ---- Brainstorm will load the available tissue in the FEM head model and the following windows appear. Select the WM anisotropy and kee all the other tissues as isotropic. The process of conversion from DWI to Conductivity tensors uses the EMA, furthermore, brainstorm proposes the option to use the adaptative EMA with the volume constraint option [ref]. In this example, we select the EMA with the VC. The process will take around 10 min, and then the FEM tensors are computed and stored in the FEM structure. [explain how it is organized and how to use it outside brainstorm ] ---- /!\ '''End of edit conflict''' ---- ---- /!\ '''Edit conflict - your version:''' ---- Brainstorm will load the available tissue in the FEM head model and the following windows appear. Select the WM anisotropy and kee all the other tissues as isotropic. The process of conversion from DWI to Conductivity tensors uses the EMA, furthermore, brainstorm proposes the option to use the adaptative EMA with the volume constraint option [ref]. In this example, we select the EMA with the VC. The process will take around 10 min, and then the FEM tensors are computed and stored in the FEM structure. [explain how it is organized and how to use it outside brainstorm ] ---- /!\ '''End of edit conflict''' ---- === Display the tensors === == Artificial/simulated conductivity tensors == ---- /!\ '''Edit conflict - other version:''' ---- ---- /!\ '''Edit conflict - other version:''' ---- In the case where the DWI is not available, or the users desire to evaluate the effect of the conductivity change on the model, the artificial conductivity can be used. ---- /!\ '''Edit conflict - your version:''' ---- In the case where the DWI is not available, or the users desire to evaluate the effect of the conductivity change on the model, the artificial conductivity can be used. ---- /!\ '''End of edit conflict''' ---- ---- /!\ '''Edit conflict - your version:''' ---- In the case where the DWI is not available, or the users desire to evaluate the effect of the conductivity change on the model, the artificial conductivity can be used. ---- /!\ '''End of edit conflict''' ---- Two approaches are integrated within Brainstorm. Either the == References == |
FEM tensors estimation with BrainSuite
Authors: Takfarinas Medani, Francois Tadel, Anand Joshi and Richard Leahy
[TUTORIAL UNDER WRITING: NOT READY FOR PUBLIC USE]
Describe brainsuite here : here we will describe the process of the brain tissues anisotrpy estimation and the different functions that brainstorm offers.
This tutorial explains how to use Brainsuite to estimate the anisotropy of the brain tissues.
refer to this page
https://neuroimage.usc.edu/brainstorm/Tutorials/SegBrainSuite?highlight=%28anand%29
The realistic tensors are estimated from the Diffusion-Weighted Images (DWI). For this purpose, Brainstorm calls internally the BrainSuite Diffusion Pipeline to compute the diffusion tensors on each brain voxel. Afterward, the Effective Medium Approach is applied to convert the diffusion tensors to the conductivity tensors. The following section shows the users how to do it from the graphical interface.
Only the NIfTI is supported. All the diffusion data, including the DWI file and direction and the value of the gradient files, respectively the *.nii, the *.bval and the *.bvec are required. Ideally, these files should have the same name and saved in the same folder.
Contents
Requirement
- You have already followed all the introduction tutorials
- You have a working copy of Brainstorm installed on your computer
Brainsuite Installation
Download the latest version of BrainSuite from http://www.brainsuite.org/download.
Install it on your computer by following the instructions in BrainSuite's quick start installation guide.
Note that you will be using BrainSuite Diffusion Pipeline(BDP), so you need to install a compatible MATLAB Compiler Runtime(last version).
Start BrainSuite to check if the installation (It's not required to open BrainSuite to run this tutorial).
The BrainSuite installation folder should be informed in the Brainstorm preferences
![/!\](/moin_static198/brainstorm1/img/alert.png "/!\"){kind=small}
Edit conflict - other version:
Edit conflict - other version:
![[ATTACH]](/brainstorm/Tutorials/FemTensors?action=AttachFile&do=upload_form&ticket=00673929b0.a8696e833e00c26ad8decd2759bc094e3d7a121f&target=brainsuiteInstalation.JPG "[ATTACH]"){kind=small}
Edit conflict - your version:
End of edit conflict
Edit conflict - your version:
End of edit conflict
Dataset
In this tutorial, we use the Brainsuite dataset example available on the Brainsuite tutorial webpage
http://brainsuite.org/tutorials/. The T1w of the subject can be download from this link BrainSuiteTutorialSVReg_Sept16.zip. The DWI from this link DWI_Feb15.zip
The first link contains the T1 MRI, with the name '2523412.nii.gz'. as in this figure,
![[ATTACH]](/brainstorm/Tutorials/FemTensors?action=AttachFile&do=upload_form&ticket=00673929b0.a8696e833e00c26ad8decd2759bc094e3d7a121f&target=vierMriFolder.JPG "[ATTACH]"){kind=small}
The second file is the DWI and should contain at least three files
![[ATTACH]](/brainstorm/Tutorials/FemTensors?action=AttachFile&do=upload_form&ticket=00673929b0.a8696e833e00c26ad8decd2759bc094e3d7a121f&target=vierDWIFolder.JPG "[ATTACH]"){kind=small}
Where the *.bval is a text file that contains the value of the gradient, and the *.bvec is also a text file that contains the orientation of the gradient. The nii.gz file is the NifTi file of the DWI where the images are stored.
Realistic condctivity tensors
First, you need to create a new subject in your protocol, let call it the 'BrainSuiteSubject'. Then import the T1 MRI of the subject and set the fiducials points as explained in the previous tutorial.
![[ATTACH]](/brainstorm/Tutorials/FemTensors?action=AttachFile&do=upload_form&ticket=00673929b0.a8696e833e00c26ad8decd2759bc094e3d7a121f&target=mri.JPG "[ATTACH]"){kind=small}
Diffusion tensor generation (DTI) from DWI
Right-click on the subject and then select the item "Convert DWI to DTI".
Edit conflict - other version:
Edit conflict - other version:
Then follow the popup windows by selecting the DWI, bval, and bval. If these files are in the same folder, Brainstorm will detect them automatically, otherwise, the user will be asked to browse the files one by one (as it's the case in this tutorial).
Edit conflict - your version:
Then follow the popup windows by selecting the DWI, bval, and bval. If these files are in the same folder, Brainstorm will detect them automatically, otherwise, the user will be asked to browse the files one by one (as it's the case in this tutorial).
End of edit conflict
Edit conflict - your version:
Then follow the popup windows by selecting the DWI, bval, and bval. If these files are in the same folder, Brainstorm will detect them automatically, otherwise, the user will be asked to browse the files one by one (as it's the case in this tutorial).
End of edit conflict
In this step, Brainstorm calls the Brainsuite internally, and the diffusion tensors are computed. At the end of this process, a new node will appear in the Brainstorm database with the name 'DTI-EIT'. This name refers to, DTI: diffusion tensors images, and EIG for eigenvalue, since the eigenvalues and eigenvectors are computed at voxel and stored in Brainstorm database.
![[ATTACH]](/brainstorm/Tutorials/FemTensors?action=AttachFile&do=upload_form&ticket=00673929b0.a8696e833e00c26ad8decd2759bc094e3d7a121f&target=mriAndDTI-EIG.JPG "[ATTACH]"){kind=small}
Conductivity tensor generation from DTI
The Effective Medium Approach is applied to convert the diffusion tensors to the conductivity tensors. This step requires the FEM mesh of the head model. You can generate the FEM head model from the MRI as explained on this page (add link).
Edit conflict - other version:
Edit conflict - other version:
For the following, we used the SimNibs FEM mesh generation. The following figure shows the FEM mesh obtained from using the T1 MRI. Note that this mesh is obtained only from the T1, the use of the T2 is highly recommended if it's available.
Edit conflict - your version:
For the following, we used the SimNibs FEM mesh generation. The following figure shows the FEM mesh obtained from using the T1 MRI. Note that this mesh is obtained only from the T1, the use of the T2 is highly recommended if it's available.
End of edit conflict
Edit conflict - your version:
For the following, we used the SimNibs FEM mesh generation. The following figure shows the FEM mesh obtained from using the T1 MRI. Note that this mesh is obtained only from the T1, the use of the T2 is highly recommended if it's available.
End of edit conflict
The FEM head model to use for tensors should be selected and highlighted with the green color (double click on the FEM mesh node to select it)
When this is done, then right-click on the subject > Convert DWI to DTI,
Edit conflict - other version:
Edit conflict - other version:
Brainstorm will load the available tissue in the FEM head model and the following windows appear.
Select the WM anisotropy and kee all the other tissues as isotropic.
The process of conversion from DWI to Conductivity tensors uses the EMA, furthermore, brainstorm proposes the option to use the adaptative EMA with the volume constraint option [ref]. In this example, we select the EMA with the VC.
The process will take around 10 min, and then the FEM tensors are computed and stored in the FEM structure. [explain how it is organized and how to use it outside brainstorm ]
Edit conflict - your version:
Brainstorm will load the available tissue in the FEM head model and the following windows appear.
Select the WM anisotropy and kee all the other tissues as isotropic.
The process of conversion from DWI to Conductivity tensors uses the EMA, furthermore, brainstorm proposes the option to use the adaptative EMA with the volume constraint option [ref]. In this example, we select the EMA with the VC.
The process will take around 10 min, and then the FEM tensors are computed and stored in the FEM structure. [explain how it is organized and how to use it outside brainstorm ]
End of edit conflict
Edit conflict - your version:
Brainstorm will load the available tissue in the FEM head model and the following windows appear.
Select the WM anisotropy and kee all the other tissues as isotropic.
The process of conversion from DWI to Conductivity tensors uses the EMA, furthermore, brainstorm proposes the option to use the adaptative EMA with the volume constraint option [ref]. In this example, we select the EMA with the VC.
The process will take around 10 min, and then the FEM tensors are computed and stored in the FEM structure. [explain how it is organized and how to use it outside brainstorm ]
End of edit conflict
Display the tensors
Artificial/simulated conductivity tensors
Edit conflict - other version:
Edit conflict - other version:
In the case where the DWI is not available, or the users desire to evaluate the effect of the conductivity change on the model, the artificial conductivity can be used.
Edit conflict - your version:
In the case where the DWI is not available, or the users desire to evaluate the effect of the conductivity change on the model, the artificial conductivity can be used.
End of edit conflict
Edit conflict - your version:
In the case where the DWI is not available, or the users desire to evaluate the effect of the conductivity change on the model, the artificial conductivity can be used.
End of edit conflict
Two approaches are integrated within Brainstorm. Either the