FFClust: Fast fiber clustering for large tractography datasets for a detailed study of brain connectivity

Fiber Microstructure Quantile (FMQ) Regression: A Novel Statistical Approach for Analyzing White Matter Bundles from Periphery to Core

The structural connections of the brain's white matter are critical for brain function. Diffusion MRI tractography enables the in-vivo reconstruction of white matter fiber bundles and the study of their relationship to covariates of interest, such as neurobehavioral or clinical factors. In this work, we introduce Fiber Microstructure Quantile (FMQ) Regression, a new statistical approach for studying the association between white matter fiber bundles and scalar factors (e.g., cognitive scores). Our approach analyzes tissue microstructure measures based on quantile-specific bundle regions. These regions are defined according to the quantiles of fractional anisotropy (FA) from the periphery to the core of a population fiber bundle, which pools all individuals' bundles. To investigate how fiber bundle tissue microstructure relates to covariates of interest, we employ the statistical technique of quantile regression. Unlike ordinary regression, which only models a conditional mean, quantile regression models the conditional quantiles of a response variable. This enables the proposed analysis, where a quantile regression is fitted for each quantile-specific bundle region. To demonstrate FMQ Regression, we perform an illustrative study in a large healthy young adult tractography dataset derived from the Human Connectome Project-Young Adult (HCP-YA), focusing on particular bundles expected to relate to particular aspects of cognition and motor function. Importantly, our analysis considers sex-specific effects in brain-behavior associations. In comparison with a traditional method, Automated Fiber Quantification (AFQ), which enables FA analysis in regions defined along the trajectory of a bundle, our results suggest that FMQ Regression is much more powerful for detecting brain-behavior associations. Importantly, FMQ Regression finds significant brain-behavior associations in multiple bundles, including findings unique to males or to females. In both males and females, language performance is significantly associated with FA in the left arcuate fasciculus, with stronger associations in the bundle's periphery. In males only, memory performance is significantly associated with FA in the left uncinate fasciculus, particularly in intermediate regions of the bundle. In females only, motor performance is significantly associated with FA in the left and right corticospinal tracts, with a slightly lower relationship at the bundle periphery and a slightly higher relationship toward the bundle core. No significant relationships are found between executive function and cingulum bundle FA. Our study demonstrates that FMQ Regression is a powerful statistical approach that can provide insight into associations from bundle periphery to bundle core. Our results also identify several brain-behavior relationships unique to males or to females, highlighting the importance of considering sex differences in future research.

Assessment of U-Net in the segmentation of short tracts: Transferring to clinical MRI routine

Accurately studying structural connectivity requires precise tract segmentation strategies. The U-Net network has been widely recognized for its exceptional capacity in image segmentation tasks and provides remarkable results in large tract segmentation when high-quality diffusion-weighted imaging (DWI) data are used. However, short tracts, which are associated with various neurological diseases, pose specific challenges, particularly when high-quality DWI data acquisition within clinical settings is concerned. Here, we aimed to evaluate the U-Net network ability to segment short tracts by using DWI data acquired in different experimental conditions. To this end, we conducted three types of training experiments involving 350 healthy subjects and 11 white matter tracts, including the anterior, posterior, and hippocampal commissure, fornix, and uncinated fasciculus. In the first experiment, the model was exclusively trained with high-quality data of the Human Connectome Project (HCP) dataset. The second experiment focused on images of healthy subjects acquired from a local hospital dataset, representing a typical clinical routine acquisition. In the third experiment, a hybrid training approach was employed, combining data of the HCP and local hospital datasets. Then, the best model was also tested in unseen DWIs of 10 epilepsy patients of the local hospital and 10 healthy subjects acquired on a scanner from another company. The outcomes of the third experiment demonstrated a notable enhancement in performance when contrasted with the preceding trials. Specifically, the short tracts within the local hospital dataset achieved Dice scores ranging between 0.60 and 0.65. Similar intervals were obtained with HCP data in the first experiment, and a substantial improvement compared to the scores between 0.37 and 0.50 obtained with the local hospital dataset at the same experiment. This improvement persisted when the method was applied to diverse scenarios, including different scanner acquisitions and epilepsy patients. These results indicate that combining datasets from different sources, coupled with resolution standardization strengthens the neural network ability to generalize predictions across a spectrum of datasets. Nevertheless, short tract segmentation performance is intricately linked to the training composition, to validation, and to testing data. Moreover, curved tracts have intricate structural nature, which adds complexities to their segmenting. Although the network training approach tested herein has provided promising results, caution must be taken when extrapolating its application to datasets acquired under distinct experimental conditions, even in the case of higher-quality data or analysis of long or short tracts.

