Algorithm for hyperfast cone-beam spiral backprojection

Cone-beam spiral backprojection is computationally highly demanding. At first sight, the backprojection requirements are similar to those of cone-beam backprojection from circular scans such as it is performed in the widely used Feldkamp algorithm. However, there is an additional complication: the illumination of each voxel, i.e. the range of angles the voxel is seen by the X-ray cone is a complex function of the voxel position. The weight function has no analytically closed form and must be numerically determined. Storage of the weights is prohibitive since the amount of memory required equals the number of voxels per spiral rotation times the number of projections a voxel receives contributions and therefore is in the order of 10(9) to 10(11) floating point values for typical spiral scans. We propose a new algorithm that combines the spiral symmetry with the ability of today's 64 bit CPUs to store large amounts of precomputed weights. Using the spiral symmetry in this way allows to exploit data-level parallelism and thereby to achieve a very high level of vectorization. An additional postprocessing step rotates these slices back to normal images. Our new backprojection algorithm achieves up to 24.6 Giga voxel updates per second (GUPS) on our systems that are equipped with two standard Intel X5570 quad core CPUs (Intel Xeon 5500 platform, 2.93 GHz, Intel Corporation). This equals the reconstruction of 410 images per second assuming each slice consists of 512 x 512 pixels, receiving contributions from 512 projections.

[1]  Frédéric Noo,et al.  Exact helical reconstruction using native cone-beam geometries. , 2003, Physics in medicine and biology.

[2]  M. Knaup,et al.  CT image reconstruction using hexagonal grids , 2007, 2007 IEEE Nuclear Science Symposium Conference Record.

[3]  M. Defrise,et al.  A solution to the long-object problem in helical cone-beam tomography. , 2000, Physics in medicine and biology.

[4]  Michael Knaup,et al.  Hyperfast parallel-beam and cone-beam backprojection using the cell general purpose hardware. , 2007, Medical physics.

[5]  R. Proksa,et al.  A quasiexact reconstruction algorithm for helical CT using a 3-Pi acquisition. , 2003, Medical physics.

[6]  Alexander Katsevich,et al.  Theoretically Exact Filtered Backprojection-Type Inversion Algorithm for Spiral CT , 2002, SIAM J. Appl. Math..

[7]  F. Makedon,et al.  High quality alias free image rotation , 1996, Conference Record of The Thirtieth Asilomar Conference on Signals, Systems and Computers.

[8]  R. W. Schafer,et al.  A comparison of rotation-based methods for iterative reconstruction algorithms , 1995 .

[9]  J. Hsieh,et al.  A three-dimensional-weighted cone beam filtered backprojection (CB-FBP) algorithm for image reconstruction in volumetric CT—helical scanning , 2006, Physics in medicine and biology.

[10]  Michael Grass,et al.  The n-PI-method for helical cone-beam CT , 2000, IEEE Transactions on Medical Imaging.

[11]  Willi A Kalender,et al.  Multithreaded cardiac CT. , 2006, Medical physics.

[12]  Katsuyuki Taguchi,et al.  A new weighting scheme for cone-beam helical CT to reduce the image noise. , 2004, Physics in medicine and biology.

[13]  Michael Knaup,et al.  GPU-based parallel-beam and cone-beam forward- and backprojection using CUDA , 2008, 2008 IEEE Nuclear Science Symposium Conference Record.

[14]  Er-Wei Bai,et al.  A Fast CT Reconstruction Scheme for a General Multi-Core PC , 2007, Int. J. Biomed. Imaging.

[15]  Willi A Kalender,et al.  Extended parallel backprojection for standard three-dimensional and phase-correlated four-dimensional axial and spiral cone-beam CT with arbitrary pitch, arbitrary cone-angle, and 100% dose usage. , 2004, Medical physics.

[16]  Alexander Katsevich,et al.  Theoretically exact FBP-type inversion algorithm for spiral CT , 2001 .

[17]  Salvatore Lanzavecchia,et al.  Image and volume data rotation with 1- and 3-pass algorithms , 1996, Comput. Appl. Biosci..

[18]  Marc Kachelrieß,et al.  Advanced single-slice rebinning in cone-beam spiral CT. , 2000 .

[19]  Michael Knaup,et al.  High performance cone-beam spiral backprojection with voxel-specific weighting. , 2009, Physics in medicine and biology.

[20]  K. Taguchi,et al.  Algorithm for image reconstruction in multi-slice helical CT. , 1998, Medical physics.

[21]  Wenyuan Bi,et al.  Accelerate helical cone-beam CT with graphics hardware , 2008, SPIE Medical Imaging.

[22]  Karl Stierstorfer,et al.  Segmented multiple plane reconstruction: a novel approximate reconstruction scheme for multi-slice spiral CT. , 2002, Physics in medicine and biology.