The ellipsoidal area ratio: an alternative anisotropy index for diffusion tensor imaging

In the processing and analysis of diffusion tensor imaging (DTI) data, certain predefined morphological features of diffusion tensors are often represented as simplified scalar indices, termed diffusion anisotropy indices (DAIs). When comparing tensor morphologies across differing voxels of an image, or across corresponding voxels in different images, DAIs are mathematically and statistically more tractable than are the full tensors, which are probabilistic ellipsoids consisting of three orthogonal vectors that each has a direction and an associated scalar magnitude. We have developed a new DAI, the "ellipsoidal area ratio" (EAR), to represent the degree of anisotropy in the morphological features of a diffusion tensor. The EAR is a normalized geometrical measure of surface curvature in the 3D diffusion ellipsoid. Monte Carlo simulations and applications to the study of in vivo human data demonstrate that, at low noise levels, EAR provides a similar contrast-to-noise ratio (CNR) but a higher signal-to-noise ratio (SNR) than does fractional anisotropy (FA), which is currently the most popular anisotropy index in active use. Moreover, at the high noise levels encountered most commonly in real-world DTI datasets, EAR compared with FA is consistently much more robust to perturbations from noise and it provides a higher CNR, features useful for the analysis of DTI data that are inherently noise sensitive.     

A unifying theoretical and algorithmic framework for least squares methods of estimation in diffusion tensor imaging

