Positive Semidefinite Generalized Diffusion Tensor Imaging Via Quadratic Semidefinite Programming.

Yannan Chen,Yuhong Dai,Deren Han,Wenyu Sun
DOI: https://doi.org/10.1137/110843526
IF: 1.938
2013-01-01
SIAM Journal on Imaging Sciences
Abstract:The positive definiteness of a diffusion tensor is importantin magnetic resonance imaging because it reflectsthe phenomenon of water molecular diffusion in complicated biological tissue environments.To preserve this property, we represent it asan explicit positive semidefinite (PSD) matrix constraintand some linear matrix equalities.The objective function is the regularized linear least squares fitting forthe log-linearized Stejskal--Tanner equation.The regularization term is the heuristic nuclear norm of the PSD matrix,since we expect it to be of low rank.In this way, we establish a convex quadratic semidefinite programming (SDP) model,whose global solution exists.The optimal solution could be solved by three efficient methods.While there are two state-of-the-art solvers---SDPT3 and QSDP---for the primal problem,we design a new augmented Lagrangian basedalternating direction method (ADM) for the dual problem.Sensitivity analyses on the coefficients of the optimal diffusion tensor andthe optimal objective function value with respect to noise-corrupted signals are presented.Experiments on synthetic data with multiple fibersshow that the new method is robust to the Rician noiseand outperforms several existing methods.Furthermore, when the fiber orientation distribution function is considered,the new method is competitive with the Q-ball imaging.Using the human brain data, we illustrate thatthe new method could capture the crossing of three nervous fiber bundles.Additionally, the new method generates positive definite generalized diffusion tensorsin all voxels, while the unconstrained least squares fitting fails.Finally, we confirm that the ADM solver is more efficientthan SDPT3 and QSDP for this special problem.
What problem does this paper attempt to address?