A comparison of sequential and GPU-accelerated implementations of B-spline signal processing operations for 2-D and 3-D images

Abstract

B-spline signal processing operations are widely used in the analysis of two and three-dimensional images. In this paper, we investigate and compare some of these basic operations (direct transformations, indirect transformations, and computation of partial derivatives) by (1) recursive filter based implementations in MATLAB and C++, and (2) GPU-accelerated implementations in CUDA. All operations are compared at a variety of resolution levels on a 2-D panoramic image as well as a 3-D magnetic resonance (MR) image. Results indicate significant improvements in efficiency for the CUDA implementations. A MATLAB toolkit implementing the various B-spline signal processing tasks as well as the C++ and CUDA implementation described here is currently publicly available.