A megavoltage scatter correction technique for cone-beam CT images acquired during VMAT delivery

Concurrent kilovoltage CBCT imaging and megavoltage beam delivery: suppression of cross-scatter with 2D antiscatter grids and grid-based scatter sampling

Objective. The concept of using kilovoltage (kV) and megavoltage (MV) beams concurrently has potential applications in cone beam computed tomography (CBCT) guided radiation therapy, such as single breath hold scans, metal artifact reduction, and simultaneous imaging during MV treatment delivery. However, MV cross-scatter generated during MV beam delivery degrades CBCT image quality. To address this, a 2D antiscatter grid and a cross-scatter correction method were investigated in the context of high dose MV treatment delivery.Approach. A 3D printed, tungsten 2D antiscatter grid prototype was utilized in kV CBCT scans to reduce MV cross-scatter fluence during concurrent MV beam delivery. Remaining cross-scatter in projections was corrected by using the 2D grid itself as a cross-scatter intensity sampling device, referred to as grid-based scatter sampling (GSS). To test this approach, kV CBCT acquisitions were performed while delivering 6 and 10 MV beams, mimicking high dose rate treatment delivery scenarios. kV and MV beam deliveries were not synchronized to eliminate MV beam delivery interruption. MV cross-scatter suppression performance of the proposed approach was evaluated in projections and CBCT images of phantoms.Main results. 2D grid reduced the intensity of MV cross-scatter in kV projections by a factor of 3 on the average, when compared to conventional antiscatter grid. Remaining cross scatter as measured by the GSS method was within 7% of measured reference intensity values, and subsequently corrected. CBCT image quality was improved substantially during concurrent kV-MV beam delivery. Median Hounsfield Unit (HU) inaccuracy was up to 182 HU without our methods, and it was reduced to a median 6.5 HU with our 2D grid and scatter correction approach. Our methods provided a factor of 2-6 improvement in contrast-to-noise ratio.Significance. This investigation demonstrates the utility of 2D antiscatter grids and grid-based scatter sampling in suppressing MV cross-scatter. Our approach successfully minimized the effects of MV cross-scatter in concurrent kV CBCT imaging and high dose MV treatment delivery scenarios. Hence, robust MV cross-scatter suppression is potentially feasible without MV beam delivery interruption or compromising kV image acquisition rate.

Megavoltage cross-scatter rejection and correction using 2D antiscatter grids in kilovoltage CBCT imaging

Simultaneous use of kilovoltage (kV) and megavoltage (MV) beams has numerous potential applications in cone beam computed tomography (CBCT)-guided radiotherapy, such as fast MV+kV CBCT for single breath-hold scan, tumor localization with kV CBCT imaging during MV therapy delivery, and metal artifact suppression. However, the introduction of MV beams results in a large MV-cross scatter fluence incident on the kV Flat Panel Detector (FPD), and thus, deteriorating the low contrast visualization and Hounsfield Unit (HU) accuracy. In this work, we introduced a novel and robust method for reducing the effects of MV cross scatter. First, we implemented a 2D antiscatter grid atop the detector which rejects a large section of MV cross scatter. This hardware-based approach, while effective, allows a fraction of MV cross scatter to be transmitted to the FPD, resulting in artifacts and degraded HU accuracy in CBCT images. We thus introduced a data correction step, which aimed to estimate and correct the remaining MV cross scatter. This approach, referred to as Grid-Based Scatter Sampling, utilized 2D antiscatter grid itself to measure and correct remaining MV cross scatter in projections. We investigated the performance of the proposed approach in experiments by simultaneously acquiring kV CBCT and delivering MV beams with a clinical linac. The results show that the proposed method can substantially reduce HU inaccuracy and increase contrast-to-noise ratio (CNR). Our method does not require synchronization of kV and MV beam pulses, reduction of kV frame acquisition rate, or MV dose rate, and therefore, it is more practical to implement in radiation therapy clinical setting.

Actual delivered dose calculation on intra-irradiation cone-beam computed tomography images: a phantom study

