Fast Multidimensional B-spline Interpolation Using Template Metaprogramming

B-spline interpolation is a widely used interpolation technique. In the field of image registration, interpolation is necessary for transforming images to obtain a measure of (dis)similarity between the images to be aligned. When gradient-based optimization methods are used, the image gradients need to be calculated as well, which also accounts for a substantial share of computation time in registration. In this paper we propose a fast multidimensional B-spline interpolation algorithm with which both image value and gradient can be computed efficiently. We present a recursive algorithm for the interpolation which is efficiently implemented with template metaprogramming (TMP). The proposed algorithm is compared with the algorithm implemented in the Insight Toolkit (ITK), for different interpolation orders and image dimensions. Also, the effect on the computation time of a typical registration problem is evaluated. The results show that the computation time of B-spline interpolation is decreased by the proposed algorithm from a factor 4.1 for a 2D image using 1st order interpolation to a factor of 19.9 for 4D using 3rd order interpolation.

[1]  Michael Unser,et al.  Optimization of mutual information for multiresolution image registration , 2000, IEEE Trans. Image Process..

[2]  Daniel Rueckert,et al.  Nonrigid registration using free-form deformations: application to breast MR images , 1999, IEEE Transactions on Medical Imaging.

[3]  Todd L. Veldhuizen,et al.  Using C++ template metaprograms , 1996 .

[4]  Max A. Viergever,et al.  elastix: A Toolbox for Intensity-Based Medical Image Registration , 2010, IEEE Transactions on Medical Imaging.

[5]  M. Unser,et al.  Interpolation revisited [medical images application] , 2000, IEEE Transactions on Medical Imaging.

[6]  Michael Unser,et al.  Splines: a perfect fit for signal and image processing , 1999, IEEE Signal Process. Mag..

[7]  Max A. Viergever,et al.  Adaptive Stochastic Gradient Descent Optimisation for Image Registration , 2009, International Journal of Computer Vision.

[8]  David Sarrut,et al.  B-LUT: Fast and low memory B-spline image interpolation , 2010, Comput. Methods Programs Biomed..