McSART: an iterative model-based, motion-compensated SART algorithm for CBCT reconstruction

PNMC: Four-dimensional conebeam CT reconstruction combining prior network and motion compensation

Four-dimensional conebeam computed tomography (4D CBCT) is an efficient technique to overcome motion artifacts caused by organ motion during breathing. 4D CBCT reconstruction in a single scan usually divides projections into different groups of sparsely sampled data based on the respiratory phases. The reconstructed images within each group present poor image quality due to the limited number of projections. To improve the image quality of 4D CBCT in a single scan, we propose a novel reconstruction scheme that combines prior knowledge with motion compensation. We apply the reconstructed images of the full projections within a single routine as prior knowledge, providing structural information for the network to enhance the restoration structure. The prior network (PN-Net) is proposed to extract features of prior knowledge and fuse them with the sparsely sampled data using an attention mechanism. The prior knowledge guides the reconstruction process to restore the approximate organ structure and alleviates severe streaking artifacts. The deformation vector field (DVF) extracted using deformable image registration among different phases is then applied in the motion-compensated ordered-subset simultaneous algebraic reconstruction algorithm to generate 4D CBCT images. Proposed method has been evaluated using simulated and clinical datasets and has shown promising results by comparative experiment. Compared with previous methods, our approach exhibits significant improvements across various evaluation metrics.

Faster and lower dose imaging: evaluating adaptive, constant gantry velocity and angular separation in fast low-dose 4D cone beam CT imaging

BACKGROUND

The adoption of four-dimensional cone beam computed tomography (4DCBCT) for image-guided lung cancer radiotherapy is increasing, especially for hypofractionated treatments. However, the drawbacks of 4DCBCT include long scan times (∼240 s), inconsistent image quality, higher imaging dose than necessary, and streaking artifacts. With the emergence of linear accelerators that can acquire 4DCBCT scans in a short period of time (9.2 s) there is a need to examine the impact that these very fast gantry rotations have on 4DCBCT image quality.

PURPOSE

This study investigates the impact of gantry velocity and angular separation between x-ray projections on image quality and its implication for fast low-dose 4DCBCT with emerging systems, such as the Varian Halcyon that provide fast gantry rotation and imaging. Large and uneven angular separation between x-ray projections is known to reduce 4DCBCT image quality through increased streaking artifacts. However, it is not known when angular separation starts degrading image quality. The study assesses the impact of constant and adaptive gantry velocity and determines the level when angular gaps impair image quality using state-of-the-art reconstruction methods.

METHODS

This study considers fast low-dose 4DCBCT acquisitions (60-80 s, 200-projection scans). To assess the impact of adaptive gantry rotations, the angular position of x-ray projections from adaptive 4DCBCT acquisitions from a 30-patient clinical trial were analyzed (referred to as patient angular gaps). To assess the impact of angular gaps, variable and static angular gaps (20°, 30°, 40°) were introduced into evenly separated 200 projections (ideal angular separation). To simulate fast gantry rotations, which are on emerging linacs, constant gantry velocity acquisitions (9.2 s, 60 s, 120 s, 240 s) were simulated by sampling x-ray projections at constant intervals using the patient breathing traces from the ADAPT clinical trial (ACTRN12618001440213). The 4D Extended Cardiac-Torso (XCAT) digital phantom was used to simulate projections to remove patient-specific image quality variables. Image reconstruction was performed using Feldkamp-Davis-Kress (FDK), McKinnon-Bates (MKB), and Motion-Compensated-MKB (MCMKB) algorithms. Image quality was assessed using Structural Similarity-Index-Measure (SSIM), Contrast-to-Noise-Ratio (CNR), Signal-to-Noise-Ratio (SNR), Tissue-Interface-Width-Diaphragm (TIW-D), and Tissue-Interface-Width-Tumor (TIW-T).

RESULTS

Patient angular gaps and variable angular gap reconstructions produced similar results to ideal angular separation reconstructions, while static angular gap reconstructions produced lower image quality metrics. For MCMKB-reconstructions, average patient angular gaps produced SSIM-0.98, CNR-13.6, SNR-34.8, TIW-D-1.5 mm, and TIW-T-2.0 mm, static angular gap 40° produced SSIM-0.92, CNR-6.8, SNR-6.7, TIW-D-5.7 mm, and TIW-T-5.9 mm and ideal produced SSIM-1.00, CNR-13.6, SNR-34.8, TIW-D-1.5 mm, and TIW-T-2.0 mm. All constant gantry velocity reconstructions produced lower image quality metrics than ideal angular separation reconstructions regardless of the acquisition time. Motion compensated reconstruction (MCMKB) produced the highest contrast images with low streaking artifacts.

