Tractography passes the test: Results from the diffusion-simulated connectivity (disco) challenge.

Estimating structural connectivity from diffusion-weighted magnetic resonance imaging is a challenging task, partly due to the presence of false-positive connections and the misestimation of connection weights. Building on previous efforts, the MICCAI-CDMRI Diffusion-Simulated Connectivity (DiSCo) challenge was carried out to evaluate state-of-the-art connectivity methods using novel large-scale numerical phantoms. The diffusion signal for the phantoms was obtained from Monte Carlo simulations. The results of the challenge suggest that methods selected by the 14 teams participating in the challenge can provide high correlations between estimated and ground-truth connectivity weights, in complex numerical environments. Additionally, the methods used by the participating teams were able to accurately identify the binary connectivity of the numerical dataset. However, specific false positive and false negative connections were consistently estimated across all methods. Although the challenge dataset doesn't capture the complexity of a real brain, it provided unique data with known macrostructure and microstructure ground-truth properties to facilitate the development of connectivity estimation methods.

Identifying subcortical connectivity during brain tumor surgery: a multimodal study.

Bipolar direct electrical stimulation (DES) of an awake patient is the reference technique for identifying brain structures to achieve maximal safe tumor resection. Unfortunately, DES cannot be performed in all cases. Alternative surgical tools are, therefore, needed to aid identification of subcortical connectivity during brain tumor removal. In this pilot study, we sought to (i) evaluate the combined use of evoked potential (EP) and tractography for identification of white matter (WM) tracts under the functional control of DES, and (ii) provide clues to the electrophysiological effects of bipolar stimulation on neural pathways. We included 12 patients (mean age of 38.4 years) who had had a dMRI-based tractography and a functional brain mapping under awake craniotomy for brain tumor removal. Electrophysiological recordings of subcortical evoked potentials (SCEPs) were acquired during bipolar low frequency (2 Hz) stimulation of the WM functional sites identified during brain mapping. SCEPs were successfully triggered in 11 out of 12 patients. The median length of the stimulated fibers was 43.24 ± 19.55 mm, belonging to tracts of median lengths of 89.84 ± 24.65 mm. The electrophysiological (delay, amplitude, and speed of propagation) and structural (number and lengths of streamlines, and mean fractional anisotropy) measures were correlated. In our experimental conditions, SCEPs were essentially limited to a subpart of the bundles, suggesting a selectivity of action of the DES on the brain networks. Correlations between functional, structural, and electrophysiological measures portend the combined use of EPs and tractography as a potential intraoperative tool to achieve maximum safe resection in brain tumor surgery.

DORIS: A diffusion MRI-based 10 tissue class deep learning segmentation algorithm tailored to improve anatomically-constrained tractography.

Modern tractography algorithms such as anatomically-constrained tractography (ACT) are based on segmentation maps of white matter (WM), gray matter (GM), and cerebrospinal fluid (CSF). These maps are generally estimated from a T1-weighted (T1w) image and then registered in diffusion weighted images (DWI) space. Registration of T1w to diffusion space and partial volume estimation are challenging and rarely voxel-perfect. Diffusion-based segmentation would, thus, potentially allow not to have higher quality anatomical priors injected in the tractography process. On the other hand, even if FA-based tractography is possible without T1 registration, the literature shows that this technique suffers from multiple issues such as holes in the tracking mask and a high proportion of generated broken and anatomically implausible streamlines. Therefore, there is an important need for a tissue segmentation algorithm that works directly in the native diffusion space. We propose DORIS, a DWI-based deep learning segmentation algorithm. DORIS outputs 10 different tissue classes including WM, GM, CSF, ventricles, and 6 other subcortical structures (putamen, pallidum, hippocampus, caudate, amygdala, and thalamus). DORIS was trained and validated on a wide range of subjects, including 1,000 individuals from 22 to 90 years old from clinical and research DWI acquisitions, from 5 public databases. In the absence of a "true" ground truth in diffusion space, DORIS used a silver standard strategy from Freesurfer output registered onto the DWI. This strategy is extensively evaluated and discussed in the current study. Segmentation maps provided by DORIS are quantitatively compared to Freesurfer and FSL-fast and the impacts on tractography are evaluated. Overall, we show that DORIS is fast, accurate, and reproducible and that DORIS-based tractograms produce bundles with a longer mean length and fewer anatomically implausible streamlines.

