A general strategy for anisotropic diffusion in MR image denoising and enhancement

Anisotropic diffusion (AD) has proven to be very effective in the denoising of magnetic resonance (MR) images. The result of AD filtering is highly dependent on several parameters, especially the conductance parameter. However, there is no automatic method to select the optimal parameter values. This paper presents a general strategy for AD filtering of MR images using an automatic parameter selection method. The basic idea is to estimate the parameters through an optimization step on a synthetic image model, which is different from traditional analytical methods. This approach can be easily applied to more sophisticated diffusion models for better denoising results. We conducted a systematic study of parameter selection for the AD filter, including the dynamic parameter decreasing rate, the parameter selection range for different noise levels and the influence of the image contrast on parameter selection. The proposed approach was validated using both simulated and real MR images. The model image generated using our approach was shown to be highly suitable for the purpose of parameter optimization. The results confirm that our method outperforms most state-of-the-art methods in both quantitative measurement and visual evaluation. By testing on real images with different noise levels, we demonstrated that our method is sufficiently general to be applied to a variety of MR images.     

A theoretical study of the effect of experimental noise on the measurement of anisotropy in diffusion imaging

Diffusion tensor imaging (DTI) is a modality known to be highly sensitive to the detrimental effects of experimental noise. Here, using Monte Carlo simulations, we compare and contrast how noise complicates the measurement of diffusion anisotropy in diffusion tensor and conventional diffusion-weighted imaging (DWI). As the signal-to-noise ratio (SNR) decreases below a value of approximately 20, the eigenvalues (λ_i) of the diffusion tensor D are found to diverge rapidly from their true values, with the result that the measured anisotropy can be significantly in error and isotropic structures falsely assigned a high level of anisotropy. The effect of noise on the rotationally variant indices, calculated from a conventional diffusion-weighted imaging experiment, is found to be much less insidious, because the apparent diffusion coefficients (ADCs) diverge only slowly as the signal-to-noise decreases. Thus, although rotationally variant indices almost always underestimate the true diffusion anisotropy, they show only a small susceptibility to experimental noise and hence, are preferred to their rotationally invariant counterparts when the signal-to-noise ratio is small.     

Evaluation of diffusional anisotropy and microscopic structure in skeletal muscles using magnetic resonance

The pulsed-gradient spin-echo (PGSE) nuclear magnetic resonance (NMR) method is used for detecting the diffusion of water molecules in biological tissues. Because tissues generally have diffusional anisotropy, their diffusion properties are denoted by a tensor. In this study, we evaluated the diffusional anisotropy and microscopic structure in atrophied skeletal muscles using the PGSE NMR method. The left sciatic nerve was severed in twelve 9-week-old rats. Neurotomy caused neurogenic muscular atrophy at the left gastrocnemius. At 2, 4 and 8 weeks after neurotomy, magnetic resonance signals were selectively acquired from a 2 x 2 x 2 mm(3) voxel, which was located on the left gastrocnemius. The diffusion tensor, the mean diffusivity (MD) and the fractional anisotropy (FA) were calculated from the signals. A theoretical model of the diffusion in muscles was derived from Tanner's equation. The muscle fiber diameter was estimated by fitting the model to the measured signals. The measurements were also performed for normal rats as controls. No significant difference was found in the MD and the estimated intracellular diffusion coefficient between the control group and the denervated group. The denervated group had significantly higher FA compared with the control group (P<.05). The estimated muscle fiber diameter of the denervated group was significantly smaller than the estimated value of the control group (P<.05). These differences were found at 8 weeks after neurotomy. The proposed method is effective for evaluating changes in the microscopic structure of skeletal muscles.     

Electromagnetic wave propagation in polycrystalline materials: effective medium approach

A theory for the propagation of electromagnetic waves in inhomogeneous solids made of anisotropic crystallites is developed in the framework of the effective medium approach. The macroscopic dielectric tensor can explicitly be expressed by characteristic integrals containing the radial distribution function and a single anisotropy parameter. The phase mismatch of the waves scattered from misoriented crystallites leads to absorption and refraction effects that are calculated using a self-consistent approach in the sense of the distorted-wave approximation. For higher frequencies resonance structures occur which can be interpreted as an interference effect between disturbed and undisturbed waves in the effective medium.     

Simultaneous magnetic resonance imaging of diffusion anisotropy and diffusion gradient