Deep multimodal saliency parcellation of cerebellar pathways: Linking microstructure and individual function through explainable multitask learning

Parcellation of human cerebellar pathways is essential for advancing our understanding of the human brain. Existing diffusion magnetic resonance imaging tractography parcellation methods have been successful in defining major cerebellar fibre tracts, while relying solely on fibre tract structure. However, each fibre tract may relay information related to multiple cognitive and motor functions of the cerebellum. Hence, it may be beneficial for parcellation to consider the potential importance of the fibre tracts for individual motor and cognitive functional performance measures. In this work, we propose a multimodal data-driven method for cerebellar pathway parcellation, which incorporates both measures of microstructure and connectivity, and measures of individual functional performance. Our method involves first training a multitask deep network to predict various cognitive and motor measures from a set of fibre tract structural features. The importance of each structural feature for predicting each functional measure is then computed, resulting in a set of structure-function saliency values that are clustered to parcellate cerebellar pathways. We refer to our method as Deep Multimodal Saliency Parcellation (DeepMSP), as it computes the saliency of structural measures for predicting cognitive and motor functional performance, with these saliencies being applied to the task of parcellation. Applying DeepMSP to a large-scale dataset from the Human Connectome Project Young Adult study (n = 1065), we found that it was feasible to identify multiple cerebellar pathway parcels with unique structure-function saliency patterns that were stable across training folds. We thoroughly experimented with all stages of the DeepMSP pipeline, including network selection, structure-function saliency representation, clustering algorithm, and cluster count. We found that a 1D convolutional neural network architecture and a transformer network architecture both performed comparably for the multitask prediction of endurance, strength, reading decoding, and vocabulary comprehension, with both architectures outperforming a fully connected network architecture. Quantitative experiments demonstrated that a proposed low-dimensional saliency representation with an explicit measure of motor versus cognitive category bias achieved the best parcellation results, while a parcel count of four was most successful according to standard cluster quality metrics. Our results suggested that motor and cognitive saliencies are distributed across the cerebellar white matter pathways. Inspection of the final k = 4 parcellation revealed that the highest-saliency parcel was most salient for the prediction of both motor and cognitive performance scores and included parts of the middle and superior cerebellar peduncles. Our proposed saliency-based parcellation framework, DeepMSP, enables multimodal, data-driven tractography parcellation. Through utilising both structural features and functional performance measures, this parcellation strategy may have the potential to enhance the study of structure-function relationships of the cerebellar pathways.

Individual differences of white matter characteristic along the anterior insula-based fiber tract circuit for pain empathy in healthy women and women with primary dysmenorrhea

Pain empathy, defined as the ability of one person to understand another person's pain, shows large individual variations. The anterior insula is the core region of the pain empathy network. However, the relationship between white matter (WM) properties of the fiber tracts connecting the anterior insula with other cortical regions and an individual's ability to modulate pain empathy remains largely unclear. In this study, we outline an automatic seed-based fiber streamline (sFS) analysis method and multivariate pattern analysis (MVPA) to predict the levels of pain empathy in healthy women and women with primary dysmenorrhoea (PDM). Using the sFS method, the anterior insula-based fiber tract network was divided into five fiber cluster groups. In healthy women, interindividual differences in pain empathy were predicted only by the WM properties of the five fiber cluster groups, suggesting that interindividual differences in pain empathy may rely on the connectivity of the anterior insula-based fiber tract network. In women with PDM, pain empathy could be predicted by a single cluster group. The mean WM properties along the anterior insular-rostroventral area of the inferior parietal lobule further mediated the effect of pain on empathy in patients with PDM. Our results suggest that chronic periodic pain may lead to maladaptive plastic changes, which could further impair empathy by making women with PDM feel more pain when they see other people experiencing pain. Our study also addresses an important gap in the analysis of the microstructural characteristics of seed-based fiber tract network.

