Fast parallel algorithms for the x-ray transform and its adjoint.

PURPOSE Iterative reconstruction methods often offer better imaging quality and allow for reconstructions with lower imaging dose than classical methods in computed tomography. However, the computational speed is a major concern for these iterative methods, for which the x-ray transform and its adjoint are two most time-consuming components. The speed issue becomes even notable for the 3D imaging such as cone beam scans or helical scans, since the x-ray transform and its adjoint are frequently computed as there is usually not enough computer memory to save the corresponding system matrix. The purpose of this paper is to optimize the algorithm for computing the x-ray transform and its adjoint, and their parallel computation. METHODS The fast and highly parallelizable algorithms for the x-ray transform and its adjoint are proposed for the infinitely narrow beam in both 2D and 3D. The extension of these fast algorithms to the finite-size beam is proposed in 2D and discussed in 3D. RESULTS The CPU and GPU codes are available at https://sites.google.com/site/fastxraytransform. The proposed algorithm is faster than Siddon's algorithm for computing the x-ray transform. In particular, the improvement for the parallel computation can be an order of magnitude. CONCLUSIONS The authors have proposed fast and highly parallelizable algorithms for the x-ray transform and its adjoint, which are extendable for the finite-size beam. The proposed algorithms are suitable for parallel computing in the sense that the computational cost per parallel thread is O(1).

[1]  Lei Zhu,et al.  Quantitative cone-beam CT imaging in radiation therapy using planning CT as a prior: first patient studies. , 2012, Medical physics.

[2]  Jiang Hsieh,et al.  Step-and-shoot data acquisition and reconstruction for cardiac x-ray computed tomography. , 2006, Medical physics.

[3]  E. Sidky,et al.  Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization , 2008, Physics in medicine and biology.

[4]  L. Xing,et al.  Overview of image-guided radiation therapy. , 2006, Medical dosimetry : official journal of the American Association of Medical Dosimetrists.

[5]  G. Christensen,et al.  A method for the reconstruction of four-dimensional synchronized CT scans acquired during free breathing. , 2003, Medical physics.

[6]  Steve B. Jiang,et al.  Quality assurance challenges for motion-adaptive radiation therapy: gating, breath holding, and four-dimensional computed tomography. , 2008, International journal of radiation oncology, biology, physics.

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

[8]  T. Guerrero,et al.  Acquiring 4D thoracic CT scans using a multislice helical method. , 2004, Physics in medicine and biology.

[9]  Hongkai Zhao,et al.  Robust principal component analysis-based four-dimensional computed tomography , 2011, Physics in medicine and biology.

[10]  Jeffrey A. Fessler,et al.  3D Forward and Back-Projection for X-Ray CT Using Separable Footprints , 2010, IEEE Transactions on Medical Imaging.

[11]  Jie Tang,et al.  Prior image constrained compressed sensing (PICCS): a method to accurately reconstruct dynamic CT images from highly undersampled projection data sets. , 2008, Medical physics.

[12]  L. Xing,et al.  Iterative image reconstruction for CBCT using edge-preserving prior. , 2008, Medical physics.

[13]  R. Mohan,et al.  Acquiring a four-dimensional computed tomography dataset using an external respiratory signal. , 2003, Physics in medicine and biology.

[14]  Huaxia Zhao,et al.  Fast ray-tracing technique to calculate line integral paths in voxel arrays , 2003, 2003 IEEE Nuclear Science Symposium. Conference Record (IEEE Cat. No.03CH37515).

[15]  S. Osher,et al.  Multi-energy CT based on a prior rank, intensity and sparsity model (PRISM) , 2011, Inverse problems.

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

[17]  Steve B. Jiang,et al.  GPU-based fast cone beam CT reconstruction from undersampled and noisy projection data via total variation , 2010 .

[18]  Tinsu Pan,et al.  Four-dimensional computed tomography: image formation and clinical protocol. , 2005, Medical physics.

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

[20]  B. De Man,et al.  Distance-driven projection and backprojection in three dimensions. , 2004, Physics in medicine and biology.

[21]  Steve B. Jiang,et al.  GPU-based iterative cone-beam CT reconstruction using tight frame regularization , 2010, Physics in medicine and biology.

[22]  Sabee Molloi,et al.  Quantification of breast density with dual energy mammography: an experimental feasibility study. , 2010, Medical physics.

[23]  Xiaochuan Pan,et al.  GPU-based 3D cone-beam CT image reconstruction: application to micro CT , 2007, 2007 IEEE Nuclear Science Symposium Conference Record.

[24]  Hengyong Yu,et al.  Compressed sensing based interior tomography , 2009, Physics in medicine and biology.

[25]  Ioannis Sechopoulos,et al.  X-ray scatter correction method for dedicated breast computed tomography. , 2012, Medical physics.

[26]  A. Boyer,et al.  Radiation dose reduction in four-dimensional computed tomography. , 2005, Medical physics.

[27]  Zhengrong Liang,et al.  Dose reduction for kilovotage cone-beam computed tomography in radiation therapy. , 2008, Physics in medicine and biology.