A Riemannian Revisiting of Structure-Function Mapping Based on Eigenmodes.

Understanding the link between brain structure and function may not only improve our knowledge of brain organization, but also lead to better quantification of pathology. To quantify this link, recent studies have attempted to predict the brain's functional connectivity from its structural connectivity. However, functional connectivity matrices live in the Riemannian manifold of the symmetric positive definite space and a specific attention must be paid to operate on this appropriate space. In this work we investigated the implications of using a distance based on an affine invariant Riemannian metric in the context of structure-function mapping. Specifically, we revisit previously proposed structure-function mappings based on eigendecomposition and test them on 100 healthy subjects from the Human Connectome Project using this adapted notion of distance. First, we show that using this Riemannian distance significantly alters the notion of similarity between subjects from a functional point of view. We also show that using this distance improves the correlation between the structural and functional similarity of different subjects. Finally, by using a distance appropriate to this manifold, we demonstrate the importance of mapping function from structure under the Riemannian manifold and show in particular that it is possible to outperform the group average and the so-called glass ceiling on the performance of mappings based on eigenmodes.

An Anisotropic 4D Filtering Approach to Recover Brain Activation From Paradigm-Free Functional MRI Data.

Context: Functional Magnetic Resonance Imaging (fMRI) is a non-invasive imaging technique that provides an indirect view into brain activity via the blood oxygen level dependent (BOLD) response. In particular, resting-state fMRI poses challenges to the recovery of brain activity without prior knowledge on the experimental paradigm, as it is the case for task fMRI. Conventional methods to infer brain activity from the fMRI signals, for example, the general linear model (GLM), require the knowledge of the experimental paradigm to define regressors and estimate the contribution of each voxel's time course to the task. To overcome this limitation, approaches to deconvolve the BOLD response and recover the underlying neural activations without a priori information on the task have been proposed. State-of-the-art techniques, and in particular the total activation (TA), formulate the deconvolution as an optimization problem with decoupled spatial and temporal regularization and an optimization strategy that alternates between the constraints. Approach: In this work, we propose a paradigm-free regularization algorithm named Anisotropic 4D-fMRI (A4D-fMRI) that is applied on the 4D fMRI image, acting simultaneously in the 3D space and 1D time dimensions. Based on the idea that large image variations should be preserved as they occur during brain activations, whereas small variations considered as noise should be removed, the A4D-fMRI applies an anisotropic regularization, thus recovering the location and the duration of brain activations. Results: Using the experimental paradigm as ground truth, the A4D-fMRI is validated on synthetic and real task-fMRI data from 51 subjects, and its performance is compared to the TA. Results show higher correlations of the recovered time courses with the ground truth compared to the TA and lower computational times. In addition, we show that the A4D-fMRI recovers activity that agrees with the GLM, without requiring or using any knowledge of the experimental paradigm.

Network alignment and similarity reveal atlas-based topological differences in structural connectomes.

The interactions between different brain regions can be modeled as a graph, called connectome, whose nodes correspond to parcels from a predefined brain atlas. The edges of the graph encode the strength of the axonal connectivity between regions of the atlas that can be estimated via diffusion magnetic resonance imaging (MRI) tractography. Herein, we aim to provide a novel perspective on the problem of choosing a suitable atlas for structural connectivity studies by assessing how robustly an atlas captures the network topology across different subjects in a homogeneous cohort. We measure this robustness by assessing the alignability of the connectomes, namely the possibility to retrieve graph matchings that provide highly similar graphs. We introduce two novel concepts. First, the graph Jaccard index (GJI), a graph similarity measure based on the well-established Jaccard index between sets; the GJI exhibits natural mathematical properties that are not satisfied by previous approaches. Second, we devise WL-align, a new technique for aligning connectomes obtained by adapting the Weisfeiler-Leman (WL) graph-isomorphism test. We validated the GJI and WL-align on data from the Human Connectome Project database, inferring a strategy for choosing a suitable parcellation for structural connectivity studies. Code and data are publicly available.

On the generalizability of diffusion MRI signal representations across acquisition parameters, sequences and tissue types: Chronicles of the MEMENTO challenge.