The theory of diffusion gradient-weighted MRI (DGWI) is presented in this paper. The Bloch-Torrey equation was modified to include the effect of intravoxel spatial-location variation of water diffusion (diffusion gradient) on MRI signal, in addition to the effect of intravoxel spatial-direction variation of water diffusion (diffusion anisotropy). An analytical solution for a diffusion-encoding spin-echo pulse sequence was derived. Unlike water diffusion which attenuates the image signal intensity, this newly derived solution relates the spatial gradient of the water diffusion with the phase of the image signal. This novel MRI technique directly measures both the water diffusion and its spatial gradient, and thus offers a noninvasive imaging tool to simultaneously investigate the intravoxel inhomogeneity and anisotropy of tissue structures. In addition, as demonstrated with our preliminary data, this new method may be utilized to delineate the interfaces of tissues with different diffusion. This method is an extension of the successful diffusion tensor MRI (DTI), but requires no additional data acquisition. In addition to the measured diffusion tensor, this new method provides measurements of the spatial derivatives of the three principal diffusivities of the tensor, thereby providing additional information for improving white matter fiber tractography.     

Rapid microscopic fractional anisotropy imaging via an optimized linear regression formulation

Water diffusion anisotropy in the human brain is affected by disease, trauma, and development. Microscopic fractional anisotropy (μFA) is a diffusion MRI (dMRI) metric that can quantify water diffusion anisotropy independent of neuron fiber orientation dispersion. However, there are several different techniques to estimate μFA and few have demonstrated full brain imaging capabilities within clinically viable scan times and resolutions. Here, we present an optimized spherical tensor encoding (STE) technique to acquire μFA directly from the 2nd order cumulant expansion of the powder averaged dMRI signal obtained from direct linear regression (i.e. diffusion kurtosis) which requires fewer powder-averaged signals than other STE fitting techniques and can be rapidly computed. We found that the optimal dMRI parameters for white matter μFA imaging were a maximum b-value of 2000 s/mm² and a ratio of STE to LTE tensor encoded acquisitions of 1.7 for our system specifications. We then compared two implementations of the direct regression approach to the well-established gamma model in 4 healthy volunteers on a 3 Tesla system. One implementation used mean diffusivity (D) obtained from a 2nd order fit of the cumulant expansion, while the other used a linear estimation of D from the low b-values. Both implementations of the direct regression approach showed strong linear correlations with the gamma model (ρ = 0.97 and ρ = 0.90) but mean biases of −0.11 and − 0.02 relative to the gamma model were also observed, respectively. All three μFA measurements showed good test-retest reliability (ρ ≥ 0.79 and bias = 0). To demonstrate the potential scan time advantage of the direct approach, 2 mm isotropic resolution μFA was demonstrated over a 10 cm slab using a subsampled data set with fewer powder-averaged signals that would correspond to a 3.3-min scan. Accordingly, our results introduce an optimization procedure that has enabled nearly full brain μFA in only several minutes.     

Covariant stringlike dynamics of scroll wave filaments in anisotropic cardiac tissue

It has been hypothesized that stationary scroll wave filaments in cardiac tissue describe a geodesic in a curved space whose metric is the inverse diffusion tensor. Several numerical studies support this hypothesis, but no analytical proof has been provided yet for general anisotropy. In this Letter, we derive dynamic equations for the filament in the case of general anisotropy. These equations are covariant under general spatial coordinate transformations and describe the motion of a stringlike object in a curved space whose metric tensor is the inverse diffusion tensor. Therefore the behavior of scroll wave filaments in excitable media with anisotropy is similar to the one of cosmic strings in a curved universe. Our dynamic equations are valid for thin filaments and for general anisotropy. We show that stationary filaments obey the geodesic equation.     

Water diffusion features as indicators of muscle structure ex vivo

The structures of two different bovine muscles, the Semitendinosus (ST) and the Triceps brachii (TB), were studied using quantitative maps obtained by diffusion tensor imaging at 4.7 T. The estimated features were: mean diffusivity, intra- and inter-voxel anisotropy and fiber tract orientation angles. Significant differences in anisotropy (fractional anisotropy and lattice index), spatial variations of anisotropy and fiber tract orientation were detected between ST and TB, and are discussed. Accumulation of free water, which diffuses more freely and isotropically than in the rest of the muscle, was detected and localized in ST. These results underline the usefulness of diffusion tensor measurements to characterize muscle structure and help understand the mechanisms of post mortem water exudation.     

Electromagnetic wave propagation in polycrystalline materials: effective medium approach

