Impact of cone-beam computed tomography artifacts on dose calculation accuracy for lung cancer.

BACKGROUND: To implement image-guided adaptive radiotherapy (IGART), many studies investigated dose calculations on cone-beam computed tomography (CBCT). A high HU accuracy is crucial for a high dose calculation accuracy and many imaging sites showed satisfactory results. It has been shown that the dose calculation accuracy for lung cancer lags behind. PURPOSE: To examine why the dose calculation accuracy for lung is insufficient, the relative effects of the field-of-view (FOV), breathing motion, and scatter on dose calculation accuracy were studied. METHODS: A framework was built to simulate CBCT scans for lung cancer patients by forward projecting repeat CT (rCT) scans for two scan geometries: small (SFOV) and medium FOV (MFOV). Breathing motion was modeled by applying a 4D deformation vector field to the mid-position rCT. Scatter was modeled by Monte-Carlo simulations with/without an anti-scatter grid (ASG). Simulated projections were reconstructed using filtered back-projection with/without scatter correction. In case of the SFOV, the CBCT images were patched with the planning CT scan in axial direction. The treatment plan was recalculated on the rCT and simulated CBCT. The mean Hounsfield unit (HU) difference (ΔHUmean), the structural similarity index measure (SSIM), and γ metrics were calculated for the CBCT datasets of various imaging settings. RESULTS: The differences in HU, SSIM and dose calculation accuracy for CBCTs with and without breathing motion were negligible (mean ΔHUmean = 6.4 vs. 13.7, mean SSIM = 0.941 vs. 0.957, mean γ (ref = MFOV) = 0.75). The SFOV resulted in a lower HU (mean ΔHUmean = -9.2 vs. 13.7) and SSIM (mean SSIM = 0.912 vs. 0.957), and therefore in dose differences compared to the MFOV (mean γ = 1.22). Scatter led to considerable discrepancies in all metrics. Adding only the ASG improved the results more than only applying a scatter correction algorithm. Combining ASG and scatter correction algorithm resulted in an even higher dose calculation accuracy. CONCLUSIONS: Scatter and FOV are the main contributors to dose inaccuracies and motion has only a minor effect on dose calculation accuracy. Therefore, utilizing an appropriate scatter correction and FOV is important to achieve sufficient dose calculation accuracy to facilitate IGART for lung.

Image quality of dual-energy cone-beam CT with total nuclear variation regularization.

Despite the improvements in image quality of cone beam computed tomography (CBCT) scans, application remains limited to patient positioning. In this study, we propose to improve image quality by dual energy (DE) imaging and iterative reconstruction using least squares fitting with total variation (TV) regularization. The generalization of TV called total nuclear variation (TNV) was used to generate DE images. We acquired single energy (SE) and DE scans of an image quality phantom (IQP) and of an anthropomorphic human male phantom (HMP). The DE scans were dual arc acquisitions of 70 kV and 130 kV with a variable dose partitioning between low energy (LE) and high energy (HE) arcs. To investigate potential benefits from a larger spectral separation between LE and HE, DE scans with an additional 2 mm copper beam filtration in the HE arc were acquired for the IQP. The DE TNV scans were compared to SE scans reconstructed with FDK and iterative TV with varying parameters. The contrast-to-noise ratio (CNR), spatial frequency, and structural similarity (SSIM) were used as image quality metrics. Results showed largely improved image quality for DE TNV over FDK for both phantoms. DE TNV with the highest dose allocation in the LE arm yielded the highest CNR. Compared to SE TV, these DE TNV results had a slightly lower CNR with similar spatial resolution for the IQP. A decrease in the dose allocated to the LE arm improved the spatial resolution with a trade-off against CNR. For the HMP, DE TNV displayed a lower CNR and/or lower spatial resolution depending on the reconstruction parameters. Regarding the SSIM, DE TNV was superior to FDK and SE TV for both phantoms. The additional beam filtration for the IQP led to improved image quality in all metrics, surpassing the SE TV results in CNR and spatial resolution.

