Fiber Orientation Distribution Function Research Articles

Reconstruction of intra- and extra-neurite conductivity tensors via conductivity at Larmor frequency and DWI data patterns

The developed magnetic resonance electrical properties tomography (MREPT) techniques visualize the internal conductivity distribution at Larmor frequency by measuring the B1 transceive phase data. In biological tissues, electrical conductivity is influenced by ion concentrations and mobility. To visualize the anisotropic conductivity tensor of biological tissues, we use the Einstein–Smoluchowski equation, which links the diffusion coefficient to particle mobility. By assuming a correlation between ion mobility and water diffusivity, we aim to decompose the internal isotropic conductivity at Larmor frequency into anisotropic conductivity tensors within the intra- and extra-neurite compartments. The multi-compartment spherical mean technique (MC-SMT), utilizing both high and low b-value diffusion-weighted imaging (DWI) data, characterizes the diffusion of water molecules within and across the intra- and extra-neurite compartments by analyzing the microstructural intricacies and the foundational architectural arrangement of the brain’s tissues. By analyzing the relationships between the measured DWI data, the microscopic diffusion signal, and the fiber orientation distribution function, we predict the DWI data for the intra- and extra-neurite compartments using spherical harmonic decomposition. Using the predicted DWI data for the intra- and extra-neurite compartments, we develop a conductivity tensor imaging method that operates specifically within the separated compartments. Human brain experiments, involving four healthy volunteers and a brain tumor patient, were performed to assess and confirm the reliability of the proposed method.

Changes in diffusion MRI and clinical motor function after physical/occupational therapies in toddler-aged children with spastic unilateral cerebral palsy.

Diffusion-weighted magnetic resonance imaging (DMRI) is a potential tool to assess changes in brain connectivity and microstructure resulting from physical and occupational therapy in young children with cerebral palsy. This works was carried out to assess whether DMRI can detect changes after 36 weeks of physical and occupational therapy in the microstructure and connectivity of the brains of children with cerebral palsy and determine whether imaging findings correlate with changes in clinical measures of motor function. Five children underwent anatomical MRI and DMRI and evaluations of motor function skills at baseline and after 36 weeks of intensive or once-weekly physical and occupational Perception-Action Approach therapies. Diffusion tensor imaging and constrained spherical deconvolution methods were used to calculate fractional anisotropy (FA) and fiber orientation distribution functions (fODFs), respectively. The fODFs were used to generate tractograms of the cerebrospinal tract (CST). After 36 weeks of physical and occupational therapy, all children showed increases in motor function. No changes were observed in anatomical MRI before and after therapy but CST tractography did show small differences indicating possible altered microstructure and connectivity in the brain. FA values along the CSTs, however, showed no significant changes. Reliable longitudinal DMRI can be employed in toddler-aged children with CP and DMRI has the potential to monitor neuroplastic changes in white matter microstructure. However, there is a high variability between subjects and clinical improvements were not always correlated with measures of FA along the CST.

Bundle-specific tractogram distribution estimation using higher-order streamline differential equation

Streamline tractography locally traces peak directions extracted from fiber orientation distribution (FOD) functions, lacking global information about the trend of the whole fiber bundle. Therefore, it is prone to producing erroneous tracks while missing true positive connections. In this work, we propose a new bundle-specific tractography (BST) method based on a bundle-specific tractogram distribution (BTD) function, which directly reconstructs the fiber trajectory from the start region to the termination region by incorporating the global information in the fiber bundle mask. A unified framework for any higher-order streamline differential equation is presented to describe the fiber bundles with disjoint streamlines defined based on the diffusion vectorial field. At the global level, the tractography process is simplified as the estimation of BTD coefficients by minimizing the energy optimization model, and is used to characterize the relations between BTD and diffusion tensor vector under the prior guidance by introducing the tractogram bundle information to provide anatomic priors. Experiments are performed on simulated Hough, Sine, Circle data, ISMRM 2015 Tractography Challenge data, FiberCup data, and in vivo data from the Human Connectome Project (HCP) for qualitative and quantitative evaluation. Results demonstrate that our approach reconstructs complex fiber geometry more accurately. BTD reduces the error deviation and accumulation at the local level and shows better results in reconstructing long-range, twisting, and large fanning tracts.

