Feature-based interpolation of diffusion tensor fields and application to human cardiac DT-MRI

A new framework for analysis of three-dimensional shape and architecture of human skeletal muscles from in vivo imaging data

A new framework is presented for comprehensive analysis of the three-dimensional shape and architecture of human skeletal muscles from magnetic resonance and diffusion tensor imaging data. The framework comprises three key features: (1) identification of points on the surface of and inside a muscle that have a correspondence to points on and inside another muscle, (2) reconstruction of average muscle shape and average muscle fibre orientations, and (3) utilization of data on between-muscle variation to visualize and make statistical inferences about changes or differences in muscle shape and architecture. The general use of the framework is demonstrated by its application to three case studies. Analysis of data obtained before and after eight weeks of strength training revealed there was little regional variation in hypertrophy of the vastus medialis and vastus lateralis, and no systematic change in pennation angle. Analysis of passive muscle lengthening revealed heterogeneous changes in shape of the medial gastrocnemius, and confirmed the ability of the methods to detect subtle changes in muscle fibre orientation. Analysis of the medial gastrocnemius of children with unilateral cerebral palsy showed that muscles in the more-affected limb were shorter, thinner and less wide than muscles in the less-affected limb, and had slightly more pennate muscle fibres in the central and proximal part of the muscle. Amongst other applications, the framework can be used to explore the mechanics of muscle contraction, investigate adaptations of muscle architecture, build anatomically realistic computational models of skeletal muscles, and compare muscle shape and architecture between species.

Spectrum-sine interpolation framework for DTI processing

Diffusion tensor imaging (DTI) data interpolation is important for DTI processing, which could affect the precision and computational complexity in the process of denoising, filtering, regularization, and DTI registration and fiber tracking. In this paper, we propose a novel DTI interpolation framework named with spectrum-sine (SS) focusing on tensor orientation variation in DTI processing. Compared with the state-of-the-art DTI interpolation method using Euler angles or quaternion to represent the orientation of DTI tensors, this method does not need to convert eigenvectors into Euler angles or quaternions, but interpolates each tensor's unit eigenvector directly. The prominent merit of this tensor interpolation method lies in tensor orientation information preservation, which is different from the existing DTI tensor interpolation methods that interpolating tensor's orientation information in a scalar way. The experimental results from both synthetic and real human brain DTI data demonstrated the proposed SS interpolation scheme not only maintains the advantages of Log-Euclidean and Riemannian interpolation frameworks, such as preserving the tensor's symmetric positive definiteness and the monotonic determinant variation, but also preserve the tensor's anisotropy property which was proposed in the spectral quaternion (SQ) method.

Accelerating cardiac diffusion tensor imaging combining local low-rank and 3D TV constraint

OBJECTIVE

Diffusion tensor magnetic resonance imaging (DT-MRI, or DTI) is a promising technique for invasively probing biological tissue structures. However, DTI is known to suffer from much longer acquisition time with respect to conventional MRI and the problem is worsened when dealing with in vivo acquisitions. Therefore, faster DTI for both ex vivo and in vivo scans is highly desired.

MATERIALS AND METHODS

This paper proposes a new compressed sensing (CS) reconstruction method that employs local low-rank (LLR) model and three-dimensional (3D) total variation (TV) constraint to reconstruct cardiac diffusion-weighted (DW) images from highly undersampled k-space data. The LLR model takes the set of DW images corresponding to different diffusion gradient directions as a 3D image volume and decomposes the latter into overlapping 3D blocks. Then, the 3D blocks are stacked as two-dimensional (2D) matrix. Finally, low-rank property is applied to each block matrix and the 3D TV constraint to the 3D image volume. The underlying constrained optimization problem is finally solved using the first-order fast method. The proposed method is evaluated on real ex vivo cardiac DTI data as a prerequisite to in vivo cardiac DTI applications.

RESULTS

The results on real human ex vivo cardiac DTI images demonstrate that the proposed method exhibits lower reconstruction errors for DTI indices, including fractional anisotropy (FA), mean diffusivities (MD), transverse angle (TA), and helix angle (HA), compared to existing CS-based DTI image reconstruction techniques.

CONCLUSION

The proposed method provides better reconstruction quality and more accurate DTI indices in comparison with the state-of-the-art CS-based DW image reconstruction methods.