Abstract

A theory for the propagation of electromagnetic waves in inhomogeneous solids made of anisotropic crystallites is developed in the framework of the effective medium approach. The macroscopic dielectric tensor can explicitly be expressed by characteristic integrals containing the radial distribution function and a single anisotropy parameter. The phase mismatch of the waves scattered from misoriented crystallites leads to absorption and refraction effects that are calculated using a self-consistent approach in the sense of the distorted-wave approximation. For higher frequencies resonance structures occur which can be interpreted as an interference effect between disturbed and undisturbed waves in the effective medium.     

How does distortion correction correlate with anisotropic indices? A diffusion tensor imaging study

PURPOSE: The purpose of this study was to determine a suitable registration algorithm for diffusion tensor imaging (DTI) using conventional preprocessing tools [statistical parametric mapping (SPM) and automated image registration (AIR)] and to investigate how anisotropic indices for clinical assessments are affected by these distortion corrections. MATERIALS AND METHODS: Brain DTI data from 15 normal healthy volunteers were used to evaluate four spatial registration schemes within subjects to correct image distortions: noncorrection, SPM-based affine registration, AIR-based affine registration and AIR-based nonlinear polynomial warping. The performance of each distortion correction was assessed using: (a) quantitative parameters: tensor-fitting error (Ef), mean dispersion index (MDI), mean fractional anisotropy (MFA) and mean variance (MV) within 11 regions of interest (ROI) defined from homogeneous fiber bundles; and (b) fiber tractography through the uncinate fasciculus and the corpus callosum. Fractional anisotropy (FA) and mean diffusivity (MD) were calculated to demonstrate the effects of distortion correction. Repeated-measures analysis of variance was used to investigate differences among the four registration paradigms. RESULTS: AIR-based nonlinear registration showed the best performance for reducing image distortions with respect to smaller Ef (P<.02), MDI (P<.01) and MV (P<.01) with larger MFA (P<.01). FA was decreased to correct distortions (P<.0001) whether the applied registration was linear or nonlinear and was lowest after nonlinear correction (P<.001). No significant differences were found in MD. CONCLUSION: In conventional DTI processing, anisotropic indices of FA can be misestimated by noncorrection or inappropriate distortion correction, which leads to an erroneous increase in FA. AIR-based nonlinear distortion correction would be required for a more accurate measurement of this diffusion parameter.     

Effect of B-value in revealing postinfarct myocardial microstructural remodeling using MR diffusion tensor imaging

Nonmonoexponential diffusion behavior has been previously reported to exist in some biological tissues, making quantification of diffusion tensor imaging (DTI) indices dependent on diffusion sensitivity of b-value. This study aims to investigate the effect of b-value in revealing postinfarct myocardial microstructural remodeling in ex vivo hearts. DTI scans were performed on heart samples 1, 3, 5, and 7 days after infarction induction as well as intact controls with b-values of 500 to 2500 s/mm². DTI indices, including fractional anisotropy (FA), and mean and directional diffusivities, were measured in infarct, adjacent and remote regions with zero and each non-zero b-values respectively using conventional DTI analysis. Experimental results showed that these DTI indices decreased gradually with b-values in all regions and groups. Optimal b-values were found to vary with targeted DTI indices, and could strengthen DTI ability in revealing myocardium degradation with using conventional DTI approach. Specifically, FA showed the most sensitive detection of fiber integrity degradation at moderate b-values (≈ 1500 to 2000 s/mm²), and the greatest ability of mean and directional diffusivities in monitoring diffusivity alteration occurred at relatively small b-values (≤ 1500 s/mm²) during the necrotic and fibrotic phases. These findings may provide useful information for DTI protocol parameter optimization in assessing heart microstructures at other pathological or in vivo states in the future.     

描述人体内水分子扩散各向异性特征的新方法

通过核磁共振扩散张量成像(DTI)得到的特定值域的扩散各向异性指数(DAI) 可用于揭示水分子扩散椭球的形态学特征, 定量反映被成像物体内部水分子扩散的优势方向和强度, 间接得到被成像物体内部的组织结构信息. DAI的可靠性直接影响对DTI数据的分析和理解. 本文基于扩散张量椭球的几何学信息, 提出利用扩散椭球几何比(EGR)定量描述水分子扩散的各向异性程度. 通过蒙特卡罗模拟实验和对人脑DTI数据进行分析, 并与当前广泛应用的水分子扩散各向异性分数(FA)和近期文献提出的扩散椭球面积比(EAR)进行对比. 实验发现EGR在不同级别噪声影响下的对比度效果和抗噪性都优于FA及EAR. 而且EGR 加入了体积修正, 增强了盘形扩散张量情况下的敏感性, 能够更好地鉴别神经纤维束交叉情况, 对于各向异性扩散程度较高的白质深层和相对均质的表层都有较好的量化区分结果. 关键词： 扩散系数 各向异性扩散 扩散张量成像 扩散椭球几何比     

