Diffusion tensor imaging: structural adaptive smoothing

Non-local means based Rician noise filtering for diffusion tensor and kurtosis imaging in human brain and spinal cord

BACKGROUND

To investigate the effect of using a Rician nonlocal means (NLM) filter on quantification of diffusion tensor (DT)- and diffusion kurtosis (DK)-derived metrics in various anatomical regions of the human brain and the spinal cord, when combined with a constrained linear least squares (CLLS) approach.

METHODS

Prospective brain data from 9 healthy subjects and retrospective spinal cord data from 5 healthy subjects from a 3 T MRI scanner were included in the study. Prior to tensor estimation, registered diffusion weighted images were denoised by an optimized blockwise NLM filter with CLLS. Mean kurtosis (MK), radial kurtosis (RK), axial kurtosis (AK), mean diffusivity (MD), radial diffusivity (RD), axial diffusivity (AD) and fractional anisotropy (FA), were determined in anatomical structures of the brain and the spinal cord. DTI and DKI metrics, signal-to-noise ratio (SNR) and Chi-square values were quantified in distinct anatomical regions for all subjects, with and without Rician denoising.

RESULTS

The averaged SNR significantly increased with Rician denoising by a factor of 2 while the averaged Chi-square values significantly decreased up to 61% in the brain and up to 43% in the spinal cord after Rician NLM filtering. In the brain, the mean MK varied from 0.70 (putamen) to 1.27 (internal capsule) while AK and RK varied from 0.58 (corpus callosum) to 0.92 (cingulum) and from 0.70 (putamen) to 1.98 (corpus callosum), respectively. In the spinal cord, FA varied from 0.78 in lateral column to 0.81 in dorsal column while MD varied from 0.91 × 10^-3 mm²/s (lateral) to 0.93 × 10^-3 mm²/s (dorsal). RD varied from 0.34 × 10^-3 mm²/s (dorsal) to 0.38 × 10^-3 mm²/s (lateral) and AD varied from 1.96 × 10^-3 mm²/s (lateral) to 2.11 × 10^-3 mm²/s (dorsal).

CONCLUSIONS

Our results show a Rician denoising NLM filter incorporated with CLLS significantly increases SNR and reduces estimation errors of DT- and KT-derived metrics, providing the reliable metrics estimation with adequate SNR levels.

Longitudinal Reproducibility of Neurite Orientation Dispersion and Density Imaging (NODDI) Derived Metrics in the White Matter

Diffusion-weighted magnetic resonance imaging (DWI) is undergoing constant evolution with the ambitious goal of developing in-vivo histology of the brain. A recent methodological advancement is Neurite Orientation Dispersion and Density Imaging (NODDI), a histologically validated multi-compartment model to yield microstructural features of brain tissue such as geometric complexity and neurite packing density, which are especially useful in imaging the white matter. Since NODDI is increasingly popular in clinical research and fields such as developmental neuroscience and neuroplasticity, it is of vast importance to characterize its reproducibility (or reliability). We acquired multi-shell DWI data in 29 healthy young subjects twice over a rescan interval of 4 weeks to assess the within-subject coefficient of variation (CV_WS), between-subject coefficient of variation (CV_BS) and the intraclass correlation coefficient (ICC), respectively. Using these metrics, we compared regional and voxel-by-voxel reproducibility of the most common image analysis approaches (tract-based spatial statistics [TBSS], voxel-based analysis with different extents of smoothing ["VBM-style"], ROI-based analysis). We observed high test-retest reproducibility for the orientation dispersion index (ODI) and slightly worse results for the neurite density index (NDI). Our findings also suggest that the choice of analysis approach might have significant consequences for the results of a study. Collectively, the voxel-based approach with Gaussian smoothing kernels of ≥4 mm FWHM and ROI-averaging yielded the highest reproducibility across NDI and ODI maps (CV_WS mostly ≤3%, ICC mostly ≥0.8), respectively, whilst smaller kernels and TBSS performed consistently worse. Furthermore, we demonstrate that image quality (signal-to-noise ratio [SNR]) is an important determinant of NODDI metric reproducibility. We discuss the implications of these results for longitudinal and cross-sectional research designs commonly employed in the neuroimaging field.

