Multi-criterial patient positioning based on dose recalculation on scatter-corrected CBCT images

Evaluation of the geometric and dosimetric accuracies of deformable image registration of targets and critical organs in prostate CBCT-guided adaptive radiotherapy

PURPOSE

Kilovoltage cone beam computed tomography (kVCBCT)-guided adaptive radiation therapy (ART) uses daily deformed CT (dCT), which is generated automatically through deformable registration methods. These registration methods may perform poorly in reproducing volumes of the target organ, rectum, and bladder during treatment. We analyzed the registration errors between the daily kVCBCTs and corresponding dCTs for these organs using the default optical flow algorithm and two registration procedures. We validated the effectiveness of these registration methods in replicating the geometry for dose calculation on kVCBCT for ART.

METHODS

We evaluated three deformable image registration (DIR) methods to assess their registration accuracy and dose calculation effeciency in mapping target and critical organs. The DIR methods include (1) default intensity-based deformable registration, (2) hybrid deformable registration, and (3) a two-step deformable registration process. Each technique was applied to a computerized imaging reference system (CIRS) phantom (Model 062 M) and to five patients who received volumetric modulated arc therapy to the prostate. Registration accuracy was assessed using the 95% Hausdorff distance (HD95) and Dice similarity coefficient (DSC), and each method was compared with the intensity-based registration method. The improvement in the dCT image quality of the CIRS phantom and five patients was assessed by comparing dCT with kVCBCT. Image quality quantitative metrics for the phantom included the signal-to-noise ratio (SNR), uniformity, and contrast-to-noise ratio (CNR), whereas those for the patients included the mean absolute error (MAE), mean error, peak signal-to-noise ratio (PSNR), and structural similarity index measure (SSIM). To determine dose metric differences, we used a dose-volume histogram (DVH) and 3.0%/0.3 mm gamma analysis to compare planning computed tomography (pCT) and kVCBCT recalculations with restimulated CT images used as a reference.

RESULTS

The dCT images generated by the hybrid (dCTH) and two-step (dCTC) registration methods resulted in significant improvements compared to kVCBCT in the phantom model. Specifically, the SNR improved by 107% and 107.2%, the uniformity improved by 90% and 75%, and the CNR improved by 212.2% and 225.6 for dCTH and dCTC methods, respectively. For the patient images, the MAEs improved by 98% and 94%, the PSNRs improved by 16.3% and 22.9%, and the SSIMs improved by 1% and 1% in the dCTH and dCTC methods, respectively. For the geometric evaluation, only the two-step registration method improved registration accuracy. The dCTH method yielded an average HD95 of 12 mm and average DSC of 0.73, whereas dCTC yielded an average HD95 of 2.9 mm and average DSC of 0.902. The DVH showed that the dCTC-based dose calculations differed by <2% from the expected results for treatment targets and volumes of organs at risk. Additionally, gamma indices for dCTC-based treatment plans were >95% at all points, whereas they were <95% for kVCBCT-based treatment plans.

CONCLUSION

The two-step registration method outperforms the intensity-based and hybrid registration methods. While the hybrid and two-step-based methods improved the image quality of kVCBCT in a linear accelerator, only the two-step method improved the registration accuracy of the corresponding structures among the pCT and kVCBCT datasets. A two-step registration process is recommended for applying kVCBCT to ART, which achieves better registration accuracy for local and global image structures. This method appears to be beneficial for radiotherapy dose calculation in patients with pelvic cancer.

Image quality evaluation of a new high-performance ring-gantry cone-beam computed tomography imager

