Formal characterization and extension of the linearized diffusion tensor model

Abstract

We analyzed the properties of the logarithm of the Rician distribution leading to a full characterization of the probability law of the errors in the linearized diffusion tensor model. An almost complete lack of bias, a simple relation between the variance and the signal-to-noise ratio in the original complex data, and a close approximation to normality facilitated estimation of the tensor components by an iterative weighted least squares algorithm. The theory of the linear model has also been used to derive the distribution of mean diffusivity, to develop an informative statistical test for relative lack of fit of the ellipsoidal (or spherical) model compared to an unrestricted linear model in which no specific shape is assumed for the diffusion process, and to estimate the signal-to-noise ratios in the original data. The false discovery rate (FDR) has been used to control thresholds for statistical significance in the context of multiple comparisons at voxel level. The methods are illustrated by application to three diffusion tensor imaging (DTI) datasets of clinical interest: a healthy volunteer, a patient with acute brain injury, and a patient with hydrocephalus. Interestingly, some salient features, such as a region normally comprising the basal ganglia and internal capsule, and areas of edema in patients with brain injury and hydrocephalus, had patterns of error largely independent from their mean diffusivities. These observations were made in brain regions with sufficiently large signal-to-noise ratios (>2) to justify the assumptions of the log Rician probability model. The combination of diffusivity and its error may provide added value in diagnostic DTI of acute pathologic expansion of the extracellular fluid compartment in brain parenchymal tissue.

Diffusion tensor imaging in schizophrenia: Relationship to symptoms

“…Susceptibility-induced geometric distortion was corrected using a 1-D correction algorithm (motion correction was also employed as part of this routine) along the phase-encoding direction with the separatelyacquired field map (Jezzard and Balaban, 1995). The software used a weighted linear least square method to calculate FA maps (Salvador et al, 2005). Intracranial voxels were separated out using a brain extraction tool-generated brain mask.…”

Probabilistic white matter fiber tracking using particle filtering and von Mises–Fisher sampling

“…It is well known that the noise in MRI can be described accurately by a Rician distribution Salvador et al (2005). Salvador et al (2005) showed that the error distribution conforms closely to a normal distribution with a zero-mean and standard deviation equal to the inverse of the signal-to-noise ratio (SNR), i.e.…”

“…Here, we simply set tr(D) to be the trace of the diffusion tensor D estimated using the linear least squares estimation (Basser et al, 1994), and set λ ⊥ = (λ 2 + λ 3 )/2. The SNR is estimated using the weighted least squares method in (Salvador et al, 2005). Panels (a) and (c) of Fig.…”

Diffusion Tensor Estimation by Maximizing Rician Likelihood

“…Jones et al [12] presented an estimation method that incorporates noise level estimation. Salvador et al [18] reviews distribution assumptions and describes a weighted least squares procedure for addressing non-Gaussianity. Recent abstracts indicate that a Rician noise model is more accurate than Gaussian estimation (e.g., [1]).…”