Denoising for Improved Parametric MRI of the Kidney: Protocol for Nonlocal Means Filtering

In order to tackle the challenges caused by the variability in estimated MRI parameters (e.g., T2* and T2) due to low SNR a number of strategies can be followed. One approach is postprocessing of the acquired data with a filter. The basic idea is that MR images possess a local spatial structure that is characterized by equal, or at least similar, noise-free signal values in vicinities of a location. Then, local averaging of the signal reduces the noise component of the signal. In contrast, nonlocal means filtering defines the weights for averaging not only within the local vicinity, bur it compares the image intensities between all voxels to define "nonlocal" weights. Furthermore, it generally compares not only single-voxel intensities but small spatial patches of the data to better account for extended similar patterns. Here we describe how to use an open source NLM filter tool to denoise 2D MR image series of the kidney used for parametric mapping of the relaxation times T2* and T2.This chapter is based upon work from the COST Action PARENCHIMA, a community-driven network funded by the European Cooperation in Science and Technology (COST) program of the European Union, which aims to improve the reproducibility and standardization of renal MRI biomarkers.

Structure Prior Constrained Estimation of Human Cardiac Diffusion Tensors

OBJECTIVE

The purpose of this paper is to increase the accuracy of human cardiac diffusion tensor (DT) estimation in diffusion magnetic resonance imaging (dMRI) with a few diffusion gradient directions.

METHODS

A structure prior constrained (SPC) method is proposed. The method consists in introducing two regularizers in the conventional nonlinear least squares estimator. The two regularizers penalize the dissimilarity between neighboring DTs and the difference between estimated and prior fiber orientations, respectively. A novel numerical solution is presented to ensure the positive definite estimation.

RESULTS

Experiments on ex vivo human cardiac data show that the SPC method is able to well estimate DTs at most voxels, and is superior to state-of-the-art methods in terms of the mean errors of principal eigenvector, second eigenvector, helix angle, transverse angle, fractional anisotropy, and mean diffusivity.

CONCLUSION

The SPC method is a practical and reliable alternative to current denoising- or regularization-based methods for the estimation of human cardiac DT.

SIGNIFICANCE

The SPC method is able to accurately estimate human cardiac DTs in dMRI with a few diffusion gradient directions.

SR-HARDI: Spatially Regularizing High Angular Resolution Diffusion Imaging

High angular resolution diffusion imaging (HARDI) has recently been of great interest in mapping the orientation of intra-voxel crossing fibers, and such orientation information allows one to infer the connectivity patterns prevalent among different brain regions and possible changes in such connectivity over time for various neurodegenerative and neuropsychiatric diseases. The aim of this paper is to propose a penalized multi-scale adaptive regression model (PMARM) framework to spatially and adaptively infer the orientation distribution function (ODF) of water diffusion in regions with complex fiber configurations. In PMARM, we reformulate the HARDI imaging reconstruction as a weighted regularized least-squares regression (WRLSR) problem. Similarity and distance weights are introduced to account for spatial smoothness of HARDI, while preserving the unknown discontinuities (e.g., edges between white matter and grey matter) of HARDI. The L₁ penalty function is introduced to ensure the sparse solutions of ODFs, while a scaled L₁ weighted estimator is calculated to correct the bias introduced by the L₁ penalty at each voxel. In PMARM, we integrate the multiscale adaptive regression models (Li et al., 2011), the propagation-separation method (Polzehl and Spokoiny, 2000), and Lasso (least absolute shrinkage and selection operator) (Tibshirani, 1996) to adaptively estimate ODFs across voxels. Experimental results indicate that PMARM can reduce the angle detection errors on fiber crossing area and provide more accurate reconstruction than standard voxel-wise methods.

POAS4SPM: a toolbox for SPM to denoise diffusion MRI data

We present an implementation of a recently developed noise reduction algorithm for dMRI data, called multi-shell position orientation adaptive smoothing (msPOAS), as a toolbox for SPM. The method intrinsically adapts to the structures of different size and shape in dMRI and hence avoids blurring typically observed in non-adaptive smoothing. We give examples for the usage of the toolbox and explain the determination of experiment-dependent parameters for an optimal performance of msPOAS.