Objective. Newer cone-beam computed tomography (CBCT) imaging systems offer reconstruction algorithms including metal artifact reduction (MAR) and extended field-of-view (eFoV) techniques to improve image quality. In this study a new CBCT imager, the new Varian HyperSight CBCT, is compared to fan-beam CT and two CBCT imagers installed in a ring-gantry and C-arm linear accelerator, respectively.Approach. The image quality was assessed for HyperSight CBCT which uses new hardware, including a large-size flat panel detector, and improved image reconstruction algorithms. The decrease of metal artifacts was quantified (structural similarity index measure (SSIM) and root-mean-squared error (RMSE)) when applying MAR reconstruction and iterative reconstruction for a dental and spine region using a head-and-neck phantom. The geometry and CT number accuracy of the eFoV reconstruction was evaluated outside the standard field-of-view (sFoV) on a large 3D-printed chest phantom. Phantom size dependency of CT numbers was evaluated on three cylindrical phantoms of increasing diameter. Signal-to-noise and contrast-to-noise were quantified on an abdominal phantom.Main results. In phantoms with streak artifacts, MAR showed comparable results for HyperSight CBCT and CT, with MAR increasing the SSIM (0.97-0.99) and decreasing the RMSE (62-55 HU) compared to iterative reconstruction without MAR. In addition, HyperSight CBCT showed better geometrical accuracy in the eFoV than CT (Jaccard Conformity Index increase of 0.02-0.03). However, the CT number accuracy outside the sFoV was lower than for CT. The maximum CT number variation between different phantom sizes was lower for the HyperSight CBCT imager (∼100 HU) compared to the two other CBCT imagers (∼200 HU), but not fully comparable to CT (∼50 HU).Significance. This study demonstrated the imaging performance of the new HyperSight CBCT imager and the potential of applying this CBCT system in more advanced scenarios by comparing the quality against fan-beam CT.

Evaluation of a cone-beam computed tomography system calibrated for accurate radiotherapy dose calculation

Background and purpose

Dose calculation on cone-beam computed tomography (CBCT) images has been less accurate than on computed tomography (CT) images due to lower image quality and discrepancies in CT numbers for CBCT. As increasing interest arises in offline and online re-planning, dose calculation accuracy was evaluated for a novel CBCT imager integrated into a ring gantry treatment machine.

Materials and methods

The new CBCT system allowed fast image acquisition (5.9 s) by using new hardware, including a large-size flat panel detector, and incorporated image-processing algorithms with iterative reconstruction techniques, leading to accurate CT numbers allowing dose calculation. In this study, CBCT- and CT-based dose calculations were compared based on three anthropomorphic phantoms, after CBCT-to-mass-density calibration was performed. Six plans were created on the CT scans covering various target locations and complexities, followed by CBCT to CT registrations, copying of contours, and re-calculation of the plans on the CBCT scans. Dose-volume histogram metrics for target volumes and organs-at-risk (OARs) were evaluated, and global gamma analyses were performed.

Results

Target coverage differences were consistently below 1.2 %, demonstrating the agreement between CT and re-calculated CBCT dose distributions. Differences in Dmean for OARs were below 0.5 Gy for all plans, except for three OARs, which were below 0.8 Gy (<1.1 %). All plans had a 3 %/1mm gamma pass rate > 97 %.

Conclusions

This study demonstrated comparable results between dose calculations performed on CBCT and CT acquisitions. The new CBCT system with enhanced image quality and CT number accuracy opens possibilities for off-line and on-line re-planning.

Unanticipated Radiation Replanning for Stage III Non-small Cell Lung Cancer

Purpose

The purpose of this study was to identify factors associated with unanticipated radiation therapy (RT) replanning in stage III non-small cell lung cancer (NSCLC).

Methods and Materials

Patients from a single institution with newly diagnosed stage III NSCLC treated with radical RT from January 1, 2016, to December 31, 2019, were retrospectively analyzed. The frequency and reasons for replanning were determined. Logistic regression analysis was used to identify factors associated with replanning.

Results

Of 144 patients included in this study, 11% (n = 16) required replanning after the start of RT. The reason for replanning in these 16 patients was changes in the target detected by cone beam computed tomography (shift in 10 patients, shrinkage in 5 patients, and growth in 1 patient). Larger planning target volume (primary and nodal) was statistically predictive of replanning (odds ratio, 2.5; 95% CI, 1.2-5.4; P = .02). The actuarial median overall survival was 33.3 months (95% CI, 10.3-43.9) for the 16 patients who were replanned and 36.3 months (95% CI, 27.4-66.5) for the remaining 128 patients (P = .96). The median time to local recurrence was 25.0 months (95% CI, 10.3-41.3) for those patients who underwent replanning, which was similar to those patients who did not undergo replanning (19.5 months; 95% CI, 11.8-23.2; P = .28).

