On removing interpolation and resampling artifacts in rigid image registration.
We show that image registration using conventional interpolation and summation approximations of continuous integrals can generally fail because of resampling artifacts. These artifacts negatively affect the accuracy of registration by producing local optima, altering the gradient, shifting the global optimum, and making rigid registration asymmetric. In this paper, after an extensive literature review, we demonstrate the causes of the artifacts by comparing inclusion and avoidance of resampling analytically. We show the sum-of-squared-differences cost function formulated as an integral to be more accurate compared with its traditional sum form in a simple case of image registration. We then discuss aliasing that occurs in rotation, which is due to the fact that an image represented in the Cartesian grid is sampled with different rates in different directions, and propose the use of oscillatory isotropic interpolation kernels, which allow better recovery of true global optima by overcoming this type of aliasing. Through our experiments on brain, fingerprint, and white noise images, we illustrate the superior performance of the integral registration cost function in both the Cartesian and spherical coordinates, and also validate the introduced radial interpolation kernel by demonstrating the improvement in registration.
Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA 02129, USA. iman@nmr.mgh.harvard.edu
Iman
I
Aganj
Boon Thye Thomas
BT
Yeo
Mert R
MR
Sabuncu
Bruce
B
Fischl
D000465:Algorithms; D016477:Artifacts; D001921:Brain / Q000033:anatomy & histology; D003878:Dermatoglyphics; D006801:Humans; D007091:Image Processing, Computer-Assisted / Q000379:methods*; D008279:Magnetic Resonance Imaging