Cone beam CT based dose calculation in the thorax region

Daily CBCT-based dose calculations for enhancing the safety of dose-escalation in lung cancer radiotherapy

PURPOSE

Dose-escalation in lung cancer comes with a high risk of severe toxicity. This study aimed to calculate the delivered dose in a Scandinavian phase-III dose-escalation trial.

METHODS

The delivered dose was evaluated for 21 locally-advanced non-small cell lung cancer (LA-NSCLC) patients treated as part of the NARLAL2 dose-escalation trial. The patients were randomized between standard and escalated heterogeneous dose-delivery. Both treatment plans were created and approved before randomization. Daily cone-beam CT (CBCT) for patient positioning, and adaptive radiotherapy were mandatory. Standard and escalated plans, including adaptive re-plans, were recalculated on each daily CBCT and accumulated on the planning CT for each patient. Dose to the clinical target volume (CTV), organs at risk (OAR), and the effects of plan adaptions were evaluated for the accumulated dose and on each treated fraction scaled to full treatment.

RESULTS

For the standard treatment, plan adaptations reduced the number of patients with CTV-T underdosage from six to one, and the total number of fractions with CTV-T underdosage from 161 to 56; while for the escalated treatment, the number of patients was reduced from five to zero and number of fractions from 81 to 11. For dose-escalation, three patients had fractions exceeding trial constraints for heart, bronchi, or esophagus, and one had an accumulated heart dose above the constraints.

CONCLUSION

Dose-escalation for LA-NSCLC patients, using daily image guidance and adaptive radiotherapy, is dosimetrically safe for the majority of patients. Dose calculation on daily CBCTs is an efficient tool to monitor target coverage and OAR doses.

Technical note: Phantom-based evaluation of CBCT dose calculation accuracy for use in adaptive radiotherapy

BACKGROUND

High-quality 3D-anatomy of the day is needed for treatment plan adaptation in radiotherapy. For online x-ray-based CBCT workflows, one approach is to create a synthetic CT or to utilize a fan-beam CT with corresponding registrations. The former potentially introduces uncertainties in the dose calculation if deformable image registration is used. The latter can introduce burden and complexity to the process, the facility, and the patient.

PURPOSE

Using the CBCT of the day, acquired on the treatment device, for direct dose calculation and plan adaptation can overcome these limitations. This study aims to assess the accuracy of the calculated dose on the CBCT scans acquired on a Halcyon linear accelerator equipped with HyperSight.

METHODS

HyperSight's new CBCT reconstruction algorithm includes improvements in scatter correction, HU calibration of the imager, and beam shape adaptation. Furthermore, HyperSight introduced a new x-ray detector. To show the effect of the implemented improvements, gamma comparisons of 2%/2 mm, 2%/1 mm, and 1%/1 mm were made between the dose distribution in phantoms calculated on the CBCT reconstructions and the simulation CT scans, considering this the standard of care. The resulting gamma passing rates were compared to those obtained with the Halcyon 3.0 reconstruction and hardware without HyperSight's technologies. Various anatomical phantoms for dosimetric evaluations on brain, head and neck, lung, breast, and prostate cases have been used in this study.

RESULTS

The overall results demonstrated that HyperSight outperformed the Halcyon 3.0 version. Based on the gamma analysis, the calculated dose using HyperSight was closer to the CT scan-based doses than the calculated dose using iCBCT Halcyon 3.0 for most cases. Over all plans and gamma criteria, Halcyon 3.0 achieved an average passing rate of 92.9%, whereas HyperSight achieved 98.1%.

CONCLUSION

Using HyperSight CBCT images for direct dose calculation, for example, in (online) plan adaptation, seems feasible for the investigated cases.

Comparative analysis of delivered and planned doses in target volumes for lung stereotactic ablative radiotherapy

BACKGROUND