Conclusions

In this study, 11% of patients treated with radical RT for NSCLC required replanning due to changes in the target detected by cone beam computed tomography. A larger planning target volume predicts a higher likelihood of requiring adaptive RT. Overall survival and local control were similar between patients who were replanned compared with those who were not replanned.

Scatter correction of 4D cone beam computed tomography to detect dosimetric effects due to anatomical changes in proton therapy for lung cancer

BACKGROUND

The treatment of moving tumor entities is expected to have superior clinical outcomes, using image-guided adaptive intensity-modulated proton therapy (IMPT).

PURPOSE

For 21 lung cancer patients, IMPT dose calculations were performed on scatter-corrected 4D cone beam CTs (4DCBCTcor ) to evaluate their potential for triggering treatment adaptation. Additional dose calculations were performed on corresponding planning 4DCTs and day-of-treatment 4D virtual CTs (4DvCTs).

METHODS

A 4DCBCT correction workflow, previously validated on a phantom, generates 4DvCT (CT-to-CBCT deformable registration) and 4DCBCTcor images (projection-based correction using 4DvCT as a prior) with 10 phase bins, using day-of-treatment free-breathing CBCT projections and planning 4DCT images as input. Using a research planning system, robust IMPT plans administering eight fractions of 7.5 Gy were created on a free-breathing planning CT (pCT) contoured by a physician. The internal target volume (ITV) was overridden with muscle tissue. Robustness settings for range and setup uncertainties were 3% and 6 mm, and a Monte Carlo dose engine was used. On every phase of planning 4DCT, day-of-treatment 4DvCT, and 4DCBCTcor , the dose was recalculated. For evaluation, image analysis as well as dose analysis were performed using mean error (ME) and mean absolute error (MAE) analysis, dose-volume histogram (DVH) parameters, and 2%/2-mm gamma pass rate analysis. Action levels (1.6% ITV D98 and 90% gamma pass rate) based on our previous phantom validation study were set to determine which patients had a loss of dosimetric coverage.

RESULTS

Quality enhancements of 4DvCT and 4DCBCTcor over 4DCBCT were observed. ITV D98% and bronchi D2% had its largest agreement for 4DCBCTcor -4DvCT, and the largest gamma pass rates (>94%, median 98%) were found for 4DCBCTcor -4DvCT. Deviations were larger and gamma pass rates were smaller for 4DvCT-4DCT and 4DCBCTcor -4DCT. For five patients, deviations were larger than the action levels, suggesting substantial anatomical changes between pCT and CBCT projections acquisition.

CONCLUSIONS

This retrospective study shows the feasibility of daily proton dose calculation on 4DCBCTcor for lung tumor patients. The applied method is of clinical interest as it generates up-to-date in-room images, accounting for breathing motion and anatomical changes. This information could be used to trigger replanning.

Technical note: Characterization of novel iterative reconstructed cone beam CT images for dose tracking and adaptive radiotherapy on L-shape linacs

BACKGROUND

Cone-beam computed tomography (CBCT) allows for patient setup and positioning, and potentially dose verification or adaptive replanning prior to each treatment delivery. Poor CBCT image quality due to scatter artifacts and patient motion has been a major limiting factor. A new image reconstruction algorithm was recently clinically implemented for improving image quality through iterative reconstruction (iCBCT).

PURPOSE

This study aims to characterize iCBCT image quality, establish image value (HU)-to-relative electron density (RED) calibration curves for dose calculation, and assess the dosimetric accuracy for different anatomical sites.

MATERIAL AND METHODS