Cone-beam computed tomography (CBCT) images acquired during volumetric modulated arc therapy (VMAT; ii-CBCT) can be used to calculate actual delivered doses (ADDs). However, such ii-CBCT images are degraded by scattered megavoltage x-rays (MV-scatters). We aimed to evaluate the dose calculation accuracy of the MV-scatter uncorrected or corrected ii-CBCT images acquired during VMAT deliveries. For MV-scatter correction on concurrent kilovoltage projections (P _MVkV), projections consisting only of MV-scatters (P _MVS) were acquired under the same MV beam parameters and gantry angles and subtracted from P _MVkV (P _MVScorr). In addition, the projections by kilovoltage beams were acquired for reference (P _kV). The corresponding CBCT images were reconstructed using the Feldkamp-Davis-Kress algorithm (CBCT_MVkV, CBCT_MVScorr, and CBCT_kV as reference). A multi-energy phantom with rods of various relative electron densities (REDs) was used to generate a CBCT-number-to-RED conversion table. First, CBCT_kV was reconstructed. Then, the mean CBCT-numbers within each rod were extracted, and a reference table was generated. Concurrent kilovoltage imaging with various field sizes was also demonstrated, and CBCT_MVkV and CBCT_MVScorr were reconstructed. The extracted CBCT-numbers of each ii-CBCT image were converted into REDs using the reference table. Next, the absolute differences of RED between the ii-CBCT image and CBCT_kV were calculated. Ten VMAT plans using a 10 MV flattening-filter-free beam were used for concurrent imaging of an anthropomorphic torso phantom. Moreover, an iterative reconstruction algorithm (IRA) was used for CBCT_MVScorr. The plans were recalculated for the corresponding CBCT_MVkV, CBCT_MVScorr, CBCT_MVScorr+IRA, and CBCT_kV with the reference table. Finally, the doses were evaluated using 3D gamma analysis (1%/1 mm). The median difference ranges between CBCT_MVkV/CBCT_MVScorr and the reference values were -0.58 to -0.10/-0.03 to 0.00. The median gamma pass rates of the doses on CBCT_MVkV, CBCT_MVScorr, and CBCT_MVScorr+IRA to the rate on CBCT_kV were 70.4, 99.5, and 98.2%, respectively. CBCT_MVScorr were comparable with CBCT_kV for calculating the ADD from VMAT.

Quantification and correction of the scattered X-rays from a megavoltage photon beam to a linac-mounted kilovoltage imaging subsystem

Objective

To quantify and correct megavoltage (MV) scattered X-rays (MV-scatter) on an image acquired using a linac-mounted kilovoltage (kV) imaging subsystem.

Methods and materials

A linac-mounted flat-panel detector (FPD) was used to acquire an image containing MV-scatter by activating the FPD only during MV beam irradiation. 6-, 10-, and 15 MV with a flattening-filter (FF; 6X-FF, 10X-FF, 15X-FF), and 6- and 10 MV without an FF (6X-FFF, 10X-FFF) were used. The maps were acquired by changing one of the irradiation parameters while the others remained fixed. The mean pixel values of the MV-scatter were normalized to the 6X-FF reference condition (MV-scatter value). An MV-scatter database was constructed using these values. An MV-scatter correction experiment with one full arc image acquisition and two square field sizes (FSs) was conducted. Measurement- and estimation-based corrections were performed using the database. The image contrast was calculated at each angle.

Results

The MV-scatter increased with a larger FS and dose rate. The MV-scatter value factor varied substantially depending on the FPD position or collimator rotation. The median relative error ranges of the contrast for the image without, and with the measurement- and estimation-based correction were -10.9 to -2.9, and -1.5 to 4.8 and -7.4 to 2.6, respectively, for an FS of 10.0 × 10.0 cm².

Conclusions

The MV-scatter was strongly dependent on the FS, dose rate, and FPD position. The MV-scatter correction improved the image contrast.

Advances in knowledge