In vivo diffusion tensor imaging of rat spinal cord at 7 T

In vivo diffusion tensor imaging of normal rat spinal cord was performed using a multi-segmented, blipped EPI sequence at 7 T field strength. At high diffusion weighting, the signal exhibited a non-monoexponential decay that was fitted to a biexponential function, associated with the fast and slow components of diffusion in the cord tissue, using a nonlinear regression analysis along with a constrained optimization procedure. From the measured tensors, the eigenvalues and the maps of invariant scalar measures (fractional anisotropy, relative anisotropy, volume ratio, and trace) were calculated and analyzed statistically. The results were combined to quantitatively characterize the anisotropic properties of the fast and slow diffusions in white- and gray matter of live spinal cords.     

Is diffusion anisotropy an accurate monitor of myelination?: Correlation of multicomponent T2 relaxation and diffusion tensor anisotropy in human brain

We compare T2-relaxation and diffusion tensor data from normal human brain. The relationships between myelin-water fraction (MWF) and various diffusion tensor measures [e.g., fractional anisotropy (FA), perpendicular diffusivity (ADC perpendicular) and mean diffusivity ] in white matter (WM) and gray matter (GM) structures in the brain were examined in 16 normal volunteers at 1.5 T and 6 normal subjects at 3.0 T and mean diffusivity. We found some degree of linear correlation between these measurements, but by using region of interest (ROI)-based analysis, we also observed several structures which seemed to deviate significantly from a linear relationship. From all investigated relationships between various diffusion tensor measures and myelin-water content, FA and ADC perpendicular yielded the highest correlation coefficients with MWF. However, diffusion anisotropy was also significantly influenced by factors other than myelin-water content. The less operator-dependent voxel-based analysis (VBA) between myelin-water and diffusional anisotropy measures is proposed as an innovative alternative to ROI-based analysis. We confirmed that WM structures, in general, have higher diffusional anisotropy than GM structures and also have higher myelin-water content. However, our findings suggest that in the highly organized fibre arrangement of compact WM structures such as the genu of the corpus callosum, elevated degrees of diffusional anisotropies are measured, which do not necessarily correspond to an elevated myelin content but more likely reflect the highly organized directionality of fibre bundles in these areas (low microscopic and macroscopic tortuosity) as well as strongly restricted diffusion in the interstitial space between the myelinated axons. Conversely, in structures with disorganized fibre bundles and multiple fibre crossings, such as the minor and major forceps, low FA values were measured, which does not necessarily reflect a decrease myelin-water content.     

Microstructural and physiological features of tissues elucidated by quantitative-diffusion-tensor MRI. 1996

Quantitative-diffusion-tensor MRI consists of deriving and displaying parameters that resemble histological or physiological stains, i.e., that characterize intrinsic features of tissue microstructure and microdynamics. Specifically, these parameters are objective, and insensitive to the choice of laboratory coordinate system. Here, these two properties are used to derive intravoxel measures of diffusion isotropy and the degree of diffusion anisotropy, as well as intervoxel measures of structural similarity, and fiber-tract organization from the effective diffusion tensor, D, which is estimated in each voxel. First, D is decomposed into its isotropic and anisotropic parts, 〈D〉 I and D – 〈D〉 I, respectively (where 〈D〉 = Trace(D)/3 is the mean diffusivity, and I is the identity tensor). Then, the tensor (dot) product operator is used to generate a family of new rotationally and translationally invariant quantities. Finally, maps of these quantitative parameters are produced from high-resolution diffusion tensor images (in which D is estimated in each voxel from a series of 2D-FT spin-echo diffusion-weighted images) in living cat brain. Due to the high inherent sensitivity of these parameters to changes in tissue architecture (i.e., macromolecular, cellular, tissue, and organ structure) and in its physiologic state, their potential applications include monitoring structural changes in development, aging, and disease.     

Influence of steady background gradients on the accuracy of molecular diffusion anisotropy measurements

