Next Article in Journal
Positive-Operator Valued Measure (POVM) Quantization
Next Article in Special Issue
Heat Kernel Embeddings, Differential Geometry and Graph Structure
Previous Article in Journal
The Yang-Baxter Equation, (Quantum) Computers and Unifying Theories
Previous Article in Special Issue
The Gromov–Wasserstein Distance: A Brief Overview
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Simplified Algorithm for Inverting Higher Order Diffusion Tensors

1
Department of Biomedical Engineering, Eindhoven University of Technology, P.O. Box 513, Eindhoven NL-5600 MB, The Netherlands
2
Department of Mathematics and Computer Science, Eindhoven University of Technology, P.O. Box 513, Eindhoven NL-5600 MB, The Netherlands
*
Author to whom correspondence should be addressed.
Axioms 2014, 3(4), 369-379; https://doi.org/10.3390/axioms3040369
Submission received: 20 April 2014 / Revised: 5 November 2014 / Accepted: 7 November 2014 / Published: 14 November 2014

Abstract

:
In Riemannian geometry, a distance function is determined by an inner product on the tangent space. In Riemann–Finsler geometry, this distance function can be determined by a norm. This gives more freedom on the form of the so-called indicatrix or the set of unit vectors. This has some interesting applications, e.g., in medical image analysis, especially in diffusion weighted imaging (DWI). An important application of DWI is in the inference of the local architecture of the tissue, typically consisting of thin elongated structures, such as axons or muscle fibers, by measuring the constrained diffusion of water within the tissue. From high angular resolution diffusion imaging (HARDI) data, one can estimate the diffusion orientation distribution function (dODF), which indicates the relative diffusivity in all directions and can be represented by a spherical polynomial. We express this dODF as an equivalent spherical monomial (higher order tensor) to directly generalize the (second order) diffusion tensor approach. To enable efficient computation of Riemann–Finslerian quantities on diffusion weighted (DW)-images, such as the metric/norm tensor, we present a simple and efficient algorithm to invert even order spherical monomials, which extends the familiar inversion of diffusion tensors, i.e., symmetric matrices.

1. Introduction