Multi-output Gaussian processes for enhancing resolution of diffusion tensor fields

Second order diffusion tensor (DT) fields are widely used in several clinical applications: brain fibers connections, diagnosis of neuro-degenerative diseases, image registration, brain conductivity models, etc. However, due to current acquisition protocols and hardware limitations in MRI machines, the diffusion magnetic resonance imaging (dMRI) data is obtained with low spatial resolution (1 or 2 mm3 for each voxel). This issue can be significant, because tissue fibers are much smaller than voxel size. Interpolation has become in a successful methodology for enhancing spatial resolution of DT fields. In this work, we present a feature-based interpolation approach through multi-output Gaussian processes (MOGP). First, we extract the logarithm of eigenvalues (direction) and the Euler angles (orientation) from diffusion tensors and we consider each feature as a separated but related output. Then, we interpolate the features along the whole DT field. In this case, the independent variables are the space coordinates (x, y, z). For this purpose, we assume that all features follow a multi-output Gaussian process with a common covariance matrix. Finally, we reconstruct new tensors from the interpolated eigenvalues and Euler angles. Accuracy of our methodology is better compared to approaches in the state of the art for performing DT interpolation, and it achieves a performance similar to the recently introduced method based on Generalized Wishart processes for interpolation of positive semidefinite matrices. We also show that MOGP preserves important properties of diffusion tensors such as fractional anisotropy.

On the averaging of cardiac diffusion tensor MRI data: the effect of distance function selection

Diffusion tensor magnetic resonance imaging (DT-MRI) allows a unique insight into the microstructure of highly-directional tissues. The selection of the most proper distance function for the space of diffusion tensors is crucial in enhancing the clinical application of this imaging modality. Both linear and nonlinear metrics have been proposed in the literature over the years. The debate on the most appropriate DT-MRI distance function is still ongoing. In this paper, we presented a framework to compare the Euclidean, affine-invariant Riemannian and log-Euclidean metrics using actual high-resolution DT-MRI rat heart data. We employed temporal averaging at the diffusion tensor level of three consecutive and identically-acquired DT-MRI datasets from each of five rat hearts as a means to rectify the background noise-induced loss of myocyte directional regularity. This procedure is applied here for the first time in the context of tensor distance function selection. When compared with previous studies that used a different concrete application to juxtapose the various DT-MRI distance functions, this work is unique in that it combined the following: (i) metrics were judged by quantitative-rather than qualitative-criteria, (ii) the comparison tools were non-biased, (iii) a longitudinal comparison operation was used on a same-voxel basis. The statistical analyses of the comparison showed that the three DT-MRI distance functions tend to provide equivalent results. Hence, we came to the conclusion that the tensor manifold for cardiac DT-MRI studies is a curved space of almost zero curvature. The signal to noise ratio dependence of the operations was investigated through simulations. Finally, the 'swelling effect' occurrence following Euclidean averaging was found to be too unimportant to be worth consideration.

Quantitative representation and description of intravoxel fiber complexity in HARDI

Diffusion tensor imaging and high angular resolution diffusion imaging are often used to analyze the fiber complexity of tissues. In these imaging techniques, the most commonly calculated metric is anisotropy, such as fractional anisotropy (FA), generalized anisotropy (GA), and generalized fractional anisotropy (GFA). The basic idea underlying these metrics is to compute the deviation from free or spherical diffusion. However, in many cases, the question is not really to know whether it concerns spherical diffusion. Instead, the main concern is to describe and quantify fiber complexity such as fiber crossing in a voxel. In this context, it would be more direct and effective to compute the deviation from a single fiber bundle instead of a sphere. We propose a new metric, called PEAM (PEAnut Metric), which is based on computing the deviation of orientation diffusion functions (ODFs) from a single fiber bundle ODF represented by a peanut. As an example, the proposed PEAM metric is used to classify intravoxel fiber configurations. The results on simulated data, physical phantom data and real brain data consistently showed that the proposed PEAM provides greater accuracy than FA, GA and GFA and enables parallel and complex fibers to be better distinguished.

Resolving intravoxel fiber architecture using nonconvex regularized blind compressed sensing