Improving linac integrated cone beam computed tomography image quality using tube current modulation.

PURPOSE: Linac integrated cone beam CT (CBCT) scanners have become widespread tool for image guidance in radiotherapy. The current implementation uses constant imaging fluence across all the projection angles, which leads to anisotropic noise properties and suboptimal image quality for noncircular symmetric objects. Tube current modulation (TCM) is widely used in conventional CT. The purpose of this work was to implement TCM on a linac integrated CBCT scanner and evaluate its impact on image quality under varying scatter conditions and scatter correction strategies. METHODS: We have implemented TCM on a nonclinical Elekta Versa HD linear accelerator with enhanced x-ray generator functionality including pulse width modulation. The pulse width was modulated using two Arduino programmable microcontrollers: one placed on the kV arm to measure the projection angle and the other connected to the kV generator control board to vary x-ray pulse width as function of gantry angle and precalculated transmission. An in-house developed phantom with a ratio of the left-right to anterior-posterior path length of 1.85:1 was scanned. Image quality was determined using the anisotropicity of the 2D noise power spectra (NPS) in the transverse plane and the contrast-to-noise ratio (CNR). In addition, to determine the impact of scatter on the applicability of the TCM method we have modified the generated scatter using three different collimators in the cranio-caudal direction as well as with and without an antiscatter grid (ASG). RESULTS: Application of the TCM led to 30-78% reduction of the angular anisotropicity of the NPS in the transverse plane. The amount of reduction depended on the scatter conditions, with lower values corresponding to higher scatter conditions. The same was true for the CNR: when scatter contribution was low (presence of an ASG or very aggressive collimation) the CNR was improved by about 30%, while in high scatter conditions the CNR was improved by about 12%. CONCLUSIONS: TCM has the potential to improve CBCT image quality, but this depends on the amount of detected x-ray scatter.

Technical Note: Long-term stability of Hounsfield unit calibration for cone beam computed tomography.

PURPOSE: The goal of this study was to investigate the stability of a phantom-based Hounsfield unit (HU) calibration for cone beam CT (CBCT). METHODS: Three consecutive scans of a large phantom configuration and a small phantom configuration were acquired and reconstructed with a uniform scatter correction method. The CBCT gray values of the phantom inserts were measured and the three values of each insert averaged. The linear calibration curve was determined and its slope and intercept were evaluated. This procedure was performed for three CBCT scanners (Elekta Synergy) over a period of 10 months with 1, 2, or 4 weeks between measurements. A dosimetric estimation of the HU fluctuations was carried out. RESULTS: The CBCT HUs were stable over time with only small variations in slope and intercept resulting in HU differences on the water level, that is, intercepts, of less than 7 HU (standard deviation). Therefore, the dosimetric influence of these HU differences was limited. The inter-scanner disparities (up to ~ 16 HU) were larger than the intra-scanner ones (up to ~ 7 HU). CONCLUSIONS: Stable HUs were observed over a period of 10 months. Due to the differences between the CBCT scanners, scanner-specific calibration curves are necessary.

Evaluating the impact of cone-beam computed tomography scatter mitigation strategies on radiotherapy dose calculation accuracy.