Deep learning microstructure estimation of developing brains from diffusion MRI: A newborn and fetal study

Diffusion-weighted magnetic resonance imaging (dMRI) is widely used to assess the brain white matter. Fiber orientation distribution functions (FODs) are a common way of representing the orientation and density of white matter fibers. However, with standard FOD computation methods, accurate estimation requires a large number of measurements that usually cannot be acquired for newborns and fetuses. We propose to overcome this limitation by using a deep learning method to map as few as six diffusion-weighted measurements to the target FOD. To train the model, we use the FODs computed using multi-shell high angular resolution measurements as target. Extensive quantitative evaluations show that the new deep learning method, using significantly fewer measurements, achieves comparable or superior results than standard methods such as Constrained Spherical Deconvolution and two state-of-the-art deep learning methods. For voxels with one and two fibers, respectively, our method shows an agreement rate in terms of the number of fibers of 77.5% and 22.2%, which is 3% and 5.4% higher than other deep learning methods, and an angular error of 10° and 20°, which is 6° and 5° lower than other deep learning methods. To determine baselines for assessing the performance of our method, we compute agreement metrics using densely sampled newborn data. Moreover, we demonstrate the generalizability of the new deep learning method across scanners, acquisition protocols, and anatomy on two clinical external datasets of newborns and fetuses. We validate fetal FODs, successfully estimated for the first time with deep learning, using post-mortem histological data. Our results show the advantage of deep learning in computing the fiber orientation density for the developing brain from in-vivo dMRI measurements that are often very limited due to constrained acquisition times. Our findings also highlight the intrinsic limitations of dMRI for probing the developing brain microstructure.

Fibre orientation distribution function mapping for short fibre polymer composite components from low resolution/large volume X-ray computed tomography

Short glass fibre injection moulded composites, used in interior and exterior automotive parts, are exposed to complex stress states, for example during a crash. As the fibre scale dominates the composite’s material properties, numerical models need to account for the local fibre orientation. In recent years, mould flow simulation results have been exploited to predict the fibre orientations for finite element models, albeit with limited accuracy. Alternatively, X-ray computed tomography can be used to directly image and analyse fibre orientations. Traditionally, achieving the necessary resolution to image individual fibres restricts the imaging to small regions of the component. However, this study takes advantage of recent advancements in imaging and image analysis to overcome this limitation. As a result, it introduces, for the first time, a reliable, fast, and automated fibre orientation mapping for a full component based on image analysis at the individual fibre level; even for cases where the pixel size is significantly larger than the fibre diameter. By scanning at lower resolutions, a drastically larger volume of interest can be achieved. The resulting fibre orientation analysis and mapping algorithm, based on X-ray computed tomography, is well matched to the level of information required for automotive crash modelling with a standard element-size of a few millimetres. The entire process, encompassing image acquisition, image analysis and fibre orientation mapping, can be directly integrated into an industrial full component application in a matter of hours.

Robust fiber orientation distribution function estimation using deep constrained spherical deconvolution for diffusion-weighted magnetic resonance imaging.