PhyberSIM: a tool for the generation of ground truth to evaluate brain fiber clustering algorithms

Diffusion Magnetic Resonance Imaging tractography is a non-invasive technique that produces a collection of streamlines representing the main white matter bundle trajectories. Methods, such as fiber clustering algorithms, are important in computational neuroscience and have been the basis of several white matter analysis methods and studies. Nevertheless, these clustering methods face the challenge of the absence of ground truth of white matter fibers, making their evaluation difficult. As an alternative solution, we present an innovative brain fiber bundle simulator that uses spline curves for fiber representation. The methodology uses a tubular model for the bundle simulation based on a bundle centroid and five radii along the bundle. The algorithm was tested by simulating 28 Deep White Matter atlas bundles, leading to low inter-bundle distances and high intersection percentages between the original and simulated bundles. To prove the utility of the simulator, we created three whole-brain datasets containing different numbers of fiber bundles to assess the quality performance of QuickBundles and Fast Fiber Clustering algorithms using five clustering metrics. Our results indicate that QuickBundles tends to split less and Fast Fiber Clustering tends to merge less, which is consistent with their expected behavior. The performance of both algorithms decreases when the number of bundles is increased due to higher bundle crossings. Additionally, the two algorithms exhibit robust behavior with input data permutation. To our knowledge, this is the first whole-brain fiber bundle simulator capable of assessing fiber clustering algorithms with realistic data.

Blurred streamlines: A novel representation to reduce redundancy in tractography

Tractography is a powerful tool to study brain connectivity in vivo, but it is well known to suffer from an intrinsic trade-off between sensitivity and specificity. A critical - but usually underrated - parameter to choose that can heavily impact the quality of the estimates is the number of streamlines to be reconstructed for a given data set. In fact, sensitivity can be improved by generating more and more streamlines, as all real anatomical connections are likely reconstructed, but lots of false positives are inevitably introduced, too. Consequently, so-called tractography filtering techniques have become increasingly popular to get rid of these false positives and improve specificity. However, increasing number of streamlines introduces redundancy in tractography reconstructions, which may negatively impact the performance of filtering algorithms, especially those based on linear formulations. To address this problem, we introduce a novel streamlines representation, called "blurred streamlines", which drastically reduces the redundancy among streamlines by (i) clustering similar trajectories and (ii) spatially blurring the corresponding signal contributions. We tested the effectiveness of the blurred streamlines both on synthetic and in vivo data. Our results clearly show that this new representation is as accurate as state-of-the-art methods despite using only 5% of the input streamlines, thus significantly decreasing the computational complexity of filtering algorithms as well as storage requirements of the resulting reconstructions.

A systematic review of automated methods to perform white matter tract segmentation

