Abstract:Purpose: 4DCT is an important technique for motion management in radiotherapy, but involves a significant amount of dose to patients. Low-dose helical 4DCT scans can be done by (a) lowering tube current and pulse duration, (b) increasing helical pitch, and (c) reducing data sampling. Conventional FDK methods may fail for these low-dose scan protocols due to noisy or incomplete data. We aim to develop an accurate and fast image reconstruction solver that enables dose reduction for helical 4DCT. Methods: For improved accuracy, the developed solver is based on a novel 4D regularization method for CT image reconstruction, which promotes both the smoothness in spatial variables for each phase and the similarity in temporal variables between phases, in contrast with the 3D regularization method with spatial regularization only. Both 3D and 4D methods here are based on total variation. The formulated optimization problem is solved by alternating direction method of multipliers (ADMM), which is chosen for its effectiveness for solving the convex problem, and more importantly its suitability for distributed computing.For improved speed, the developed solver utilizes an efficient parallel algorithm for X-ray transform and its adjoint with optimized computational cost per parallel thread O (1) that takes flying focal spots into account, and multi-GPU platforms for solving such a massive parallel computing problem in helical 4DCT. Results: The solver was validated using both simulations (4D NCAT phantoms) and experimental sinogram data (Siemens, Somatom 64-slices) with low-dose, high-pitch and sparse-view helical 4DCT scans. Both the reconstructed images and quantitative results demonstrate that the 4D method improved image quality compared with the 3D method, which is in turn better than FDK. Conclusion: We have developed a fast, accurate solver that synergizes state-of-art iterative reconstruction methods and parallel computing techniques that can potentially enable low-dose, high-pitch and sparse-view helical 4DCT scans. Minghao Guo and Hao Gao were partially supported by the NSFC (#11405105), the 973 Program (#2015CB856000) and the Shanghai Pujiang Talent Program (#14PJ1404500).

Fast 3-D Reconstruction of CT Images Using a Geometric Parameter Table

Fast and Quantitative Reconstruction of 4D SPECT for Fan-beam Geometry Based on Novikov's Inverse Formula

A fast iteration approach to undersampled cone-beam CT reconstruction.

A Fast Regularized Iterative Algorithm for Fan-Beam CT Reconstruction.

Accelerated 3-D T-FDK Reconstruction Algorithm Using the Commordity Graphics Hardware

CT image reconstruction via industrial CT fast scanning

Accelerate Helical Cone-Beam CT with Graphics Hardware

3D Algebraic Iterative Reconstruction for Cone-Beam X-Ray Differential Phase-Contrast Computed Tomography

Accelerating 3D cone beam T-FDK algorithm on commodity PC graphics hardware

Fast Imaging By A Single-Slice-Detector Helical Ct

An Exact and Fast CBCT Reconstruction via Pseudo-Polar Fourier Transform based Discrete Grangeat's Formula

GPU-based 3D cone-beam CT image reconstruction: application to micro CT

iFDK: A Scalable Framework for Instant High-resolution Image Reconstruction

FPGA Acceleration for 3-D Low-Dose Tomographic Reconstruction.

GPU-based iterative cone-beam CT reconstruction using tight frame regularization

TH-CD-303-02: A GPU-Based Iterative Image Reconstruction Solver with 4D Regularization for Low-Dose Helical 4DCT

FPGA Acceleration of Ray-Based Iterative Algorithm for 3D Low-Dose CT Reconstruction.

Parallel processing system for fast 3D-CT image reconstruction by circular shifting float memory

A Free-Geometry Cone Beam Ct and Its Fdk-Type Reconstruction

GPU Accelerated Voxel-Driven Forward Projection for Iterative Reconstruction of Cone-Beam CT

Generalized-Equiangular Geometry CT: Concept and Shift-Invariant FBP Algorithms