Adaptive therapy has been enormously improved based on the art of generating adaptive computed tomography (ACT) from planning CT (PCT) and the on-board image used for the patient setup. Exploiting the ACT, this study evaluated the dose delivered to patients with non-small-cell lung cancer (NSCLC) patients treated with stereotactic ablative radiotherapy (SABR) and derived relationship between the delivered dose and the parameters obtained through the evaluation procedure.

METHODS

SABR treatment records of 72 patients with NSCLC who were prescribed a dose of 60 Gy (Dprescribed) to the 95% volume of the planning target volume (PTV) in four fractions were analysed in this retrospective study; 288 ACTs were generated by rigid and deformable registration of a PCT to a cone-beam computed tomography (CBCT) per fraction. Each ACT was sent to the treatment planning system (TPS) and treated as an individual PCT to calculate the dose. Delivered dose to a patient was estimated by averaging four doses calculated from four ACTs per treatment. Through the process, each ACT provided the geometric parameters, such as mean displacement of the deformed PTV voxels (Warpmean) and Dice similarity coefficient (DSC) from deformation vector field, and dosimetric parameters, e.g. difference of homogeneity index (ΔHI, HI defined as (D2%-D98%)/Dprescribed*100) and mean delivered dose to the PTV (Dmean), obtained from the dose statistics in the TPS. Those parameters were analyzed using multiple linear regression and one-way-ANOVA of SPSS® (version 27).

RESULTS

The prescribed dose was confirmed to be fully delivered to internal target volume (ITV) within maximum difference of 1%, and the difference between the planned and delivered doses to the PTV was agreed within 6% for more than 95% of the ACT cases. Volume changes of the ITV during the treatment course were observed to be minor in comparison of their standard deviations. Multiple linear regression analysis between the obtained parameters and the dose delivered to 95% volume of the PTV (D95%) revealed four PTV parameters [Warpmean, DSC, ΔHI between the PCT and ACT, Dmean] and the PTV D95% to be significantly related with P-values < 0.05. The ACT cases of high ΔHI were caused by higher values of the Warpmean and DSC from the deformable image registration, resulting in lower PTV D95% delivered. The mean values of PTV D95% and Warpmean showed significant differences depending on the lung lobe where the tumour was located.

CONCLUSIONS

Evaluation of the dose delivered to patients with NSCLC treated with SABR using ACTs confirmed that the prescribed dose was accurately delivered to the ITV. However, for the PTV, certain ACT cases characterised by high HI deviations from the original plan demonstrated variations in the delivered dose. These variations may potentially arise from factors such as patient setup during treatment, as suggested by the statistical analyses of the parameters obtained from the dose evaluation process.

Virtual cone-beam computed tomography simulator with human phantom library and its application to the elemental material decomposition

PURPOSE

The purpose of this study is to develop a virtual CBCT simulator with a head and neck (HN) human phantom library and to demonstrate the feasibility of elemental material decomposition (EMD) for quantitative CBCT imaging using this virtual simulator.

METHODS

The library of 36 HN human phantoms were developed by extending the ICRP 110 adult phantoms based on human age, height, and weight statistics. To create the CBCT database for the library, a virtual CBCT simulator that simulated the direct and scattered X-ray on a flat panel detector using ray-tracing and deep-learning (DL) models was used. Gaussian distributed noise was also included on the flat panel detector, which was evaluated using a real CBCT system. The usefulness of the virtual CBCT system was demonstrated through the application of the developed DL-based EMD model for case involving virtual phantom and real patient.

RESULTS

The virtual simulator could generate various virtual CBCT images based on the human phantom library, and the prediction of the EMD could be successfully performed by preparing the CBCT database from the proposed virtual system, even for a real patient. The CBCT image degradation owing to the scattered X-ray and the statistical noise affected the prediction accuracy, although these effects were minimal. Furthermore, the elemental distribution using the real CBCT image was also predictable.

CONCLUSIONS

This study demonstrated the potential of using computer vision for medical data preparation and analysis, which could have important implications for improving patient outcomes, especially in adaptive radiation therapy.

