Potential of iterative reconstruction for maxillofacial cone beam CT imaging: technical note.

Iterative reconstruction has been proven to be an effective tool for low-dose computed tomography imaging. However, this technology is currently not available in commercial diagnostic maxillofacial cone beam CT. For this technical note, an iterative reconstruction technique was applied to cone beam CT raw data of two maxillofacial clinical cases to explore its potential for dose reduction and metal artifact reduction. Low-dose imaging was emulated by using only fractions of the clinical projection dataset. The reconstruction algorithms tested were filtered backprojection (FBP) as a reference method, and a total variation minimization (TV) regularized ordered subsets convex (OSC-TV) method as the iterative technique. Upon qualitative examination, the OSC-TV technique was found to conserve most diagnostic information using half the projections. Test images have also shown that at 1/4 of the projections, OSC-TV was more robust than FBP with respect to streaking and metal artifacts.

Iterative reconstruction for image enhancement and dose reduction in diagnostic cone beam CT imaging.

BACKGROUND: Iterative reconstruction is well-established in diagnostic multidetector computed tomography (MDCT) for dose reduction and image quality enhancement. Its application to diagnostic cone beam computed tomography (CBCT) is only emerging and warrants a quantitative evaluation. METHODS: Several phantoms and a canine head specimen were imaged using a commercially available small-field CBCT scanner. Raw projection data were reconstructed using the Feldkamp-Davis-Kress (FDK) method with different filters, including denoising via total variation (TV) minimization (FDK-TV). Iterative reconstruction was carried out using the TV-regularized ordered subsets convex technique (OSC-TV). Signal-to-noise ratio (SNR), noise power spectrum (NPS) and spatial resolution of images were estimated. Dose levels were measured via the weighted computed tomography dose index, while low-dose image quality degradation was estimated via structural similarity (SSIM). RESULTS: OSC-TV and FDK-TV were shown to significantly improve image signal-to-noise ratio (SNR) compared to FDK with a standard filter, 5.8 and 4.0 times, respectively. Spatial resolution attained with different algorithms varied moderately across different experiments. For low-dose acquisitions, image quality decreased dramatically for FDK but not for FDK-TV nor OSC-TV. For low-dose canine head images acquired using about 1/5 of the dose compared to a reference image, SSIM dropped to about 0.3 for FDK, while remaining at 0.92 for FDK-TV and 0.96 for OSC-TV. CONCLUSION: OSC-TV was shown to improve image quality compared to FDK and FDK-TV. Moreover, this iterative approach allowed for significant dose reduction while maintaining image quality.

A fast 4D cone beam CT reconstruction method based on the OSC-TV algorithm.

BACKGROUND: Four-dimensional cone beam computed tomography allows for temporally resolved imaging with useful applications in radiotherapy, but raises particular challenges in terms of image quality and computation time. OBJECTIVE: The purpose of this work is to develop a fast and accurate 4D algorithm by adapting a GPU-accelerated ordered subsets convex algorithm (OSC), combined with the total variation minimization regularization technique (TV). METHODS: Different initialization schemes were studied to adapt the OSC-TV algorithm to 4D reconstruction: each respiratory phase was initialized either with a 3D reconstruction or a blank image. Reconstruction algorithms were tested on a dynamic numerical phantom and on a clinical dataset. 4D iterations were implemented for a cluster of 8 GPUs. RESULTS: All developed methods allowed for an adequate visualization of the respiratory movement and compared favorably to the McKinnon-Bates and adaptive steepest descent projection onto convex sets algorithms, while the 4D reconstructions initialized from a prior 3D reconstruction led to better overall image quality. CONCLUSION: The most suitable adaptation of OSC-TV to 4D CBCT was found to be a combination of a prior FDK reconstruction and a 4D OSC-TV reconstruction with a reconstruction time of 4.5 minutes. This relatively short reconstruction time could facilitate a clinical use.

System matrix computation vs storage on GPU: A comparative study in cone beam CT.