Diffusion-weighted magnetic resonance imaging (DW-MRI) is a critical imaging method for capturing and modeling tissue microarchitecture at a millimeter scale. A common practice to model the measured DW-MRI signal is via fiber orientation distribution function (fODF). This function is the essential first step for the downstream tractography and connectivity analyses. With recent advantages in data sharing, large-scale multisite DW-MRI datasets are being made available for multisite studies. However, measurement variabilities (e.g., inter- and intrasite variability, hardware performance, and sequence design) are inevitable during the acquisition of DW-MRI. Most existing model-based methods [e.g., constrained spherical deconvolution (CSD)] and learning-based methods (e.g., deep learning) do not explicitly consider such variabilities in fODF modeling, which consequently leads to inferior performance on multisite and/or longitudinal diffusion studies. In this paper, we propose a data-driven deep CSD method to explicitly constrain the scan-rescan variabilities for a more reproducible and robust estimation of brain microstructure from repeated DW-MRI scans. Specifically, the proposed method introduces a three-dimensional volumetric scanner-invariant regularization scheme during the fODF estimation. We study the Human Connectome Project (HCP) young adults test-retest group as well as the MASiVar dataset (with inter- and intrasite scan/rescan data). The Baltimore Longitudinal Study of Aging dataset is employed for external validation. From the experimental results, the proposed data-driven framework outperforms the existing benchmarks in repeated fODF estimation. By introducing the contrastive loss with scan/rescan data, the proposed method achieved a higher consistency while maintaining higher angular correlation coefficients with the CSD modeling. The proposed method is assessing the downstream connectivity analysis and shows increased performance in distinguishing subjects with different biomarkers. We propose a deep CSD method to explicitly reduce the scan-rescan variabilities, so as to model a more reproducible and robust brain microstructure from repeated DW-MRI scans. The plug-and-play design of the proposed approach is potentially applicable to a wider range of data harmonization problems in neuroimaging.

Multimodal anatomical mapping of subcortical regions in marmoset monkeys using high-resolution MRI and matched histology with multiple stains

Subcortical nuclei and other deep brain structures play essential roles in regulating the central and peripheral nervous systems. However, many of these nuclei and their subregions are challenging to identify and delineate in conventional MRI due to their small size, hidden location, and often subtle contrasts compared to neighboring regions. To address these limitations, we scanned the whole brain of the marmoset monkeys in ex vivo using a clinically feasible diffusion MRI method, called the mean apparent propagator (MAP)-MRI, along with T2W and MTR (T1-like contrast) images acquired at 7 Tesla. Additionally, we registered these multimodal MRI volumes to the high-resolution images of matched whole-brain histology sections with seven different stains obtained from the same brain specimens. At high spatial resolution, the microstructural parameters and fiber orientation distribution functions derived with MAP-MRI can distinguish the subregions of many subcortical and deep brain structures, including fiber tracts of different sizes and orientations. The good correlation with multiple but distinct histological stains from the same brain serves as a thorough validation of the structures identified with MAP-MRI and other MRI parameters. Moreover, the anatomical details of deep brain structures found in the volumes of MAP-MRI parameters are not visible in conventional T1W or T2W images. The high-resolution mapping using novel MRI contrasts, combined and correlated with histology, can elucidate structures that were previously invisible radiologically. Thus, this multimodal approach offers a roadmap toward identifying salient brain areas in vivo in future neuroradiological studies. It also provides a useful anatomical standard reference for the region definition of subcortical targets and the generation of a 3D digital template atlas for the marmoset brain research (Saleem et al., 2023). Additionally, we conducted a cross-species comparison between marmoset and macaque monkeys using results from our previous studies (Saleem et al., 2021). We found that the two species had distinct patterns of iron distribution in subregions of the basal ganglia, red nucleus, and deep cerebellar nuclei, confirmed with T2W MRI and histology.

TrGANet: Transforming 3T to 7T dMRI using Trapezoidal Rule and Graph based Attention Modules.

Diffusion MRI (dMRI) is a non-invasive tool for assessing the white matter region of the brain by approximating the fiber streamlines, structural connectivity, and estimation of microstructure. This modality can yield useful information for diagnosing several mental diseases as well as for surgical planning. The higher angular resolution diffusion imaging (HARDI) technique is helpful in obtaining more robust fiber tracts by getting a good approximation of regions where fibers cross. Moreover, HARDI is more sensitive to tissue changes and can accurately represent anatomical details in the human brain at higher magnetic strengths. In other words, magnetic strengths affect the quality of the image, and hence high magnetic strength has good tissue contrast with better spatial resolution. However, a higher magnetic strength scanner (like 7T) is costly and unaffordable to most hospitals. Hence, in this work, we have proposed a novel CNN architecture for the transformation of 3T to 7T dMRI. Additionally, we have also reconstructed the multi-shell multi-tissue fiber orientation distribution function (MSMT fODF) at 7T from single-shell 3T. The proposed architecture consists of a CNN-based ODE solver utilizing the Trapezoidal rule and graph-based attention layer alongwith L1 and total variation loss. Finally, the model has been validated on the HCP data set quantitatively and qualitatively.

