Diffusion tensor imaging is widely used in brain connectivity research. As more and more studies recruit large numbers of subjects, it is important to design registration methods which are not only theoretically rigorous, but also computationally efficient. However, the requirement of reorienting diffusion tensors complicates and considerably slows down registration procedures, due to the correlated impacts of registration forces at adjacent voxel locations. Based on the diffeomorphic Demons algorithm (Vercauteren , 2009), we propose a fast local trust region algorithm for handling inseparable registration forces for quadratic energy functions. The method guarantees that, at any time and at any voxel location, the velocity is always within its local trust region. This local regularization allows efficient calculation of the transformation update with numeric integration instead of completely solving a large linear system at every iteration. It is able to incorporate exact reorientation and regularization into the velocity optimization, and preserve the linear complexity of the diffeomorphic Demons algorithm. In an experiment with 84 diffusion tensor images involving both pair-wise and group-wise registrations, the proposed algorithm achieves better registration in comparison with other methods solving large linear systems (Yeo , 2009). At the same time, this algorithm reduces the computation time and memory demand tenfold.

Human brain connectivity can be studied using graph theory. Many connectivity studies parcellate the brain into regions and count fibres extracted between them. The resulting network analyses require validation of the tractography, as well as region and parameter selection. Here we investigate whole brain connectivity from a different perspective. We propose a mathematical formulation based on studying the eigenvalues of the Laplacian matrix of the diffusion tensor field at the voxel level. This voxelwise matrix has over a million parameters, but we derive the Kirchhoff complexity and eigen-spectrum through elegant mathematical theorems, without heavy computation. We use these novel measures to accurately estimate the voxelwise connectivity in multiple biomedical applications such as Alzheimer's disease and intelligence prediction.