Fiber-driven resolution enhancement of diffusion-weighted images

Diffusion-weighted imaging (DWI), while giving rich information about brain circuitry, is often limited by insufficient spatial resolution and low signal-to-noise ratio (SNR). This paper describes an algorithm that will increase the resolution of DW images beyond the scan resolution, allowing for a closer investigation of fiber structures and more accurate assessment of brain connectivity. The algorithm is capable of generating a dense vector-valued field, consisting of diffusion data associated with the full set of diffusion-sensitizing gradients. The fundamental premise is that, to best preserve information, interpolation should always be performed along axonal fibers. To achieve this, at each spatial location, we probe neighboring voxels in various directions to gather diffusion information for data interpolation. Based on the fiber orientation distribution function (ODF), directions that are more likely to be traversed by fibers will be given greater weights during interpolation and vice versa. This ensures that data interpolation is only contributed by diffusion data coming from fibers that are aligned with a specific direction. This approach respects local fiber structures and prevents blurring resulting from averaging of data from significantly misaligned fibers. Evaluations suggest that this algorithm yields results with significantly less blocking artifacts, greater smoothness in anatomical structures, and markedly improved structural visibility

Improving Fiber Alignment in HARDI by Combining Contextual PDE Flow with Constrained Spherical Deconvolution

We propose two strategies to improve the quality of tractography results computed from diffusion weighted magnetic resonance imaging (DW-MRI) data. Both methods are based on the same PDE framework, defined in the coupled space of positions and orientations, associated with a stochastic process describing the enhancement of elongated structures while preserving crossing structures. In the first method we use the enhancement PDE for contextual regularization of a fiber orientation distribution (FOD) that is obtained on individual voxels from high angular resolution diffusion imaging (HARDI) data via constrained spherical deconvolution (CSD). Thereby we improve the FOD as input for subsequent tractography. Secondly, we introduce the fiber to bundle coherence (FBC), a measure for quantification of fiber alignment. The FBC is computed from a tractography result using the same PDE framework and provides a criterion for removing the spurious fibers. We validate the proposed combination of CSD and enhancement on phantom data and on human data, acquired with different scanning protocols. On the phantom data we find that PDE enhancements improve both local metrics and global metrics of tractography results, compared to CSD without enhancements. On the human data we show that the enhancements allow for a better reconstruction of crossing fiber bundles and they reduce the variability of the tractography output with respect to the acquisition parameters. Finally, we show that both the enhancement of the FODs and the use of the FBC measure on the tractography improve the stability with respect to different stochastic realizations of probabilistic tractography. This is shown in a clinical application: the reconstruction of the optic radiation for epilepsy surgery planning

Left-Invariant Diffusion on the Motion Group in terms of the Irreducible Representations of SO(3)

In this work we study the formulation of convection/diffusion equations on the 3D motion group SE(3) in terms of the irreducible representations of SO(3). Therefore, the left-invariant vector-fields on SE(3) are expressed as linear operators, that are differential forms in the translation coordinate and algebraic in the rotation. In the context of 3D image processing this approach avoids the explicit discretization of SO(3) or $S_2$, respectively. This is particular important for SO(3), where a direct discretization is infeasible due to the enormous memory consumption. We show two applications of the framework: one in the context of diffusion-weighted magnetic resonance imaging and one in the context of object detection

Quantitative Perfusion-Sensitive Mri Phantoms

