Multi-scale regularization approaches of non-parametric deformable registrations

Most deformation algorithms use a single-value smoother during optimization. We investigate multi-scale regularizations (smoothers) during the multi-resolution iteration of two non-parametric deformable registrations (demons and diffeomorphic algorithms) and compare them to a conventional single-value smoother. Our results show that as smoothers increase, their convergence rate decreases; however, smaller smoothers also have a large negative value of the Jacobian determinant suggesting that the one-to-one mapping has been lost; i.e., image morphology is not preserved. A better one-to-one mapping of the multiscale scheme has also been established by the residual vector field measures. In the demons method, the multiscale smoother calculates faster than the large single-value smoother (Gaussian kernel width larger than 0.5) and is equivalent to the smallest single-value smoother (Gaussian kernel width equals to 0.5 in this study). For the diffeomorphic algorithm, since our multi-scale smoothers were implemented at the deformation field and the update field, calculation times are longer. For the deformed images in this study, the similarity measured by mean square error, normal correlation, and visual comparisons show that the multi-scale implementation has better results than large single-value smoothers, and better or equivalent for smallest single-value smoother. Between the two deformable registrations, diffeormophic method constructs better coherence space of the deformation field while the deformation is large between images.