320 likes | 485 Views
Multi-valued Geodesic based Fiber Tracking for Diffusion Tensor Imaging. Neda Sepasian. Supervised by Prof. Bart ter Haar Romeny, Dr. Anna Vilanova Bartoli Dr. J.H.M. ten Thije Boonkkamp. Overview. Diffusion tensor imaging(DTI) Fiber tracking Results Conclusion. Fiber Tracking.
E N D
Multi-valued Geodesic based Fiber Tracking for Diffusion Tensor Imaging Neda Sepasian Supervised by Prof. Bart ter Haar Romeny, Dr. Anna Vilanova Bartoli Dr. J.H.M. ten Thije Boonkkamp
Overview • Diffusion tensor imaging(DTI) • Fiber tracking • Results • Conclusion
Fiber Tracking Results Conclusion DTI MRI can be used to obtain local chemical and physical properties of water. Molecular diffusion Flow
Fiber Tracking Results Conclusion DTI Diffusion Tensor Imaging Measuring the diffusion of water molecules gives us the shape and orientation of the diffusion ellipsoid.
Fiber Tracking Results Conclusion DTI Low anisotropy Suitable for understanding the structure locally. Clutter in 3D Difficult to understand global structure High anisotropy
Fiber Tracking Results Conclusion DTI Fiber Tracking: Provides a potential method for exploring a connectivity network of the brain.
Fiber Tracking Results Conclusion DTI Streamline • Using only the dominant eigenvalue. • deviations in the eigenvectors caused the accumulate error. • In an isotropic region • We are locally maximizing the diffusion.
Fiber Tracking Results Conclusion DTI Streamline
Fiber Tracking Results Conclusion DTI Geodesics • The shortest path between points on the space. • Geodesics can be reconstructed using: • PDE based algorithms(eg. Eikonaleq.) • ODE based algorithms(Euler Lagrange eq.) Euler-Lagrange(EL) solution Correct solution Eikonal Solution
Fiber Tracking Results Conclusion DTI Eikonal equation
Fiber Tracking Results Conclusion DTI Eikonal equation • Solve the Eikonal equation using the numerical approximation: • Charpit’s system to reconstruct the fibers:
Fiber Tracking Results Conclusion DTI Eikonal equation Fibers are selected using connectivity measure:
Fiber Tracking Results Conclusion DTI Eikonal equation
Fiber Tracking Results Conclusion DTI Eikonal equation • It is globally minimizing the geodesics using the inverse of the diffusion tensors. • Therefore it is more robust to noise but at the same time less sensitive to local orientations. • Only the first arrival time (unique solution) is computed at each grid point.
Fiber Tracking Results Conclusion DTI Euler-Lagrange Equation
Fiber Tracking Results Conclusion DTI Euler-Lagrange Equation • Solve the geodesic ODEs using well-known ODE solver like RK4.
Fiber Tracking Results Conclusion DTI Euler-Lagrange Equation • Shoot rays in different initial direction with the same initial position. • Apply ray-tracing algorithm for finding the geodesic connecting two given points.
Fiber Tracking Results Conclusion DTI Euler-Lagrange Equation
Fiber Tracking Results Conclusion DTI Euler-Lagrange Equation
Fiber Tracking Results Conclusion DTI Eikonal EL
Fiber Tracking Results Conclusion DTI Classic fiber-tracking PDE based fiber-tracking
Fiber Tracking Results Conclusion DTI EL based fiber-tracking
Fiber Tracking Results Conclusion DTI HJ EL
Fiber Tracking Results Conclusion DTI i iii Corpus Callosum (CC) trackts based on atlas Gray’s anatomy CC tracts using EL based algorithm ii
Fiber Tracking Results Conclusion DTI EL based method
Fiber Tracking Results Conclusion DTI (a) Arcuate fasciculus (ARC) ( f ) Uncinate fasciculus (UNC) EL based fiber-tracking
Fiber Tracking Conclusion Results DTI Eikonal solution EL solution • Global minimization • Robust to noise • Accuracy for quantitative analysis • Algorithm efficiency • Only the first arrival time • Global minimization • Robust to noise • Accuracy for quantitative analysis • Algorithm efficiency • Multi-valued solution. • Less information is deduced from the computation
Fiber Tracking Results Conclusion DTI What could be an ideal algorithm ???
Fiber Tracking Conclusion Results DTI Other Challenges • Single tensor models are not sufficient • Fiber-tracking algorithms are still imperfect
Fiber Tracking Conclusion Results DTI Work in progress!!! • Multi-valued HARDI fiber-tracking in single processor DTI HARDI • Multi-valued HARDI fiber-tracking in GPU (using CUDA)
Acknowledgements • Dr. Olof Runborg (KTH, Stockholm, Sweden) • Luc Florack • Laura Astola • Evert Aart