Versatile, robust, and efficient tractography with constrained higher-order tensor fODFs

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.

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.

Comparison of basis functions and q-space sampling schemes for robust compressed sensing reconstruction accelerating diffusion spectrum imaging

Time constraints placed on magnetic resonance imaging often restrict the application of advanced diffusion MRI (dMRI) protocols in clinical practice and in high throughput research studies. Therefore, acquisition strategies for accelerated dMRI have been investigated to allow for the collection of versatile and high quality imaging data, even if stringent scan time limits are imposed. Diffusion spectrum imaging (DSI), an advanced acquisition strategy that allows for a high resolution of intra-voxel microstructure, can be sufficiently accelerated by means of compressed sensing (CS) theory. CS theory describes a framework for the efficient collection of fewer samples of a data set than conventionally required followed by robust reconstruction to recover the full data set from sparse measurements. For an accurate recovery of DSI data, a suitable acquisition scheme for sparse q-space sampling and the sensing and sparsifying bases for CS reconstruction need to be selected. In this work we explore three different types of q-space undersampling schemes and two frameworks for CS reconstruction based on either Fourier or SHORE basis functions. After CS recovery, diffusion and microstructural parameters and orientational information are estimated from the reconstructed data by means of state-of-the-art processing techniques for dMRI analysis. By means of simulation, diffusion phantom and in vivo DSI data, an isotropic distribution of q-space samples was found to be optimal for sparse DSI. The CS reconstruction results indicate superior performance of Fourier-based CS-DSI compared to the SHORE-based approach. Based on these findings we outline an experimental design for accelerated DSI and robust CS reconstruction of the sparse measurements that is suitable for the application within time-limited studies.

Compressed Sensing Diffusion Spectrum Imaging for Accelerated Diffusion Microstructure MRI in Long-Term Population Imaging

Mapping non-invasively the complex microstructural architecture of the living human brain, diffusion magnetic resonance imaging (dMRI) is one of the core imaging modalities in current population studies. For the application in longitudinal population imaging, the dMRI protocol should deliver reliable data with maximum potential for future analysis. With the recent introduction of novel MRI hardware, advanced dMRI acquisition strategies can be applied within reasonable scan time. In this work we conducted a pilot study based on the requirements for high resolution dMRI in a long-term and high throughput population study. The key question was: can diffusion spectrum imaging accelerated by compressed sensing theory (CS-DSI) be used as an advanced imaging protocol for microstructure dMRI in a long-term population imaging study? As a minimum requirement we expected a high level of agreement of several diffusion metrics derived from both CS-DSI and a 3-shell high angular resolution diffusion imaging (HARDI) acquisition, an established imaging strategy used in other population studies. A wide spectrum of state-of-the-art diffusion processing and analysis techniques was applied to the pilot study data including quantitative diffusion and microstructural parameter mapping, fiber orientation estimation and white matter fiber tracking. When considering diffusion weighted images up to the same maximum diffusion weighting for both protocols, group analysis across 20 subjects indicates that CS-DSI performs comparable to 3-shell HARDI in the estimation of diffusion and microstructural parameters. Further, both protocols provide similar results in the estimation of fiber orientations and for local fiber tracking. CS-DSI provides high radial resolution while maintaining high angular resolution and it is well-suited for analysis strategies that require high b-value acquisitions, such as CHARMED modeling and biomarkers from the diffusion propagator.

Differentiation of High-Grade Glioma and Intracranial Metastasis Using Volumetric Diffusion Tensor Imaging Tractography

OBJECTIVE

A reliable, noninvasive method to differentiate high-grade glioma (HGG) and intracranial metastasis (IM) has remained elusive. The aim of this study was to differentiate between HGG and IM using tumoral and peritumoral diffusion tensor imaging characteristics.

METHODS

A semiautomated script generated volumetric regions of interest (ROIs) for the tumor and a peritumoral shell at a predetermined voxel thickness. ROI differences in diffusion tensor imaging-related metrics between HGG and IM groups were estimated, including fractional anisotropy, mean diffusivity, total fiber tract counts, and tract density.

RESULTS

The HGG group (n = 46) had a significantly higher tumor-to-brain volume ratio than the IM group (n = 35) (P < 0.001). The HGG group exhibited significantly higher mean fractional anisotropy and significantly lower mean diffusivity within peritumoral ROI than the IM group (P < 0.05). The HGG group exhibited significantly higher total tract count and higher tract density in tumoral and peritumoral ROIs than the IM group (P < 0.05). Tumoral tract count and peritumoral tract density were the most optimal metrics to differentiate the groups based on receiver operating characteristic curve analysis. Predictive analysis using receiver operating characteristic curve thresholds was performed on 13 additional participants. Compared with correct clinical diagnoses, the 2 thresholds exhibited equal specificities (66.7%), but the tumoral tract count (85.7%) seemed more sensitive in differentiating the 2 groups.

CONCLUSIONS

Tract count and tract density were significantly different in tumoral and peritumoral regions between HGG and IM. Differences in microenvironmental interactions between the tumor types may cause these tract differences.

DT-MRI Streamsurfaces Revisited

DT-MRI streamsurfaces, defined as surfaces that are everywhere tangential to the major and medium eigenvector fields, have been proposed as a tool for visualizing regions of predominantly planar behavior in diffusion tensor MRI. Even though it has long been known that their construction assumes that the involved eigenvector fields satisfy an integrability condition, it has never been tested systematically whether this condition is met in real-world data. We introduce a suitable and efficiently computable test to the visualization literature, demonstrate that it can be used to distinguish integrable from nonintegrable configurations in simulations, and apply it to whole-brain datasets of 15 healthy subjects. We conclude that streamsurface integrability is approximately satisfied in a substantial part of the brain, but not everywhere, including some regions of planarity. As a consequence, algorithms for streamsurface extraction should explicitly test local integrability. Finally, we propose a novel patch-based approch to streamsurface visualization that reduces visual artifacts, and is shown to more fully sample the extent of streamsurfaces.