In diffusion magnetic resonance imaging, accurate and reliable estimation of intravoxel fiber architectures is a major prerequisite for tractography algorithms or any other derived statistical analysis. Several methods have been proposed that estimate intravoxel fiber architectures using low angular resolution acquisitions owing to their shorter acquisition time and relatively low b-values. But these methods are highly sensitive to noise. In this work, we propose a nonconvex regularized blind compressed sensing approach to estimate intravoxel fiber architectures in low angular resolution acquisitions. The method models diffusion-weighted (DW) signals as a sparse linear combination of unfixed reconstruction basis functions and introduces a nonconvex regularizer to enhance the noise immunity. We present a general solving framework to simultaneously estimate the sparse coefficients and the reconstruction basis. Experiments on synthetic, phantom, and real human brain DW images demonstrate the superiority of the proposed approach.

Free-breathing diffusion tensor imaging and tractography of the human heart in healthy volunteers using wavelet-based image fusion

Free-breathing cardiac diffusion tensor imaging (DTI) is a promising but challenging technique for the study of fiber structures of the human heart in vivo. This work proposes a clinically compatible and robust technique to provide three-dimensional (3-D) fiber architecture properties of the human heart. To this end, 10 short-axis slices were acquired across the entire heart using a multiple shifted trigger delay (TD) strategy under free breathing conditions. Interscan motion was first corrected automatically using a nonrigid registration method. Then, two post-processing schemes were optimized and compared: an algorithm based on principal component analysis (PCA) filtering and temporal maximum intensity projection (TMIP), and an algorithm that uses the wavelet-based image fusion (WIF) method. The two methods were applied to the registered diffusion-weighted (DW) images to cope with intrascan motion-induced signal loss. The tensor fields were finally calculated, from which fractional anisotropy (FA), mean diffusivity (MD), and 3-D fiber tracts were derived and compared. The results show that the comparison of the FA values (FA(PCATMIP) = 0.45 ±0.10, FA(WIF) = 0.42 ±0.05, P=0.06) showed no significant difference, while the MD values ( MD(PCATMIP)=0.83 ±0.12×10(-3) mm (2)/s, MD(WIF)=0.74±0.05×10(-3) mm (2)/s, P=0.028) were significantly different. Improved helix angle variations through the myocardium wall reflecting the rotation characteristic of cardiac fibers were observed with WIF. This study demonstrates that the combination of multiple shifted TD acquisitions and dedicated post-processing makes it feasible to retrieve in vivo cardiac tractographies from free-breathing DTI acquisitions. The substantial improvements were observed using the WIF method instead of the previously published PCATMIP technique.

Cardiac fiber unfolding by semidefinite programming

Diffusion-tensor imaging allows noninvasive assessment of the myocardial fiber architecture, which is fundamental in understanding the mechanics of the heart. In this context, tractography techniques are often used for representing and visualizing cardiac fibers, but their output is only qualitative. We introduce here a new framework toward a more quantitative description of the cardiac fiber architecture from tractography results. The proposed approach consists in taking three-dimensional (3-D) fiber tracts as inputs, and then unfolding these fibers in the Euclidean plane under local isometry constraints using semidefinite programming. The solution of the unfolding problem takes the form of a Gram matrix which defines the two-dimensional (2-D) embedding of the fibers and whose spectrum provides quantitative information on their organization. Experiments on synthetic and real data show that unfolding makes it easier to observe and to study the cardiac fiber architecture. Our conclusion is that 2-D embedding of cardiac fibers is a promising approach to supplement 3-D rendering for understanding the functioning of the heart.

A comparative study of different level interpolations for improving spatial resolution in diffusion tensor imaging

This paper studies and evaluates the feasibility and the performance of different level interpolations for improving spatial resolution of diffusion tensor magnetic resonance imaging (DT-MRI or DTI). In particular, the following techniques are investigated: anisotropic interpolation operating on scalar gray-level images, log-Euclidean interpolation method, and the quaternion interpolation method, which operate on diffusion tensor fields. The performance is evaluated both qualitatively and quantitatively using criteria such as tensor determinant, fractional anisotropy (FA), mean diffusivity (MD), fiber length, etc. We conclude that tensor field interpolations allow avoiding undesirable swelling effect in DTI, which is not the case with scalar gray-level interpolation, and that scalar gray-level image interpolation and log-Euclidean tensor field interpolation suffer from decrease in FA and MD, which may mislead the interpretation of the clinical parameters FA and MD. In contrast, the quaternion tensor field interpolation avoids such FA and MD decrease, which suggests its use for clinical applications.