Evaluation of non-local means based denoising filters for diffusion kurtosis imaging using a new phantom

Image denoising has a profound impact on the precision of estimated parameters in diffusion kurtosis imaging (DKI). This work first proposes an approach to constructing a DKI phantom that can be used to evaluate the performance of denoising algorithms in regard to their abilities of improving the reliability of DKI parameter estimation. The phantom was constructed from a real DKI dataset of a human brain, and the pipeline used to construct the phantom consists of diffusion-weighted (DW) image filtering, diffusion and kurtosis tensor regularization, and DW image reconstruction. The phantom preserves the image structure while minimizing image noise, and thus can be used as ground truth in the evaluation. Second, we used the phantom to evaluate three representative algorithms of non-local means (NLM). Results showed that one scheme of vector-based NLM, which uses DWI data with redundant information acquired at different b-values, produced the most reliable estimation of DKI parameters in terms of Mean Square Error (MSE), Bias and standard deviation (Std). The result of the comparison based on the phantom was consistent with those based on real datasets.

Spatially Varying Coefficient Model for Neuroimaging Data with Jump Discontinuities

Motivated by recent work on studying massive imaging data in various neuroimaging studies, we propose a novel spatially varying coefficient model (SVCM) to capture the varying association between imaging measures in a three-dimensional (3D) volume (or 2D surface) with a set of covariates. Two stylized features of neuorimaging data are the presence of multiple piecewise smooth regions with unknown edges and jumps and substantial spatial correlations. To specifically account for these two features, SVCM includes a measurement model with multiple varying coefficient functions, a jumping surface model for each varying coefficient function, and a functional principal component model. We develop a three-stage estimation procedure to simultaneously estimate the varying coefficient functions and the spatial correlations. The estimation procedure includes a fast multiscale adaptive estimation and testing procedure to independently estimate each varying coefficient function, while preserving its edges among different piecewise-smooth regions. We systematically investigate the asymptotic properties (e.g., consistency and asymptotic normality) of the multiscale adaptive parameter estimates. We also establish the uniform convergence rate of the estimated spatial covariance function and its associated eigenvalues and eigenfunctions. Our Monte Carlo simulation and real data analysis have confirmed the excellent performance of SVCM.

Brain imaging of mild cognitive impairment and Alzheimer's disease

The rapidly increasing prevalence of cognitive impairment and Alzheimer's disease has the potential to create a major worldwide healthcare crisis. Structural MRI studies in patients with Alzheimer's disease and mild cognitive impairment are currently attracting considerable interest. It is extremely important to study early structural and metabolic changes, such as those in the hippocampus, entorhinal cortex, and gray matter structures in the medial temporal lobe, to allow the early detection of mild cognitive impairment and Alzheimer's disease. The microstructural integrity of white matter can be studied with diffusion tensor imaging. Increased mean diffusivity and decreased fractional anisotropy are found in subjects with white matter damage. Functional imaging studies with positron emission tomography tracer compounds enable detection of amyloid plaques in the living brain in patients with Alzheimer's disease. In this review, we will focus on key findings from brain imaging studies in mild cognitive impairment and Alzheimer's disease, including structural brain changes studied with MRI and white matter changes seen with diffusion tensor imaging, and other specific imaging methodologies will also be discussed.

Non-local mean denoising in diffusion tensor space

The aim of the present study was to present a novel non-local mean (NLM) method to denoise diffusion tensor imaging (DTI) data in the tensor space. Compared with the original NLM method, which uses intensity similarity to weigh the voxel, the proposed method weighs the voxel using tensor similarity measures in the diffusion tensor space. Euclidean distance with rotational invariance, and Riemannian distance and Log-Euclidean distance with affine invariance were implemented to compare the geometric and orientation features of the diffusion tensor comprehensively. The accuracy and efficacy of the proposed novel NLM method using these three similarity measures in DTI space, along with unbiased novel NLM in diffusion-weighted image space, were compared quantitatively and qualitatively in the present study.

SGPP: spatial Gaussian predictive process models for neuroimaging data