Diffusion tensor imaging is a technique based on nuclear magnetic resonance imaging, where one can non-invasively measure the orientation-dependent diffusivity of water molecules in human tissue, such as brain white matter [1,2]. As the behavior of water molecules inside tissue is influenced by the local geometry, this allows the indirect study of the structure and orientation of tissue. One of the goals in brain diffusion imaging is to use this information to trace the axonal fiber bundles that carry information between neurons in different parts of the brain. There are numerous models that relate the diffusion weighted imaging (DWI) signal to local physiological parameters [3]. One of the earlier models, and essentially the only one that is currently used in the clinical setting, is diffusion tensor imaging (DTI). In this model, one assumes anisotropic Gaussian diffusion.
By imaging the signal attenuation, first without an additional magnetic gradient and then comparing this with the signal when a directed magnetic field is applied, one can measure the so-called apparent diffusion coefficient (ADC) in each voxel. Measurements are repeated in several, e.g., 35 directions, and a second order model is fitted to describe this radially-symmetric diffusion profile. The diffusion length, typically measured at a scale of square micrometers per milliseconds, in direction v = ( cos ϕ sin θ , sin ϕ sin θ , cos θ ) , can be computed as:
D ( v ) = i , j = 1 3 D i j v i v j
where D i j is the element at the i-th row and j-th column of a 3 × 3 matrix D (diffusion tensor) and v i is the i-th component of vector v .
Throughout this paper, the following notation is used: we express scalars with regular lower case fonts s (except for the diffusivity D), vectors with bold face lower cases v and tensors with bold face upper case fonts T . Similar to the vectors, the components of the tensors are denoted with indices in the subscript T i j .
In the DTI case, two popular ways to visualize the diffusion profile is to either plot the set of points { v · D 1 2 · v · v | | v | = 1 } , i.e., a glyph or { D · v | | v | = 1 } , that is, an ellipsoid (see Figure 1). Note that in the D-induced metric, the vectors that constitute the ellipsoid have a unit length. The semi-axes of these ellipsoids are λ i v i , where λ i are the eigenvalues and v i the eigenvectors of the diffusion tensor.
Later, O’Donnell et al. connected the inverse of the diffusion tensor to a metric tensor [4,5]. Whereas the diffusion tensor tells how far a particle can diffuse in a given time, say a millisecond, a metric tensor tells how long the particle has to diffuse to reach, say, a micrometer distance from the starting point. Together with the assumption that fiber bundles coincide with the shortest curves relative to this metric, this interpretation of the diffusion tensor led to a fiber tracking algorithm, called geodesic tractography. Of course, given the complex architecture of the white matter, such an assumption is not always valid. Nevertheless, whenever axon tracking can be cast into a shortest path problem, efficient methods for extracting the geodesics become relevant and are valuable as such.
Figure 1. (A) The ellipsoid, i.e., the unit level set of the inverse of a diffusion tensor; and (B) the glyph determined by the same diffusion tensor.
Figure 1. (A) The ellipsoid, i.e., the unit level set of the inverse of a diffusion tensor; and (B) the glyph determined by the same diffusion tensor.
Axioms 03 00369 g001
High resolution diffusion magnetic resonance imaging (HARDI) brought axon tracking to another level by scanning far more imaging directions than is necessary for DTI. The diffusion profiles (glyphs) could now have a more complex form than a symmetric peanut. The most obvious disadvantage of the DTI model is the fact that it works very well when a voxel only contains a single fiber bundle, but fails completely when there are two or more fiber populations with distinct orientations. A number of models have been suggested to overcome this, but none of them has the intuitive Riemannian interpretation as in DTI. In an effort to extend the Riemannian framework into HARDI, approaches employing Riemann–Finsler geometry were suggested [6,7,8,9,10]. The essential difference with respect to Riemannian geometry is that the distance function is determined by a position- and direction-dependent norm function instead of a position-dependent inner product. Axons that intersect each other in a single voxel, which were impossible to distinguish in the DTI setting, could now be resolved to some extent using the Riemann–Finsler geometry [6]. An algorithm extending the regular streamline tracking to the Finslerian framework, where the metric is determined by the Hessian of the local Finsler-norm, could apparently successfully extract crossing bundles of fibers in human brain, such as typical of corpus callosum, corona radiata and cingulum [8]. However, the approach in [8], while being a straightforward higher order generalization of DTI, has a shortcoming that we wish to remove.
The Finsler-norm function was computed as the (multiplicative) Möbius-inverse [11] of the diffusion profile, which is typically not a spherical tensor (monomial) itself.
We would like to treat diffusion and norm profiles on an equal footing, as is the case for the diffusion and metric tensor in the Riemannian framework [4,12]. This requires that the (higher order) diffusion and metric tensors can be related to each other via a unique inverse operation.
In this contribution, we describe a simple heuristic algorithm for computing the inverse of a higher order tensor. In Section 2, we briefly provide the motivation for our approach, although, in essence, equivalent to that of other authors [13,14], besides aesthetic, but also practical and efficient. In Section 3, we describe our simple inversion algorithm in detail, as well as by means of visual illustration and show some examples. We finish with conclusions in Section 5.

2. The Parsimony of Monomials