A unifying theoretical and algorithmic framework for diffusion tensor estimation is presented. Theoretical connections among the least squares (LS) methods, (linear least squares (LLS), weighted linear least squares (WLLS), nonlinear least squares (NLS) and their constrained counterparts), are established through their respective objective functions, and higher order derivatives of these objective functions, i.e., Hessian matrices. These theoretical connections provide new insights in designing efficient algorithms for NLS and constrained NLS (CNLS) estimation. Here, we propose novel algorithms of full Newton-type for the NLS and CNLS estimations, which are evaluated with Monte Carlo simulations and compared with the commonly used Levenberg-Marquardt method. The proposed methods have a lower percent of relative error in estimating the trace and lower reduced chi2 value than those of the Levenberg-Marquardt method. These results also demonstrate that the accuracy of an estimate, particularly in a nonlinear estimation problem, is greatly affected by the Hessian matrix. In other words, the accuracy of a nonlinear estimation is algorithm-dependent. Further, this study shows that the noise variance in diffusion weighted signals is orientation dependent when signal-to-noise ratio (SNR) is low (相似文献   

Diffusion tensor imaging at low SNR: nonmonotonic behaviors of tensor contrasts

Diffusion tensor imaging (DTI) provides measurements of directional diffusivities and has been widely used to characterize changes in the tissue microarchitecture of the brain. DTI is gaining prominence in applications outside of the brain, where resolution, motion and short T2 values often limit the achievable signal-to-noise ratio (SNR). Consequently, it is important to revisit the topic of tensor estimation in low-SNR regimes. A theoretical framework is developed to model noise in DTI, and by using simulations based on this theory, the degree to which the noise, tensor estimation method and acquisition protocol affect tensor-derived quantities, such as fractional anisotropy and apparent diffusion coefficient, is clarified. These results are then validated against clinical data. It is shown that reliability of tensor contrasts depends on the noise level, estimation method, diffusion-weighting scheme and underlying anatomy. The propensity for bias and errors does not monotonically increase with noise. Comparative results are shown in both graphical and tabular forms, so that decisions about suitable acquisition protocols and processing methods can be made on a case-by-case basis without exhaustive experimentation.     

A statistical framework for the classification of tensor morphologies in diffusion tensor images

Tractography algorithms for diffusion tensor (DT) images consecutively connect directions of maximal diffusion across neighboring DTs in order to reconstruct the 3-dimensional trajectories of white matter tracts in vivo in the human brain. The performance of these algorithms, however, is strongly influenced by the amount of noise in the images and by the presence of degenerate tensors-- i.e., tensors in which the direction of maximal diffusion is poorly defined. We propose a simple procedure for the classification of tensor morphologies that uses test statistics based on invariant measures of DTs, such as fractional anisotropy, while accounting for the effects of noise on tensor estimates. Examining DT images from seven human subjects, we demonstrate that this procedure validly classifies DTs at each voxel into standard types (nondegenerate DTs, as well as degenerate oblate, prolate or isotropic DTs), and we provide preliminary estimates for the prevalence and spatial distribution of degenerate tensors in these brains. We also show that the P values for test statistics are more sensitive tools for classifying tensor morphologies than are invariant measures of anisotropy alone.     

Condition number as a measure of noise performance of diffusion tensor data acquisition schemes with MRI

Diffusion tensor mapping with MRI can noninvasively track neural connectivity and has great potential for neural scientific research and clinical applications. For each diffusion tensor imaging (DTI) data acquisition scheme, the diffusion tensor is related to the measured apparent diffusion coefficients (ADC) by a transformation matrix. With theoretical analysis we demonstrate that the noise performance of a DTI scheme is dependent on the condition number of the transformation matrix. To test the theoretical framework, we compared the noise performances of different DTI schemes using Monte-Carlo computer simulations and experimental DTI measurements. Both the simulation and the experimental results confirmed that the noise performances of different DTI schemes are significantly correlated with the condition number of the associated transformation matrices. We therefore applied numerical algorithms to optimize a DTI scheme by minimizing the condition number, hence improving the robustness to experimental noise. In the determination of anisotropic diffusion tensors with different orientations, MRI data acquisitions using a single optimum b value based on the mean diffusivity can produce ADC maps with regional differences in noise level. This will give rise to rotational variances of eigenvalues and anisotropy when diffusion tensor mapping is performed using a DTI scheme with a limited number of diffusion-weighting gradient directions. To reduce this type of artifact, a DTI scheme with not only a small condition number but also a large number of evenly distributed diffusion-weighting gradients in 3D is preferable.     

Diffusion tensor encoding schemes optimized for white matter fibers with selected orientations

A method to produce gradient encoding schemes that minimize the noise of diffusion tensor imaging (DTI) indices for selected fiber orientations has been developed. The accuracy of DTI measurements depends on the gradient encoding scheme used. Most current acquisition schemes contain diffusion directions uniformly distributed in 3D space in order to provide equal noise levels for fibers in any orientation. However, when considering specific fiber bundles such as the corticospinal tract (CST) or parts of fiber bundles, the range of fiber orientations of interest may be limited. We hypothesized that, when studying fiber tracts with a limited range of orientations, measuring diffusion in directions that are uniformly distributed in 3D space may be suboptimal for the noise levels of various DTI indices. Therefore, we first used simulations to determine six diffusion directions that minimize the noise of DTI measurements for selected fiber orientations. The resulting optimized set of directions was then tested on the right CST of a healthy human subject, and its performance was compared with that of conventional acquisition strategies. Both the simulations and the experiments on the human subject demonstrated that the new scheme significantly reduced the standard deviation of DTI indices for tensors with primary eigenvectors within a selected range of orientations.     

Error bounds in diffusion tensor estimation using multiple-coil acquisition systems

We extend the diffusion tensor (DT) signal model for multiple-coil acquisition systems. Considering the sum-of-squares reconstruction method, we compute the Cramér–Rao bound (CRB) assuming the widely accepted noncentral chi distribution. Within this framework, we assess the effect of noise in DT estimation and other measures derived from it, as a function of the number of acquisition coils, as well as other system parameters. We show the applications of CRB in many actual problems related to DT estimation: we compare different gradient field setup schemes proposed in the literature and show how the CRB can be used to choose a convenient one; we show that for fiber-type anisotropy tensors the ellipsoidal area ratio (EAR) can be estimated with less error than other scalar factors such as the fractional anisotropy (FA) or the relative anisotropy (RA), and that for this type of anisotropy tensors, increasing the number of coils is equivalent to increasing the signal-to-noise ratio, i.e., the information of the different coils can be regarded as independent. Also, we present results showing the CRB of several parameters for actual DT-MRI data. We conclude that the CRB is a valuable tool to optimal experiment design in DT-related studies.     

A continuous tensor field approximation of discrete DT-MRI data for extracting microstructural and architectural features of tissue

The effective diffusion tensor of water, D, measured by diffusion tensor MRI (DT-MRI), is inherently a discrete, noisy, voxel-averaged sample of an underlying macroscopic effective diffusion tensor field, D(x). Within fibrous tissues this field is presumed to be continuous and smooth at a gross anatomical length scale. Here a new, general mathematical framework is proposed that uses measured DT-MRI data to produce a continuous approximation to D(x). One essential finding is that the continuous tensor field representation can be constructed by repeatedly performing one-dimensional B-spline transforms of the DT-MRI data. The fidelity and noise-immunity of this approximation are tested using a set of synthetically generated tensor fields to which background noise is added via Monte Carlo methods. Generally, these tensor field templates are reproduced faithfully except at boundaries where diffusion properties change discontinuously or where the tensor field is not microscopically homogeneous. Away from such regions, the tensor field approximation does not introduce bias in useful DT-MRI parameters, such as Trace(D(x)). It also facilitates the calculation of several new parameters, particularly differential quantities obtained from the tensor of spatial gradients of D(x). As an example, we show that they can identify tissue boundaries across which diffusion properties change rapidly using in vivo human brain data. One important application of this methodology is to improve the reliability and robustness of DT-MRI fiber tractography.     

Noise and nonlinear estimation with optimal schemes in DTI

In general, the estimation of the diffusion properties for diffusion tensor experiments (DTI) is accomplished via least squares estimation (LSE). The technique requires applying the logarithm to the measurements, which causes bad propagation of errors. Moreover, the way noise is injected to the equations invalidates the least squares estimate as the best linear unbiased estimate. Nonlinear estimation (NE), despite its longer computation time, does not possess any of these problems. However, all of the conditions and optimization methods developed in the past are based on the coefficient matrix obtained in a LSE setup. In this article, NE for DTI is analyzed to demonstrate that any result obtained relatively easily in a linear algebra setup about the coefficient matrix can be applied to the more complicated NE framework. The data, obtained using non-optimal and optimized diffusion gradient schemes, are processed with NE. In comparison with LSE, the results show significant improvements, especially for the optimization criterion. However, NE does not resolve the existing conflicts and ambiguities displayed with LSE methods.     

A new gravitational energy tensor

The Bel-Robinson tensor is the most used gravitational energy tensor; however, it has the dimensions of energy squared. How to construct tensors with the dimensions of energy by using Lancoz tensors is shown here. The resulting tensors have a large number of arbitrary parameters, frequently have spacelike currents, and frequently do not reduce to familiar pseudo-energy tensors in the weak field limit. Two particular examples of interest are one with well-behaved currents and one which reduces to an energy pseudo-tensor in the weak field limit.     

Diffusion tensor image up-sampling: a registration-based approach

Diffusion weighted images (DWI), from which the corresponding diffusion tensor images (DTI) are estimated, are commonly acquired with anisotropic discretizations. Traditional methods to up-sample diffusion weighted images generally rely on scene-based interpolation and do not exploit structural information from the images. In this study, a DTI up-sampling framework is presented that incorporates the underlying anatomical shape information by means of non-rigid inter-slice registration. A strategy is proposed to reorient the interpolated tensor in order to maintain its proper orientation. Tests on phantom as well as on real data sets show that the proposed method is able to produce better results compared to scene based interpolation methods in terms of the accuracy of DWI/DTI interpolation, especially when diffusion tensor orientation is taken into account.     

Retrospective correction of bias in diffusion tensor imaging arising from coil combination mode

PurposeTo quantify and retrospectively correct for systematic differences in diffusion tensor imaging (DTI) measurements due to differences in coil combination mode.BackgroundMulti-channel coils are now standard among MRI systems. There are several options for combining signal from multiple coils during image reconstruction, including sum-of-squares (SOS) and adaptive combine (AC). This contribution examines the bias between SOS- and AC-derived measures of tissue microstructure and a strategy for limiting that bias.MethodsFive healthy subjects were scanned under an institutional review board-approved protocol. Each set of raw image data was reconstructed twice—once with SOS and once with AC. The diffusion tensor was calculated from SOS- and AC-derived data by two algorithms—standard log-linear least squares and an approach that accounts for the impact of coil combination on signal statistics. Systematic differences between SOS and AC in terms of tissue microstructure (axial diffusivity, radial diffusivity, mean diffusivity and fractional anisotropy) were evaluated on a voxel-by-voxel basis.ResultsSOS-based tissue microstructure values are systematically lower than AC-based measures throughout the brain in each subject when using the standard tensor calculation method. The difference between SOS and AC can be virtually eliminated by taking into account the signal statistics associated with coil combination.ConclusionsThe impact of coil combination mode on diffusion tensor-based measures of tissue microstructure is statistically significant but can be corrected retrospectively. The ability to do so is expected to facilitate pooling of data among imaging protocols.     

Translation invariant tensor product states in a finite lattice system

We show that the matrix (or more generally tensor) product states in a finite translation invariant system can be accurately constructed from a same set of local matrices (or tensors) that are determined from an infinite lattice system in one or higher dimensions. This provides an efficient approach for studying translation invariant tensor product states in finite lattice systems. Two methods are introduced to determine the size-independent local tensors.     

Efficient anisotropic filtering of diffusion tensor images

To improve the accuracy of structural and architectural characterization of living tissue with diffusion tensor imaging, an efficient smoothing algorithm is presented for reducing noise in diffusion tensor images. The algorithm is based on anisotropic diffusion filtering, which allows both image detail preservation and noise reduction. However, traditional numerical schemes for anisotropic filtering have the drawback of inefficiency and inaccuracy due to their poor stability and first order time accuracy. To address this, an unconditionally stable and second order time accuracy semi-implicit Craig-Sneyd scheme is adapted in our anisotropic filtering. By using large step size, unconditional stability allows this scheme to take much fewer iterations and thus less computation time than the explicit scheme to achieve a certain degree of smoothing. Second-order time accuracy makes the algorithm reduce noise more effectively than a first order scheme with the same total iteration time. Both the efficiency and effectiveness are quantitatively evaluated based on synthetic and in vivo human brain diffusion tensor images, and these tests demonstrate that our algorithm is an efficient and effective tool for denoising diffusion tensor images.     

基于快速层次交替最小二乘非负张量Tucker分解的干涉高光谱图像光谱信息压缩方法

提出一种基于快速层次交替最小二乘非负张量Tucker分解的高光谱图像光谱信息压缩算法。首先,将干涉高光谱图像光程差方向的三维信息采用三维光程差方向提升小波变换(3D OPT-LDWT)进行分解,将三维小波子带系数看作三阶非负张量,采用快速层次交替最小二乘非负张量Tucker分解(FHALS-NTD)算法对进行分解,得到核心张量和模式矩阵。对每个模式矩阵进行量化,对核心张量采用比特平面重要系数编码算法进行编码,得出最终的压缩码流。结果表明,此压缩算法可以稳定可靠地工作。与传统压缩算法比较,平均信噪比提高了1.23 dB。有效的提高了干涉高光谱图像压缩性能。     

Optimal real-time estimation in diffusion tensor imaging

Diffusion tensor imaging (DTI) constitutes the most used paradigm among the diffusion-weighted magnetic resonance imaging (DW-MRI) techniques due to its simplicity and application potential. Recently, real-time estimation in DW-MRI has deserved special attention, with several proposals aiming at the estimation of meaningful diffusion parameters during the repetition time of the acquisition sequence. Specifically focusing on DTI, the underlying model of the noise present in the acquired data is not taken into account, leading to a suboptimal estimation of the diffusion tensor. In this paper, we propose an optimal real-time estimation framework for DTI reconstruction in single-coil acquisitions. By including an online estimation of the time-changing noise variance associated to the acquisition process, the proposed method achieves the sequential best linear unbiased estimator. Results on both synthetic and real data show that our method outperforms those so far proposed, reaching the best performance of the existing proposals by processing a substantially lower number of diffusion images.     

A Monte Carlo simulation of image misalignment effects in diffusion tensor imaging

Diffusion tensor imaging (DTI) and tractography are noninvasive MRI methods, providing an insight on microscopic structural information of anisotropic tissues in vivo. The success of this technique stems on a watchful choice of imaging parameters and post-acquisition reconstruction. In the present work, we have focused on the problem of residual linear image misalignment in the DTI data and its effects on the parameters of the diffusion tensor and fiber tracking in human brain. We demonstrate substantial sensitivity of the reconstructed diffusion tensor and fiber tractography on increasing amplitude of artificially induced random image misalignment in the DTI. We show that already a submillimeter image misalignment in the DTI is an important source of error, which may potentially mask pathological presentations of the diseases and may partially explain variations in the results obtained from the DTI. Finally, we evaluated four implementations of image registrations and demonstrate their variable performance. This further supports the fact that a robust image registration must be performed to ensure reliable and reproducible diffusion tensor mapping and reconstruction of white matter (WM) fibers.     

Lovelock tensor as generalized Einstein tensor

We show that the splitting feature of the Einstein tensor, as the first term of the Lovelock tensor, into two parts, namely the Ricci tensor and the term proportional to the curvature scalar, with the trace relation between them is a common feature of any other homogeneous terms in the Lovelock tensor. Motivated by the principle of general invariance, we find that this property can be generalized, with the aid of a generalized trace operator which we define, for any inhomogeneous Euler–Lagrange expression that can be spanned linearly in terms of homogeneous tensors. Then, through an application of this generalized trace operator, we demonstrate that the Lovelock tensor analogizes the mathematical form of the Einstein tensor, hence, it represents a generalized Einstein tensor. Finally, we apply this technique to the scalar Gauss–Bonnet gravity as an another version of string–inspired gravity. This work was partially supported by a grant from the MSRT/Iran.     

Estimation and application of spatially variable noise fields in diffusion tensor imaging

Optimal interpretation of magnetic resonance image content often requires an estimate of the underlying image noise, which is typically realized as a spatially invariant estimate of the noise distribution. This is not an ideal practice in diffusion tensor imaging because the noise distribution is usually spatially varying due to the use of fast imaging and noise suppression techniques. A new estimation approach for spatially varying noise fields (NFs) is proposed in this article. The approach is based on a noise invariance property in scenarios in which more than one image, each with potentially different signal levels, is acquired on each slice, as in diffusion-weighted MRI. This technique leads to improved NF estimates in simulations, phantom experiments and in vivo studies when compared to traditional NF estimators that use regional variability or background intensity histograms. The proposed method reduces the NF estimation error by a factor of 100 in simulations, shows a strong linear correlation (R²=0.99) between theoretical and estimated noise changes in phantoms and demonstrates consistent (<5% variability) NF estimates in vivo. The advantages of spatially varying NF estimation are demonstrated for power analysis, outlier detection and tensor estimation.