The aim of this paper is to develop a spatial Gaussian predictive process (SGPP) framework for accurately predicting neuroimaging data by using a set of covariates of interest, such as age and diagnostic status, and an existing neuroimaging data set. To achieve a better prediction, we not only delineate spatial association between neuroimaging data and covariates, but also explicitly model spatial dependence in neuroimaging data. The SGPP model uses a functional principal component model to capture medium-to-long-range (or global) spatial dependence, while SGPP uses a multivariate simultaneous autoregressive model to capture short-range (or local) spatial dependence as well as cross-correlations of different imaging modalities. We propose a three-stage estimation procedure to simultaneously estimate varying regression coefficients across voxels and the global and local spatial dependence structures. Furthermore, we develop a predictive method to use the spatial correlations as well as the cross-correlations by employing a cokriging technique, which can be useful for the imputation of missing imaging data. Simulation studies and real data analysis are used to evaluate the prediction accuracy of SGPP and show that SGPP significantly outperforms several competing methods, such as voxel-wise linear model, in prediction. Although we focus on the morphometric variation of lateral ventricle surfaces in a clinical study of neurodevelopment, it is expected that SGPP is applicable to other imaging modalities and features.

Smoothing fields of weighted collections with applications to diffusion MRI processing

Using modern diffusion weighted magnetic resonance imaging protocols, the orientations of multiple neuronal fiber tracts within each voxel can be estimated. Further analysis of these populations, including application of fiber tracking and tract segmentation methods, is often hindered by lack of spatial smoothness of the estimated orientations. For example, a single noisy voxel can cause a fiber tracking method to switch tracts in a simple crossing tract geometry. In this work, a generalized spatial smoothing framework that handles multiple orientations as well as their fractional contributions within each voxel is proposed. The approach estimates an optimal fuzzy correspondence of orientations and fractional contributions between voxels and smooths only between these correspondences. Avoiding a requirement to obtain exact correspondences of orientations reduces smoothing anomalies due to propagation of erroneous correspondences around noisy voxels. Phantom experiments are used to demonstrate both visual and quantitative improvements in postprocessing steps. Improvement over smoothing in the measurement domain is also demonstrated using both phantoms and in vivo human data.

Feature-preserving smoothing of diffusion weighted images using nonstationarity adaptive filtering

Although promising for studying the microstructure of in vivo tissues, the performance and the potentiality of diffusion tensor magnetic resonance imaging are hampered by the presence of high-level noise in diffusion weighted (DW) images. This paper proposes a novel smoothing approach, called the nonstationarity adaptive filtering, which estimates the intensity of a pixel by averaging intensities in its adaptive homogeneous neighborhood. The latter is determined according to five constraints and spatiodirectional nonstationarity measure maps. The proposed approach is compared with an anisotropic diffusion method used in DW image smoothing. Experimental results on both synthetic and real human DW images show that the proposed method achieves a better compromise between the smoothness of homogeneous regions and the preservation of desirable features such as boundaries, even for highly noisy data, thus leading to homogeneously consistent tensor fields and consequently more coherent fibers.

Local Tests for Identifying Anisotropic Diffusion Areas in Human Brain with DTI

Diffusion tensor imaging (DTI) plays a key role in analyzing the physical structures of biological tissues, particularly in reconstructing fiber tracts of the human brain in vivo. On the one hand, eigenvalues of diffusion tensors (DTs) estimated from diffusion weighted imaging (DWI) data usually contain systematic bias, which subsequently biases the diffusivity measurements popularly adopted in fiber tracking algorithms. On the other hand, correctly accounting for the spatial information is important in the construction of these diffusivity measurements since the fiber tracts are typically spatially structured. This paper aims to establish test-based approaches to identify anisotropic water diffusion areas in the human brain. These areas in turn indicate the areas passed by fiber tracts. Our proposed test statistic not only takes into account the bias components in eigenvalue estimates, but also incorporates the spatial information of neighboring voxels. Under mild regularity conditions, we demonstrate that the proposed test statistic asymptotically follows a χ² distribution under the null hypothesis. Simulation and real DTI data examples are provided to illustrate the efficacy of our proposed methods.

