
Study with the several resources on Docsity
Earn points by helping other students or get them with a premium plan
Prepare for your exams
Study with the several resources on Docsity
Earn points to download
Earn points by helping other students or get them with a premium plan
Community
Ask the community for help and clear up your study doubts
Discover the best universities in your country according to Docsity users
Free resources
Download our free guides on studying techniques, anxiety management strategies, and thesis advice from Docsity tutors
The cone of uncertainty in diffusion tensor imaging (dti) tractography and how the uncertainty in fiber directions is not described by a right circular cone but a right elliptical cone. The authors present three methods to demonstrate this finding and discuss its implications for dti tractography. The document also mentions previous studies on the topic and provides references.
Typology: Lecture notes
1 / 1
This page cannot be seen from the preview
Don't miss anything!
C. Koay^1 , L-C. Chang^1 , and P. J. Basser 1 (^1) NICHD, NIH, Bethesda, MD, United States
The diffusion tensor is always measured in the presence of background noise. Within the framework of recently proposed methods of error propagation and of 1st^ -order matrix perturbation for DTI [1,2], one can easily characterize the noise induced variability of the estimated tensor itself and of tensor- derived quantities, scalar or vector. Of particular interest in tractography is the uncertainty of the eigenvectors or principal directions. With this end in view, the cone of uncertainty was proposed as a local measure of tract dispersion by Basser [3] and later used by Basser & Pajevic [4] and Jones [5] as a tool for visualizing and quantifying uncertainty along the fiber tract. The cone of uncertainty is not the only measure for quantifying tract dispersion, there is another approach that based on the perturbation method, which was proposed by Anderson [6] and used by Lazar [7,8]. None of these studies provided the key geometric object relating to the local measure of tract dispersion, which is the covariance matrix of the major eigenvector of the diffusion tensor. Here we present three different methods based on the error propagation framework [1], the 1st^ -order matrix perturbation method [2], and the dyadic formalism [3,9], to show that the local variability in the purported fiber direction is, in general, not described by a right circular cone of uncertainty but a right elliptical cone. Consequences of this finding are discussed. METHODS Error Propagation Method Based on the framework presented in [1], propagation of errors amounts to finding an appropriate ìcoordinateî transformation between covariance structures of interest. In DTI, the diffusion tensor can be reparametrized by the eigenvalues and the Euler angles in order to construct the covariance matrix of the major eigenvector of the diffusion tensor. The eigenvectors of this covariance matrix whose eigenvalues are nonzero are then the major and minor axes of the ellipse of the cone of uncertainty. 1 st^ -order Matrix Perturbation Method The pioneering work of Hext [10] provided a method to construct the covariance structure of the major eigenvector based on the 1st^ -order matrix perturbation, see Eq[4.17] of [10], but the importance of this formula has not been pointed out in previous DTI studies. Here, we suggest Hextís approach as an alternative to the other two methods proposed here for constructing the elliptical cone of uncertainty. Dyadic Product Formalism This formalism was proposed by Basser et al.[3] to overcome the problem of bias in computing sample mean of the eigenvector. Here, we show how the cone of uncertainty can be constructed from the average dyadics. Let {^ q^11^^ ,L, q^1 N } be a sample of major eigenvectors and let =^ ∑^ = iN i Ti N q (^) 1 q 1 T^^11 q 1 q 1 be the average dyadic [3]. Further, let the eigenvalue decomposition of the average dyadics be
The analysis presented above shows that the cone of uncertainty has an elliptical not circular cross section. Particularly, the 1st^ -order matrix perturbation method shows that the cone of uncertainty is elliptical when the two smaller eigenvalues of the diffusion tensor are not equal to each other. Only when they are equal will the cone of uncertainty have the shape of a right circular cone. Similar observations can be made based on the other two methods discussed here. Figure 1B shows the cones of uncertainty constructed from a simulated human brain data based on the error propagation method [1] and the nonlinear least squares method [11]; Figure 1A is the FA map. The implications for DTI tractography are clear. Since the medium and minor eigenvalues of most tensor estimates in the brain are generally different [12], the axial symmetry of the tract dispersion is usually broken. In streamline based or probabilistic methods [13,14], the chance of deviating away from a tract in the direction along which uncertainty is low is less likely than from deviating away from the tract in the direction along which uncertainty is high. Thus, we would predict an asymmetric spread of tracts trajectories normal to the mean trajectory in which paths tending to meander along the direction with lower variance will remain closer to the mean tract whereas tracts meandering along the direction with higher variance will exhibit more dispersion and lie farther from the mean. Generalizing the approach given in [15], we can view the tractography process as a biased random walk in which there is anisotropic diffusion of the tract trajectory normal to the mean tract. In fact, such behavior appears to have been observed in bootstrap simulations of tracts generated using a streamline following method. In Lazar et al., white matter ìaxial asymmetryî (i.e., inequality of λ 2 and λ 3 ) caused anisotropic dispersion patterns in the estimated tract trajectories [8] but a theoretical explanation was not provided for this finding.
Figure 1A Figure 1B REFERENCES [1] Koay CG et al. Proc. ISMRM. 14 (2006). [2] Chang LC et al. Proc. ISMRM. 14 (2006). [3] Basser PJ. Proc. ISMRM. 5 (1997). [4] Basser PJ et al. MRM. 2000; 44: 41-50. [5] Jones DK. MRM 2003;49(1):7-12. [6] Anderson AW. MRM 2001; 46:1174-1188. [7] Lazar M et al. NeuroImage 2003; 20:1140-1153. [8] Lazar M et al. MRM. 2005; 54: 860-867. [9] Wu YC et al. MRM. 2004; 52:1146-1155. [10] Hext GR. Biometrika 1963; 50: p 353. [11] Koay CG et al. J Magn Reson. 2006; 182: 115-125. [12] Pierpaoli C et al. Radiology 1996; 201:637-648. [13] Mori S et al.Ann Neurol 1999; 45:265-269. [14] Behrens TEJ et al. MRM. 2003; 50: 1077-1088. [15] Boguna M et al. New Journal of Physics 2005;7.
CONCLUSION We have presented three independent and complementary methods to show that the cone of uncertainty is generally a right elliptical cone rather than a right circular cone. We showed that the key geometric object for constructing the local measure of tract dispersion is the covariance matrix of the major eigenvector. We have also provided the theoretical basis on which the observation of anistropic dispersion patterns in the estimated tract trajectories caused by axial asymmetric diffusion of water diffusion in white matter can now be explained.