ID: 216730 Comparing Tractography Methods for Preoperative Deep Brain Stimulation Planning in Essential Tremor Patients

Accurate lead placement is essential to successful treatment with deep brain stimulation (DBS). The Ventralis Intermediate Nucleus (VIM) of the thalamus has been identified as a surgical target for DBS neuromodulation of motor input with successful outcomes [1-3]. To aid in this, prospective pre-operative tractography can be implemented for targeting the VIM in ET patients [4,5]. For this, FSL’s ball-and-stick model based probabilistic tracking” (Probtrackx) has been validated to be an effective tractography method (TM) [5]. However, with many other tractography widely available and embedded in commercial surgical platforms, there is benefit in evaluating the relative performance of tractography methods (TMs). Here, we assessed the relative performance of five TMs from the MRTRIX toolbox: three of these techniques were based on fiber orientation distribution functions (ODFs) represented by spherical harmonic series: iFOD1, iFOD2, iFOD-rk4, and two tensor-based models: Tensor-Deterministic, and Tensor-Probabilistic.

CNTSeg: A multimodal deep-learning-based network for cranial nerves tract segmentation.

The segmentation of cranial nerves (CNs) tracts based on diffusion magnetic resonance imaging (dMRI) provides a valuable quantitative tool for the analysis of the morphology and course of individual CNs. Tractography-based approaches can describe and analyze the anatomical area of CNs by selecting the reference streamlines in combination with ROIs-based (regions-of-interests) or clustering-based. However, due to the slender structure of CNs and the complex anatomical environment, single-modality data based on dMRI cannot provide a complete and accurate description, resulting in low accuracy or even failure of current algorithms in performing individualized CNs segmentation. In this work, we propose a novel multimodal deep-learning-based multi-class network for automated cranial nerves tract segmentation without using tractography, ROI placement or clustering, called CNTSeg. Specifically, we introduced T1w images, fractional anisotropy (FA) images, and fiber orientation distribution function (fODF) peaks into the training data set, and design the back-end fusion module which uses the complementary information of the interphase feature fusion to improve the segmentation performance. CNTSeg has achieved the segmentation of 5 pairs of CNs (i.e. optic nerve CN II, oculomotor nerve CN III, trigeminal nerve CN V, and facial-vestibulocochlear nerve CN VII/VIII). Extensive comparisons and ablation experiments show promising results and are anatomically convincing even for difficult tracts. The code will be openly available athttps://github.com/IPIS-XieLei/CNTSeg.

Deep Constrained Spherical Deconvolution for Robust Harmonization.

Diffusion weighted magnetic resonance imaging (DW-MRI) captures tissue microarchitecture at millimeter scale. With recent advantages in data sharing, large-scale multi-site DW-MRI datasets are being made available for multi-site studies. However, DW-MRI suffers from measurement variability (e.g., inter- and intra-site variability, hardware performance, and sequence design), which consequently yields inferior performance on multi-site and/or longitudinal diffusion studies. In this study, we propose a novel, deep learning-based method to harmonize DW-MRI signals for a more reproducible and robust estimation of microstructure. Our method introduces a data-driven scanner-invariant regularization scheme to model a more robust fiber orientation distribution function (FODF) estimation. We study the Human Connectome Project (HCP) young adults test-retest group as well as the MASiVar dataset (with inter- and intra-site scan/rescan data). The 8th order spherical harmonics coefficients are employed as data representation. The results show that the proposed harmonization approach maintains higher angular correlation coefficients (ACC) with the ground truth signals (0.954 versus 0.942), while achieves higher consistency of FODF signals for intra-scanner data (0.891 versus 0.826), as compared with the baseline supervised deep learning scheme. Furthermore, the proposed data-driven framework is flexible and potentially applicable to a wider range of data harmonization problems in neuroimaging.

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.

