Non-rigid brain image registration using a statistical deformation model

In this article, we propose a new registration method, based on a statistical analysis of deformation fields. At first, a set of MRI brain images was registered using a viscous fluid algorithm.(1) The obtained deformation fields are then used to calculate a Principal Component Analysis (PCA) based decomposition. Since PCA models the deformations as a linear combination of statistically uncorrelated principal components, new deformations can be created by changing the coefficients in the linear combination. We then use the PCA representation of the deformation fields to non-rigidly align new sets of images. We use a gradient descent method to adjust the coefficients of the principal components, such that the resulting deformation maximizes the mutual information between the deformed image and an atlas image. The results of our method are promising. Viscous fluid registrations of new images can be recovered with an accuracy of about half a voxel. Better results can be obtained by using a more extensive database of learning images (we only used 84). Also, the optimization method used here can be improved, especially to shorten computation time.