A constrained variational principle for direct estimation and smoothing of the diffusion tensor field from complex DWI

In this paper, we present a novel constrained variational principle for simultaneous smoothing and estimation of the diffusion tensor field from complex valued diffusion-weighted images (DWI). The constrained variational principle involves the minimization of a regularization term of L^p norms, subject to a nonlinear inequality constraint on the data. The data term we employ is the original Stejskal-Tanner equation instead of the linearized version usually employed in literature. The complex valued nonlinear form leads to a more accurate (when compared to the linearized version) estimate of the tensor field. The inequality constraint requires that the nonlinear least squares data term be bounded from above by a known tolerance factor. Finally, in order to accommodate the positive definite constraint on the diffusion tensor, it is expressed in terms of Cholesky factors and estimated. The constrained variational principle is solved using the augmented Lagrangian technique in conjunction with the limited memory quasi-Newton method. Experiments with complex-valued synthetic and real data are shown to depict the performance of our tensor field estimation and smoothing algorithm.