Spatial susceptibility variations of body components lead to local gradients of the static magnetic field. Effects of such background gradients on fractional diffusion anisotropy (FA) measurements on whole-body magnetic resonance units operating at 1.5, 3.0 and 7.0 T were analyzed theoretically and experimentally. Analytical expressions were derived for the cases of diffusion occurring in isotropic media and in tissues with cylindrical symmetry (e.g., white matter tracts or skeletal musculature). Typical magnitudes of background gradient strengths were estimated from in vivo and in vitro measurements with B0 field mapping sequences. Additionally, numerical simulations of magnetic field distributions and resulting field gradients were performed considering tissue-air interfaces in simplified geometrical arrangements. For media with isotropic diffusion, both measurements and analytical calculations showed increasing FA inaccuracy with stronger coupling between diffusion-encoding and background gradients. For cylindrical symmetry, FA values were estimated for a standard diffusion tensor imaging protocol in a realistic scenario. At 1 mm distance from a water-air interface, susceptibility-related background gradients amount to approximately 9 mT/m at 7 T and lead to a relative error of the measured FA of up to 48%. The error in the anisotropy assessment rises considerably with increasing field strength and must be taken into account especially for experimental and clinical studies on modern high-field systems.     

A geometric theory for scroll wave filaments in anisotropic excitable media

Scroll waves are an important example of self-organisation in excitable media. In cardiac tissue, scroll waves of electrical activity underlie lethal ventricular arrhythmias and fibrillation. They rotate around a topological line defect which has been termed the filament. Numerical investigation has shown that anisotropy can substantially affect the dynamics of scroll waves. It has recently been hypothesised that stationary scroll wave filaments in cardiac tissue describe geodesics in a space whose metric is the inverse diffusion tensor. Several computational studies have validated this hypothesis, but until now no quantitative theory has been provided to study the effects of anisotropy on scroll wave filaments. Here, we review in detail the recently developed covariant formalism for scroll wave dynamics in general anisotropy and derive the equations of motion of filaments. These equations are fully covariant under general spatial coordinate transformations and describe the motion of filaments in a curved space whose metric tensor is the inverse diffusion tensor. Our dynamic equations are valid for thin filaments and for general anisotropy and we show that stationary filaments obey the geodesic equation. We extend previous work by allowing spatial variations in the determinant of the diffusion tensor and the reaction parameters, leading to drift of the filament.     

Determination of the rotational diffusion tensor of macromolecules in solution from nmr relaxation data with a combination of exact and approximate methods--application to the determination of interdomain orientation in multidomain proteins

In this paper we present a method for determining the rotational diffusion tensor from NMR relaxation data using a combination of approximate and exact methods. The approximate method, which is computationally less intensive, computes values of the principal components of the diffusion tensor and estimates the Euler angles, which relate the principal axis frame of the diffusion tensor to the molecular frame. The approximate values of the principal components are then used as starting points for an exact calculation by a downhill simplex search for the principal components of the tensor over a grid of the space of Euler angles relating the diffusion tensor frame to the molecular frame. The search space of Euler angles is restricted using the tensor orientations calculated using the approximate method. The utility of this approach is demonstrated using both simulated and experimental relaxation data. A quality factor that determines the extent of the agreement between the measured and predicted relaxation data is provided. This approach is then used to estimate the relative orientation of SH3 and SH2 domains in the SH(32) dual-domain construct of Abelson kinase complexed with a consolidated ligand.     

A diffusion gradient optimization framework for spinal cord diffusion tensor imaging

The uncertainty in the estimation of diffusion model parameters in diffusion tensor imaging (DTI) can be reduced by optimally selecting the diffusion gradient directions utilizing some prior structural information. This is beneficial for spinal cord DTI, where the magnetic resonance images have low signal-to-noise ratio and thus high uncertainty in diffusion model parameter estimation. Presented is a gradient optimization scheme based on D-optimality, which reduces the overall estimation uncertainty by minimizing the Rician Cramer-Rao lower bound of the variance of the model parameter estimates. The tensor-based diffusion model for DTI is simplified to a four-parameter axisymmetric DTI model where diffusion transverse to the principal eigenvector of the tensor is assumed isotropic. Through simulations and experimental validation, we demonstrate that an optimized gradient scheme based on D-optimality is able to reduce the overall uncertainty in the estimation of diffusion model parameters for the cervical spinal cord and brain stem white matter tracts.