Assessment of cardiac motion effects on the fiber architecture of the human heart in vivo

The use of diffusion tensor imaging (DTI) for studying the human heart in vivo is very challenging due to cardiac motion. This paper assesses the effects of cardiac motion on the human myocardial fiber architecture. To this end, a model for analyzing the effects of cardiac motion on signal intensity is presented. A Monte-Carlo simulation based on polarized light imaging data is then performed to calculate the diffusion signals obtained by the displacement of water molecules, which generate diffusion weighted (DW) images. Rician noise and in vivo motion data obtained from DENSE acquisition are added to the simulated cardiac DW images to produce motion-induced datasets. An algorithm based on principal components analysis filtering and temporal maximum intensity projection (PCATMIP) is used to compensate for motion-induced signal loss. Diffusion tensor parameters derived from motion-reduced DW images are compared to those derived from the original simulated DW images. Finally, to assess cardiac motion effects on in vivo fiber architecture, in vivo cardiac DTI data processed by PCATMIP are compared to those obtained from one trigger delay (TD) or one single phase acquisition. The results showed that cardiac motion produced overestimated fractional anisotropy and mean diffusivity as well as a narrower range of fiber angles. The combined use of shifted TD acquisitions and postprocessing based on image registration and PCATMIP effectively improved the quality of in vivo DW images and subsequently, the measurement accuracy of fiber architecture properties. This suggests new solutions to the problems associated with obtaining in vivo human myocardial fiber architecture properties in clinical conditions.

Generating fibre orientation maps in human heart models using Poisson interpolation

Smoothly varying muscle fibre orientations in the heart are critical to its electrical and mechanical function. From detailed histological studies and diffusion tensor imaging, we now know that fibre orientations in humans vary gradually from approximately -70° in the outer wall to +80° in the inner wall. However, the creation of fibre orientation maps for computational analyses remains one of the most challenging problems in cardiac electrophysiology and cardiac mechanics. Here, we show that Poisson interpolation generates smoothly varying vector fields that satisfy a set of user-defined constraints in arbitrary domains. Specifically, we enforce the Poisson interpolation in the weak sense using a standard linear finite element algorithm for scalar-valued second-order boundary value problems and introduce the feature to be interpolated as a global unknown. User-defined constraints are then simply enforced in the strong sense as Dirichlet boundary conditions. We demonstrate that the proposed concept is capable of generating smoothly varying fibre orientations, quickly, efficiently and robustly, both in a generic bi-ventricular model and in a patient-specific human heart. Sensitivity analyses demonstrate that the underlying algorithm is conceptually able to handle both arbitrarily and uniformly distributed user-defined constraints; however, the quality of the interpolation is best for uniformly distributed constraints. We anticipate our algorithm to be immediately transformative to experimental and clinical settings, in which it will allow us to quickly and reliably create smooth interpolations of arbitrary fields from data-sets, which are sparse but uniformly distributed.

Linear invariant tensor interpolation applied to cardiac diffusion tensor MRI

PURPOSE

Various methods exist for interpolating diffusion tensor fields, but none of them linearly interpolate tensor shape attributes. Linear interpolation is expected not to introduce spurious changes in tensor shape.

METHODS

Herein we define a new linear invariant (LI) tensor interpolation method that linearly interpolates components of tensor shape (tensor invariants) and recapitulates the interpolated tensor from the linearly interpolated tensor invariants and the eigenvectors of a linearly interpolated tensor. The LI tensor interpolation method is compared to the Euclidean (EU), affine-invariant Riemannian (AI), log-Euclidean (LE) and geodesic-loxodrome (GL) interpolation methods using both a synthetic tensor field and three experimentally measured cardiac DT-MRI datasets.

RESULTS

EU, AI, and LE introduce significant microstructural bias, which can be avoided through the use of GL or LI.

CONCLUSION

GL introduces the least microstructural bias, but LI tensor interpolation performs very similarly and at substantially reduced computational cost.