Improving the brain image resolution of generalized q-sampling MRI revealed by a three-dimensional CNN-based method.

Understanding neural connections facilitates the neuroscience and cognitive behavioral research. There are many nerve fiber intersections in the brain that need to be observed, and the size is between 30 and 50 nanometers. Improving image resolution has become an important issue for mapping the neural connections non-invasively. Generalized q-sampling imaging (GQI) was used to reveal the fiber geometry of straight and crossing. In this work, we attempted to achieve super-resolution with a deep learning method on diffusion weighted imaging (DWI). A three-dimensional super-resolution convolutional neural network (3D SRCNN) was utilized to achieve super-resolution on DWI. Then, generalized fractional anisotropy (GFA), normalized quantitative anisotropy (NQA), and the isotropic value of the orientation distribution function (ISO) mapping were reconstructed using GQI with super-resolution DWI. We also reconstructed the orientation distribution function (ODF) of brain fibers using GQI. With the proposed super-resolution method, the reconstructed DWI was closer to the target image than the interpolation method. The peak signal-to-noise ratio (PSNR) and structural similarity index measure (SSIM) were also significantly improved. The diffusion index mapping reconstructed by GQI also had higher performance. The ventricles and white matter regions were much clearer. This super-resolution method can assist in postprocessing low-resolution images. With SRCNN, high-resolution images can be effectively and accurately generated. The method can clearly reconstruct the intersection structure in the brain connectome and has the potential to accurately describe the fiber geometry on a subvoxel scale.

Undersampled single-shell to MSMT fODF reconstruction using CNN-based ODE solver

Background and Objective: Diffusion MRI (dMRI) has been considered one of the most popular non-invasive techniques for studying the human brain’s white matter (WM). dMRI is used to delineate the brain’s microstructure by approximating the WM region’s fiber tracts. The achieved fiber tracts can be utilized to assess mental diseases like Multiple sclerosis, ADHD, Seizures, Intellectual disability, and others. New techniques such as high angular resolution diffusion-weighted imaging (HARDI) have been developed, providing precise fiber directions, and overcoming the limitation of traditional DTI. Unlike Single-shell, Multi-shell HARDI provides tissue fractions for white matter, gray matter, and cerebrospinal fluid, resulting in a Multi-shell Multi-tissue fiber orientation distribution function (MSMT fODF). This MSMT fODF comes up with more precise fiber directions than a Single-shell, which helps to get correct fiber tracts. In addition, various multi-compartment diffusion models, including as CHARMED and NODDI, have been developed to describe the brain tissue microstructural information. This type of model requires multi-shell data to obtain more specific tissue microstructural information. However, a major concern with multi-shell is that it takes a longer scanning time restricting its use in clinical applications. In addition, most of the existing dMRI scanners with low gradient strengths commonly acquire a single b-value (shell) upto b=1000s/mm2 due to SNR (Signal-to-noise ratio) reasons and severe imaging artifacts.Methods: To address this issue, we propose a CNN-based ordinary differential equations solver for the reconstruction of MSMT fODF from under-sampled and fully sampled Single-shell (b=1000s/mm2) dMRI. The proposed architecture consists of CNN-based Adams-Bash-forth and Runge-Kutta modules along with two loss functions, including L1 and total variation.Results: We have shown quantitative results and visualization of fODF, fiber tracts, and structural connectivity for several brain regions on the publicly available HCP dataset. In addition, the obtained angular correlation coefficients for white matter and full brain are high, showing the proposed network’s utility.Finally, we have also demonstrated the effect of noise by adjusting SNR from 5 to 50 and observed the network robustness.Conclusion: We can conclude that our model can accurately predict MSMT fODF from under-sampled or fully sampled Single-shell dMRI volumes.

Reconstruction for 7T high-resolution whole-brain diffusion MRI using two-stage N/2 ghost correction and L1-SPIRiT without single-band reference.