By a (spherical) monomial, we mean a function of the following type:
i 1 , , i n = 1 3 A i 1 i n v i 1 v i n
where A i n R , | v | = 1 and n is a fixed even positive number. This in contrast with a (spherical) polynomial:
n = 0 , 2 , i 1 , , i n = 1 3 A i 1 i 2 i n v i 1 v i n
In both cases, running through all unit vectors v on the two-sphere, one obtains a spherical surface homeomorphic to the two-sphere.
The general protocol of converting raw HARDI data to a diffusion profile (a spherical polynomial or an equivalent monomial) is described elsewhere [7,13,14], and we do not repeat it here. The number of directions along which the ADC can be measured and the resulting angular resolution is, in practice, limited. This implies that it makes no sense to use a very high order model to capture the diffusion profile [15]. Therefore, a spherical polynomial of order from 2 to 6 is typically sufficient to estimate the diffusion profiles. The conventional formula that translates the signal S D obtained with a diffusion-sensitized gradient to an “apparent diffusion coefficient” D is as follows:
D = ln ( S 0 / S D ) b
Here, S 0 is the signal without diffusion weighting magnetic gradient and:
b = ( γ G Δ ) 2 ( Δ - δ / 3 )
where δ is the gradient pulse duration, Δ is the time between gradient pulses, γ is the gyromagnetic ratio and G is the magnitude of the magnetic gradient. In DWI, one cannot distinguish between the two possible directions (forward/backward) along a specified orientation. For this reason, at least in the conventional HARDI, the diffusion profile is always radially symmetric. As a consequence, any function describing such a profile is also radially symmetric.
For the reasons mentioned above, we model diffusion profiles with fourth order monomials fitted to the signal data as in [7,14]. It can be shown that a monomial expression is equivalent to spherical harmonics or a polynomial expression of equal order [7,16,17]. Therefore, why use monomials? One reason is that they are a compact expression that directly generalizes the DTI tensor concept. Another reason is that in exploiting the Riemann–Finsler concepts of norms and metric tensors, one computes from a position ( x ) -dependent norm a position and direction (x, y) -dependent metric tensor by calculating the Hessian of the squared norm function with respect to y as follows:
G i j ( x , y ) = 1 2 2 F 2 ( x ) y i y j
where G is the metric tensor and F is the norm function. Indeed, monomials are extremely simple/fast to differentiate twice. The third reason is that, as we will show in the next section, an even order symmetric monomial can be inverted, so that the inverse is also (an equal order) monomial, and moreover, this inversion can be done in a simple setting of matrix multiplication.
The main motivation for computing these position- and orientation-dependent glyphs and (Finsler) metric tensors is that this gives a straightforward extension of DTI streamline tracking to HARDI, allowing the extraction of multiple fiber orientations in a single voxel.

3. Inversion of Even Order Fully-Symmetric Tensors