Diffusion MRI (dMRI) has become an invaluable tool to assess the microstructural organization of brain tissue. Depending on the specific acquisition settings, the dMRI signal encodes specific properties of the underlying diffusion process. In the last two decades, several signal representations have been proposed to fit the dMRI signal and decode such properties. Most methods, however, are tested and developed on a limited amount of data, and their applicability to other acquisition schemes remains unknown. With this work, we aimed to shed light on the generalizability of existing dMRI signal representations to different diffusion encoding parameters and brain tissue types. To this end, we organized a community challenge - named MEMENTO, making available the same datasets for fair comparisons across algorithms and techniques. We considered two state-of-the-art diffusion datasets, including single-diffusion-encoding (SDE) spin-echo data from a human brain with over 3820 unique diffusion weightings (the MASSIVE dataset), and double (oscillating) diffusion encoding data (DDE/DODE) of a mouse brain including over 2520 unique data points. A subset of the data sampled in 5 different voxels was openly distributed, and the challenge participants were asked to predict the remaining part of the data. After one year, eight participant teams submitted a total of 80 signal fits. For each submission, we evaluated the mean squared error, the variance of the prediction error and the Bayesian information criteria. The received submissions predicted either multi-shell SDE data (37%) or DODE data (22%), followed by cartesian SDE data (19%) and DDE (18%). Most submissions predicted the signals measured with SDE remarkably well, with the exception of low and very strong diffusion weightings. The prediction of DDE and DODE data seemed more challenging, likely because none of the submissions explicitly accounted for diffusion time and frequency. Next to the choice of the model, decisions on fit procedure and hyperparameters play a major role in the prediction performance, highlighting the importance of optimizing and reporting such choices. This work is a community effort to highlight strength and limitations of the field at representing dMRI acquired with trending encoding schemes, gaining insights into how different models generalize to different tissue types and fiber configurations over a large range of diffusion encodings.

Towards linking diffusion MRI based macro- and microstructure measures with cortico-cortical transmission in brain tumor patients.

We aimed to link macro- and microstructure measures of brain white matter obtained from diffusion MRI with effective connectivity measures based on a propagation of cortico-cortical evoked potentials induced with intrasurgical direct electrical stimulation. For this, we compared streamline lengths and log-transformed ratios of streamlines computed from presurgical diffusion-weighted images, and the delays and amplitudes of N1 peaks recorded intrasurgically with electrocorticography electrodes in a pilot study of 9 brain tumor patients. Our results showed positive correlation between these two modalities in the vicinity of the stimulation sites (Pearson coefficient 0.54±0.13 for N1 delays, and 0.47±0.23 for N1 amplitudes), which could correspond to the neural propagation via U-fibers. In addition, we reached high sensitivities (0.78±0.07) and very high specificities (0.93±0.03) in a binary variant of our comparison. Finally, we used the structural connectivity measures to predict the effective connectivity using a multiple linear regression model, and showed a significant role of brain microstructure-related indices in this relation.

Diffusion MRI tractography filtering techniques change the topology of structural connectomes.

Objective.The use of non-invasive techniques for the estimation of structural brain networks (i.e. connectomes) opened the door to large-scale investigations on the functioning and the architecture of the brain, unveiling the link between neurological disorders and topological changes of the brain network. This study aims at assessing if and how the topology of structural connectomes estimated non-invasively with diffusion MRI is affected by the employment of tractography filtering techniques in structural connectomic pipelines. Additionally, this work investigates the robustness of topological descriptors of filtered connectomes to the common practice of density-based thresholding.Approach.We investigate the changes in global efficiency, characteristic path length, modularity and clustering coefficient on filtered connectomes obtained with the spherical deconvolution informed filtering of tractograms and using the convex optimization modelling for microstructure informed tractography. The analysis is performed on both healthy subjects and patients affected by traumatic brain injury and with an assessment of the robustness of the computed graph-theoretical measures with respect to density-based thresholding of the connectome.Main results.Our results demonstrate that tractography filtering techniques change the topology of brain networks, and thus alter network metrics both in the pathological and the healthy cases. Moreover, the measures are shown to be robust to density-based thresholding.Significance.The present work highlights how the inclusion of tractography filtering techniques in connectomic pipelines requires extra caution as they systematically change the network topology both in healthy subjects and patients affected by traumatic brain injury. Finally, the practice of low-to-moderate density-based thresholding of the connectomes is confirmed to have negligible effects on the topological analysis.

A unified framework for multimodal structure-function mapping based on eigenmodes.