The MV-scatters on the TrueBeam linac kV imaging subsystem were quantified with various MV beam parameters, and strongly depended on the fieldsize, dose rate, and flat panel detector position. The MV-scatter correction using the constructed database improved the image quality.

Image quality evaluation of intra-irradiation cone-beam computed tomography acquired during one- and two-arc prostate volumetric-modulated arc therapy delivery: A phantom study

Purpose

To evaluate (a) the effects of megavoltage (MV)‐scatter on concurrent kilovoltage (kV) projections (P_MVkV) acquired during rotational delivery, and (b) the image quality of intra‐irradiation cone‐beam computed tomography (ii‐CBCT) images acquired during prostate volumetric‐modulated arc therapy (VMAT) delivery.

Methods

Experiment (1): P_MVkVs were acquired with various MV beam parameters using a cylindrical phantom: field size (FS), MV energy (6 or 15 MV), dose rate (DR), and gantry speed. The average pixel values were calculated in a region on each P_MVkV which were extracted at eight equally spaced gantry angles. Experiment (2): 11 one‐arc and seven two‐arc 15 MV prostate VMAT plans were used along with a pelvis phantom. One plan was selected from each of arc plans and its MV energy was changed to 6 MV. After P_MVkVs were acquired, projections consisting of MV‐scatter only (P_MVS) were acquired with closing kV blades and subtracted from P_MVkV (P_MVScorr). Projections by kV beams only were acquired (P_kV). The corresponding CBCT images were reconstructed (CBCT_MVkV, CBCT_MVScorr, and CBCT_kV). The root‐mean‐square errors (RMSEs) were calculated in prostate region and 3D gamma analysis was conducted, in which the CBCT‐number was used instead of doses between ii‐CBCT images and CBCT_kV (30 HU/1 mm).

Results

Experiment (1): The MV‐scatters were dependent on the FSs, MV energies, and DRs. Experiment (2): The median RMSEs for CBCT_MVScorr were decreased by 107.5 HU (1‐arc) and 42.9 HU (2‐arc) compared to those for CBCT_MVkV. The median GPRs for CBCT_MVScorr were 94.7% (1‐arc) and 93.4% (2‐arc), while those for CBCT_MVkV were 61.1% and 79.9%, respectively. GPRs for 6 MV plans were smaller than those for 15 MV plans.

Conclusions

The number of MV‐scatters increased with larger FSs and DRs, and smaller MV energy. The MV‐scatters were corrected on the CBCT_MVScorr regardless of the number of arcs.

Image quality of 4D in-treatment CBCT acquired during lung SBRT using FFF beam: a phantom study

Background

Rotational beam delivery enables concurrent acquisition of cone-beam CT (CBCT), thereby facilitating further geometric verification of patient setup during radiation treatment. However, it is challenging to acquire CBCT during stereotactic body radiation therapy (SBRT) using flattening-filter free X-ray beams, in which a high radiation dose is delivered. This study presents quantitative evaluation results of the image quality in four-dimensional (4D) in-treatment CBCT acquired during SBRT delivery.

Methods

The impact of megavoltage (MV) scatter and acquisition parameters on the image quality was evaluated using Catphan 503 and XSight lung tracking phantoms. The in-treatment CBCT images of the phantoms were acquired while delivering 16 SBRT plans. The uniformity, contrast, and contrast-to-noise ratio (CNR) of the in-treatment CBCT images were calculated and compared to those of CBCT images acquired without SBRT delivery. Furthermore, the localizing accuracy of the moving target in the XSight lung phantom was evaluated for 10 respiratory phases.

Results

The CNR of the 3D-reconstucted Catphan CBCT images was reduced from 6.3 to 2.6 due to the effect of MV treatment scatter. Both for the Catphan and XSight phantoms, the CBCT image quality was affected by the tube current and monitor units (MUs) of the treatment plan. The lung target in the XSight tracking phantom was most visible for extreme phases; the mean CNRs of the lung target in the in-treatment CBCT images (with 40 mA tube current) across the SBRT plans were 3.2 for the end-of-exhalation phase and 3.0 for the end-of-inhalation phase. The lung target was localized with sub-millimeter accuracy for the extreme respiratory phases.