BACKGROUND AND PURPOSE: The scatter induced image quality degradation of cone-beam computed tomography (CBCT) prevents more advanced applications in radiotherapy. We evaluated the dose calculation accuracy on CBCT of various disease sites using different scatter mitigation strategies. MATERIALS AND METHODS: CBCT scans of two patient cohorts (C1, C2) were reconstructed using a uniform (USC) and an iterative scatter correction (ISC) method, combined with an anti-scatter grid (ASG). Head and neck (H&N), lung, pelvic region, and prostate patients were included. To achieve a high accuracy Hounsfield unit and physical density calibrations were performed. The dose distributions of the original treatment plans were analyzed with the γ evaluation method using criteria of 1%/2 mm using the planning CT as the reference. The investigated parameters were the mean γ (γmean), the points in agreement (Pγ≤1) and the 99th percentile (γ1%). RESULTS: Significant differences between USC and ISC in C1 were found for the lung and prostate, where the latter using the ISC produced the best results with medians of 0.38, 98%, and 1.1 for γmean, Pγ≤1 and γ1%, respectively. For C2 the ISC with ASG showed an improvement for all imaging sites. The lung demonstrated the largest relative increase in accuracy with improvements between 48% and 54% for the medians of γmean, Pγ≤1 and γ1%. CONCLUSIONS: The introduced method demonstrated high dosimetric accuracy for H&N, prostate and pelvic region if an ASG is applied. A significantly lower accuracy was seen for lung. The ISC yielded a higher robustness against scatter variations than the USC.

Reducing 4DCBCT imaging time and dose: the first implementation of variable gantry speed 4DCBCT on a linear accelerator.

Four dimensional cone beam computed tomography (4DCBCT) uses a constant gantry speed and imaging frequency that are independent of the patient's breathing rate. Using a technique called respiratory motion guided 4DCBCT (RMG-4DCBCT), we have previously demonstrated that by varying the gantry speed and imaging frequency, in response to changes in the patient's real-time respiratory signal, the imaging dose can be reduced by 50-70%. RMG-4DCBCT optimally computes a patient specific gantry trajectory to eliminate streaking artefacts and projection clustering that is inherent in 4DCBCT imaging. The gantry trajectory is continuously updated as projection data is acquired and the patient's breathing changes. The aim of this study was to realise RMG-4DCBCT for the first time on a linear accelerator. To change the gantry speed in real-time a potentiometer under microcontroller control was used to adjust the current supplied to an Elekta Synergy's gantry motor. A real-time feedback loop was developed on the microcontroller to modulate the gantry speed and projection acquisition in response to the real-time respiratory signal so that either 40, RMG-4DCBCT40, or 60, RMG-4DCBCT60, uniformly spaced projections were acquired in 10 phase bins. Images of the CIRS dynamic Thorax phantom were acquired with sinusoidal breathing periods ranging from 2 s to 8 s together with two breathing traces from lung cancer patients. Image quality was assessed using the contrast to noise ratio (CNR) and edge response width (ERW). For the average patient, with a 3.8 s breathing period, the imaging time and image dose were reduced by 37% and 70% respectively. Across all respiratory rates, RMG-4DCBCT40 had a CNR in the range of 6.5 to 7.5, and RMG-4DCBCT60 had a CNR between 8.7 and 9.7, indicating that RMG-4DCBCT allows consistent and controllable CNR. In comparison, the CNR for conventional 4DCBCT drops from 20.4 to 6.2 as the breathing rate increases from 2 s to 8 s. With RMG-4DCBCT, the ERW in the direction of motion of the imaging insert decreases from 2.1 mm to 1.1 mm as the breathing rate increases from 2 s to 8 s while for conventional 4DCBCT the ERW increases from 1.9 mm to 2.5 mm. Image quality can be controlled during 4DCBCT acquisition by varying the gantry speed and the projection acquisition in response to the patient's real-time respiratory signal. However, although the image sharpness, i.e. ERW, is improved with RMG-4DCBCT, the ERW depends on the patient's breathing rate and breathing regularity.

Optimal combination of anti-scatter grids and software correction for CBCT imaging.