The geometric and dosimetric accuracy of kilovoltage cone beam computed tomography images for adaptive treatment: a systematic review

Objectives

To provide an overview and meta-analysis of different techniques adopted to accomplish kVCBCT for dose calculation and automated segmentation.

Methods

A systematic review and meta-analysis were performed on eligible studies demonstrating kVCBCT-based dose calculation and automated contouring of different tumor features. Meta-analysis of the performance was accomplished on the reported γ analysis and dice similarity coefficient (DSC) score of both collected results as three subgroups (head and neck, chest, and abdomen).

Results

After the literature scrutinization (n = 1008), 52 papers were recognized for the systematic review. Nine studies of dosimtric studies and eleven studies of geometric analysis were suitable for inclusion in meta-analysis. Using kVCBCT for treatment replanning depends on a method used. Deformable Image Registration (DIR) methods yielded small dosimetric error (≤2%), γ pass rate (≥90%) and DSC (≥0.8). Hounsfield Unit (HU) override and calibration curve-based methods also achieved satisfactory yielded small dosimetric error (≤2%) and γ pass rate ((≥90%), but they are prone to error due to their sensitivity to a vendor-specific variation in kVCBCT image quality.

Conclusions

Large cohorts of patients ought to be undertaken to validate methods achieving low levels of dosimetric and geometric errors. Quality guidelines should be established when reporting on kVCBCT, which include agreed metrics for reporting on the quality of corrected kVCBCT and defines protocols of new site-specific standardized imaging used when obtaining kVCBCT images for adaptive radiotherapy.

Advances in knowledge

This review gives useful knowledge about methods making kVCBCT feasible for kVCBCT-based adaptive radiotherapy, simplifying patient pathway and reducing concomitant imaging dose to the patient.

Dosimetric benefits of adaptive radiation therapy for patients with stage III non-small cell lung cancer

BACKGROUND

Daily adaptive radiation therapy (ART) of patients with non-small cell lung cancer (NSCLC) lowers organs at risk exposure while maintaining the planning target volume (PTV) coverage. Thus, ART allows an isotoxic approach with increased doses to the PTV that could improve local tumor control. Herein we evaluate daily online ART strategies regarding their impact on relevant dose-volume metrics.

METHODS

Daily cone-beam CTs (1 × n = 28, 1 × n = 29, 11 × n = 30) of 13 stage III NSCLC patients were converted into synthetic CTs (sCTs). Treatment plans (TPs) were created retrospectively on the first-fraction sCTs (sCT₁) and subsequently transferred unaltered to the sCTs of the remaining fractions of each patient (sCT_2-n) (IGRT scenario). Two additional TPs were generated on sCT_2-n: one minimizing the lung-dose while preserving the D_95%(PTV) (isoeffective scenario), the other escalating the D_95%(PTV) with a constant V_20Gy(lung_ipsilateral) (isotoxic scenario).

RESULTS

Compared to the original TPs predicted dose, the median D_95%(PTV) in the IGRT scenario decreased by 1.6 Gy ± 4.2 Gy while the V_20Gy(lung_ipsilateral) increased in median by 1.1% ± 4.4%. The isoeffective scenario preserved the PTV coverage and reduced the median V_20Gy(lung_ipsilateral) by 3.1% ± 3.6%. Furthermore, the median V_5%(heart) decreased by 2.9% ± 6.4%. With an isotoxic prescription, a median dose-escalation to the gross target volume of 10.0 Gy ± 8.1 Gy without increasing the V_20Gy(lung_ipsilateral) and V_5%(heart) was feasible.

CONCLUSIONS

We demonstrated that even without reducing safety margins, ART can reduce lung-doses, while still reaching adequate target coverage or escalate target doses without increasing ipsilateral lung exposure. Clinical benefits by means of toxicity and local control of both strategies should be evaluated in prospective clinical trials.

Evaluating on-board kVCT- and MVCT-based dose calculation accuracy using a thorax phantom for helical tomotherapy treatments

Purpose.To evaluate the impact of CT number calibration and imaging parameter selection on dose calculation accuracy relative to the CT planning process in thoracic treatments for on-board helical CT imaging systems used in helical tomotherapy.Methods and Materials.Direct CT number calibrations were performed with appropriate protocols for each imaging system using an electron density phantom. Large volume and SBRT treatment plans were simulated and optimized for planning CT scans of an anthropomorphic thorax phantom and transferred to registered kVCT and MVCT scans of the phantom as appropriate. Relevant DVH metrics and dose-difference maps were used to evaluate and compare dose calculation accuracy relative to the planning CT based on a variation in imaging parameters applied for the on-board systems.Results.For helical kVCT scans of the thorax phantom, median differences in DVH parameters for the large volume treatment plan were less than ±1% with dose to the target volume either over- or underestimated depending on the imaging parameters utilized for CT number calibration and thorax phantom acquisition. For the lung SBRT plan calculated on helical kVCT scans, median dose differences were up to -2.7% with a more noticeable dependence on parameter selection. For MVCT scans, median dose differences for the large volume plan were within +2% with dose to the target overestimated regardless of the imaging protocol.Conclusion.Accurate dose calculations (median errors of <±1%) using a thorax phantom simulating realistic patient geometry and scatter conditions can be achieved with images acquired with a helical kVCT system on a helical tomotherapy unit. This accuracy is considerably improved relative to that achieved with the MV-based approach. In a clinical setting, careful consideration should be made when selecting appropriate kVCT imaging parameters for this process as dose calculation accuracy was observed to vary with both parameter selection and treatment type.

Evaluation of Daily CT for EPID-Based Transit In Vivo Dosimetry

Purpose

The difference in anatomical structure and positioning between planning and treatment may lead to bias in electronic portal image device (EPID)-based in vivo dosimetry calculations. The purpose of this study was to use daily CT instead of planning CT as a reference for EPID-based in vivo dosimetry calculations and to analyze the necessity of using daily CT for EPID-based in vivo dosimetry calculations in terms of patient quality assurance.

Materials and Methods

Twenty patients were enrolled in this study. The study design included eight different sites (the cervical, nasopharyngeal, and oral cavities, rectum, prostate, bladder, lung, and esophagus). All treatments were delivered with a CT-linac 506c (UIH, Shanghai) using 6 MV photon beams. This machine is equipped with diagnosis-level fan-beam CT and an amorphous silicon EPID XRD1642 (Varex Imaging Corporation, UT, USA). A Monte Carlo algorithm was developed to calculate the transmit EPID image. A pretreatment measurement was performed to assess system accuracy by delivering based on a homogeneous phantom (RW3 slab, PTW, Freiburg). During treatment, each patient underwent CT scanning before delivery either once or twice for a total of 268 fractions obtained daily CT images. Patients may have had a position correction that followed our image-guided radiation therapy (IGRT) procedure. Meanwhile, transmit EPID images were acquired for each field during delivery. After treatment, all patient CTs were reviewed to ensure that there was no large anatomical change between planning and treatment. The reference of transmit EPID images was calculated based on both planning and daily CTs, and the IGRT correction was corrected for the EPID calculation. The gamma passing rate (3 mm 3%, 2 mm 3%, and 2 mm 2%) was calculated and compared between the planning CT and daily CT. Mechanical errors [ ± 1 mm, ± 2 mm, ± 5 mm multileaf collimator (MLC) systematic shift and 3%, 5% monitor unit (MU) scaling] were also introduced in this study for comparing detectability between both types of CT.

Result

The average (standard deviation) gamma passing rate (3 mm 3%, 2 mm 3%, and 2 mm 2%) in the RW3 slab phantom was 99.6% ± 1.0%, 98.9% ± 2.1%, and 97.2% ± 3.9%. For patient measurement, the average (standard deviation) gamma passing rates were 87.8% ± 14.0%, 82.2% ± 16.9%, and 74.2% ± 18.9% for using planning CTs as reference and 93.6% ± 8.2%, 89.7% ± 11.0%, and 82.8% ± 14.7% for using daily CTs as reference. There were significant differences between the planning CT and daily CT results. All p-values (Mann–Whitney test) were less than 0.001. In terms of error simulation, nonparametric test shows that there were significant differences between practical daily results and error simulation results (p < 0.001). The receiver operating characteristic (ROC) analysis indicated that the detectability of mechanical delivery error using daily CT was better than that of planning CT. AUC_{Daily CT} = 0.63–0.96 and AUC_{Planning CT} = 0.49–0.93 in MLC systematic shift and AUC_{Daily CT} = 0.56–0.82 and AUC_{Planning CT} = 0.45–0.73 in MU scaling.

Conclusion

This study shows the feasibility and effectiveness of using two-dimensional (2D) EPID portal image and daily CT-based in vivo dosimetry for intensity-modulated radiation therapy (IMRT) verification during treatment. The daily CT-based in vivo dosimetry has better sensitivity and specificity to identify the variation of IMRT in MLC-related and dose-related errors than planning CT-based.

Clinical suitability of deep learning based synthetic CTs for adaptive proton therapy of lung cancer

Purpose

Adaptive proton therapy (APT) of lung cancer patients requires frequent volumetric imaging of diagnostic quality. Cone‐beam CT (CBCT) can provide these daily images, but x‐ray scattering limits CBCT‐image quality and hampers dose calculation accuracy. The purpose of this study was to generate CBCT‐based synthetic CTs using a deep convolutional neural network (DCNN) and investigate image quality and clinical suitability for proton dose calculations in lung cancer patients.

Methods

A dataset of 33 thoracic cancer patients, containing CBCTs, same‐day repeat CTs (rCT), planning‐CTs (pCTs), and clinical proton treatment plans, was used to train and evaluate a DCNN with and without a pCT‐based correction method. Mean absolute error (MAE), mean error (ME), peak signal‐to‐noise ratio, and structural similarity were used to quantify image quality. The evaluation of clinical suitability was based on recalculation of clinical proton treatment plans. Gamma pass ratios, mean dose to target volumes and organs at risk, and normal tissue complication probabilities (NTCP) were calculated. Furthermore, proton radiography simulations were performed to assess the HU‐accuracy of sCTs in terms of range errors.

Results

On average, sCTs without correction resulted in a MAE of 34 ± 6 HU and ME of 4 ± 8 HU. The correction reduced the MAE to 31 ± 4HU (ME to 2 ± 4HU). Average 3%/3 mm gamma pass ratios increased from 93.7% to 96.8%, when the correction was applied. The patient specific correction reduced mean proton range errors from 1.5 to 1.1 mm. Relative mean target dose differences between sCTs and rCT were below ± 0.5% for all patients and both synthetic CTs (with/without correction). NTCP values showed high agreement between sCTs and rCT (<2%).

Conclusion

CBCT‐based sCTs can enable accurate proton dose calculations for APT of lung cancer patients. The patient specific correction method increased the image quality and dosimetric accuracy but had only a limited influence on clinically relevant parameters.

Møller D, Assenholt MS, Hansen R, Nyvang L, Ravkilde T, Thomsen MS, Hoffmann L. Density calibrated cone beam CT as a tool for adaptive radiotherapy

BACKGROUND

Visual inspections of anatomical changes observed on daily cone-beam CT (CBCT) images are often used as triggers for radiotherapy plan adaptation to avoid unacceptable dose levels to the target or OARs. Direct CBCT dose calculations would improve the ability to adapt only those plans where dosimetric changes are observed. This study investigates the accuracy of dose calculations on CBCTs.

MATERIALS AND METHODS

Calibration curves were obtained for CBCT imagers at nine identical accelerators. CBCT scans of a phantom with different density inserts were recorded for two scan modes (Head-Neck and Pelvis) and mean calibration curves were calculated. Subsequently, CBCT scans of the phantom with six different density inserts were recorded, the dose distributions on the CBCTs were calculated and compared to dose on the planning CT (pCT). The uncertainty was quantified by the dosimetric difference between the pCT and the CBCT. The two mean calibration curves were used to calculate the daily delivered CBCT dose for ten Head-Neck-, eleven Lung-, and ten pelvic patients. Additional patient calculations were performed using low-HU empirically corrected calibration curves. Patient doses were compared on target coverage and mean dose, and D1cc for OARs.

RESULTS

The dose differences between pCT and CBCT for phantom data were small for all DVH parameters, with mean deviations below ±0.6% for both CBCT modes. For patient data, it was found that low-HU corrected calibration curves performed the best. The mean deviations for the mean dose and coverage of the target were 0.2%±0.7% and 0.1%±0.6%, across all patient groups.

CONCLUSION

Dose calculation on CBCT images results in target coverage and mean dose with an accuracy of the order of 1%, which makes this acceptable for clinical use. The CBCT mode specific calibration curves can be used at all identical imaging devices and for all patient groups.

The effect of the choice of patient model on the performance of in vivo 3D EPID dosimetry to detect variations in patient position and anatomy

PURPOSE

In vivo EPID dosimetry is meant to trigger on relevant differences between delivered and planned dose distributions and should therefore be sensitive to changes in patient position and patient anatomy. Three-dimensional (3D) EPID back-projection algorithms can use either the planning computed tomography (CT) or the daily patient anatomy as patient model for dose reconstruction. The purpose of this study is to quantify the effect of the choice of patient model on the performance of in vivo 3D EPID dosimetry to detect patient-related variations.

METHODS

Variations in patient position and patient anatomy were simulated by transforming the reference planning CT images (pCT) into synthetic daily CT images (dCT) representing a variation of a given magnitude in patient position or in patient anatomy. For each variation, synthetic in vivo EPID data were also generated to simulate the reconstruction of in vivo EPID dose distributions. Both the planning CT images and the synthetic daily CT images could be used as patient model in the reconstructions yielding e D pCT and e D dCT EPID reconstructed dose distributions respectively. The accuracy of e D pCT and e D dCT reconstructions was evaluated against absolute dose measurements made in different phantom setups, and against dose distributions calculated by the treatment planning system (TPS). The comparison was performed by γ-analysis (3% local dose/2 mm). The difference in sensitivity between e D pCT and e D dCT reconstructions to detect variations in patient position and in patient anatomy was investigated using receiver operating characteristic analysis and the number of triggered alerts for 100 volumetric modulated arc therapy plans and 12 variations.

RESULTS

e D dCT showed good agreement with both absolute point dose measurements (<0.5%) and TPS data (γ-mean = 0.52 ± 0.11). The agreement degraded with e D pCT , with the magnitude of the deviation varying with each specific case. e D dCT readily detected combined 3 mm translation setup errors in all directions (AUC = 1.0) and combined 3° rotation setup errors around all axes (AUC = 0.86) whereas e D pCT showed good detectability only for 12 mm translations (AUC = 0.85) and 9° rotations (AUC = 0.80). Conversely, e D pCT manifested a higher sensitivity to patient anatomical changes resulting in AUC values of 0.92/0.95 for a 6 mm patient contour expansion/contraction compared to 0.70/0.64 with e D dCT . Using |ΔPTV_D50 | > 3% as clinical tolerance level, the percentage of alerts for 6 mm changes in patient contour were 85%/27% with e D pCT / e D dCT .

CONCLUSIONS

With planning CT images as patient model, EPID dose reconstructions underestimate the dosimetric effects caused by errors in patient positioning and overestimate the dosimetric effects caused by changes in patient anatomy. The use of the daily patient position and anatomy as patient model for in vivo 3D EPID transit dosimetry improves the ability of the system to detect uncorrected errors in patient position and it reduces the likelihood of false positives due to patient anatomical changes.

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.