Brazell et al. [18] show how even order tensors can be inverted in the sense of the so-called Einstein contracted product. In our context, the tensors are totally symmetric, which allows some essential simplifications to the more general formulation of inversion in [18].
To keep the notation simple, we define the Einstein contracted product (ECP) only for fourth order tensors. For a general definition, see [18]. Let A and B be fourth order tensors. The ECP of these, C , is a fourth order tensor with elements:
C i j k l = m , n = 1 3 A i j m n · B m n k l
Recall that even order tensors can be regarded as nested matrices. In particular, a fourth order tensor can be denoted by a matrix of matrices. For example, a fourth order tensor A can be explicitly written down as follows:
A = A 1111 A 1112 A 1113 A 1121 A 1122 A 1123 A 1131 A 1132 A 1133 A 1211 A 1212 A 1213 A 1221 A 1222 A 1223 A 1231 A 1232 A 1233 A 1311 A 1312 A 1313 A 1321 A 1322 A 1323 A 1331 A 1332 A 1333 A 2111 A 2112 A 2113 A 2121 A 2122 A 2123 A 2131 A 2132 A 2133 A 2211 A 2212 A 2213 A 2221 A 2222 A 2223 A 2231 A 2232 A 2233 A 2311 A 2312 A 2313 A 2321 A 2322 A 2323 A 2331 A 2332 A 2333 A 3111 A 3112 A 3113 A 3121 A 3122 A 3123 A 3131 A 3132 A 3133 A 3211 A 3212 A 3213 A 3221 A 3222 A 3223 A 3231 A 3232 A 3233 A 3311 A 3312 A 3313 A 3321 A 3322 A 3323 A 3331 A 3332 A 3333
Further, recalling that in our case, the tensors are totally symmetric, we have, in particular, the following equivalence:
A i j k l = A k l i j
This means that the ECP of A and B can be seen as componentwise products of all “inner” matrices. Furthermore, if we “reshape” the inner matrices, ECP is actually a matrix product (see Figure 2).
Figure 2. (A) An illustration of how the nested matrix should be reshaped to replace the Einstein contracted product (ECP) with an ordinary matrix product; (B) the correspondence between ECP and the matrix dot product after manipulating the tensor shape.
Figure 2. (A) An illustration of how the nested matrix should be reshaped to replace the Einstein contracted product (ECP) with an ordinary matrix product; (B) the correspondence between ECP and the matrix dot product after manipulating the tensor shape.
Axioms 03 00369 g002
To find the inverse tensor A - 1 of A , we need to solve this from the following equation:
s y m m , n = 1 3 A i j m n · A m n k l - 1 = I i j k l
where the left-hand side is the symmetrized ECP of A and A - 1 , and I is the fourth order identity tensor. The fully symmetric identity tensor I is defined as follows:
I i j k l = 1 3 δ i k δ j l + δ i l δ j k + δ i j δ k l
where:
δ i j = 1 , i = j 0 , i j
The symmetric identity tensor I in matrix notation is then:
1 0 0 0 1 3 0 0 0 1 3 0 1 3 0 1 3 0 0 0 0 0 0 0 1 3 0 0 0 1 3 0 0 0 1 3 0 1 3 0 0 0 0 0 1 3 0 0 0 1 0 0 0 1 3 0 0 0 0 0 1 3 0 1 3 0 0 0 1 3 1 3 0 0 1 3 0 0 0 0 0 0 0 1 3 0 1 3 0 1 3 0 0 0 1 3 0 0 0 1
This inversion clearly generalizes the matrix inversion used to convert the familiar second order diffusion tensor to a metric tensor. Due to the symmetries between the tensor elements, the inverse operation can be performed in the three simple steps of Algorithm 1.
Algorithm 1: Computing the inverse of a symmetric fourth order tensor A .
1.
Compute the Einstein product P of the chosen tensor A with a symmetric variable tensor V ,  with components v 1111 , v 1112 , , e.g., as a matrix product, as in Figure 2;
2.
reshape the result to a four tensor and symmetrize with formula P i j k l s y m = 1 4 ! σ S 4 P σ ( i ) σ ( j ) I σ ( k ) σ ( l ) ;
3.
Solve the 15 equations P 1111 s y m = I 1111 , P 1112 s y m = I 1112 , , P 3333 s y m = I 3333 .
Here, S 4 refers to the set of all permutations of a set of four elements. For a visually intuitive explanation on symmetrizing a fourth order tensor, see e.g., [7]. One should not directly solve the tensor with the matrix inverse, since the highly symmetric tensor in matrix form is singular. By solving only the 15 necessary equations, we avoid a system of over-determined equations and get reliable results, which are different from that obtained by the plain pseudo-inverse approach.
It is interesting to see that this also perfectly matches with the intuition of inverting the so-called glyph, which is the set:
i , j , k , l = 1 3 A i j k l y i y k y l | | y | = 1
The glyph (and the unit level set) of I is the unit two-sphere. For an illustration of dODF-based glyphs and their ECP-inverses, see Figure 3. Clearly, in the orientations that correspond to larger diffusivity, the ECP-inverse is small, just as the case with inverting the second order tensor. Using this inverted tensor, one can further compute Finslerian quantities, such as geodesics, that can be exploited, e.g., in the context of diffusion-based fiber tracking.
Figure 3. (Left) The glyph represents the diffusion orientation distribution function. (middle) The glyph represents the ECP-inverse of this; and (right) another tensor glyph with its inverse is plotted together. For visual comparison, lines of longitudes are indicated in red and latitudes in blue.
Figure 3. (Left) The glyph represents the diffusion orientation distribution function. (middle) The glyph represents the ECP-inverse of this; and (right) another tensor glyph with its inverse is plotted together. For visual comparison, lines of longitudes are indicated in red and latitudes in blue.
Axioms 03 00369 g003

