A rebinned backprojection-filtration algorithm for image reconstruction in helical cone-beam CT

A weighted rebinned backprojection-filtration algorithm from partially beam-blocked data for a single-scan cone-beam CT with hybrid type scatter correction

PURPOSE

Scatter contamination constitutes a dominant source of degradation of image quality in cone-beam computed tomography (CBCT). We have recently developed an analytic image reconstruction method with a scatter correction capability from the partially blocked cone-beam data out of a single scan. Despite its easy implementation and its computational efficiency, the developed method may result in additional image artifacts for a large cone angle geometry due to data inconsistency. To improve the image quality at a large cone angle, we propose a weighted rebinned backprojection-filtration (wrBPF) algorithm in conjunction with a hybrid type scatter correction approach.

METHODS

The proposed method uses a beam-blocker array that provides partial data for scatter correction and image reconstruction and that only blocks the beam within a limited cone angle. This design allows a chance to keep the image quality at larger cone angles by use of data redundancy since the projection data corresponding to larger cone angles are not blocked. However, the scatter correction would not be straightforward. In order to correct for the scatter in the projections at larger cone angles, we propose a novel scatter correction method combining a measurement-based and a convolution-based method. We first estimated the scatter signal using a measurement-based method in the partially beam-blocked regions, and then optimized the fitting parameters of a convolution-kernel that can be used for scatter correction in the projections at larger cone angles. For image reconstruction, we developed a wrBPF with butterfly filtering. We have conducted an experimental study to validate the proposed algorithm for image reconstruction and scatter correction.

RESULTS

The experimental results revealed that the developed reconstruction method makes full use of the benefits of partial beam-blocking for scatter correction and image reconstruction and at the same time enhances image quality at larger cone angles by use of an optimized convolution-based scatter correction.

CONCLUSIONS

The proposed method that enjoys the advantages of both measurement-based and convolution-based methods for scatter correction has successfully demonstrated its capability of reconstructing accurate images out of a single scan in circular CBCT.

Half-Fan-Based Intensity-Weighted Region-of-Interest Imaging for Low-Dose Cone-Beam CT in Image-Guided Radiation Therapy

Objectives

With the increased use of computed tomography (CT) in clinics, dose reduction is the most important feature people seek when considering new CT techniques or applications. We developed an intensity-weighted region-of-interest (IWROI) imaging method in an exact half-fan geometry to reduce the imaging radiation dose to patients in cone-beam CT (CBCT) for image-guided radiation therapy (IGRT). While dose reduction is highly desirable, preserving the high-quality images of the ROI is also important for target localization in IGRT.

Methods

An intensity-weighting (IW) filter made of copper was mounted in place of a bowtie filter on the X-ray tube unit of an on-board imager (OBI) system such that the filter can substantially reduce radiation exposure to the outer ROI. In addition to mounting the IW filter, the lead-blade collimation of the OBI was adjusted to produce an exact half-fan scanning geometry for a further reduction of the radiation dose. The chord-based rebinned backprojection-filtration (BPF) algorithm in circular CBCT was implemented for image reconstruction, and a humanoid pelvis phantom was used for the IWROI imaging experiment.

Results

The IWROI image of the phantom was successfully reconstructed after beam-quality correction, and it was registered to the reference image within an acceptable level of tolerance. Dosimetric measurements revealed that the dose is reduced by approximately 61% in the inner ROI and by 73% in the outer ROI compared to the conventional bowtie filter-based half-fan scan.

Conclusions

The IWROI method substantially reduces the imaging radiation dose and provides reconstructed images with an acceptable level of quality for patient setup and target localization. The proposed half-fan-based IWROI imaging technique can add a valuable option to CBCT in IGRT applications.

Reconstructing cone-beam CT with spatially varying qualities for adaptive radiotherapy: a proof-of-principle study

With the aim of maximally reducing imaging dose while meeting requirements for adaptive radiation therapy (ART), we propose in this paper a new cone beam CT (CBCT) acquisition and reconstruction method that delivers images with a low noise level inside a region of interest (ROI) and a relatively high noise level outside the ROI. The acquired projection images include two groups: densely sampled projections at a low exposure with a large field of view (FOV) and sparsely sampled projections at a high exposure with a small FOV corresponding to the ROI. A new algorithm combining the conventional filtered back-projection algorithm and the tight-frame iterative reconstruction algorithm is also designed to reconstruct the CBCT based on these projection data. We have validated our method on a simulated head-and-neck (HN) patient case, a semi-real experiment conducted on a HN cancer patient under a full-fan scan mode, as well as a Catphan phantom under a half-fan scan mode. Relative root-mean-square errors (RRMSEs) of less than 3% for the entire image and ~1% within the ROI compared to the ground truth have been observed. These numbers demonstrate the ability of our proposed method to reconstruct high-quality images inside the ROI. As for the part outside ROI, although the images are relatively noisy, it can still provide sufficient information for radiation dose calculations in ART. Dose distributions calculated on our CBCT image and on a standard CBCT image are in agreement, with a mean relative difference of 0.082% inside the ROI and 0.038% outside the ROI. Compared with the standard clinical CBCT scheme, an imaging dose reduction of approximately 3-6 times inside the ROI was achieved, as well as an 8 times outside the ROI. Regarding computational efficiency, it takes 1-3 min to reconstruct a CBCT image depending on the number of projections used. These results indicate that the proposed method has the potential for application in ART.

Region of Interest Imaging for a General Trajectory with the Rebinned BPF Algorithm

The back-projection-filtration (BPF) algorithm has been applied to image reconstruction for cone-beam configurations with general source trajectories. The BPF algorithm can reconstruct 3-D region-of-interest (ROI) images from data containing truncations. However, like many other existing algorithms for cone-beam configurations, the BPF algorithm involves a back-projection with a spatially varying weighting factor, which can result in the non-uniform noise levels in reconstructed images and increased computation time. In this work, we propose a BPF algorithm to eliminate the spatially varying weighting factor by using a rebinned geometry for a general scanning trajectory. This proposed BPF algorithm has an improved noise property, while retaining the advantages of the original BPF algorithm such as minimum data requirement.

Accurate helical cone-beam CT reconstruction with redundant data

We present a new image reconstruction algorithm for helical cone-beam computed tomography (CT). This algorithm is designed for data collected at or near maximum pitch, and provides a theoretically exact and stable reconstruction while beneficially using all measured data. The main operations involved are a differentiated backprojection and a finite-support Hilbert transform inversion. These operations are applied onto M-lines, and the beneficial use of all measured data is gained from averaging three volumes reconstructed each with a different choice of M-lines. The technique is overall similar to that presented by one of the authors in a previous publication, but operates volume-wise, instead of voxel-wise, which yields a significantly more efficient reconstruction procedure. The algorithm is presented in detail. Also, preliminary results from computer-simulated data are provided to demonstrate the numerical stability of the algorithm, the beneficial use of redundant data and the ability to process data collected with an angular flying focal spot.