Multiscale adaptive generalized estimating equations for longitudinal neuroimaging data

Many large-scale longitudinal imaging studies have been or are being widely conducted to better understand the progress of neuropsychiatric and neurodegenerative disorders and normal brain development. The goal of this article is to develop a multiscale adaptive generalized estimation equation (MAGEE) method for spatial and adaptive analysis of neuroimaging data from longitudinal studies. MAGEE is applicable to making statistical inference on regression coefficients in both balanced and unbalanced longitudinal designs and even in twin and familial studies, whereas standard software platforms have several major limitations in handling these complex studies. Specifically, conventional voxel-based analyses in these software platforms involve Gaussian smoothing imaging data and then independently fitting a statistical model at each voxel. However, the conventional smoothing methods suffer from the lack of spatial adaptivity to the shape and spatial extent of region of interest and the arbitrary choice of smoothing extent, while independently fitting statistical models across voxels does not account for the spatial properties of imaging observations and noise distribution. To address such drawbacks, we adapt a powerful propagation-separation (PS) procedure to sequentially incorporate the neighboring information of each voxel and develop a new novel strategy to solely update a set of parameters of interest, while fixing other nuisance parameters at their initial estimators. Simulation studies and real data analysis show that MAGEE significantly outperforms voxel-based analysis.

Diffusion tensor smoothing through weighted Karcher means

Diffusion tensor magnetic resonance imaging (MRI) quantifies the spatial distribution of water Diffusion at each voxel on a regular grid of locations in a biological specimen by Diffusion tensors- 3 × 3 positive definite matrices. Removal of noise from DTI is an important problem due to the high scientific relevance of DTI and relatively low signal to noise ratio it provides. Leading approaches to this problem amount to estimation of weighted Karcher means of Diffusion tensors within spatial neighborhoods, under various metrics imposed on the space of tensors. However, it is unclear how the behavior of these estimators varies with the magnitude of DTI sensor noise (the noise resulting from the thermal e!ects of MRI scanning) as well as the geometric structure of the underlying Diffusion tensor neighborhoods. In this paper, we combine theoretical analysis, empirical analysis of simulated DTI data, and empirical analysis of real DTI scans to compare the noise removal performance of three kernel-based DTI smoothers that are based on Euclidean, log-Euclidean, and affine-invariant metrics. The results suggest, contrary to conventional wisdom, that imposing a simplistic Euclidean metric may in fact provide comparable or superior noise removal, especially in relatively unstructured regions and/or in the presence of moderate to high levels of sensor noise. On the contrary, log-Euclidean and affine-invariant metrics may lead to better noise removal in highly structured anatomical regions, especially when the sensor noise is of low magnitude. These findings emphasize the importance of considering the interplay of sensor noise magnitude and tensor field geometric structure when assessing Diffusion tensor smoothing options. They also point to the necessity for continued development of smoothing methods that perform well across a large range of scenarios.

TwinMARM: two-stage multiscale adaptive regression methods for twin neuroimaging data

Twin imaging studies have been valuable for understanding the relative contribution of the environment and genes on brain structures and their functions. Conventional analyses of twin imaging data include three sequential steps: spatially smoothing imaging data, independently fitting a structural equation model at each voxel, and finally correcting for multiple comparisons. However, conventional analyses are limited due to the same amount of smoothing throughout the whole image, the arbitrary choice of smoothing extent, and the decreased power in detecting environmental and genetic effects introduced by smoothing raw images. The goal of this paper is to develop a two-stage multiscale adaptive regression method (TwinMARM) for spatial and adaptive analysis of twin neuroimaging and behavioral data. The first stage is to establish the relationship between twin imaging data and a set of covariates of interest, such as age and gender. The second stage is to disentangle the environmental and genetic influences on brain structures and their functions. In each stage, TwinMARM employs hierarchically nested spheres with increasing radii at each location and then captures spatial dependence among imaging observations via consecutively connected spheres across all voxels. Simulation studies show that our TwinMARM significantly outperforms conventional analyses of twin imaging data. Finally, we use our method to detect statistically significant effects of genetic and environmental variations on white matter structures in a neonatal twin study.

Progression of microstructural putamen alterations in a case of symptomatic recurrent seizures using diffusion tensor imaging