PURPOSE: Iterative reconstruction algorithms in computed tomography (CT) require a fast method for computing the intersection distances between the trajectories of photons and the object, also called ray tracing or system matrix computation. This work focused on the thin-ray model is aimed at comparing different system matrix handling strategies using graphical processing units (GPUs). METHODS: In this work, the system matrix is modeled by thin rays intersecting a regular grid of box-shaped voxels, known to be an accurate representation of the forward projection operator in CT. However, an uncompressed system matrix exceeds the random access memory (RAM) capacities of typical computers by one order of magnitude or more. Considering the RAM limitations of GPU hardware, several system matrix handling methods were compared: full storage of a compressed system matrix, on-the-fly computation of its coefficients, and partial storage of the system matrix with partial on-the-fly computation. These methods were tested on geometries mimicking a cone beam CT (CBCT) acquisition of a human head. Execution times of three routines of interest were compared: forward projection, backprojection, and ordered-subsets convex (OSC) iteration. RESULTS: A fully stored system matrix yielded the shortest backprojection and OSC iteration times, with a 1.52× acceleration for OSC when compared to the on-the-fly approach. Nevertheless, the maximum problem size was bound by the available GPU RAM and geometrical symmetries. On-the-fly coefficient computation did not require symmetries and was shown to be the fastest for forward projection. It also offered reasonable execution times of about 176.4 ms per view per OSC iteration for a detector of 512 × 448 pixels and a volume of 3843 voxels, using commodity GPU hardware. Partial system matrix storage has shown a performance similar to the on-the-fly approach, while still relying on symmetries. CONCLUSION: Partial system matrix storage was shown to yield the lowest relative performance. On-the-fly ray tracing was shown to be the most flexible method, yielding reasonable execution times. A fully stored system matrix allowed for the lowest backprojection and OSC iteration times and may be of interest for certain performance-oriented applications.

Evaluation of the OSC-TV iterative reconstruction algorithm for cone-beam optical CT.

PURPOSE: The present work evaluates an iterative reconstruction approach, namely, the ordered subsets convex (OSC) algorithm with regularization via total variation (TV) minimization in the field of cone-beam optical computed tomography (optical CT). One of the uses of optical CT is gel-based 3D dosimetry for radiation therapy, where it is employed to map dose distributions in radiosensitive gels. Model-based iterative reconstruction may improve optical CT image quality and contribute to a wider use of optical CT in clinical gel dosimetry. METHODS: This algorithm was evaluated using experimental data acquired by a cone-beam optical CT system, as well as complementary numerical simulations. A fast GPU implementation of OSC-TV was used to achieve reconstruction times comparable to those of conventional filtered backprojection. Images obtained via OSC-TV were compared with the corresponding filtered backprojections. Spatial resolution and uniformity phantoms were scanned and respective reconstructions were subject to evaluation of the modulation transfer function, image uniformity, and accuracy. The artifacts due to refraction and total signal loss from opaque objects were also studied. RESULTS: The cone-beam optical CT data reconstructions showed that OSC-TV outperforms filtered backprojection in terms of image quality, thanks to a model-based simulation of the photon attenuation process. It was shown to significantly improve the image spatial resolution and reduce image noise. The accuracy of the estimation of linear attenuation coefficients remained similar to that obtained via filtered backprojection. Certain image artifacts due to opaque objects were reduced. Nevertheless, the common artifact due to the gel container walls could not be eliminated. CONCLUSIONS: The use of iterative reconstruction improves cone-beam optical CT image quality in many ways. The comparisons between OSC-TV and filtered backprojection presented in this paper demonstrate that OSC-TV can potentially improve the rendering of spatial features and reduce cone-beam optical CT artifacts.

GPU-accelerated regularized iterative reconstruction for few-view cone beam CT.

PURPOSE: The present work proposes an iterative reconstruction technique designed for x-ray transmission computed tomography (CT). The main objective is to provide a model-based solution to the cone-beam CT reconstruction problem, yielding accurate low-dose images via few-views acquisitions in clinically acceptable time frames. METHODS: The proposed technique combines a modified ordered subsets convex (OSC) algorithm and the total variation minimization (TV) regularization technique and is called OSC-TV. The number of subsets of each OSC iteration follows a reduction pattern in order to ensure the best performance of the regularization method. Considering the high computational cost of the algorithm, it is implemented on a graphics processing unit, using parallelization to accelerate computations. RESULTS: The reconstructions were performed on computer-simulated as well as human pelvic cone-beam CT projection data and image quality was assessed. In terms of convergence and image quality, OSC-TV performs well in reconstruction of low-dose cone-beam CT data obtained via a few-view acquisition protocol. It compares favorably to the few-view TV-regularized projections onto convex sets (POCS-TV) algorithm. It also appears to be a viable alternative to full-dataset filtered backprojection. Execution times are of 1-2 min and are compatible with the typical clinical workflow for nonreal-time applications. CONCLUSIONS: Considering the image quality and execution times, this method may be useful for reconstruction of low-dose clinical acquisitions. It may be of particular benefit to patients who undergo multiple acquisitions by reducing the overall imaging radiation dose and associated risks.