Accelerating simultaneous algebraic reconstruction technique with motion compensation using CUDA-enabled GPU

PurposeTo accelerate the simultaneous algebraic reconstruction technique (SART) with motion compensation for speedy and quality computed tomography reconstruction by exploiting CUDA-enabled GPU.MethodsTwo core techniques are proposed to fit SART into the CUDA architecture: (1) a ray-driven projection along with hardware trilinear interpolation, and (2) a voxel-driven back-projection that can avoid redundant computation by combining CUDA shared memory. We utilize the independence of each ray and voxel on both techniques to design CUDA kernel to represent a ray in the projection and a voxel in the back-projection respectively. Thus, significant parallelization and performance boost can be achieved. For motion compensation, we rectify each ray’s direction during the projection and back-projection stages based on a known motion vector field.ResultsExtensive experiments demonstrate the proposed techniques can provide faster reconstruction without compromising image quality. The process rate is nearly 100 projections s−1, and it is about 150 times faster than a CPU-based SART. The reconstructed image is compared against ground truth visually and quantitatively by peak signal-to-noise ratio (PSNR) and line profiles. We further evaluate the reconstruction quality using quantitative metrics such as signal-to-noise ratio (SNR) and mean-square-error (MSE). All these reveal that satisfactory results are achieved. The effects of major parameters such as ray sampling interval and relaxation parameter are also investigated by a series of experiments. A simulated dataset is used for testing the effectiveness of our motion compensation technique. The results demonstrate our reconstructed volume can eliminate undesirable artifacts like blurring.ConclusionOur proposed method has potential to realize instantaneous presentation of 3D CT volume to physicians once the projection data are acquired.

[1]  Klaus Mueller,et al.  Rapid 3-D cone-beam reconstruction with the simultaneous algebraic reconstruction technique (SART) using 2-D texture mapping hardware , 2000, IEEE Transactions on Medical Imaging.

[2]  A. Kak,et al.  Simultaneous Algebraic Reconstruction Technique (SART): A Superior Implementation of the Art Algorithm , 1984, Ultrasonic imaging.

[3]  Seung-Oh Jin,et al.  Fast image reconstruction from fan beam projections using parallel digital signal processors and special purpose processors , 1999, Proceedings of IEEE. IEEE Region 10 Conference. TENCON 99. 'Multimedia Technology for Asia-Pacific Information Infrastructure' (Cat. No.99CH37030).

[4]  C. Ling,et al.  Effect of respiratory gating on quantifying PET images of lung cancer. , 2002, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[5]  Brian Cabral,et al.  Accelerated volume rendering and tomographic reconstruction using texture mapping hardware , 1994, VVS '94.

[6]  Gabor T. Herman.  Algebraic Reconstruction Techniques , 2009 .

[7]  M. Knaup,et al.  Hyperfast Perspective Cone--Beam Backprojection , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

[8]  E. Sidky,et al.  Accurate image reconstruction from few-views and limited-angle data in divergent-beam CT , 2009, 0904.4495.

[9]  S. Rit,et al.  Cone-beam projection of a deformable volume for motion compensated algebraic reconstruction , 2007, 2007 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[10]  Torsten Möller,et al.  Rapid Emission Tomography Reconstruction , 2003, VG.

[11]  J. Hornegger,et al.  Fast GPU-Based CT Reconstruction using the Common Unified Device Architecture (CUDA) , 2007, 2007 IEEE Nuclear Science Symposium Conference Record.

[12]  Fang Xu,et al.  Accelerating popular tomographic reconstruction algorithms on commodity PC graphics hardware , 2005, IEEE Transactions on Nuclear Science.

[13]  Klaus Mueller,et al.  IOP PUBLISHING PHYSICS IN MEDICINE AND BIOLOGY , 2007 .

[14]  L. Feldkamp,et al.  Practical cone-beam algorithm , 1984 .

[15]  M. V. van Herk,et al.  Precise and real-time measurement of 3D tumor motion in lung due to breathing and heartbeat, measured during radiotherapy. , 2002, International journal of radiation oncology, biology, physics.

[16]  Rüdiger Westermann,et al.  MR image reconstruction using the GPU , 2006, SPIE Medical Imaging.

[17]  Dmitri Riabkov,et al.  Accelerated cone-beam backprojection using GPU-CPU hardware , 2022 .

[18]  Michel Desvignes,et al.  Hardware/software 2D-3D backprojection on a SoPC platform , 2006, SAC.

[19]  G Malandain,et al.  Model-based respiratory motion compensation for emission tomography image reconstruction , 2007, Physics in medicine and biology.

[20]  Avinash C. Kak,et al.  Principles of computerized tomographic imaging , 2001, Classics in applied mathematics.

[21]  R. Siddon.  Fast calculation of the exact radiological path for a three-dimensional CT array. , 1985, Medical physics.

[22]  Fumihiko Ino,et al.  Accelerating cone beam reconstruction using the CUDA-enabled GPU , 2008, HiPC'08.

[23]  K. Lam,et al.  Uncertainties in CT-based radiation therapy treatment planning associated with patient breathing. , 1996, International journal of radiation oncology, biology, physics.

[24]  Myoung-Hee Kim,et al.  A fast algorithm for voxel-based deterministic simulation of X-ray imaging , 2008, Comput. Phys. Commun..

[25]  Pierre Grangeat,et al.  Dynamic X-ray computed tomography , 2003, Proc. IEEE.