PURPOSE: Cone beam computed tomography (CBCT) has been widely adopted in clinical practice for image-guided radiotherapy. Soft tissue contrast and Hounsfield units are impaired to the presence of scattered radiation. In our previous work, we proposed a high selectivity anti-scatter grid (ASG) as a possible solution to the problem. An alternative approach is the application of iterative scatter correction using deconvolution with scatter point spread function (PSF). The purpose of this work was to compare the performance of ASGs with different selectivity with and without the iterative and uniform scatter corrections in terms of CBCT image quality. A secondary objective of this study was to develop a novel measurement approach to measure the scatter point spread functions. METHODS: The scatter PSF was modeled as a sum of two bivariate Gaussian functions. The PSF parameters were estimated from a series of transmission measurements through polystyrene slabs of varying thickness with lead partial beam-blocker for three different ASG designs ranging from low (5.6), medium (9), and high (11) selectivity. The scatter correction scheme is based on iterative convolution of the current estimate of the primary with the scatter PSF until the root mean square deviation (RMSD) of the measured projection and the sum of the estimate of primary and scatter falls below a predefined threshold. The image quality was evaluated with the CIRS CBCT Image Quality and Electron Density phantom in a head and neck and pelvis configuration and the CIRS Virtual Male Human Patient. The image quality was quantified by the contrast-to-noise ratio (CNR) relative to the uncorrected scans and the root mean square deviation of the average gray values for different regions with respect to the nominal Hounsfield units and the mean difference of the reconstructed HU between the planning CT and CBCTs of the virtual human phantom. RESULTS: For the head and neck phantom, the CNR increased with more advanced scatter correction algorithm and the ASG selectivity, reaching 3.9, 3.7, 3.5, and 3.1 for the high, medium, light, and with no grid configuration, respectively, combined with the iterative software correction. The same is true for the pelvis phantom with CNR improvement reaching 1.5 for the heavy and medium grid, 1.3 for the light grid, and 1.1 on its own. The HU RMSD for the head and neck phantom was 22 HU, 13 HU, 12 HU, and 6 HU for iterative correction without the grid, with the light grid, medium grid and the heavy grid, respectively. For same correction strategies, the values for the pelvis phantom where 170, 120, 34, and 27 HU. The average difference with the PCT of the virtual human phantom was 59 ± 48 HU and 63 ± 59 HU with scans reconstructed with the iterative correction and two higher selectivity grids. Visual inspection revealed similar trends for a head-and-neck and prostate cancer patient. CONCLUSIONS: The best scatter mitigation strategy was found to be a combination of a grid with selectivity larger than 9, combined with iterative scatter estimation. None of the investigated grids required increasing the imaging dose. The PSF determined using proposed method leads to image quality improvements results for all but one of the investigated scenarios.

Mitigating differential baseline shifts in locally advanced lung cancer patients using an average anatomy model.

BACKGROUND AND PURPOSE: Differential baseline shifts between primary tumor and involved lymph nodes in locally advanced lung cancer patients compromise the accuracy of radiotherapy. The purpose of this study was to evaluate the performance of an average anatomy model (AAM) derived from repeat imaging and deformable registration to reduce these geometrical uncertainties. METHODS AND MATERIALS: An in-house implementation of a B-Spline deformable image registration (DIR) algorithm was first validated using three different validation approaches: (a) a circle method to test the consistency of the DIR, (b) fiducial marker target registration error, and (c) the recovery of a known deformation vector field (DVF). Subsequently, AAM was generated by first averaging five DVFs resulting from cone beam CT (CBCT) to planning CT (pCT) DIR and second by applying the inverse of the average DVF to the pCT. The proposed method was evaluated on 15 locally advanced lung cancer patients receiving daily motion compensated CBCT and a repeat CT (rCT) for adaptive radiotherapy. Reduction of systematic baseline shifts of the primary tumor were quantified for the fractions used to build the AAM as well as over the whole treatment and compared to the performance of the rCT. RESULTS: The deformable registration accuracy was ≤ 2 mm vector length for the first two validation methods and about 3 mm for the third method. The systematic baseline shifts over the five fractions prior to the rCT used to build the AAM reduced from 5.9 mm vector length relative to the pCT to 2.3 and 4.2 mm relative to the AAM and rCT, respectively. The overall systematic errors in the left-right, cranio-caudal, and anterior-posterior directions were [3.4, 3.8, 3.3] mm, [2.3, 2.9, 2.6] mm, and [2.3, 3.1, 2.7] mm for the pCT, AAM, and rCT, respectively. CONCLUSIONS: The AAM mitigates systematic errors occurring during treatment due to differential baseline shifts between the primary tumor and involved lymph nodes similar to (or even better than) rCT. The superior performance of the AAM in terms of the systematic error derived from the initial fractions indicates that further analysis of the optimum intervention time is required. This model has the potential to be used as an efficient and accurate alternative for rCT in adaptive radiotherapy of locally advanced lung cancer patients, obviating the need for rescanning and recontouring.