Both conventional CBCT and iCBCT scans were acquired from a Varian TrueBeam On-Board Imager system. A Catphan 604 phantom was scanned to compare image quality between the traditional Feldkamp-Davis-Kress (FDK) and novel iterative reconstruction techniques. Computerized Imaging Reference Systems (CIRS) electron density phantom was used to construct site-specific HU-RED curves corresponding to various scan settings. The CIRS Dynamic Thorax phantom, Rando pelvis phantom, and BrainLab head phantom were used for assessing dosimetric accuracy calculated on iCBCT images, compared to that on traditional FDK-based CBCT images. All phantoms were scanned on a computed tomography (CT) to obtain baseline HU values for comparison.

RESULTS

Test results obtained from Catphan showed statistically significant improvement with iCBCT, compared to FDK CBCT. Average HU differences from the baseline CT values were improved to within ±30 HU for iCBCT, compared to FDK CBCT for phantom studies. Dose calculated on iCBCT for both phantoms and patient cases directly using baseline HU-RED calibration from CT showed 0.5%-2.0% accuracy from the baseline dose calculated on CT, which is comparable to doses calculated using site-specific HU-RED calibration curves.

CONCLUSION

iCBCT provides improved image quality, improved HU accuracy compared to CT baseline, and has potential to provide online dose verification as part of the adaptive radiotherapy workflow directly using the baseline HU-RED calibration curve from CT.

Real-time dose-guidance in radiotherapy: Proof of principle

PURPOSE

The outcome of radiotherapy is a direct consequence of the dose delivered to the patient. Yet image-guidance and motion management to date focus on geometrical considerations as a practical surrogate for dose. Here, we propose real-time dose-guidance realized through continuous motion-including dose reconstructions and demonstrate this new concept in simulated liver stereotactic body radiotherapy (SBRT) delivery.

MATERIALS AND METHODS

During simulated liver SBRT delivery, in-house developed software performed real-time motion-including reconstruction of the tumor dose delivered so far and continuously predicted the remaining fraction tumor dose. The total fraction dose was estimated as the sum of the delivered and predicted doses, both with and without the emulated couch correction that maximized the predicted final CTV D95% (minimum dose to 95% of the clinical target volume). Dose-guided treatments were simulated for 15 liver SBRT patients previously treated with tumor motion monitoring, using both sinusoidal tumor motion and the actual patient-measured motion. A dose-guided couch correction was triggered if it improved the predicted final CTV D95% with 3, 4 or 5 %-points. The final CTV D95% of the dose-guidance strategy was compared with simulated treatments using geometry guided couch corrections (Wilcoxon signed-rank test).

RESULTS

Compared to geometry guidance, dose-guided couch corrections improved the median CTV D95% with 0.5-1.5 %-points (p < 0.01) for sinusoidal motions and with 0.9 %-points (p < 0.01, 3 %-points trigger threshold), 0.5 %-points (p = 0.03, 4 %-points threshold) and 1.2 %-points (p = 0.09, 5 %-points threshold) for patient-measured tumor motion.

CONCLUSION

Real-time dose-guidance was proposed and demonstrated to be superior to geometrical adaptation in liver SBRT simulations.

Validation of proton dose calculation on scatter corrected 4D cone beam computed tomography using a porcine lung phantom

Proton therapy treatment for lungs remains challenging as images enabling the detection of inter- and intra-fractional motion, which could be used for proton dose adaptation, are not readily available. 4D computed tomography (4DCT) provides high image quality but is rarely available in-room, while in-room 4D cone beam computed tomography (4DCBCT) suffers from image quality limitations stemming mostly from scatter detection. This study investigated the feasibility of using virtual 4D computed tomography (4DvCT) as a prior for a phase-per-phase scatter correction algorithm yielding a 4D scatter corrected cone beam computed tomography image (4DCBCT_cor), which can be used for proton dose calculation. 4DCT and 4DCBCT scans of a porcine lung phantom, which generated reproducible ventilation, were acquired with matching breathing patterns. Diffeomorphic Morphons, a deformable image registration algorithm, was used to register the mid-position 4DCT to the mid-position 4DCBCT and yield a 4DvCT. The 4DCBCT was reconstructed using motion-aware reconstruction based on spatial and temporal regularization (MA-ROOSTER). Successively for each phase, digitally reconstructed radiographs of the 4DvCT, simulated without scatter, were exploited to correct scatter in the corresponding CBCT projections. The 4DCBCT_corwas then reconstructed with MA-ROOSTER using the corrected CBCT projections and the same settings and deformation vector fields as those already used for reconstructing the 4DCBCT. The 4DCBCT_corand the 4DvCT were evaluated phase-by-phase, performing proton dose calculations and comparison to those of a ground truth 4DCT by means of dose-volume-histograms (DVH) and gamma pass-rates (PR). For accumulated doses, DVH parameters deviated by at most 1.7% in the 4DvCT and 2.0% in the 4DCBCT_corcase. The gamma PR for a (2%, 2 mm) criterion with 10% threshold were at least 93.2% (4DvCT) and 94.2% (4DCBCT_cor), respectively. The 4DCBCT_cortechnique enabled accurate proton dose calculation, which indicates the potential for applicability to clinical 4DCBCT scans.