Characterizing the connection between brain structure and brain function is essential for understanding how behaviour emerges from the underlying anatomy. A number of studies have shown that the network structure of the white matter shapes functional connectivity. Therefore, it should be possible to predict, at least partially, functional connectivity given the structural network. Many structure-function mappings have been proposed in the literature, including several direct mappings between the structural and functional connectivity matrices. However, the current literature is fragmented and does not provide a uniform treatment of current methods based on eigendecompositions. In particular, existing methods have never been compared to each other and their relationship explicitly derived in the context of brain structure-function mapping. In this work, we propose a unified computational framework that generalizes recently proposed structure-function mappings based on eigenmodes. Using this unified framework, we highlight the link between existing models and show how they can be obtained by specific choices of the parameters of our framework. By applying our framework to 50 subjects of the Human Connectome Project, we reproduce 6 recently published results, devise two new models and provide a direct comparison between all mappings. Finally, we show that a glass ceiling on the performance of mappings based on eigenmodes seems to be reached and conclude with possible approaches to break this performance limit.

Analytical and fast Fiber Orientation Distribution reconstruction in 3D-Polarized Light Imaging.

Three dimensional Polarized Light Imaging (3D-PLI) is an optical technique which allows mapping the spatial fiber architecture of fibrous postmortem tissues, at sub-millimeter resolutions. Here, we propose an analytical and fast approach to compute the fiber orientation distribution (FOD) from high-resolution vector data provided by 3D-PLI. The FOD is modeled as a sum of K orientations/Diracs on the unit sphere, described on a spherical harmonics basis and analytically computed using the spherical Fourier transform. Experiments are performed on rich synthetic data which simulate the geometry of the neuronal fibers and on human brain data. Results indicate the analytical FOD is computationally efficient and very fast, and has high angular precision and angular resolution. Furthermore, investigations on the right occipital lobe illustrate that our strategy of FOD computation enables the bridging of spatial scales from microscopic 3D-PLI information to macro- or mesoscopic dimensions of diffusion Magnetic Resonance Imaging (MRI), while being a means to evaluate prospective resolution limits for diffusion MRI to reconstruct region-specific white matter tracts. These results demonstrate the interest and great potential of our analytical approach.

Non-invasive inference of information flow using diffusion MRI, functional MRI, and MEG.

OBJECTIVE: To infer information flow in the white matter of the brain and recover cortical activity using functional MRI, diffusion MRI, and MEG without a manual selection of the white matter connections of interest. APPROACH: A Bayesian network which encodes the priors knowledge of possible brain states is built from imaging data. Diffusion MRI is used to enumerate all possible connections between cortical regions. Functional MRI is used to prune connections without manual intervention and increase the likelihood of specific regions being active. MEG data is used as evidence into this network to obtain a posterior distribution on cortical regions and connections. MAIN RESULTS: We show that our proposed method is able to identify connections associated with the a sensory-motor task. This allows us to build the Bayesian network with no manual selection of connections of interest. Using sensory-motor MEG evoked response as evidence into this network, our method identified areas known to be involved in a visuomotor task. In addition, information flow along white matter fiber bundles connecting those regions was also recovered. SIGNIFICANCE: Current methods to estimate white matter information flow are extremely invasive, therefore limiting our understanding of the interaction between cortical regions. The proposed method makes use of functional MRI, diffusion MRI, and M/EEG to infer communication between cortical regions, therefore opening the door to the non-invasive exploration of information flow in the white matter.

Adaptive phase correction of diffusion-weighted images.

Phase correction (PC) is a preprocessing technique that exploits the phase of images acquired in Magnetic Resonance Imaging (MRI) to obtain real-valued images containing tissue contrast with additive Gaussian noise, as opposed to magnitude images which follow a non-Gaussian distribution, e.g. Rician. PC finds its natural application to diffusion-weighted images (DWIs) due to their inherent low signal-to-noise ratio and consequent non-Gaussianity that induces a signal overestimation bias that propagates to the calculated diffusion indices. PC effectiveness depends upon the quality of the phase estimation, which is often performed via a regularization procedure. We show that a suboptimal regularization can produce alterations of the true image contrast in the real-valued phase-corrected images. We propose adaptive phase correction (APC), a method where the phase is estimated by using MRI noise information to perform a complex-valued image regularization that accounts for the local variance of the noise. We show, on synthetic and acquired data, that APC leads to phase-corrected real-valued DWIs that present a reduced number of alterations and a reduced bias. The substantial absence of parameters for which human input is required favors a straightforward integration of APC in MRI processing pipelines.