4. Applications in Fiber Tracking

Once a proper diffusion-induced metric has been defined for a DWI image, one may apply tractography in order to reconstruct the fibrous tissue structure of the underlying neural axons of the white matter. One prominent approach in deterministic tractography is to model the fibers as minimum length paths, i.e., the geodesics in the given metric space. The main assumption in this technique is that the fibers tend to follow the most efficient diffusion propagation, thereby corresponding to the paths of shortest distance [19,20]. A recently-introduced geodesic tractography technique obtains the pathways as solutions of Euler–Lagrange (EL) equations in Riemannian or Finslerian space [21,22]. This tractography approach can capture multi-valued geodesics connecting two given points or regions by considering the geodesics as a function of position and direction. A manifold is defined using as a metric the adjugate of the diffusion tensor or the inverted dODF. Please note that we use the dODF obtained by fitting the HARDI data with a procedure described by Descoteaux et al. [23]. By applying the inversion technique proposed here to the dODF-tensors and subsequently feeding them as input for the tractography, we have obtained some preliminary results. The tractography was done in a 10 × 10 × 30 cube cut from the healthy human dataset obtained by a Philips scanner with resolution 2 × 2 × 2 mm, b-value of 1000 s/mm 2 and 128 gradient directions.
Given the region of interests, geodesics are computed from all discrete coordinates to a large number of directions until they meet one of the boundaries. To determine which of the numerous geodesics corresponds to the most likely fiber, we have used the line-plane intersection method, as in [22,24]. In this way, the geodesics are cut once they enter one of the user-defined regions of interest. Figure 4 indicates the ODFs and the inverted ODFs for a part of the centrum semiovale. In Figure 5, few fiber tracts for corpus callosum (CC) have been computed, and their corresponding ODFs and the inverted ODFs per fiber voxel have been indicated. Figure 6 shows the tractography result for the area of corona radiata (CR), CC and cingulum (CG).
Figure 4. Inverted fourth order orientation distribution function (ODF) field (top) and its corresponding ODF field (bottom) from part of the centrum semiovale.
Figure 4. Inverted fourth order orientation distribution function (ODF) field (top) and its corresponding ODF field (bottom) from part of the centrum semiovale.
Axioms 03 00369 g004
Figure 5. Corpus callosum fibers: (left) the corresponding ODFs per fiber voxel and (right) the corresponding inverted ODFs per fiber voxel.
Figure 5. Corpus callosum fibers: (left) the corresponding ODFs per fiber voxel and (right) the corresponding inverted ODFs per fiber voxel.
Axioms 03 00369 g005
Figure 6. Fiber tracts for corpus callosum (CC), corona radiata (CR) and cingulum (CG).
Figure 6. Fiber tracts for corpus callosum (CC), corona radiata (CR) and cingulum (CG).
Axioms 03 00369 g006

5. Conclusions

We have proposed to use a symmetrized contracted Einstein product in the inversion of higher even order tensors in the context of HRDW-image analysis. We have shown that this approach generalizes the approach of O’Donnell et al., which applies tools from Riemannian geometry in the DTI setting. The unit level set of the ECP-inverted monomial generalizes the metric ellipsoid used in the DTI-setting. By employing monomials in modeling the diffusion profiles, we obtain an intuitive and consistent framework for the analysis and tractography on DW-images. We also illustrated how this inversion can be efficiently computed using a simple Algorithm 1. This is useful, since, in the case a tensor is of high order, the computation of indexed sums can easily become time consuming.

Acknowledgments

Tom Dela Haije gratefully acknowledges the Netherlands Organization for Scientific Research for financial support.

Author Contributions