White matter tract segmentation is a pivotal research area that leverages diffusion-weighted magnetic resonance imaging (dMRI) for the identification and mapping of individual white matter tracts and their trajectories. This study aims to provide a comprehensive systematic literature review on automated methods for white matter tract segmentation in brain dMRI scans. Articles on PubMed, ScienceDirect [NeuroImage, NeuroImage (Clinical), Medical Image Analysis], Scopus and IEEEXplore databases and Conference proceedings of Medical Imaging Computing and Computer Assisted Intervention Society (MICCAI) and International Symposium on Biomedical Imaging (ISBI), were searched in the range from January 2013 until September 2023. This systematic search and review identified 619 articles. Adhering to the specified search criteria using the query, "white matter tract segmentation OR fiber tract identification OR fiber bundle segmentation OR tractography dissection OR white matter parcellation OR tract segmentation," 59 published studies were selected. Among these, 27% employed direct voxel-based methods, 25% applied streamline-based clustering methods, 20% used streamline-based classification methods, 14% implemented atlas-based methods, and 14% utilized hybrid approaches. The paper delves into the research gaps and challenges associated with each of these categories. Additionally, this review paper illuminates the most frequently utilized public datasets for tract segmentation along with their specific characteristics. Furthermore, it presents evaluation strategies and their key attributes. The review concludes with a detailed discussion of the challenges and future directions in this field.

A, Hernandez Larzabal H, Román C, Poupon C, Mangin JF, Hernández C, Guevara P. Phybers: a package for brain tractography analysis

We present a Python library (Phybers) for analyzing brain tractography data. Tractography datasets contain streamlines (also called fibers) composed of 3D points representing the main white matter pathways. Several algorithms have been proposed to analyze this data, including clustering, segmentation, and visualization methods. The manipulation of tractography data is not straightforward due to the geometrical complexity of the streamlines, the file format, and the size of the datasets, which may contain millions of fibers. Hence, we collected and structured state-of-the-art methods for the analysis of tractography and packed them into a Python library, to integrate and share tools for tractography analysis. Due to the high computational requirements, the most demanding modules were implemented in C/C++. Available functions include brain Bundle Segmentation (FiberSeg), Hierarchical Fiber Clustering (HClust), Fast Fiber Clustering (FFClust), normalization to a reference coordinate system, fiber sampling, calculation of intersection between sets of brain fibers, tools for cluster filtering, calculation of measures from clusters, and fiber visualization. The library tools were structured into four principal modules: Segmentation, Clustering, Utils, and Visualization (Fibervis). Phybers is freely available on a GitHub repository under the GNU public license for non-commercial use and open-source development, which provides sample data and extensive documentation. In addition, the library can be easily installed on both Windows and Ubuntu operating systems through the pip library.

Deep Learning Methods for Identification of White Matter Fiber Tracts: Review of State-of-the-Art and Future Prospective

Quantitative analysis of white matter fiber tracts from diffusion Magnetic Resonance Imaging (dMRI) data is of great significance in health and disease. For example, analysis of fiber tracts related to anatomically meaningful fiber bundles is highly demanded in pre-surgical and treatment planning, and the surgery outcome depends on accurate segmentation of the desired tracts. Currently, this process is mainly done through time-consuming manual identification performed by neuro-anatomical experts. However, there is a broad interest in automating the pipeline such that it is fast, accurate, and easy to apply in clinical settings and also eliminates the intra-reader variabilities. Following the advancements in medical image analysis using deep learning techniques, there has been a growing interest in using these techniques for the task of tract identification as well. Recent reports on this application show that deep learning-based tract identification approaches outperform existing state-of-the-art methods. This paper presents a review of current tract identification approaches based on deep neural networks. First, we review the recent deep learning methods for tract identification. Next, we compare them with respect to their performance, training process, and network properties. Finally, we end with a critical discussion of open challenges and possible directions for future works.

Deep fiber clustering: Anatomically informed fiber clustering with self-supervised deep learning for fast and effective tractography parcellation