A computational Framework for generating rotation invariant features and its application in diffusion MRI.

In this work, we present a novel computational framework for analytically generating a complete set of algebraically independent Rotation Invariant Features (RIF) given the Laplace-series expansion of a spherical function. Our computational framework provides a closed-form solution for these new invariants, which are the natural expansion of the well known spherical mean, power-spectrum and bispectrum invariants. We highlight the maximal number of algebraically independent invariants which can be obtained from a truncated Spherical Harmonic (SH) representation of a spherical function and show that most of these new invariants can be linked to statistical and geometrical measures of spherical functions, such as the mean, the variance and the volume of the spherical signal. Moreover, we demonstrate their application to dMRI signal modeling including the Apparent Diffusion Coefficient (ADC), the diffusion signal and the fiber Orientation Distribution Function (fODF). In addition, using both synthetic and real data, we test the ability of our invariants to estimate brain tissue microstructure in healthy subjects and show that our framework provides more flexibility and open up new opportunities for innovative development in the domain of microstructure recovery from diffusion MRI.

The Dmipy Toolbox: Diffusion MRI Multi-Compartment Modeling and Microstructure Recovery Made Easy.

Non-invasive estimation of brain microstructure features using diffusion MRI (dMRI)-known as Microstructure Imaging-has become an increasingly diverse and complicated field over the last decades. Multi-compartment (MC)-models, representing the measured diffusion signal as a linear combination of signal models of distinct tissue types, have been developed in many forms to estimate these features. However, a generalized implementation of MC-modeling as a whole, providing deeper insights in its capabilities, remains missing. To address this fact, we present Diffusion Microstructure Imaging in Python (Dmipy), an open-source toolbox implementing PGSE-based MC-modeling in its most general form. Dmipy allows on-the-fly implementation, signal modeling, and optimization of any user-defined MC-model, for any PGSE acquisition scheme. Dmipy follows a "building block"-based philosophy to Microstructure Imaging, meaning MC-models are modularly constructed to include any number and type of tissue models, allowing simultaneous representation of a tissue's diffusivity, orientation, volume fractions, axon orientation dispersion, and axon diameter distribution. In particular, Dmipy is geared toward facilitating reproducible, reliable MC-modeling pipelines, often allowing the whole process from model construction to parameter map recovery in fewer than 10 lines of code. To demonstrate Dmipy's ease of use and potential, we implement a wide range of well-known MC-models, including IVIM, AxCaliber, NODDI(x), Bingham-NODDI, the spherical mean-based SMT and MC-MDI, and spherical convolution-based single- and multi-tissue CSD. By allowing parameter cascading between MC-models, Dmipy also facilitates implementation of advanced approaches like CSD with voxel-varying kernels and single-shell 3-tissue CSD. By providing a well-tested, user-friendly toolbox that simplifies the interaction with the otherwise complicated field of dMRI-based Microstructure Imaging, Dmipy contributes to more reproducible, high-quality research.

White matter information flow mapping from diffusion MRI and EEG.

The human brain can be described as a network of specialized and spatially distributed regions. The activity of individual regions can be estimated using electroencephalography and the structure of the network can be measured using diffusion magnetic resonance imaging. However, the communication between the different cortical regions occurring through the white matter, coined information flow, cannot be observed by either modalities independently. Here, we present a new method to infer information flow in the white matter of the brain from joint diffusion MRI and EEG measurements. This is made possible by the millisecond resolution of EEG which makes the transfer of information from one region to another observable. A subject specific Bayesian network is built which captures the possible interactions between brain regions at different times. This network encodes the connections between brain regions detected using diffusion MRI tractography derived white matter bundles and their associated delays. By injecting the EEG measurements as evidence into this model, we are able to estimate the directed dynamical functional connectivity whose delays are supported by the diffusion MRI derived structural connectivity. We present our results in the form of information flow diagrams that trace transient communication between cortical regions over a functional data window. The performance of our algorithm under different noise levels is assessed using receiver operating characteristic curves on simulated data. In addition, using the well-characterized visual motor network as grounds to test our model, we present the information flow obtained during a reaching task following left or right visual stimuli. These promising results present the transfer of information from the eyes to the primary motor cortex. The information flow obtained using our technique can also be projected back to the anatomy and animated to produce videos of the information path through the white matter, opening a new window into multi-modal dynamic brain connectivity.