To combine a new two-stage N/2 ghost correction and an adapted L1-SPIRiT method for reconstruction of 7T highly accelerated whole-brain diffusion MRI (dMRI) using only autocalibration scans (ACS) without the need of additional single-band reference (SBref) scans. The proposed ghost correction consisted of a 3-line reference approach in stage 1 and the reference-free entropy method in stage 2. The adapted L1-SPIRiT method was formulated within the 3D k-space framework. Its efficacy was examined by acquiring two dMRI data sets at 1.05-mm isotropic resolutions with a total acceleration of 6 or 9 (i.e., 2-fold or 3-fold slice and 3-fold in-plane acceleration). Diffusion analysis was performed to derive DTI metrics and estimate fiber orientation distribution functions (fODFs). The results were compared with those of 3D k-space GRAPPA using only ACS, all in reference to 3D k-space GRAPPA using both ACS and SBref (serving as a reference). The proposed ghost correction eliminated artifacts more robustly than conventional approaches. Our adapted L1-SPIRiT method outperformed 3D k-space GRAPPA when using only ACS, improving image quality to what was achievable with 3D k-space GRAPPA using both ACS and SBref scans. The improvement in image quality further resulted in an improvement in estimation performances for DTI and fODFs. The combination of our new ghost correction and adapted L1-SPIRiT method can reliably reconstruct 7T highly accelerated whole-brain dMRI without the need of SBref scans, increasing acquisition efficiency and reducing motion sensitivity.

Effect of fibre orientation on the light scattering during laser transmission welding

During laser welding of thermoplastic (LTW) composites reinforced with short fibres, a divergence of the laser beam is observed due to internal refraction of the beam at each matrix–fibre interface. This phenomenon leads to scattering of the laser beam in this heterogeneous medium, resulting in reduction of energy reaching the weld interface. This work presents a numerical study of the effect of fibre orientation in fibre-reinforced composites during the laser transmission welding simulation. A three dimensional (3D) numerical structure is generated to take into account the real microstructure of a composite material reinforced with short fibres. The information of fibre volume fraction, fibre length distribution and fibre orientation distribution is extracted from injection moulding simulations. Fibre orientation distribution function is extracted from fibre orientation tensors. An algorithm is implemented to trace rays propagating through composite material at two scales: microscopic and macroscopic. This algorithm uses optical properties of the composite material to simulate laser beam reflection and refraction in a complex structure. Laser beam scatter at the weld interface is simulated in 3D geometries with various fibre orientation distributions. The effect of fibre orientation on the light scattering phenomenon of laser is compared for different cases of fibre orientation. The simulation results are validated experimentally.

Investigating Tissue-Specific Abnormalities in Alzheimer's Disease with Multi-Shell Diffusion MRI.

Most studies using diffusion-weighted MRI (DW-MRI) in Alzheimer's disease (AD) have focused their analyses on white matter (WM) microstructural changes using the diffusion (kurtosis) tensor model. Although recent works have addressed some limitations of the tensor model, such as the representation of crossing fibers and partial volume effects with cerebrospinal fluid (CSF), the focus remains in modeling and analyzing the WM. In this work, we present a brain analysis approach for DW-MRI that disentangles multiple tissue compartments as well as micro- and macroscopic effects to investigate differences between groups of subjects in the AD continuum and controls. By means of the multi-tissue constrained spherical deconvolution of multi-shell DW-MRI, underlying brain tissue is modeled with a WM fiber orientation distribution function along with the contributions of gray matter (GM) and CSF to the diffusion signal. From this multi-tissue model, a set of measures capturing tissue diffusivity properties and morphology are extracted. Group differences were interrogated following fixel-, voxel-, and tensor-based morphometry approaches while including strong FWE control across multiple comparisons. Abnormalities related to AD stages were detected in WM tracts including the splenium, cingulum, longitudinal fasciculi, and corticospinal tract. Changes in tissue composition were identified, particularly in the medial temporal lobe and superior longitudinal fasciculus. This analysis framework constitutes a comprehensive approach allowing simultaneous macro and microscopic assessment of WM, GM, and CSF, from a single DW-MRI dataset.