Laura Astola conceived the study and prepared the manuscript. Neda Sepasian implemented and performed the fiber tractography. Tom Dela Haije, Andrea Fuster, Luc Florack and Neda Sepasian edited the manuscript, discussed the results and implications and commented on the manuscript at all stages.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Stejskal, E.; Tanner, J. Spin Diffusion Measurements: Spin Echoes in the Presence of a Time-Dependent Field Gradient. J. Chem. Phys. 1965, 42. [Google Scholar] [CrossRef]
  2. Mori, S. Introduction to Diffusion Tensor Imaging; Elsevier: Amsterdam, The Netherlands, 2007. [Google Scholar]
  3. Johansen-Berg, H.; Behrens, T. Diffusion MRI, 2nd ed.; Academic Press: London, UK, 2013. [Google Scholar]
  4. O’Donnell, L.; Haker, S.; Westin, C.F. New Approaches to Estimation of White Matter Connectivity in Diffusion Tensor MRI: Elliptic PDEs and Geodesics in a Tensor-Warped Space. Lect. Notes Comput. Sci. 2002, 2488, 459–466. [Google Scholar]
  5. Lenglet, C.; Deriche, R.; Faugeras, O. Inferring White Matter Geometry from Diffusion Tensor MRI: Application to Connectivity Mapping. In Proceedings of the Eighth European Conference on Computer Vision, Prague, Czech Republic, 11–14 May 2004.
  6. Melonakos, J.; Pichon, E.; Angenent, S.; Tannenbaum, A. Finsler Active Contours. IEEE Trans. Pattern Anal. Mach. Intell. 2008, 30, 412–423. [Google Scholar] [CrossRef] [PubMed]
  7. Astola, L. Multi-Scale Riemann-Finsler Geometry, Applications to Diffusion Tensor Imaging and High Angular Resolution Diffusion Imaging. Ph.D. Thesis, Technische Universiteit Eindhoven, Eindhoven, The Netherlands, January 2010. [Google Scholar]
  8. Astola, L.; Jalba, A.; Balmashnova, E.; Florack, L. Finsler Streamline Tracking with Single Tensor Orientation Distribution Function for High Angular Resolution Diffusion Imaging. J. Math. Imaging Vis. 2011, 41, 170–181. [Google Scholar] [CrossRef]
  9. Florack, L.; Fuster, A.; Dela Haije, T. Riemann-Finsler Geometry and its Applications to Diffusion Magnetic Resonance Imaging. In Proceedings of the 2013 International BASP Frontiers Workshop, Villars-sur-Ollon, Switzerland, 27 January–1 February 2013.
  10. Florack, L.; Fuster, A. Visualization and Processing of Tensors and Higher Order Descriptors for Multi-Valued Data; Mathematics and Visualization; Westin, C.F., Vilanova, A., Burgeth, B., Eds.; Springer: Berlin, Germany, 2014. [Google Scholar]
  11. Vuorinen, M. Conformal Geometry and Quasiregular Mappings; Springer-Verlag: New York, NY, USA, 1980. [Google Scholar]
  12. Astola, L.; Florack, L.; ter Haar Romeny, B. Measures for Pathway Analysis in Brain White Matter using Diffusion Tensor Images. Inf. Process. Med. Imaging 2007, 20, 642–649. [Google Scholar] [PubMed]
  13. Tuch, D.; Reese, T.; Wiegell, M.; Makris, N.; Belliveau, J.; van Wedeen, J. High angular resolution diffusion imaging reveals intravoxel white matter fiber heterogeneity. Magn. Reson. Med. 2004, 48, 1358–1372. [Google Scholar] [CrossRef] [PubMed]
  14. Özarslan, E.; Mareci, T. Generalized Diffusion tensor Imaging and Analytical Relationships between Diffusion Tensor Imaging and High Angular Resolution Diffusion Imaging. Magn. Reson. Med. 2003, 50, 955–965. [Google Scholar] [CrossRef] [PubMed]
  15. Chung, H.W.; Chou, M.C.; Chen, C.Y. Principles and Limitations of Computational Algorithms in Clinical Diffusion Tensor MR Tractography. MAJNR Am. J. Neuroradiol. 2011, 32, 3–13. [Google Scholar] [CrossRef] [PubMed]
  16. Müller, C. Analysis of Spherical Symmetries in Euclidean Spaces; Springer: New York, NY, USA, 1998; Volume 129. [Google Scholar]
  17. Balmashnova, E.; Fuster, A.; Florack, L. Decomposition of Higher-Order Homogeneous Tensors and Applications to HARDI. In Proceedings of the Computational Diffusion MRI 2012 (CDMRI’12), Nice, France, 5 October 2012.
  18. Brazell, M.; Li, N.; Navasca, C.; Tamon, C. Solving Multilinear Systems via Tensor Inversion. SIAM J. Matrix Anal. Appl. 2013, 34, 542–570. [Google Scholar] [CrossRef]
  19. Hao, X.; Zygmunt, K.; Whitaker, R.; Fletcher, P.T. Improved Segmentation of White Matter Tracts with Adaptive Riemannian Metrics. Med. Image Anal. 2014, 18, 161–175. [Google Scholar] [CrossRef] [PubMed]
  20. Hao, X.; Whitaker, R.; Fletcher, P.T. Adaptive Riemannian metrics for improved geodesic tracking of white matter. In Information Processing in Medical Imaging; Springer: Berlin, Germany, 2011. [Google Scholar]
  21. Sepasian, N.; ten Thije Boonkkamp, J.H.M.; Florack, L.; Romeny, B.M.T.H.; Vilanova, A. Visualization and Processing of Tensors and Higher Order Descriptors for Multi-Valued Data; Westin, C.F., Vilanova, A., Burgeth, B., Eds.; Springer: Berlin, Germany, 2014. [Google Scholar]
  22. Sepasian, N.; ten Thije Boonkkamp, J.H.M.; Romeny, B.M.T.H.; Vilanova, A. Multivalued Geodesic Ray-Tracing for Computing Brain Connections Using Diffusion Tensor Imaging. SIAM J. Imaging Sci. 2012, 5, 483–504. [Google Scholar] [CrossRef]
  23. Descoteaux, M.; Angelino, E.; Fitzgibbons, S.; Deriche, R. Regularized, Fast and Robust Analytical Q-Ball Imaging. Magn. Reson. Med. 2007, 58, 497–510. [Google Scholar] [CrossRef] [PubMed]
  24. Van Aart, E.; Sepasian, N.; Jalba, A.; Vilanova, A. CUDA-Accelerated Geodesic Ray-Tracing for Fiber Tracking. Int. J. Biomed. Imaging 2011, 2011. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Astola, L.; Sepasian, N.; Haije, T.D.; Fuster, A.; Florack, L. A Simplified Algorithm for Inverting Higher Order Diffusion Tensors. Axioms 2014, 3, 369-379. https://doi.org/10.3390/axioms3040369

AMA Style

Astola L, Sepasian N, Haije TD, Fuster A, Florack L. A Simplified Algorithm for Inverting Higher Order Diffusion Tensors. Axioms. 2014; 3(4):369-379. https://doi.org/10.3390/axioms3040369

Chicago/Turabian Style

Astola, Laura, Neda Sepasian, Tom Dela Haije, Andrea Fuster, and Luc Florack. 2014. "A Simplified Algorithm for Inverting Higher Order Diffusion Tensors" Axioms 3, no. 4: 369-379. https://doi.org/10.3390/axioms3040369

APA Style

Astola, L., Sepasian, N., Haije, T. D., Fuster, A., & Florack, L. (2014). A Simplified Algorithm for Inverting Higher Order Diffusion Tensors. Axioms, 3(4), 369-379. https://doi.org/10.3390/axioms3040369

Article Metrics

Back to TopTop