White matter fiber clustering is an important strategy for white matter parcellation, which enables quantitative analysis of brain connections in health and disease. In combination with expert neuroanatomical labeling, data-driven white matter fiber clustering is a powerful tool for creating atlases that can model white matter anatomy across individuals. While widely used fiber clustering approaches have shown good performance using classical unsupervised machine learning techniques, recent advances in deep learning reveal a promising direction toward fast and effective fiber clustering. In this work, we propose a novel deep learning framework for white matter fiber clustering, Deep Fiber Clustering (DFC), which solves the unsupervised clustering problem as a self-supervised learning task with a domain-specific pretext task to predict pairwise fiber distances. This process learns a high-dimensional embedding feature representation for each fiber, regardless of the order of fiber points reconstructed during tractography. We design a novel network architecture that represents input fibers as point clouds and allows the incorporation of additional sources of input information from gray matter parcellation. Thus, DFC makes use of combined information about white matter fiber geometry and gray matter anatomy to improve the anatomical coherence of fiber clusters. In addition, DFC conducts outlier removal naturally by rejecting fibers with low cluster assignment probability. We evaluate DFC on three independently acquired cohorts, including data from 220 individuals across genders, ages (young and elderly adults), and different health conditions (healthy control and multiple neuropsychiatric disorders). We compare DFC to several state-of-the-art white matter fiber clustering algorithms. Experimental results demonstrate superior performance of DFC in terms of cluster compactness, generalization ability, anatomical coherence, and computational efficiency.

Accurate corresponding fiber tract segmentation via FiberGeoMap learner with application to autism

Fiber tract segmentation is a prerequisite for tract-based statistical analysis. Brain fiber streamlines obtained by diffusion magnetic resonance imaging and tractography technology are usually difficult to be leveraged directly, thus need to be segmented into fiber tracts. Previous research mainly consists of two steps: defining and computing the similarity features of fiber streamlines, then adopting machine learning algorithms for fiber clustering or classification. Defining the similarity feature is the basic premise and determines its potential reliability and application. In this study, we adopt geometric features for fiber tract segmentation and develop a novel descriptor (FiberGeoMap) for the corresponding representation, which can effectively depict fiber streamlines' shapes and positions. FiberGeoMap can differentiate fiber tracts within the same subject, meanwhile preserving the shape and position consistency across subjects, thus can identify common fiber tracts across brains. We also proposed a Transformer-based encoder network called FiberGeoMap Learner, to perform segmentation based on the geometric features. Experimental results showed that the proposed method can differentiate the 103 various fiber tracts, which outperformed the existing methods in both the number of categories and segmentation accuracy. Furthermore, the proposed method identified some fiber tracts that were statistically different on fractional anisotropy (FA), mean diffusion (MD), and fiber number ration in autism.

Spatially regularized low-rank tensor approximation for accurate and fast tractography

Tractography based on diffusion Magnetic Resonance Imaging (dMRI) is the prevalent approach to the in vivo delineation of white matter tracts in the human brain. Many tractography methods rely on models of multiple fiber compartments, but the local dMRI information is not always sufficient to reliably estimate the directions of secondary fibers. Therefore, we introduce two novel approaches that use spatial regularization to make multi-fiber tractography more stable. Both represent the fiber Orientation Distribution Function (fODF) as a symmetric fourth-order tensor, and recover multiple fiber orientations via low-rank approximation. Our first approach computes a joint approximation over suitably weighted local neighborhoods with an efficient alternating optimization. The second approach integrates the low-rank approximation into a current state-of-the-art tractography algorithm based on the unscented Kalman filter (UKF). These methods were applied in three different scenarios. First, we demonstrate that they improve tractography even in high-quality data from the Human Connectome Project, and that they maintain useful results with a small fraction of the measurements. Second, on the 2015 ISMRM tractography challenge, they increase overlap, while reducing overreach, compared to low-rank approximation without joint optimization or the traditional UKF, respectively. Finally, our methods permit a more comprehensive reconstruction of tracts surrounding a tumor in a clinical dataset. Overall, both approaches improve reconstruction quality. At the same time, our modified UKF significantly reduces the computational effort compared to its traditional counterpart, and to our joint approximation. However, when used with ROI-based seeding, joint approximation more fully recovers fiber spread.

Disentangling the variability of the superficial white matter organization using regional-tractogram-based population stratification