Fiber orientation distribution function from non-negative sparse recovery with quantitative analysis of local fiber orientations and tractography using DW-MRI datasets.

Diffusion weighted MRI (DW-MRI) is the unique non-invasive imaging modality capable of estimating in vivo the structure of the white matter. In this paper, we propose, evaluate and validate a new DW-MRI method to model and recover high quality tractogram even with multiple fiber populations in a voxel and from a limited number of acquisitions. Our method relies on the estimation of the Fiber Orientation Distribution (FOD) function, parameterized as a non-negative sum of rank-1 tensors and the use of a non-negative sparse recovery scheme to efficiently recover the tensors, and their number. Each fiber population of a voxel is characterized by the orientation and the weight of a rank-1 tensor. Using both deterministic and probabilistic tractography algorithms, we show that our method is able to accurately reconstruct narrow crossing fibers and obtain a high quality connectivity reconstruction even from a limited number of acquisitions. To this end, a validation scheme based on the connectivity recovered from tractography is developed to quantitatively evaluate and analyze the performance of our method. The tractometer tool is used to quantify the tractography obtained from a simulated DW-MRI dataset including a high angular resolution dataset of 60 gradient directions and a dataset of 30 gradient directions, each of them corrupted with Rician noise of SNR 10 and 20. The performance of our FOD model and its impact on the tractography results are also demonstrated and illustrated on in vivo DW-MRI datasets with high and low angular resolutions.

Reducing the number of samples in spatiotemporal dMRI acquisition design.

PURPOSE: Acquisition time is a major limitation in recovering brain white matter microstructure with diffusion magnetic resonance imaging. The aim of this paper is to bridge the gap between growing demands on spatiotemporal resolution of diffusion signal and the real-world time limitations. The authors introduce an acquisition scheme that reduces the number of samples under adjustable quality loss. METHODS: Finding a sampling scheme that maximizes signal quality and satisfies given time constraints is NP-hard. Therefore, a heuristic method based on genetic algorithm is proposed in order to find suboptimal solutions in acceptable time. The analyzed diffusion signal representation is defined in the qτ space, so that it captures both spacial and temporal phenomena. RESULTS: The experiments on synthetic data and in vivo diffusion images of the C57Bl6 wild-type mouse corpus callosum reveal superiority of the proposed approach over random sampling and even distribution in the qτ space. CONCLUSIONS: The use of genetic algorithm allows to find acquisition parameters that guarantee high signal reconstruction accuracy under given time constraints. In practice, the proposed approach helps to accelerate the acquisition for the use of qτ-dMRI signal representation.

A Second Order Multi-Stencil Fast Marching Method with a Non-Constant Local Cost Model.

The Fast Marching method is widely employed in several fields of image processing. Some years ago a Multi-Stencil version (MSFM) was introduced to improve its accuracy by solving the equation for a set of stencils and choosing the best solution at each considered node. The following work proposes a modified numerical scheme for MSFM to take into account the variation of the local cost, which has proven to be second order. The influence of the stencil set choice on the algorithm outcome with respect to stencil orthogonality and axis swapping is also explored, where stencils are taken from neighborhoods of varying radius. The experimental results show that the proposed schemes improve the accuracy of their original counterparts, and that the use of permutation-invariant stencil sets provides robustness against shifted vector coordinates in the stencil set.

Groupwise structural parcellation of the whole cortex: A logistic random effects model based approach.

Current theories hold that brain function is highly related to long-range physical connections through axonal bundles, namely extrinsic connectivity. However, obtaining a groupwise cortical parcellation based on extrinsic connectivity remains challenging. Current parcellation methods are computationally expensive; need tuning of several parameters or rely on ad-hoc constraints. Furthermore, none of these methods present a model for the cortical extrinsic connectivity of the cortex. To tackle these problems, we propose a parsimonious model for the extrinsic connectivity and an efficient parceling technique based on clustering of tractograms. Our technique allows the creation of single subject and groupwise parcellations of the whole cortex. The parcellations obtained with our technique are in agreement with structural and functional parcellations in the literature. In particular, the motor and sensory cortex are subdivided in agreement with the human homunculus of Penfield. We illustrate this by comparing our resulting parcels with the motor strip mapping included in the Human Connectome Project data.