Comparison of CBCT conversion methods for dose calculation in the head and neck region

The purpose of this study was to compare different methods of CBCT conversion respect to dose calculation accuracy. Twelve head and neck cancer patients treated with VMAT using simultaneous integrated boost technique were selected for the study. For each patient a planning CT (pCT), a control. CT acquired in the fourth week of treatment and a CBCT scan acquired on the closest day with the control CT were used. In order to re-calculate dose directly on CBCT image sets, a population based approach (CBCTPop) and a Histogram Matching (HM) approach based on rigid (CBCTHM-R) and deformable registration (CBCTHM-D) were used. Additionally, virtual CTs (vCTs) were generated using two deformable image registration algorithms (CTELX and CTANC) of the planning CT to the CBCT by using two different deformable image registration (DIR) algorithms. The corresponding control CTs were selected as ground truth and dose distributions on CBCT were analyzed using 3D global gamma index analysis applying a threshold of 10% with respect to the prescribed dose. Using the 2%/2mm gamma criterion, the results were 89.9%(±8.3%), 94.1%(±5.0%), 94.3%(±5.7%), 96.1%(±3.9%), 93.4%(±6.3%) for the CBCTPop, CBCTHM-R, CBCTHM-D, CTELX and CTANC, respectively. On average, the HM and DIR techniques showed a higher accuracy compared to the population based approach, but Kruskal-Wallis test did not show significant difference among the investigated dose calculation techniques assuming p<0.05. More sophisticated CBCT dose calculation methods seem to improve the dose calculation accuracy, but statistical significance remains to be demonstrated.

"Dose of the day" based on cone beam computed tomography and deformable image registration for lung cancer radiotherapy

Purpose

Adaptive radiotherapy (ART) has potential to reduce toxicity and facilitate safe dose escalation. Dose calculations with the planning CT deformed to cone beam CT (CBCT) have shown promise for estimating the “dose of the day”. The purpose of this study is to investigate the “dose of the day” calculation accuracy based on CBCT and deformable image registration (DIR) for lung cancer radiotherapy.

Methods

A total of 12 lung cancer patients were identified, for which daily CBCT imaging was performed for treatment positioning. A re‐planning CT (rCT) was acquired after 20 Gy for all patients. A virtual CT (vCT) was created by deforming initial planning CT (pCT) to the simulated CBCT that was generated from deforming CBCT to rCT acquired on the same day. Treatment beams from the initial plan were copied to the vCT and rCT for dose calculation. Dosimetric agreement between vCT‐based and rCT‐based accumulated doses was evaluated using the Bland‐Altman analysis.

Results

Mean differences in dose‐volume metrics between vCT and rCT were smaller than 1.5%, and most discrepancies fell within the range of ± 5% for the target volume, lung, esophagus, and heart. For spinal cord D_max, a large mean difference of −5.55% was observed, which was largely attributed to very limited CBCT image quality (e.g., truncation artifacts).

Conclusion

This study demonstrated a reasonable agreement in dose‐volume metrics between dose accumulation based on vCT and rCT, with the exception for cases with poor CBCT image quality. These findings suggest potential utility of vCT for providing a reasonable estimate of the “dose of the day”, and thus facilitating the process of ART for lung cancer.