Microstructural alterations of the putamen were recently reported in patients with partial and generalized epilepsy disorders. However, it is unknown whether these alterations pre-exist or are secondary to recurrent seizures. Here we investigated the progression of putamen fractional anisotropy (FA) alterations in a case of recurrent psychomotor seizures using longitudinal diffusion tensor imaging (DTI) shortly before (DTI-1) and after a psychomotor seizure (DTI-2). We obtained FA values of a hypothesis-guided putamen region-of-interest (ROI) and seven exploratory ROIs. FA values from both DTIs were compared with reference values from 19 controls. Relative to controls, the patient's putamen FA was increased at DTI-1 (13% left putamen, 7% right putamen), an effect that was exacerbated at DTI-2 (24% left putamen (p < 0.05), 20% right putamen). In the exploratory ROIs we found FA reductions in the corticospinal tract, temporal lobe, and occipital lobe (p < 0.05) relative to controls at DTI-1 and DTI-2. In contrast to the putamen, all exploratory ROIs showed no relevant FA change between DTI-1 and DTI-2. These results suggest that recurrent seizures may lead to progressive microstructural putamen alterations.

Multiscale Adaptive Regression Models for Neuroimaging Data

Neuroimaging studies aim to analyze imaging data with complex spatial patterns in a large number of locations (called voxels) on a two-dimensional (2D) surface or in a 3D volume. Conventional analyses of imaging data include two sequential steps: spatially smoothing imaging data and then independently fitting a statistical model at each voxel. However, conventional analyses suffer from the same amount of smoothing throughout the whole image, the arbitrary choice of smoothing extent, and low statistical power in detecting spatial patterns. We propose a multiscale adaptive regression model (MARM) to integrate the propagation-separation (PS) approach (Polzehl and Spokoiny, 2000, 2006) with statistical modeling at each voxel for spatial and adaptive analysis of neuroimaging data from multiple subjects. MARM has three features: being spatial, being hierarchical, and being adaptive. We use a multiscale adaptive estimation and testing procedure (MAET) to utilize imaging observations from the neighboring voxels of the current voxel to adaptively calculate parameter estimates and test statistics. Theoretically, we establish consistency and asymptotic normality of the adaptive parameter estimates and the asymptotic distribution of the adaptive test statistics. Our simulation studies and real data analysis confirm that MARM significantly outperforms conventional analyses of imaging data.

Image analysis and statistical inference in neuroimaging with R

R is a language and environment for statistical computing and graphics. It can be considered an alternative implementation of the S language developed in the 1970s and 1980s for data analysis and graphics (Becker and Chambers, 1984; Becker et al., 1988). The R language is part of the GNU project and offers versions that compile and run on almost every major operating system currently available. We highlight several R packages built specifically for the analysis of neuroimaging data in the context of functional MRI, diffusion tensor imaging, and dynamic contrast-enhanced MRI. We review their methodology and give an overview of their capabilities for neuroimaging. In addition we summarize some of the current activities in the area of neuroimaging software development in R.

Smoothing that does not blur: effects of the anisotropic approach for evaluating diffusion tensor imaging data in the clinic

PURPOSE

To compare the effects of anisotropic and Gaussian smoothing on the outcomes of diffusion tensor imaging (DTI) voxel-based (VB) analyses in the clinic, in terms of signal-to-noise ratio (SNR) enhancement and directional information and boundary structures preservation.

MATERIALS AND METHODS

DTI data of 30 Alzheimer's disease (AD) patients and 30 matched control subjects were obtained at 3T. Fractional anisotropy (FA) maps with variable degrees and quality (Gaussian and anisotropic) of smoothing were created and compared with an unsmoothed dataset. The two smoothing approaches were evaluated in terms of SNR improvements, capability to separate differential effects between patients and controls by a standard VB analysis, and level of artifacts introduced by the preprocessing.

RESULTS

Gaussian smoothing regionally biased the FA values and introduced a high variability of results in clinical analysis, greatly dependent on the kernel size. On the contrary, anisotropic smoothing proved itself capable of enhancing the SNR of images and maintaining boundary structures, with only moderate dependence of results on smoothing parameters.