CONCLUSION

Very fast 4DCBCT scans can be acquired provided that the entire scan range is adaptively sampled, and motion-compensated reconstruction is performed. Importantly, the angular separation between x-ray projections within each individual respiratory bin had minimal effect on the image quality of fast low-dose 4DCBCT imaging. The results will assist the development of future 4DCBCT acquisition protocols that can now be achieved in very short time frames with emerging linear accelerators.

Surrogate-driven respiratory motion model for projection-resolved motion estimation and motion compensated cone-beam CT reconstruction from unsorted projection data

Objective.As the most common solution to motion artefact for cone-beam CT (CBCT) in radiotherapy, 4DCBCT suffers from long acquisition time and phase sorting error. This issue could be addressed if the motion at each projection could be known, which is a severely ill-posed problem. This study aims to obtain the motion at each time point and motion-free image simultaneously from unsorted projection data of a standard 3DCBCT scan.Approach.Respiration surrogate signals were extracted by the Intensity Analysis method. A general framework was then deployed to fit a surrogate-driven motion model that characterized the relation between the motion and surrogate signals at each time point. Motion model fitting and motion compensated reconstruction were alternatively and iteratively performed. Stochastic subset gradient based method was used to significantly reduce the computation time. The performance of our method was comprehensively evaluated through digital phantom simulation and also validated on clinical scans from six patients.Results.For digital phantom experiments, motion models fitted with ground-truth or extracted surrogate signals both achieved a much lower motion estimation error and higher image quality, compared with non motion-compensated results.For the public SPARE Challenge datasets, more clear lung tissues and less blurry diaphragm could be seen in the motion compensated reconstruction, comparable to the benchmark 4DCBCT images but with a higher temporal resolution. Similar results were observed for two real clinical 3DCBCT scans.Significance.The motion compensated reconstructions and motion models produced by our method will have direct clinical benefit by providing more accurate estimates of the delivered dose and ultimately facilitating more accurate radiotherapy treatments for lung cancer patients.

Four-Dimensional Cone Beam CT Imaging Using a Single Routine Scan via Deep Learning

A novel method is proposed to obtain four-dimensional (4D) cone-beam computed tomography (CBCT) images from a routine scan in patients with upper abdominal cancer. The projections are sorted according to the location of the lung diaphragm before being reconstructed to phase-sorted data. A multiscale-discriminator generative adversarial network (MSD-GAN) is proposed to alleviate the severe streaking artifacts in the original images. The MSD-GAN is trained using simulated CBCT datasets from patient planning CT images. The enhanced images are further used to estimate the deformable vector field (DVF) among breathing phases using a deformable image registration method. The estimated DVF is then applied in the motion-compensated ordered-subset simultaneous algebraic reconstruction approach to generate 4D CBCT images. The proposed MSD-GAN is compared with U-Net on the performance of image enhancement. Results show that the proposed method significantly outperforms the total variation regularization-based iterative reconstruction approach and the method using only MSD-GAN to enhance original phase-sorted images in simulation and patient studies on 4D reconstruction quality. The MSD-GAN also shows higher accuracy than the U-Net. The proposed method enables a practical way for 4D-CBCT imaging from a single routine scan in upper abdominal cancer treatment including liver and pancreatic tumors.

Reducing 4DCBCT imaging dose and time: exploring the limits of adaptive acquisition and motion compensated reconstruction

This study investigates the dose and time limits of adaptive 4DCBCT acquisitions (adaptive-acquisition) compared with current conventional 4DCBCT acquisition (conventional-acquisition). We investigate adaptive-acquisitions as low as 60 projections (∼25 s scan, 6 projections per respiratory phase) in conjunction with emerging image reconstruction methods. 4DCBCT images from 20 patients recruited into the adaptive CT acquisition for personalized thoracic imaging clinical study (NCT04070586) were resampled to simulate faster and lower imaging dose acquisitions. All acquisitions were reconstructed using Feldkamp–Davis–Kress (FDK), McKinnon–Bates (MKB), motion compensated FDK (MCFDK), motion compensated MKB (MCMKB) and simultaneous motion estimation and image reconstruction (SMEIR) algorithms. All reconstructions were compared against conventional-acquisition 4DFDK-reconstruction using Structural SIMilarity Index (SSIM), signal-to-noise ratio (SNR), contrast-to-noise-ratio (CNR), tissue interface sharpness diaphragm (TIS-D), tissue interface sharpness tumor (TIS-T) and center of mass trajectory (COMT) for difference in diaphragm and tumor motion. All reconstruction methods using 110-projection adaptive-acquisition (11 projections per respiratory phase) had a SSIM of greater than 0.92 relative to conventional-acquisition 4DFDK-reconstruction. Relative to conventional-acquisition 4DFDK-reconstruction, 110-projection adaptive-acquisition MCFDK-reconstructions images had 60% higher SNR, 10% higher CNR, 30% higher TIS-T and 45% higher TIS-D on average. The 110-projection adaptive-acquisition SMEIR-reconstruction images had 123% higher SNR, 90% higher CNR, 96% higher TIS-T and 60% higher TIS-D on average. The difference in diaphragm and tumor motion compared to conventional-acquisition 4DFDK-reconstruction was within submillimeter accuracy for all acquisition reconstruction methods. Adaptive-acquisitions resulted in faster scans with lower imaging dose and equivalent or improved image quality compared to conventional-acquisition. Adaptive-acquisition with motion compensated-reconstruction enabled scans with as low as 110 projections to deliver acceptable image quality. This translates into 92% lower imaging dose and 80% less scan time than conventional-acquisition.