CTtrack: A CNN+Transformer-based framework for fiber orientation estimation & tractography

In recent years, multiple data-driven fiber orientation distribution function (fODF) estimation algorithms and automatic tractography pipelines have been proposed to address the limitations of traditional methods. However, these approaches lack precision and generalizability. To tackle these shortcomings, we introduce CTtrack, a CNN+Transformer-based pipeline to estimate fODFs and perform tractography. In this approach, a convolutional neural network (CNN) module is employed to project the resampled diffusion-weighted magnetic resonance imaging (DW-MRI) data to a lower dimension. Then, a transformer model estimates the fiber orientation distribution functions using the projected data within a local block around each voxel. The proposed model represents the extracted fODFs by spherical harmonics coefficients. The predicted fiber ODFs can be used for both deterministic and probabilistic tractography. Our pipeline was tested in terms of the precision and robustness in estimating fODFs and performing tractography using both simulated and real diffusion data. The Tractometer tool was employed to compare our method with the classical and data-driven tractography approaches. The qualitative and quantitative assessments illustrate the competitive performance of our framework compared to other available algorithms.

A discrete fiber network finite element model of arterial elastin network considering inter-fiber crosslinking property and density

Inter-fiber crosslinks within the extracellular matrix (ECM) play important roles in determining the mechanical properties of the fibrous network. Discrete fiber network (DFN) models have been used to study fibrous biological material, however the contribution of inter-fiber crosslinks to the mechanics of the ECM network is not well understood. In this study, a DFN model of arterial elastin network was developed based on measured structural features to study the contribution of inter-fiber crosslinking properties and density to the mechanics and fiber kinematics of the network. The DFN was generated by randomly placing line segments into a given domain following a fiber orientation distribution function obtained from multiphoton microscopy until a desired fiber areal fraction was reached. Intersections between the line segments were treated as crosslinks. The generated DFN model was then incorporated into an ABAQUS finite element model to simulate the network under equi- and nonequi-biaxial deformation. The inter-fiber crosslinks were modeled using connector elements with either zero (pin joint) or infinite (weld joint) rotational stiffness. Furthermore, inter-fiber crosslinking density was systematically reduced and its effect on both network- and fiber-level mechanics was studied. The DFN model showed good fitting and predicting capabilities of the stress-strain behavior of the elastin network. While the pin and weld joints do not seem to have noticeable effect on the network stress-strain behavior, the crosslinking properties can affect the local fiber mechanics and kinematics. Overall, our study suggests that inter-fiber crosslinking properties are important to the multiscale mechanics and fiber kinematics of the ECM network.

Bridging the gap between constrained spherical deconvolution and diffusional variance decomposition via tensor-valued diffusion MRI.

Diffusion tensor imaging (DTI) is widely used to extract valuable tissue measurements and white matter (WM) fiber orientations, even though its lack of specificity is now well-known, especially for WM fiber crossings. Models such as constrained spherical deconvolution (CSD) take advantage of high angular resolution diffusion imaging (HARDI) data to compute fiber orientation distribution functions (fODF) and tackle the orientational part of the DTI limitations. Furthermore, the recent introduction of tensor-valued diffusion MRI allows for diffusional variance decomposition (DIVIDE), enabling the computation of measures more specific to microstructure than DTI measures, such as microscopic fractional anisotropy (μFA). Recent work on making CSD compatible with tensor-valued diffusion MRI data opens the door for methods combining CSD and DIVIDE to get both fODFs and microstructure measures. However, the impacts of such atypical data on fODF reconstruction with CSD are yet to be fully known and understood. In this work, we use simulated data to explore the effects of various combinations of diffusion encodings on the angular resolution of extracted fOFDs and on the versatility of CSD in multiple realistic situations. We also compare the combinations with regards to their performance at producing accurate and precise μFA with DIVIDE, and present an optimized 10min protocol combining linear and spherical b-tensor encodings for both methods. We show that our proposed protocol enables the reconstruction of both fODFs and μFA on in vivo data.