Conclusions

Full-arc acquisition with an increased tube current (e.g. 40 mA) is recommended to compensate for degradation in the CBCT image quality due to unflattened MV beam scatter. Acquiring in-treatment CBCT with a high-MU treatment beam is also suggested to improve the resulting CBCT image quality.

Artificial intelligence in radiotherapy: a technological review

Radiation therapy (RT) is widely used to treat cancer. Technological advances in RT have occurred in the past 30 years. These advances, such as three-dimensional image guidance, intensity modulation, and robotics, created challenges and opportunities for the next breakthrough, in which artificial intelligence (AI) will possibly play important roles. AI will replace certain repetitive and labor-intensive tasks and improve the accuracy and consistency of others, particularly those with increased complexity because of technological advances. The improvement in efficiency and consistency is important to manage the increasing cancer patient burden to the society. Furthermore, AI may provide new functionalities that facilitate satisfactory RT. The functionalities include superior images for real-time intervention and adaptive and personalized RT. AI may effectively synthesize and analyze big data for such purposes. This review describes the RT workflow and identifies areas, including imaging, treatment planning, quality assurance, and outcome prediction, that benefit from AI. This review primarily focuses on deep-learning techniques, although conventional machine-learning techniques are also mentioned.

Intrafraction 4D-cone beam CT acquired during volumetric arc radiotherapy delivery: kV parameter optimization and 4D motion accuracy for lung stereotactic body radiotherapy (SBRT) patients

Purpose

Elekta XVI 5.0 allows for four‐dimensional cone beam computed tomography (4D CBCT) image acquisition during treatment delivery to monitor intrafraction motion. These images can have poorer image quality due to undersampling of kV projections and treatment beam MV scatter effects. We determine if a universal intrafraction preset can be used for stereotactic body radiotherapy (SBRT) lung patients and validate the accuracy of target motion characterized by XVI intrafraction 4D CBCT.

Methods

The most critical parameter within the intrafraction preset is the nominal AcquisitionInterval, which controls kV imaging acquisition frequency. An optimal value was determined by maximizing the kV frame number acquired up to 1000 frames, typical of pretreatment 4D CBCT. CIRS motion phantom intrafraction phase images for 16 SBRT beams were obtained. Mean target position, time‐weighted standard deviation, and amplitude for these images as well as target motion for three SBRT lung patients were compared to respective pretreatment 4D CBCTs. Evaluation of intrafraction 4D CBCT reconstruction revealed inclusion of MV only images acquired to remove MV scatter effects. A workaround to remove these images was developed.

Results

AcquisitionInterval of 0.1°/frame was optimal. The number of kV frames acquired was 567–1116 and showed strong linear correlation with beam monitor unit (MUs). Phantom target motion accuracy was excellent with average differences in target position, standard deviation and amplitude range of ≤0.5 mm. Target tracking for SBRT patients also showed good agreement. Evaluation of phase sorting wave forms showed that inclusion of MV only images significantly impacts intrafraction image reconstruction for patients and use of workaround is necessary.

Conclusions

A universal intrafraction imaging preset can be used safely for SBRT lung patients. The number of kV projections with MV delivery parameters varies; however images with fewer kV projections still provided accurate target position information. Impact of the reconstruction workaround was significant and is mandated for all 4D CBCT intrafraction imaging performed at our institution.

A Bayesian approach for three-dimensional markerless tumor tracking using kV imaging during lung radiotherapy