Clinical introduction of image lag correction for a cone beam CT system.

PURPOSE: Image lag in the flat-panel detector used for Linac integrated cone beam computed tomography (CBCT) has a degrading effect on CBCT image quality. The most prominent visible artifact is the presence of bright semicircular structure in the transverse view of the scans, known also as radar artifact. Several correction strategies have been proposed, but until now the clinical introduction of such corrections remains unreported. In November 2013, the authors have clinically implemented a previously proposed image lag correction on all of their machines at their main site in Amsterdam. The purpose of this study was to retrospectively evaluate the effect of the correction on the quality of CBCT images and evaluate the required calibration frequency. METHODS: Image lag was measured in five clinical CBCT systems (Elekta Synergy 4.6) using an in-house developed beam interrupting device that stops the x-ray beam midway through the data acquisition of an unattenuated beam for calibration. A triple exponential falling edge response was fitted to the measured data and used to correct image lag from projection images with an infinite response. This filter, including an extrapolation for saturated pixels, was incorporated in the authors' in-house developed clinical cbct reconstruction software. To investigate the short-term stability of the lag and associated parameters, a series of five image lag measurement over a period of three months was performed. For quantitative analysis, the authors have retrospectively selected ten patients treated in the pelvic region. The apparent contrast was quantified in polar coordinates for scans reconstructed using the parameters obtained from different dates with and without saturation handling. RESULTS: Visually, the radar artifact was minimal in scans reconstructed using image lag correction especially when saturation handling was used. In patient imaging, there was a significant reduction of the apparent contrast from 43 ± 16.7 to 15.5 ± 11.9 HU without the saturation handling and to 9.6 ± 12.1 HU with the saturation handling, depending on the date of the calibration. The image lag correction parameters were stable over a period of 3 months. The computational load was increased by approximately 10%, not endangering the fast in-line reconstruction. CONCLUSIONS: The lag correction was successfully implemented clinically and removed most image lag artifacts thus improving the image quality. Image lag correction parameters were stable for 3 months indicating low frequency of calibration requirements.

The first implementation of respiratory triggered 4DCBCT on a linear accelerator.

Four dimensional cone beam computed tomography (4DCBCT) is an image guidance strategy used for patient positioning in radiotherapy. In conventional implementations of 4DCBCT, a constant gantry speed and a constant projection pulse rate are used. Unfortunately, this leads to higher imaging doses than are necessary because a large number of redundant projections are acquired. In theoretical studies, we have previously demonstrated that by suppressing redundant projections the imaging dose can be reduced by 40-50% for a majority of patients with little reduction in image quality. The aim of this study was to experimentally realise the projection suppression technique, which we have called Respiratory Triggered 4DCBCT (RT-4DCBCT). A real-time control system was developed that takes the respiratory signal as input and computes whether to acquire, or suppress, the next projection trigger during 4DCBCT acquisition. The CIRS dynamic thorax phantom was programmed with a 2 cm peak-to-peak motion and periods ranging from 2 to 8 s. Image quality was assessed by computing the edge response width of a 3 cm imaging insert placed in the phantom as well as the signal to noise ratio of the phantoms tissue and the contrast to noise ratio between the phantoms lung and tissue. The standard deviation in the superior-inferior direction of the 3 cm imaging insert was used to assess intra-phase bin displacement variations with a higher standard deviation implying more motion blur. The 4DCBCT imaging dose was reduced by 8.6%, 41%, 54%, 70% and 77% for patients with 2, 3, 4, 6 and 8 s breathing periods respectively when compared to conventional 4DCBCT. The standard deviation of the intra-phase bin displacement variation of the 3 cm imaging insert was reduced by between 13% and 43% indicating a more consistent position for the projections within respiratory phases. For the 4 s breathing period, the edge response width was reduced by 39% (0.8 mm) with only a 6-7% decrease in the signal to noise and contrast to noise ratios. RT-4DCBCT has been experimentally realised and reduced to practice on a linear accelerator with a measurable imaging dose reductions over conventional 4DCBCT and little degradation in image quality.