Validation of a CT-based motion model with in-situ fluoroscopy for lung surface deformation estimation

Many surrogate-based motion models (SMMs), proposed to guide motion management in radiotherapy, are constructed by correlating motion of an external surrogate and internal anatomy during CT-simulation. Changes in this correlation define model break down. We validate a methodology that incorporates fluoroscopic images (FL) acquired during treatment for SMM construction and update. Under a prospective IRB, 4DCT scans, VisionRT surfaces, and orthogonal FLs were collected from five lung cancer patients. VisionRT surfaces and two FL time-series were acquired pre- and post-treatment. A simulated annealing optimization scheme was used to estimate optimal lung deformations by maximizing the mutual information between digitally reconstructed radiographs (DRRs) of the SMM-estimated 3D images and FLs. Our SMM used partial-least-regression and was trained using the optimal deformations and VisionRT surfaces from the first breathing-cycle. SMM performance was evaluated using the mutual information score between reference FLs and the corresponding SMM or phase-assigned 4DCT DRRs. The Hausdorff distance for contoured landmarks was used to evaluate target position estimation error. For four out of five patients, two principal components approximated lung surface deformations with submillimeter accuracy. Analysis of the mutual information score between more than 4,000 pairs of FL and DRR demonstrated that our model led to more similarity between the FL and DRR images compared to 4DCT and DRR images from a model based on an a priori correlation model. Our SMM consistently displayed lower mean and 95^th percentile Hausdorff distances. For one patient, 95^th percentile Hausdorff distance was reduced by 11mm. Patient-averaged reductions in mean and 95^th percentile Hausdorff distances were 3.6mm and 7mm for right-lung, and 3.1mm and 4mm for left-lung targets. FL data were used to evaluate model performance and investigate the feasibility of model update. Despite variability in breathing, use of post-treatment FL preserved model fidelity and consistently outperformed 4DCT for position estimation.

Total variation with modified group sparsity for CT reconstruction under low SNR

BACKGROUND AND OBJECTIVE

Since the stair artifacts may affect non-destructive testing (NDT) and diagnosis in the later stage, an applicable model is desperately needed, which can deal with the stair artifacts and preserve the edges. However, the classical total variation (TV) algorithm only considers the sparsity of the gradient transformed image. The objective of this study is to introduce and test a new method based on group sparsity to address the low signal-to-noise ratio (SNR) problem.

METHODS

This study proposes a weighted total variation with overlapping group sparsity model. This model combines the Gaussian kernel and overlapping group sparsity into TV model denoted as GOGS-TV, which considers the structure sparsity of the image to be reconstructed to deal with the stair artifacts. On one hand, TV is the accepted commercial algorithm, and it can work well in many situations. On the other hand, the Gaussian kernel can associate the points around each pixel. Quantitative assessments are implemented to verify this merit.

RESULTS

Numerical simulations are performed to validate the presented method, compared with the classical simultaneous algebraic reconstruction technique (SART) and the state-of-the-art TV algorithm. It confirms the significantly improved SNR of the reconstruction images both in suppressing the noise and preserving the edges using new GOGS-TV model.

CONCLUSIONS

The proposed GOGS-TV model demonstrates its advantages to reduce stair artifacts especially in low SNR reconstruction because this new model considers both the sparsity of the gradient image and the structured sparsity. Meanwhile, the Gaussian kernel is utilized as a weighted factor that can be adapted to the global distribution.

Technical Note: Comparison of the internal target volume (ITV) contours and dose calculations on 4DCT, average CBCT, and 4DCBCT imaging for lung stereotactic body radiation therapy (SBRT)

