-
Notifications
You must be signed in to change notification settings - Fork 76
DWI_Connectome
The dwi-connectome
pipeline computes a weighted graph encoding anatomical connections between a set of brain regions from corrected DWI datasets.
To aim that, it relies on the MRtrix3 [Tournier et al., 2019] software to compute the constrained spherical deconvolution diffusion model, perform probabilistic tractography and computes a connectome using the Desikan & Destrieux atlases from FreeSurfer.
You need to preprocess your DWI data and run the t1-freesurfer
pipeline on your T1-weighted MRI images prior to running this pipeline.
If you only installed the core of Clinica, this pipeline needs the installation of FSL 6.0 and MRtrix3 on your computer.
You can find how to install these software packages on the third-party page.
The pipeline can be run with the following command line:
clinica run dwi-connectome <caps_directory>
where:
-
caps_directory
is the input/output folder containing the results in a CAPS hierarchy.
If you want to run the pipeline on a subset of your CAPS dataset, you can use the -tsv
flag to specify in a TSV file the participants belonging to your subset.
!!! note "Number of streamlines (--n_tracks
parameters)"
The quality of the tractography and, as a result, the connectome mainly depends on the number of streamlines you can generate (the more the better). However, increasing the number of streamlines increases the need for computational resources and space to store the results. On default, 1 million streamlines are generated which represents 1 hour of computation time.
Results are stored in the following folder of the CAPS hierarchy: subjects/sub-<participant_label>/ses-<session_label>/dwi/connectome_based_processing/
.
The main output files are:
-
<source_file>_space-{b0|T1w}_model-CSD_diffmodel.nii.gz
: Constrained spherical deconvolution (CSD) diffusion model. -
<source_file>_space-{b0|T1w}_model-CSD_tractography.tck
: The whole-brain tractography. -
<source_file>_space-{b0|T1w}_model-CSD_parcellation-{desikan|destrieux}_connectivity.tsv
: The connectivity matrix based on the Desikan or Destrieux parcellation.
!!! note "Atlases available for the Connectome-based processing pipeline:" - Desikan [Desikan et al., 2006]: This atlas is a subdivision of the cerebral cortex into gyri and contains 34 regions per hemisphere. It was built using a dataset of 40 MRI scans from which 34 cortical ROIs were manually identified in each of the individual hemispheres. - Destrieux [Destrieux et al., 2010]: This atlas is a subdivision of the cerebral cortex into gyri and sulci, and contains 74 regions per hemisphere. It was built on anatomical MRI of 24 healthy subjects from which 74 cortical ROIs were manually identified in each of the individual hemispheres.
!!! note The full list of output files can be found in the The ClinicA Processed Structure (CAPS) specifications.
!!! cite "Example of paragraph for the dwi-connectome
pipeline"
These results have been obtained using the dwi-connectome
pipeline of Clinica [Routier et al] relying on the MRtrix3 [Tournier et al., 2019] software package. Fiber orientation distributions (FOD) at highly anisotropic voxels (FA >0.7) was computed to determine the response function, which was used for constrained spherical deconvolution to accurately estimate the FOD [Tournier et al., 2007]. Then, <n_tracks>
fibers with a probabilistic tracking algorithm [Tournier et al., 2010] were generated. Default parameters included minimum length 20 mm, a step size of 0.2 mm, minimum radius of curvature of 1 mm and FOD cutoff of 0.1. All voxels in the 1-mm dilated white-matter mask were used as seeds and the tracking procedure was stopped if a fiber reached a voxel outside the mask or if a stopping criterion was met (high fiber curvature or low FOD). Finally, the connectome is estimated by counting the number of tracks connecting each pair of nodes according to the [Desikan|Destrieux] parcellation.
!!! tip Easily access the papers cited on this page on Zotero.
- You can use the Clinica Google Group to ask for help!
- Report an issue on GitHub.