The ability to monitor tumor motion without implanted markers can potentially enable broad access to more accurate and precise lung radiotherapy. A major challenge is that kilovoltage (kV) imaging based methods are rarely able to continuously track the tumor due to the inferior tumor visibility on 2D kV images. Another challenge is the estimation of 3D tumor position based on only 2D imaging information. The aim of this work is to address both challenges by proposing a Bayesian approach for markerless tumor tracking for the first time. The proposed approach adopts the framework of the extended Kalman filter, which combines a prediction and measurement steps to make the optimal tumor position update. For each imaging frame, the tumor position is first predicted by a respiratory-correlated model. The 2D tumor position on the kV image is then measured by template matching. Finally, the prediction and 2D measurement are combined based on the 3D distribution of tumor positions in the past 10 s and the estimated uncertainty of template matching. To investigate the clinical feasibility of the proposed method, a total of 13 lung cancer patient datasets were used for retrospective validation, including 11 cone-beam CT scan pairs and two stereotactic ablative body radiotherapy cases. The ground truths for tumor motion were generated from the the 3D trajectories of implanted markers or beacons. The mean, standard deviation, and 95th percentile of the 3D tracking error were found to range from 1.6-2.9 mm, 0.6-1.5 mm, and 2.6-5.8 mm, respectively. Markerless tumor tracking always resulted in smaller errors compared to the standard of care. The improvement was the most pronounced in the superior-inferior (SI) direction, with up to 9.5 mm reduction in the 95th-percentile SI error for patients with  >10 mm 5th-to-95th percentile SI tumor motion. The percentage of errors with 3D magnitude  <5 mm was 96.5% for markerless tumor tracking and 84.1% for the standard of care. The feasibility of 3D markerless tumor tracking has been demonstrated on realistic clinical scenarios for the first time. The clinical implementation of the proposed method will enable more accurate and precise lung radiotherapy using existing hardware and workflow. Future work is focused on the clinical and real-time implementation of this method.

Evaluating the image quality of cone beam CT acquired during rotational delivery

OBJECTIVE

The aim of this work was to evaluate the quality of kilovoltage (kV) cone beam CT (CBCT) images acquired during arc delivery.

METHODS

Arc plans were delivered on a Catphan(®) 600 phantom (The Phantom Laboratory Inc., Salem, NY), and kV CBCT images were acquired during the treatment. The megavoltage (MV) scatter effect on kV CBCT image quality was evaluated using parameters such as Hounsfield unit (HU) accuracy, spatial resolution, contrast-to-noise ratio (CNR) and spatial non-uniformity (SNU). These CBCT images were compared with reference scans acquired with the same acquisition parameters without MV "beam on". This evaluation was carried out for different photon beams (6 and 15 MV), arc types (half vs full arc), static field sizes (10 × 10 and 25 × 25 cm(2)) and source-to-imager distances (SID) (150 and 170 cm).

RESULTS AND CONCLUSION

HU accuracy, CNR and SNU were considerably affected by MV scatter, and this effect was increased with increasing field size and decreasing photon energy, whereas the spatial resolution was almost unchanged. The MV scatter effect was observed to be more for full-rotation arc delivery than for half-arc delivery. In addition, increasing the SID resulted in decreased MV scatter effect and improved the image quality.

ADVANCES IN KNOWLEDGE

Nowadays, volumetric modulated arc therapy (VMAT) is increasingly used in clinics, and this arc therapy enables us to acquire CBCT imaging simultaneously. But, the main issue of concurrent imaging is the "MV scatter" effect on CBCT imaging. This study aims to experimentally quantify the effect of MV scatter on CBCT image quality.

A moving blocker-based strategy for simultaneous megavoltage and kilovoltage scatter correction in cone-beam computed tomography image acquired during volumetric modulated arc therapy

PURPOSE

To evaluate a moving blocker-based approach in estimating and correcting megavoltage (MV) and kilovoltage (kV) scatter contamination in kV cone-beam computed tomography (CBCT) acquired during volumetric modulated arc therapy (VMAT).

METHODS AND MATERIALS