Each variation of the cortical folding pattern implies a particular rearrangement of the geometry of the fibers of the underlying white matter. While this rearrangement only impacts the ends of the long pathways, it may affect most of the trajectory of the short bundles. Therefore, mapping the short fibers of the human brain using diffusion-based tractography requires a dedicated strategy to overcome the variability of the folding patterns. In this paper, we propose a fiber-based stratification strategy splitting the population into homogeneous groups for disentangling the superficial white matter bundle organization. This strategy introduces a new refined fiber distance which includes angular considerations for inferring fine-grained atlases of the short bundles surrounding a specific sulcus and a subtractogram distance that quantifies the similitude between fiber sets of two different subjects. The stratification splits the population into groups with similar regional fiber organization using manifold learning. We first successfully test the hypothesis that the main source of variability of the regional fiber organization is the variability of the regional folding pattern. Then, in each group, we proceed with the automatic identification of the most stable bundles, at a higher granularity level than what can be achieved with the non-stratified whole population, enabling the disentanglement of the very variable configuration of the short fibers. Finally, the method searches for bundle correspondence across groups to build a population level atlas. As a proof of concept, the atlas refinement achieved by this strategy is illustrated for the fibers that surround the central sulcus and the superior temporal sulcus using the HCP dataset.

Fast Streamline Search: An Exact Technique for Diffusion MRI Tractography

In this work, a hierarchical search algorithm is proposed to efficiently compute the distance between similar tractography streamlines. This hierarchical framework offers an upper bound and a lower bound for the point-wise distance between two streamlines, which guarantees the validity of a proximity search. The proposed streamline representation enables the use of space-partitioning search trees to increase the tractography clustering speed without reducing its accuracy. The resulting approach enables a fast reconstruction a sparse distance matrix between two sets of streamlines, for all similar streamlines within a given radius. Alongside a white matter atlas, this fast streamline search can be used for accurate and reproducible tractogram clustering.

White matter variability, cognition, and disorders: a systematic review

Inter-individual differences can inform treatment procedures and—if accounted for—have the potential to significantly improve patient outcomes. However, when studying brain anatomy, these inter-individual variations are commonly unaccounted for, despite reports of differences in gross anatomical features, cross-sectional, and connectional anatomy. Brain connections are essential to facilitate functional organization and, when severed, cause impairments or complete loss of function. Hence, the study of cerebral white matter may be an ideal compromise to capture inter-individual variability in structure and function. We reviewed the wealth of studies that associate cognitive functions and clinical symptoms with individual tracts using diffusion tractography. Our systematic review indicates that tractography has proven to be a sensitive method in neurology, psychiatry, and healthy populations to identify variability and its functional correlates. However, the literature may be biased, as the most commonly studied tracts are not necessarily those with the highest sensitivity to cognitive functions and pathologies. Additionally, the hemisphere of the studied tract is often unreported, thus neglecting functional laterality and asymmetries. Finally, we demonstrate that tracts, as we define them, are not correlated with one, but multiple cognitive domains or pathologies. While our systematic review identified some methodological caveats, it also suggests that tract–function correlations might still be a promising tool in identifying biomarkers for precision medicine. They can characterize variations in brain anatomy, differences in functional organization, and predicts resilience and recovery in patients.

Enhanced Automatic Segmentation for Superficial White Matter Fiber Bundles for Probabilistic Tractography Datasets

This paper presents an enhanced algorithm for automatic segmentation of superficial white matter (SWM) bundles from probabilistic dMRI tractography datasets, based on a multi-subject bundle atlas. Previous segmentation methods use the maximum Euclidean distance between corresponding points of the subject fibers and the atlas centroids. However, this scheme might include noisy fibers. Here, we propose a three step approach to discard noisy fibers improving the identification of fibers. The first step applies a fiber clustering and the segmentation is performed between the centroids of the clusters and the atlas centroids. This step removes outliers and enables a better identification of fibers with similar shapes. The second step applies a fiber filter based on two different fiber similarities. One is the Symmetrized Segment-Path Distance (SSPD) over 2D ISOMAP and the other is an adapted version of SSPD for 3D space. The last step eliminates noisy fibers by removing those that connect regions that are far from the main atlas bundle connections. We perform an experimental evaluation using ten subjects of the Human Connectome (HCP) database. The evaluation only considers the bundles connecting precentral and postcentral gyri, with a total of seven bundles per hemisphere. For comparison, the bundles of the ten subjects were manually segmented. Bundles segmented with our method were evaluated in terms of similarity to manually segmented bundles and the final number of fibers. The results show that our approach obtains bundles with a higher similarity score than the state-of-the-art method and maintains a similar number of fibers.Clinical relevance-Many brain pathologies or disorders can occur in specific regions of the SWM automatic segmentation of reliable SWM bundles would help applications to clinical research.