Comparing Unet training with three different datasets to correct CBCT images for prostate radiotherapy dose calculations

Image intensity correction is crucial to enable cone beam computed tomography (CBCT) based radiotherapy dose calculations. This study evaluated three different deep learning based correction methods using a U-shaped convolutional neural network architecture (Unet) in terms of their photon and proton dose calculation accuracy. CT and CBCT imaging data of 42 prostate cancer patients were included. For target ground truth data generation, a CBCT correction method based on CT to CBCT deformable image registration (DIR) was used. The method yields a deformed CT called (i) virtual CT (vCT) which is used to generate (ii) corrected CBCT projections allowing the reconstruction of (iii) a final corrected CBCT image. The single Unet architecture was trained using these three different datasets: (Unet1) raw and corrected CBCT projections, (Unet2) raw CBCT and vCT image slices and (Unet3) raw and reference corrected CBCT image slices. Volumetric arc therapy (VMAT) and proton pencil beam scanning (PBS) single field uniform dose (SFUD) plans were optimized on the reference corrected image and recalculated on the obtained Unet-corrected CBCT images. The mean error (ME) and mean absolute error (MAE) for Unet1/2/3 were [Formula: see text] Hounsfield units (HU) and [Formula: see text] HU. The 1% dose difference pass rates were better than 98.4% for VMAT for 8 test patients not seen during training, with little difference between Unets. Gamma evaluation results were even better. For protons a gamma evaluation was employed to account for small range shifts, and [Formula: see text] mm pass rates for Unet1/2/3 were better than [Formula: see text] and 91%. A 3 mm range difference threshold was established. Only for Unet3 the 5th and 95th percentiles of the range difference distributions over all fields, test patients and dose profiles were within this threshold. A single Unet architecture was successfully trained using both CBCT projections and CBCT image slices. Since the results of the other Unets were poorer than Unet3, we conclude that training using corrected CBCT image slices as target data is optimal for PBS SFUD proton dose calculations, while for VMAT all Unets provided sufficient accuracy.

Dose-guided patient positioning in proton radiotherapy using multicriteria-optimization

Proton radiotherapy (PT) requires accurate target alignment before each treatment fraction, ideally utilizing 3D in-room X-ray computed tomography (CT) imaging. Typically, the optimal patient position is determined based on anatomical landmarks or implanted markers. In the presence of non-rigid anatomical changes, however, the planning scenario cannot be exactly reproduced and positioning should rather aim at finding the optimal position in terms of the actually applied dose. In this work, dose-guided patient alignment, implemented as multicriterial optimization (MCO) problem, was investigated in the scope of intensity-modulated and double-scattered PT (IMPT and DSPT) for the first time. A method for automatically determining the optimal patient position with respect to pre-defined clinical goals was implemented. Linear dose interpolation was used to access a continuous space of potential patient shifts. Fourteen head and neck (H&N) and eight prostate cancer patients with up to five repeated CTs were included. Dose interpolation accuracy was evaluated and the potential dosimetric advantages of dose-guided over bony-anatomy-based patient alignment investigated by comparison of clinically relevant target and organ-at-risk (OAR) dose-volume histogram (DVH) parameters. Dose interpolation was found sufficiently accurate with average pass-rates of 90% and 99% for an exemplary H&N and prostate patient, respectively, using a 2% dose-difference criterion. Compared to bony-anatomy-based alignment, the main impact of automated MCO-based dose-guided positioning was a reduced dose to the serial OARs (spinal cord and brain stem) for the H&N cohort. For the prostate cohort, under-dosage of the target structures could be efficiently diminished. Limitations of dose-guided positioning were mainly found in reducing target over-dosage due to weight loss for H&N patients, which might require adaptation of the treatment plan. Since labor-intense online quality-assurance is not required for dose-guided patient positioning, it might, nevertheless, be considered an interesting alternative to full online re-planning for initially mitigating the effects of anatomical changes.