CONCLUSION

Our study suggests that anisotropic smoothing is more suitable in DTI studies; however, regardless of technique, a moderate level of smoothing seems to be preferable considering the artifacts introduced by this manipulation.

Diffusion tensor imaging in Alzheimer's disease and mild cognitive impairment

Structural magnetic resonance imaging (MRI) studies of Alzheimer’s disease and mild cognitive impairment (MCI) have focused on the hippocampus and entorhinal cortex; gray matter structures in the medial temporal lobe. Few studies have investigated the integrity of white matter in patients with AD or MCI. Diffusion tensor imaging (DTI) is a MRI technique that allows for the interrogation of the microstructural integrity of white matter. Based on increases in translational diffusion (mean diffusivity: MD) and decreases directional diffusion (fractional anisotropy: FA) damage to white matter can be assessed. Studies have identified regions of increased MD and decreased FA in patients with AD and MCI in all lobes of the brain, as well as medial temporal lobe structures including the hippocampus, entorhinal cortex and parahippocampal white matter. The pattern of white matter integrity disruption tends to follow an anterior to posterior gradient with greater damage noted in posterior regions in AD and MCI. Recent studies have exploited inter-voxel directional similarities to develop models of white matter pathways, and have used these models to assess the integrity of inter-cerebral connections. Particular focus has been applied to the parahippocampal white matter (including the perforant path) and the posterior cingulum. Although many studies have found DTI indicators of impaired white matter in AD and MCI, other studies have failed to detect any differences in MD or FA between the groups, demonstrating the need for large replicative studies. DTI is an evolving technique and advances in its application ought to provide new insights into AD and MCI.

Using Perturbation theory to reduce noise in diffusion tensor fields

We propose the use of Perturbation theory to reduce noise in Diffusion Tensor (DT) fields. Diffusion Tensor Imaging (DTI) encodes the diffusion of water molecules along different spatial directions in a positive definite, 3 x 3 symmetric tensor. Eigenvectors and eigenvalues of DTs allow the in vivo visualization and quantitative analysis of white matter fiber bundles across the brain. The validity and reliability of these analyses are limited, however, by the low spatial resolution and low Signal-to-Noise Ratio (SNR) in DTI datasets. Our procedures can be applied to improve the validity and reliability of these quantitative analyses by reducing noise in the tensor fields. We model a tensor field as a three-dimensional Markov Random Field and then compute the likelihood and the prior terms of this model using Perturbation theory. The prior term constrains the tensor field to be smooth, whereas the likelihood term constrains the smoothed tensor field to be similar to the original field. Thus, the proposed method generates a smoothed field that is close in structure to the original tensor field. We evaluate the performance of our method both visually and quantitatively using synthetic and real-world datasets. We quantitatively assess the performance of our method by computing the SNR for eigenvalues and the coherence measures for eigenvectors of DTs across tensor fields. In addition, we quantitatively compare the performance of our procedures with the performance of one method that uses a Riemannian distance to compute the similarity between two tensors, and with another method that reduces noise in tensor fields by anisotropically filtering the diffusion weighted images that are used to estimate diffusion tensors. These experiments demonstrate that our method significantly increases the coherence of the eigenvectors and the SNR of the eigenvalues, while simultaneously preserving the fine structure and boundaries between homogeneous regions, in the smoothed tensor field.

Voxelwise regularisation of high angular resolution diffusion imaging data

The problem of noise suppression in high angular resolution diffusion MRI data is approached through direct regularisation of the apparent diffusion coefficient profiles. The proposed algorithm is derived in a Bayesian framework in the style of the traditional techniques for image restoration using Markov random field models. In a novel departure from the classical approach, a Markov random field model is applied within each voxel across gradient directions, thus smoothing the image data without inducing additional spatial dependencies that would render region-of-interest statistical testing of diffusion characteristics invalid. The anisotropic smoothing algorithm exploits the heterogeneous distribution of gradient directions and their antipodal pairs on the sphere and, in application to both simulated and experimental high angular resolution imaging datasets, is demonstrated to be superior to the isotropic Markov random field variant and the maximum likelihood estimator.