Perfusion-sensitive MR methods are increasingly utilized in preclinical and clinical MR research studies with the promise of providing quantitative estimates of parameters that describe in vivo microvasculature. One of these techniques, dynamic contrast enhanced: DCE) MRI, has found particularly common use in oncology for the detection, staging, and monitoring of highly vascularized tumors. DCE-MRI has been qualitatively validated by various studies that show a high correlation between modeled parameters from DCE and histologically measured microvascular density: MVD). However, in the absence of a matching gold-standard technique, DCE-MRI has not yet been quantitatively validated: i.e., the accuracy of the estimated parameters is unknown). Partly because of this inability to determine the accuracy of the measured parameters, there remains debate in the literature about which DCE signal model(s) best reflect(s) experimental data. In order to address these scientific challenges, realistic DCE tissue phantoms have been constructed. These phantoms implement semi-permeable hollow fibers, found commonly in commercial hemodialysis cartridges, to simulate leaky vasculature. Their design and construction are cataloged in this thesis. In addition, the phantoms have been experimentally characterized. In conjunction with these experiments, an interesting example of diffusion driven longitudinal relaxation was observed and is described herein. Lastly, the permeability of the fiber wall with respect to Gd-based contrast agents has been measured independently and compared with values derived from a mock-DCE experiment performed on the phantoms. In general, the results of these experiments support current DCE-MRI methods

Axon diameters and myelin content modulate microscopic fractional anisotropy at short diffusion times in fixed rat spinal cord

Mapping tissue microstructure accurately and noninvasively is one of the frontiers of biomedical imaging. Diffusion Magnetic Resonance Imaging (MRI) is at the forefront of such efforts, as it is capable of reporting on microscopic structures orders of magnitude smaller than the voxel size by probing restricted diffusion. Double Diffusion Encoding (DDE) and Double Oscillating Diffusion Encoding (DODE) in particular, are highly promising for their ability to report on microscopic fractional anisotropy ({\mu}FA), a measure of the pore anisotropy in its own eigenframe, irrespective of orientation distribution. However, the underlying correlates of {\mu}FA have insofar not been studied. Here, we extract {\mu}FA from DDE and DODE measurements at ultrahigh magnetic field of 16.4T in the aim to probe fixed rat spinal cord microstructure. We further endeavor to correlate {\mu}FA with Myelin Water Fraction (MWF) derived from multiexponential T2 relaxometry, as well as with literature-based spatially varying axonal diameters. In addition, a simple new method is presented for extracting unbiased {\mu}FA from three measurements at different b-values. Our findings reveal strong anticorrelations between {\mu}FA (derived from DODE) and axon diameter in the distinct spinal cord tracts; a moderate correlation was also observed between {\mu}FA derived from DODE and MWF. These findings suggest that axonal membranes strongly modulate {\mu}FA, which - owing to its robustness towards orientation dispersion effects - reflects axon diameter much better than its typical FA counterpart. The {\mu}FA exhibited modulations when measured via oscillating or blocked gradients, suggesting selective probing of different parallel path lengths and providing insight into how those modulate {\mu}FA metrics. Our findings thus shed light into the underlying microstructural correlates of {\mu}FA and are (...

Multiple indices of diffusion identifies white matter damage in mild cognitive impairment and Alzheimer's disease

The study of multiple indices of diffusion, including axial (DA), radial (DR) and mean diffusion (MD), as well as fractional anisotropy (FA), enables WM damage in Alzheimer's disease (AD) to be assessed in detail. Here, tract-based spatial statistics (TBSS) were performed on scans of 40 healthy elders, 19 non-amnestic MCI (MCIna) subjects, 14 amnestic MCI (MCIa) subjects and 9 AD patients. Significantly higher DA was found in MCIna subjects compared to healthy elders in the right posterior cingulum/precuneus. Significantly higher DA was also found in MCIa subjects compared to healthy elders in the left prefrontal cortex, particularly in the forceps minor and uncinate fasciculus. In the MCIa versus MCIna comparison, significantly higher DA was found in large areas of the left prefrontal cortex. For AD patients, the overlap of FA and DR changes and the overlap of FA and MD changes were seen in temporal, parietal and frontal lobes, as well as the corpus callosum and fornix. Analysis of differences between the AD versus MCIna, and AD versus MCIa contrasts, highlighted regions that are increasingly compromised in more severe disease stages. Microstructural damage independent of gross tissue loss was widespread in later disease stages. Our findings suggest a scheme where WM damage begins in the core memory network of the temporal lobe, cingulum and prefrontal regions, and spreads beyond these regions in later stages. DA and MD indices were most sensitive at detecting early changes in MCIa