Group-Wise Cortical Surface Parcellation Based on Inter-Subject Fiber Clustering. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2021;2021:2655-2659. [PMID: 34891798 DOI: 10.1109/embc46164.2021.9631099] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

We present an automatic algorithm for the group-wise parcellation of the cortical surface. The method is based on the structural connectivity obtained from representative brain fiber clusters, calculated via an inter-subject clustering scheme. Preliminary regions were defined from cluster-cortical mesh intersection points. The final parcellation was obtained using parcel probability maps to model and integrate the connectivity information of all subjects, and graphs to represent the overlap between parcels. Two inter-subject clustering schemes were tested, generating a total of 171 and 109 parcels, respectively. The resulting parcels were quantitatively compared with three state-of-the-art atlases. The best parcellation returned 69 parcels with a Dice similarity coefficient greater than 0.5. To the best of our knowledge, this is the first diffusion-based cortex parcellation method based on whole-brain inter-subject fiber clustering.

Fiber Clustering Acceleration With a Modified Kmeans++ Algorithm Using Data Parallelism

Fiber clustering methods are typically used in brain research to study the organization of white matter bundles from large diffusion MRI tractography datasets. These methods enable exploratory bundle inspection using visualization and other methods that require identifying brain white matter structures in individuals or a population. Some applications, such as real-time visualization and inter-subject clustering, need fast and high-quality intra-subject clustering algorithms. This work proposes a parallel algorithm using a General Purpose Graphics Processing Unit (GPGPU) for fiber clustering based on the FFClust algorithm. The proposed GPGPU implementation exploits data parallelism using both multicore and GPU fine-grained parallelism present in commodity architectures, including current laptops and desktop computers. Our approach implements all FFClust steps in parallel, improving execution times in all of them. In addition, our parallel approach includes a parallel Kmeans++ algorithm implementation and defines a new variant of Kmeans++ to reduce the impact of choosing outliers as initial centroids. The results show that our approach provides clustering quality results very similar to FFClust, and it requires an execution time of 3.5 s for processing about a million fibers, achieving a speedup of 11.5 times compared to FFClust.

ABrainVis: an android brain image visualization tool

Background

The visualization and analysis of brain data such as white matter diffusion tractography and magnetic resonance imaging (MRI) volumes is commonly used by neuro-specialist and researchers to help the understanding of brain structure, functionality and connectivity. As mobile devices are widely used among users and their technology shows a continuous improvement in performance, different types of applications have been designed to help users in different work areas.

Results

We present, ABrainVis, an Android mobile tool that allows users to visualize different types of brain images, such as white matter diffusion tractographies, represented as fibers in 3D, segmented fiber bundles, MRI 3D images as rendered volumes and slices, and meshes. The tool enables users to choose and combine different types of brain imaging data to provide visual anatomical context for specific visualization needs. ABrainVis provides high performance over a wide range of Android devices, including tablets and cell phones using medium and large tractography datasets. Interesting visualizations including brain tumors and arteries, along with fiber, are given as examples of case studies using ABrainVis.

Conclusions

The functionality, flexibility and performance of ABrainVis tool introduce an improvement in user experience enabling neurophysicians and neuroscientists fast visualization of large tractography datasets, as well as the ability to incorporate other brain imaging data such as MRI volumes and meshes, adding anatomical contextual information.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12938-021-00909-0.