During the concurrent CBCT/VMAT acquisition, a physical attenuator (i.e., "blocker") consisting of equally spaced lead strips was mounted and moved constantly between the CBCT source and patient. Both kV and MV scatter signals were estimated from the blocked region of the imaging panel, and interpolated into the unblocked region. A scatter corrected CBCT was then reconstructed from the unblocked projections after scatter subtraction using an iterative image reconstruction algorithm based on constraint optimization. Experimental studies were performed on a Catphan® phantom and an anthropomorphic pelvis phantom to demonstrate the feasibility of using a moving blocker for kV-MV scatter correction.

RESULTS

Scatter induced cupping artifacts were substantially reduced in the moving blocker corrected CBCT images. Quantitatively, the root mean square error of Hounsfield units (HU) in seven density inserts of the Catphan phantom was reduced from 395 to 40.

CONCLUSIONS

The proposed moving blocker strategy greatly improves the image quality of CBCT acquired with concurrent VMAT by reducing the kV-MV scatter induced HU inaccuracy and cupping artifacts.

Digital tomosynthesis (DTS) for verification of target position in early stage lung cancer patients

PURPOSE

The ability to verify intrafraction tumor position is clinically useful for hypofractionated treatments. Short arc kV digital tomosynthesis (DTS) could facilitate more frequent target verification. The authors used DTS combined with triangulation to determine the mean temporal position of small-volume lung tumor targets treated with stereotactic radiotherapy. DTS registration results were benchmarked against online clinical localization using registration between free-breathing cone-beam computed tomography (CBCT) and the average intensity projection (AvIP) of the planning 4DCT.

METHODS

In this retrospective study, 76 sets of kV-projection images from online CBCT scans of 13 patients were used to generate DTS image slices (CB-DTS) with nonclinical research software (DTS Toolkit, Varian Medical Systems). Three-dimensional tumor motion was 1.3-4 mm in six patients and 6.1-25.4 mm in seven patients on 4DCT (significant difference in the mean of the groups, P < 0.01). The 4DCT AvIP was used to digitally reconstruct the Reference-DTS. DTS registration and DTS registration combined with triangulation were investigated. Progressive shortening of total DTS arc lengths from 95° to 35° around 0° gantry position was evaluated for different scenarios: DTS registration using the entire arc; DTS registration plus triangulation using two nonoverlapping arcs; and for 55° and 45° total gantry rotation, DTS registration plus triangulation using two overlapping arcs. Finally, DTS registration plus triangulation performed at eight gantry angles, each separated by 45° was evaluated using full fan kV projection data for one patient with an immobile tumor and five patients with mobile tumors.

RESULTS

For DTS registration alone, shortening arc length did not influence accuracy in X- and Y-directions, but in Z-direction, mean deviations from online CBCT localization systematically increased for shorter arc length (P < 0.05). For example, using a 95° arc mean DTS-CBCT difference was 0.8 mm (1 SD = 0.6 mm) and for a 35° arc the mean was 2.4 mm (1 SD = 1.7 mm). DTS plus triangulation using nonoverlapping-arcs increased accuracy in Z-direction for tested arc lengths ≤55° (P < 0.01). Overlapping arcs increased accuracy in Y-direction for tumors with motion >4 mm (P < 0.02) but increased Z-direction accuracy was only observed with 55° total gantry rotation. The 95th percentile deviations with this overlapping technique in X-, Y-, and Z-directions were 1.3, 2.0, and 2.5 mm, respectively. For the five patients with mobile tumors where DTS + triangulation was performed with 45° intervals, the pooled deviation from online CBCT correction showed, for X-, Y-, and Z-directions, mean of 1.1 mm, standard deviations (SD) of 0.9, 1.0, and 0.9 mm, respectively. The mean + 2 SD was <3 mm for each direction.

CONCLUSIONS

Short-arc DTS verification of time averaged lung tumor position is feasible using free-breathing kV projection data and the AvIP of the 4DCT as a reference. Observed differences between DTS and online CBCT registration with AvIP were ≤3 mm (mean + 2 SD), however, the increased temporal resolution of DTS + triangulation also identified short period deviations from the average target position on the CBCT. Short-arc DTS appears promising for intrafraction tumor position monitoring during stereotactic lung radiotherapy delivered with a rotational technique.