Improved image quality of cone beam CT scans for radiotherapy image guidance using fiber-interspaced antiscatter grid.

PURPOSE: Medical linear accelerator mounted cone beam CT (CBCT) scanner provides useful soft tissue contrast for purposes of image guidance in radiotherapy. The presence of extensive scattered radiation has a negative effect on soft tissue visibility and uniformity of CBCT scans. Antiscatter grids (ASG) are used in the field of diagnostic radiography to mitigate the scatter. They usually do increase the contrast of the scan, but simultaneously increase the noise. Therefore, and considering other scatter mitigation mechanisms present in a CBCT scanner, the applicability of ASGs with aluminum interspacing for a wide range of imaging conditions has been inconclusive in previous studies. In recent years, grids using fiber interspacers have appeared, providing grids with higher scatter rejection while maintaining reasonable transmission of primary radiation. The purpose of this study was to evaluate the impact of one such grid on CBCT image quality. METHODS: The grid used (Philips Medical Systems) had ratio of 21:1, frequency 36 lp/cm, and nominal selectivity of 11.9. It was mounted on the kV flat panel detector of an Elekta Synergy linear accelerator and tested in a phantom and a clinical study. Due to the flex of the linac and presence of gridline artifacts an angle dependent gain correction algorithm was devised to mitigate resulting artifacts. Scan reconstruction was performed using XVI4.5 augmented with inhouse developed image lag correction and Hounsfield unit calibration. To determine the necessary parameters for Hounsfield unit calibration and software scatter correction parameters, the Catphan 600 (The Phantom Laboratory) phantom was used. Image quality parameters were evaluated using CIRS CBCT Image Quality and Electron Density Phantom (CIRS) in two different geometries: one modeling head and neck and other pelvic region. Phantoms were acquired with and without the grid and reconstructed with and without software correction which was adapted for the different acquisition scenarios. Parameters used in the phantom study were t(cup) for nonuniformity and contrast-to-noise ratio (CNR) for soft tissue visibility. Clinical scans were evaluated in an observer study in which four experienced radiotherapy technologists rated soft tissue visibility and uniformity of scans with and without the grid. RESULTS: The proposed angle dependent gain correction algorithm suppressed the visible ring artifacts. Grid had a beneficial impact on nonuniformity, contrast to noise ratio, and Hounsfield unit accuracy for both scanning geometries. The nonuniformity reduced by 90% for head sized object and 91% for pelvic-sized object. CNR improved compared to no corrections on average by a factor 2.8 for the head sized object, and 2.2 for the pelvic sized phantom. Grid outperformed software correction alone, but adding additional software correction to the grid was overall the best strategy. In the observer study, a significant improvement was found in both soft tissue visibility and nonuniformity of scans when grid is used. CONCLUSIONS: The evaluated fiber-interspaced grid improved the image quality of the CBCT system for broad range of imaging conditions. Clinical scans show significant improvement in soft tissue visibility and uniformity without the need to increase the imaging dose.