EPFL | Biomedical Imaging Group

Discretization of the Radon Transform and of Its Inverse by Spline Convolutions

S. Horbelt, M. Liebling, M. Unser

IEEE Transactions on Medical Imaging, vol. 21, no. 4, pp. 363–376, April 2002.

We present an explicit formula for B-spline convolution kernels; these are defined as the convolution of several B-splines of variable widths h_i and degrees n_i. We apply our results to derive spline-convolution-based algorithms for two closely related problems: the computation of the Radon transform and of its inverse. First, we present an efficient discrete implementation of the Radon transform that is optimal in the least-squares sense. We then consider the reverse problem and introduce a new spline-convolution version of the filtered back-projection algorithm for tomographic reconstruction. In both cases, our explicit kernel formula allows for the use of high-degree splines; these offer better approximation performance than the conventional lower-degree formulations (e.g., piecewise constant or piecewise linear models). We present multiple experiments to validate our approach and to find the parameters that give the best tradeoff between image quality and computational complexity. In particular, we find that it can be computationally more efficient to increase the approximation degree than to increase the sampling rate.

@ARTICLE(http://bigwww.epfl.ch/publications/horbelt0201.html,
AUTHOR="Horbelt, S. and Liebling, M. and Unser, M.",
TITLE="Discretization of the {R}adon Transform and of Its Inverse by
	Spline Convolutions",
JOURNAL="{IEEE} Transactions on Medical Imaging",
YEAR="2002",
volume="21",
number="4",
pages="363--376",
month="April",
note="")

© 2002 IEEE. Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or redistribution to servers or lists, or to reuse any copyrighted component of this work in other works must be obtained from IEEE. This material is presented to ensure timely dissemination of scholarly and technical work. Copyright and all rights therein are retained by authors or by other copyright holders. All persons copying this information are expected to adhere to the terms and constraints invoked by each author's copyright. In most cases, these works may not be reposted without the explicit permission of the copyright holder.