PURPOSE

To investigate the differences between internal target volumes (ITVs) contoured on the simulation 4DCT and daily 4DCBCT images for lung cancer patients treated with stereotactic body radiotherapy (SBRT) and determine the dose delivered on 4D planning technique.

METHODS

For nine patients, 4DCBCTs were acquired before each fraction to assess tumor motion. An ITV was contoured on each phase of the 4DCBCT and a union of the 10 ITVs was used to create a composite ITV. Another ITV was drawn on the average 3DCBCT (avgCBCT) to compare with current clinical practice. The Dice coefficient, Hausdorff distance, and center of mass (COM) were averaged over four fractions to compare the ITVs contoured on the 4DCT, avgCBCT, and 4DCBCT for each patient. Planning was done on the average CT, and using the online registration, plans were calculated on each phase of the 4DCBCT and on the avgCBCT. Plan dose calculations were tested by measuring ion chamber dose in the CIRS lung phantom.

RESULTS

The Dice coefficients were similar for all three comparisons: avgCBCT-to-4DCBCT (0.7 ± 0.1), 4DCT-to-avgCBCT (0.7 ± 0.1), and 4DCT-to-4DCBCT (0.7 ± 0.1); while the mean COM differences were also comparable (2.6 ± 2.2mm, 2.3 ± 1.4mm, and 3.1 ± 1.1mm, respectively). The Hausdorff distances for the comparisons with 4DCBCT (8.2 ± 2.9mm and 8.1 ± 3.2mm) were larger than the comparison without (6.5 ± 2.5mm). The differences in ITV D95% between the treatment plan and avgCBCT calculations were 4.3 ± 3.0% and -0.5 ± 4.6%, between treatment plan and 4DCBCT plans, respectively, while the ITV V100% coverages were 99.0 ± 1.9% and 93.1 ± 8.0% for avgCBCT and 4DCBCT, respectively.

CONCLUSION

There is great potential for 4DCBCT to evaluate the extent of tumor motion before treatment, but image quality challenges the clinician to consistently delineate lung target volumes.

Rigid and non-rigid motion artifact reduction in X-ray CT using attention module

Motion artifacts are a major factor that can degrade the diagnostic performance of computed tomography (CT) images. In particular, the motion artifacts become considerably more severe when an imaging system requires a long scan time such as in dental CT or cone-beam CT (CBCT) applications, where patients generate rigid and non-rigid motions. To address this problem, we proposed a new real-time technique for motion artifacts reduction that utilizes a deep residual network with an attention module. Our attention module was designed to increase the model capacity by amplifying or attenuating the residual features according to their importance. We trained and evaluated the network by creating four benchmark datasets with rigid motions or with both rigid and non-rigid motions under a step-and-shoot fan-beam CT (FBCT) or a CBCT. Each dataset provided a set of motion-corrupted CT images and their ground-truth CT image pairs. The strong modeling power of the proposed network model allowed us to successfully handle motion artifacts from the two CT systems under various motion scenarios in real-time. As a result, the proposed model demonstrated clear performance benefits. In addition, we compared our model with Wasserstein generative adversarial network (WGAN)-based models and a deep residual network (DRN)-based model, which are one of the most powerful techniques for CT denoising and natural RGB image deblurring, respectively. Based on the extensive analysis and comparisons using four benchmark datasets, we confirmed that our model outperformed the aforementioned competitors. Our benchmark datasets and implementation code are available at https://github.com/youngjun-ko/ct_mar_attention.

A modern review of the uncertainties in volumetric imaging of respiratory-induced target motion in lung radiotherapy

Radiotherapy has become a critical component for the treatment of all stages and types of lung cancer, often times being the primary gateway to a cure. However, given that radiation can cause harmful side effects depending on how much surrounding healthy tissue is exposed, treatment of the lung can be particularly challenging due to the presence of moving targets. Careful implementation of every step in the radiotherapy process is absolutely integral for attaining optimal clinical outcomes. With the advent and now widespread use of stereotactic body radiation therapy (SBRT), where extremely large doses are delivered, accurate, and precise dose targeting is especially vital to achieve an optimal risk to benefit ratio. This has largely become possible due to the rapid development of image-guided technology. Although imaging is critical to the success of radiotherapy, it can often be plagued with uncertainties due to respiratory-induced target motion. There has and continues to be an immense research effort aimed at acknowledging and addressing these uncertainties to further our abilities to more precisely target radiation treatment. Thus, the goal of this article is to provide a detailed review of the prevailing uncertainties that remain to be investigated across the different imaging modalities, as well as to highlight the more modern solutions to imaging motion and their role in addressing the current challenges.