Quantitative use of cone-beam computed tomography in proton therapy: challenges and opportunities

The fundamental goal in radiation therapy (RT) is to simultaneously maximize tumor cell killing and healthy tissue sparing. Reducing uncertainty margins improves normal tissue sparing, but generally requires advanced techniques. Adaptive RT (ART) is a compelling technique that leverages daily imaging and anatomical information to support reduced margins and to optimize plan quality for each treatment fraction. An especially exciting avenue for ART is proton therapy (PT), which aims to combine daily plan re-optimization with the unique advantages provided by protons, including reduced integral dose and near-zero dose deposition distal to the target along the beam direction. A core component for ART is onboard image guidance, and currently two options are available on proton systems, including cone-beam computed tomography (CBCT) and CT-on-rail (CToR) imaging. While CBCT suffers from poorer image quality compared to CToR imaging, CBCT platforms can be more easily integrated with PT systems and thus may support more streamlined adaptive proton therapy (APT). In this review, we present current status of CBCT application to proton therapy dose evaluation and plan adaptation, including progress, challenges and future directions.

Cone-beam computed tomography (CBCT) image-quality improvement using a denoising diffusion probabilistic model conditioned by pseudo-CBCT of pelvic regions

Cone-beam computed tomography (CBCT) is widely used in radiotherapy to image patient configuration before treatment but its image quality is lower than planning CT due to scattering, motion, and reconstruction methods. This reduces the accuracy of Hounsfield units (HU) and limits its use in adaptive radiation therapy (ART). However, synthetic CT (sCT) generation using deep learning methods for CBCT intensity correction faces challenges due to deformation. To address these issues, we propose enhancing CBCT quality using a conditional denoising diffusion probability model (CDDPM), which is trained on pseudo-CBCT created by adding pseudo-scatter to planning CT. The CDDPM transforms CBCT into high-quality sCT, improving HU accuracy while preserving anatomical configuration. The performance evaluation of the proposed sCT showed a reduction in mean absolute error (MAE) from 81.19 HU for CBCT to 24.89 HU for the sCT. Peak signal-to-noise ratio (PSNR) improved from 31.20 dB for CBCT to 33.81 dB for the sCT. The Dice and Jaccard coefficients between CBCT and sCT for the colon, prostate, and bladder ranged from 0.69 to 0.91. When compared to other deep learning models, the proposed sCT outperformed them in terms of accuracy and anatomical preservation. The dosimetry analysis for prostate cancer revealed a dose error of over 10% with CBCT but nearly 0% with the sCT. Gamma pass rates for the proposed sCT exceeded 90% for all dose criteria, indicating high agreement with CT-based dose distributions. These results show that the proposed sCT improves image quality, dosimetry accuracy, and treatment planning, advancing ART for pelvic cancer.

Deep learning based super-resolution for CBCT dose reduction in radiotherapy

BACKGROUND

Cone-beam computed tomography (CBCT) is a crucial daily imaging modality in image-guided and adaptive radiotherapy. However, the use of ionizing radiation in CBCT imaging increases the risk of secondary cancers, which is particularly concerning for pediatric patients. Deep learning super-resolution has shown promising results in enhancing the resolution of photographic and medical images but has not yet been explored in the context of CBCT dose reduction.

PURPOSE

To facilitate CBCT imaging dose reduction, we propose using an enhanced super-resolution generative adversarial network (ESRGAN) in both the projection and image domains to restore the image quality of low-dose CBCT.

METHODS

An extensive projection database, containing 2997 CBCT scans from head and neck cancer patients, was used to train two different ESRGAN models to generate super-resolution CBCTs. One model operated in the projection domain, using pairs of simulated low-resolution (low-dose) and original high-resolution (high-dose) projections and yielded CBCTSRpro. The other model operated in the image domain, using pairs of axial slices from reconstructed low-resolution and high-resolution CBCTs (CBCTLR and CBCTHR) and resulted in CBCTSRimg. Super-resolution CBCTs were evaluated in terms of image similarity (MAE, ME, PSNR, and SSIM), noise characteristics, spatial resolution, and registration accuracy, using the original CBCT as a reference. To test the perceptual difference between the original and super-resolution CBCT, we performed a visual Turing test.

RESULTS

Visually, both super-resolution approaches in the projection and image domains improved the image quality of low-dose CBCTs. This was confirmed by the visual Turing test, that showed low accuracy, sensitivity, and specificity, indicating almost no perceptual difference between CBCTHR and the super-resolution CBCTs. CBCTSRimg (accuracy: 0.55, sensitivity: 0.59, specificity: 0.50) performed slightly better than CBCTSRpro (accuracy: 0.59, sensitivity: 0.61, specificity: 0.57). Image similarity metrics were affected by varying noise levels and did not reflect the visual improvements, with MAE/ME/PSNR/SSIM values of 110.4 HU/2.9 HU/40.4 dB/0.82 for CBCTLR, 136.6 HU/-0.4 HU/38.6 dB/0.77 for CBCTSRpro, and 128.2 HU/1.9 HU/39.0 dB/0.80 for CBCTSRimg. In terms of spatial resolution, quantified by calculating 10% levels of the task transfer function, both CBCTSRpro and CBCTSRimg outperformed CBCTLR and nearly matched the reference CBCTHR (CBCTLR: 0.66 lp/mm, CBCTSRpro: 0.88 lp/mm, CBCTSRimg: 0.95 lp/mm, CBCTHR: 1.01 lp/mm). Noise characteristics of CBCTSRimg and CBCTSRpro were comparable to the reference CBCTHR. Registration parameters showed negligible differences for all CBCTs (CBCTLR, CBCTSRpro, CBCTSRimg), with average absolute differences in registration parameters being below 0.4° for rotations and below 0.06 mm for translations (CBCTHR as reference).

CONCLUSIONS

This study demonstrates that deep learning can be a valuable tool for CBCT dose reduction in CBCT-guided radiotherapy by acquiring low-dose CBCTs and restoring the image quality using deep learning super-resolution. The results suggest that higher quality images can be generated when super-resolution is performed in the image domain compared to the projection domain.

Investigation of 2D anti-scatter grid implementation in a gantry-mounted cone beam computed tomography system for proton therapy

Background and purpose

Robust scatter mitigation by 2D anti-scatter grids (2D-ASG) in proton therapy cone beam computed tomography (CBCT) may improve target visualization and computed tomography (CT) number fidelity, allowing online dose verifications and plan adaptations. However, grid artifact-free implementation of 2D-ASG depends on the CBCT system characteristics. Thus, we investigated the feasibility of 2D-ASG implementation in a proton therapy gantry-mounted CBCT system and evaluated its impact on image quality.

Materials and methods

A prototype 2D-ASG and a grid support platform were developed for a proton therapy CBCT system with a 340 cm source to imager distance. The effect of gantry flex on 2D-ASG's wall shadows and scan-to-scan reproducibility of 2D-ASG's wall shadows were evaluated. Experiments were conducted to assess 2D-ASG's wall shadow suppression and the effect of 2D-ASG on image quality.

Results

While maximum displacement in 2D-ASG wall shadows was 103 µm during gantry rotation, the drift from baseline over 3 months was 8 µm and 1 µm in the transverse and axial directions. 2D-ASG shadows were successfully suppressed in CBCT images. With 2D-ASG, maximum Hounsfield Unit (HU) nonuniformity decreased from 134 to 45 HU, contrast-to-noise ratio (CNR) increased by a factor of 2.5, and HU errors were reduced from 34 % to 5 %.

Conclusions

Proton therapy gantry flex was highly reproducible and did not noticeably affect 2D-ASG wall shadow suppression in CBCT images, supporting its feasibility in proton therapy CBCT system. Improved CT accuracy and artifact reduction with 2D-ASG could enhance CBCT-based proton therapy dose calculations.

Cone Beam Computed Tomography Image-Quality Improvement Using "One-Shot" Super-resolution

Cone beam computed tomography (CBCT) images are convenient representations for obtaining information about patients' internal organs, but their lower image quality than those of treatment planning CT images constitutes an important shortcoming. Several proposed CBCT image-quality improvement methods based on deep learning require large amounts of training data. Our newly developed model using a super-resolution method, "one-shot" super-resolution (OSSR) based on the "zero-shot" super-resolution method, requires only small amounts of training data to improve CBCT image quality using only the target CBCT image and the paired treatment planning CT image. For this study, pelvic CBCT images and treatment planning CT images of 30 prostate cancer patients were used. We calculated the root mean squared error (RMSE), peak signal-to-noise ratio (PSNR), and structural similarity (SSIM) to evaluate image-quality improvement and normalized mutual information (NMI) as a quantitative evaluation of positional accuracy. Our proposed method can improve CBCT image quality without requiring large amounts of training data. After applying our proposed method, the resulting RMSE, PSNR, SSIM, and NMI between the CBCT images and the treatment planning CT images were as much as 0.86, 1.05, 1.03, and 1.31 times better than those obtained without using our proposed method. By comparison, CycleGAN exhibited values of 0.91, 1.03, 1.02, and 1.16. The proposed method achieved performance equivalent to that of CycleGAN, which requires images from approximately 30 patients for training. Findings demonstrated improvement of CBCT image quality using only the target CBCT images and the paired treatment planning CT images.

Feasibility of proton dosimetry overriding planning CT with daily CBCT elaborated through generative artificial intelligence tools

Radiotherapy commonly utilizes cone beam computed tomography (CBCT) for patient positioning and treatment monitoring. CBCT is deemed to be secure for patients, making it suitable for the delivery of fractional doses. However, limitations such as a narrow field of view, beam hardening, scattered radiation artifacts, and variability in pixel intensity hinder the direct use of raw CBCT for dose recalculation during treatment. To address this issue, reliable correction techniques are necessary to remove artifacts and remap pixel intensity into Hounsfield Units (HU) values. This study proposes a deep-learning framework for calibrating CBCT images acquired with narrow field of view (FOV) systems and demonstrates its potential use in proton treatment planning updates. Cycle-consistent generative adversarial networks (cGAN) processes raw CBCT to reduce scatter and remap HU. Monte Carlo simulation is used to generate CBCT scans, enabling the possibility to focus solely on the algorithm's ability to reduce artifacts and cupping effects without considering intra-patient longitudinal variability and producing a fair comparison between planning CT (pCT) and calibrated CBCT dosimetry. To showcase the viability of the approach using real-world data, experiments were also conducted using real CBCT. Tests were performed on a publicly available dataset of 40 patients who received ablative radiation therapy for pancreatic cancer. The simulated CBCT calibration led to a difference in proton dosimetry of less than 2%, compared to the planning CT. The potential toxicity effect on the organs at risk decreased from about 50% (uncalibrated) up the 2% (calibrated). The gamma pass rate at 3%/2 mm produced an improvement of about 37% in replicating the prescribed dose before and after calibration (53.78% vs 90.26%). Real data also confirmed this with slightly inferior performances for the same criteria (65.36% vs 87.20%). These results may confirm that generative artificial intelligence brings the use of narrow FOV CBCT scans incrementally closer to clinical translation in proton therapy planning updates.

Use of Virtual CT and On-Treatment MRI to Reduce Radiation Dose and Anesthesia Exposure Associated With the Adaptive Workflow in Pediatric Patients Treated With Intensity Modulated Proton Therapy

Purpose

The purpose of this study was to determine whether virtual computed tomography (vCT) derived from daily cone beam computed tomography (CBCT), or on-treatment magnetic resonance imaging (MRItx) can replace quality assurance computed tomography (qCT) in our clinical workflow to minimize imaging dose and potentially anesthesia exposure in patients requiring plan adaptation.

Methods and Materials

Pediatric patients (age <24 years) treated from 2020 to 2023 with intensity modulated proton therapy with at least 1 qCT during proton therapy were eligible. For cases that required plan adaptation, the dose was recalculated on vCT and compared with same-day qCT as well as the original planning computed tomography (pCT). Anatomic changes triggering plan adaptation were grouped into categories. Two pediatric radiation oncologists verified whether these changes could be detected using CBCT, qCT, and/or MRItx. A new adaptive imaging workflow was proposed to limit imaging dose and anesthesia exposure.

Results

One hundred sixty-eight pediatric patients were treated from 2020 to 2023. Across all patients, there were 517 qCT scans and 61 MRItx acquired. The median number of qCT scans per patient was 3 (range, 1-5). The treatment plans for 20 patients (12%) were adapted. In all patients requiring plan adaptation, there was a correlation between dose differences in target coverage and maximum body dose when comparing vCT with pCT and qCT with pCT (n = 20, r2 = 0.79, P < .01, and r2 = 0.32 P = .01, respectively). The most common reason for adaptation was tissue change (eg, inflammation, changes in abdominal gas, or diaphragmatic variability) in the beam path (10/20) and changes in tumor volume (6/20). All cases of weight change, tissue change in beam path, and unreproducible setup could be detected on CBCT. All cases of change in tumor volume within the brain were detected on MRItx. Replacing the qCT with the vCT was associated with an estimated median reduction of imaging dose by 50% and anesthesia exposure by 1.5 hours.

Conclusions

vCT derived from daily CBCT only or MRItx can safely replace qCT for monitoring dosimetric changes to trigger a new pCT in our clinical workflow. This change would potentially reduce imaging dose and anesthesia exposure.

Post-radiotherapy stage III/IV non-small cell lung cancer radiomics research: a systematic review and comparison of CLEAR and RQS frameworks

BACKGROUND

Lung cancer, the second most common cancer, presents persistently dismal prognoses. Radiomics, a promising field, aims to provide novel imaging biomarkers to improve outcomes. However, clinical translation faces reproducibility challenges, despite efforts to address them with quality scoring tools.

OBJECTIVE

This study had two objectives: 1) identify radiomics biomarkers in post-radiotherapy stage III/IV nonsmall cell lung cancer (NSCLC) patients, 2) evaluate research quality using the CLEAR (CheckList_for_EvaluAtion_of_Radiomics_research), RQS (Radiomics_Quality_Score) frameworks, and formulate an amalgamated CLEAR-RQS tool to enhance scientific rigor.

MATERIALS AND METHODS

A systematic literature review (Jun-Aug 2023, MEDLINE/PubMed/SCOPUS) was conducted concerning stage III/IV NSCLC, radiotherapy, and radiomic features (RF). Extracted data included study design particulars, such as sample size, radiotherapy/CT technique, selected RFs, and endpoints. CLEAR and RQS were merged into a CLEAR-RQS checklist. Three readers appraised articles utilizing CLEAR, RQS, and CLEAR-RQS metrics.

RESULTS

Out of 871 articles, 11 met the inclusion/exclusion criteria. The Median cohort size was 91 (range: 10-337) with 9 studies being single-center. No common RF were identified. The merged CLEAR-RQS checklist comprised 61 items. Most unreported items were within CLEAR's "methods" and "open-source," and within RQS's "phantom-calibration," "registry-enrolled prospective-trial-design," and "cost-effective-analysis" sections. No study scored above 50% on RQS. Median CLEAR scores were 55.74% (32.33/58 points), and for RQS, 17.59% (6.3/36 points). CLEAR-RQS article ranking fell between CLEAR and RQS and aligned with CLEAR.

CONCLUSION

Radiomics research in post-radiotherapy stage III/IV NSCLC exhibits variability and frequently low-quality reporting. The formulated CLEAR-RQS checklist may facilitate education and holds promise for enhancing radiomics research quality.

CLINICAL RELEVANCE STATEMENT

Current radiomics research in the field of stage III/IV postradiotherapy NSCLC is heterogenous, lacking reproducibility, with no identified imaging biomarker. Radiomics research quality assessment tools may enhance scientific rigor and thereby facilitate radiomics translation into clinical practice.

KEY POINTS

There is heterogenous and low radiomics research quality in postradiotherapy stage III/IV nonsmall cell lung cancer. Barriers to reproducibility are small cohort size, nonvalidated studies, missing technical parameters, and lack of data, code, and model sharing. CLEAR (CheckList_for_EvaluAtion_of_Radiomics_research), RQS (Radiomics_Quality_Score), and the amalgamated CLEAR-RQS tool are useful frameworks for assessing radiomics research quality and may provide a valuable resource for educational purposes in the field of radiomics.

Supervised deep learning-based synthetic computed tomography from kilovoltage cone-beam computed tomography images for adaptive radiation therapy in head and neck cancer

PURPOSE

To generate and investigate a supervised deep learning algorithm for creating synthetic computed tomography (sCT) images from kilovoltage cone-beam computed tomography (kV-CBCT) images for adaptive radiation therapy (ART) in head and neck cancer (HNC).

MATERIALS AND METHODS

This study generated the supervised U-Net deep learning model using 3,491 image pairs from planning computed tomography (pCT) and kV-CBCT datasets obtained from 40 HNC patients. The dataset was split into 80% for training and 20% for testing. The evaluation of the sCT images compared to pCT images focused on three aspects: Hounsfield units accuracy, assessed using mean absolute error (MAE) and root mean square error (RMSE); image quality, evaluated using the peak signal-to-noise ratio (PSNR) and structural similarity index (SSIM) between sCT and pCT images; and dosimetric accuracy, encompassing 3D gamma passing rates for dose distribution and percentage dose difference.

RESULTS

MAE, RMSE, PSNR, and SSIM showed improvements from their initial values of 53.15 ± 40.09, 153.99 ± 79.78, 47.91 ± 4.98 dB, and 0.97 ± 0.02 to 41.47 ± 30.59, 130.39 ± 78.06, 49.93 ± 6.00 dB, and 0.98 ± 0.02, respectively. Regarding dose evaluation, 3D gamma passing rates for dose distribution within sCT images under 2%/2 mm, 3%/2 mm, and 3%/3 mm criteria, yielded passing rates of 92.1% ± 3.8%, 93.8% ± 3.0%, and 96.9% ± 2.0%, respectively. The sCT images exhibited minor variations in the percentage dose distribution of the investigated target and structure volumes. However, it is worth noting that the sCT images exhibited anatomical variations when compared to the pCT images.

CONCLUSION

These findings highlight the potential of the supervised U-Net deep learningmodel in generating kV-CBCT-based sCT images for ART in patients with HNC.

2D antiscatter grid and scatter sampling based CBCT method for online dose calculations during CBCT guided radiation therapy of pelvis

BACKGROUND

Online dose calculations before the delivery of radiation treatments have applications in dose delivery verification, online adaptation of treatment plans, and simulation-free treatment planning. While dose calculations by directly utilizing CBCT images are desired, dosimetric accuracy can be compromised due to relatively lower HU accuracy in CBCT images.

PURPOSE

In this work, we propose a novel CBCT imaging pipeline to enhance the accuracy of CBCT-based dose calculations in the pelvis region. Our approach aims to improve the HU accuracy in CBCT images, thereby improving the overall accuracy of CBCT-based dose calculations prior to radiation treatment delivery.

METHODS

An in-house developed quantitative CBCT pipeline was implemented to address the CBCT raw data contamination problem. The pipeline combines algorithmic data correction strategies and 2D antiscatter grid-based scatter rejection to achieve high CT number accuracy. To evaluate the effect of the quantitative CBCT pipeline on CBCT-based dose calculations, phantoms mimicking pelvis anatomy were scanned using a linac-mounted CBCT system, and a gold standard multidetector CT used for treatment planning (pCT). A total of 20 intensity-modulated treatment plans were generated for five targets, using 6 and 10 MV flattening filter-free beams, and utilizing small and large pelvis phantom images. For each treatment plan, four different dose calculations were performed using pCT images and three CBCT imaging configurations: quantitative CBCT, clinical CBCT protocol, and a high-performance 1D antiscatter grid (1D ASG). Subsequently, dosimetric accuracy was evaluated for both targets and organs at risk as a function of patient size, target location, beam energy, and CBCT imaging configuration.

RESULTS

When compared to the gold-standard pCT, dosimetric errors in quantitative CBCT-based dose calculations were not significant across all phantom sizes, beam energies, and treatment sites. The largest error observed was 0.6% among all dose volume histogram metrics and evaluated dose calculations. In contrast, dosimetric errors reached up to 7% and 97% in clinical CBCT and high-performance ASG CBCT-based treatment plans, respectively. The largest dosimetric errors were observed in bony targets in the large phantom treated with 6 MV beams. The trends of dosimetric errors in organs at risk were similar to those observed in the targets.

CONCLUSIONS

The proposed quantitative CBCT pipeline has the potential to provide comparable dose calculation accuracy to the gold-standard planning CT in photon radiation therapy for the abdomen and pelvis. These robust dose calculations could eliminate the need for density overrides in CBCT images and enable direct utilization of CBCT images for dose delivery monitoring or online treatment plan adaptations before the delivery of radiation treatments.

Generation and evaluation of anatomy-preserving virtual CT for online adaptive proton therapy

BACKGROUND

Daily CTs generated by CBCT correction are required for daily replanning in online-adaptive proton therapy (APT) to effectively deal with inter-fractional changes. Out of the currently available methods, the suitability of a daily CT generation method for proton dose calculation also depends on the anatomical site.

PURPOSE

We propose an anatomy-preserving virtual CT (APvCT) method as a hybrid method of CBCT correction, which is especially suitable for large anatomy deformations. The accuracy of the hybrid method was assessed by comparison with the corrected CBCT (cCBCT) and virtual CT (vCT) methods in the context of online APT.

METHODS

Seventy-one daily CBCTs of four prostate cancer patients treated with intensity modulated proton therapy (IMPT) were converted to daily CTs using cCBCT, vCT, and the newly proposed APvCT method. In APvCT, planning CT (pCT) were mapped to CBCT geometry using deformable image registration with boundary conditions on controlling regions of interest (ROIs) created with deep learning segmentation on cCBCT. The relative frequency distribution (RFD) of HU, mass density and stopping power ratio (SPR) values were assessed and compared with the pCT. The ROIs in the APvCT and vCT were compared with cCBCT in terms of Dice similarity coefficient (DSC) and mean distance-to-agreement (mDTA). For each patient, a robustly optimized IMPT plan was created on the pCT and subsequent daily adaptive plans on daily CTs. For dose distribution comparison on the same anatomy, the daily adaptive plans on cCBCT and vCT were recalculated on the corresponding APvCT. The dose distributions were compared in terms of isodose volumes and 3D global gamma-index passing rate (GPR) at γ(2%, 2 mm) criterion.

RESULTS

For all patients, no noticeable difference in RFDs was observed amongst APvCT, vCT, and pCT except in cCBCT, which showed a noticeable difference. The minimum DSC value was 0.96 and 0.39 for contours in APvCT and vCT respectively. The average value of mDTA for APvCT was 0.01 cm for clinical target volume and ≤0.01 cm for organs at risk, which increased to 0.18 cm and ≤0.52 cm for vCT. The mean GPR value was 90.9%, 64.5%, and 67.0% for APvCT versus cCBCT, vCT versus cCBCT, and APvCT versus vCT, respectively. When recalculated on APvCT, the adaptive cCBCT and vCT plans resulted in mean GPRs of 89.5 ± 5.1% and 65.9 ± 19.1%, respectively. The mean DSC values for 80.0%, 90.0%, 95.0%, 98.0%, and 100.0% isodose volumes were 0.97, 0.97, 0.97, 0.95, and 0.91 for recalculated cCBCT plans, and 0.89, 0.88, 0.87, 0.85, and 0.81 for recalculated vCT plans. Hausdorff distance for the 100.0% isodose volume in some cases of recalculated cCBCT plans on APvCT exceeded 1.00 cm.

CONCLUSIONS

APvCT contours showed good agreement with reference contours of cCBCT which indicates anatomy preservation in APvCT. A vCT with erroneous anatomy can result in an incorrect adaptive plan. Further, slightly lower values of GPR between the APvCT and cCBCT-based adaptive plans can be explained by the difference in the cCBCT's SPR RFD from the pCT.

Cone beam CT-based adaptive intensity modulated proton therapy assessment using automated planning for head-and-neck cancer

BACKGROUND

To assess the feasibility of CBCT-based adaptive intensity modulated proton therapy (IMPT) using automated planning for treatment of head and neck (HN) cancers.

METHODS

Twenty HN cancer patients who received radiotherapy and had pretreatment CBCTs were included in this study. Initial IMPT plans were created using automated planning software for all patients. Synthetic CTs (sCT) were then created by deforming the planning CT (pCT) to the pretreatment CBCTs. To assess dose calculation accuracy on sCTs, repeat CTs (rCTs) were deformed to the pretreatment CBCT obtained on the same day to create deformed rCT (rCTdef), serving as gold standard. The dose recalculated on sCT and on rCTdef were compared by using Gamma analysis. The accuracy of DIR generated contours was also assessed. To explore the potential benefits of adaptive IMPT, two sets of plans were created for each patient, a non-adapted IMPT plan and an adapted IMPT plan calculated on weekly sCT images. The weekly doses for non-adaptive and adaptive IMPT plans were accumulated on the pCT, and the accumulated dosimetric parameters of two sets were compared.

RESULTS

Gamma analysis of the dose recalculated on sCT and rCTdef resulted in a passing rate of 97.9% ± 1.7% using 3 mm/3% criteria. With the physician-corrected contours on the sCT, the dose deviation range of using sCT to estimate mean dose for the most organ at risk (OARs) can be reduced to (- 2.37%, 2.19%) as compared to rCTdef, while for V95 of primary or secondary CTVs, the deviation can be controlled within (- 1.09%, 0.29%). Comparison of the accumulated doses from the adaptive planning against the non-adaptive plans reduced mean dose to constrictors (- 1.42 Gy ± 2.79 Gy) and larynx (- 2.58 Gy ± 3.09 Gy). The reductions result in statistically significant reductions in the normal tissue complication probability (NTCP) of larynx edema by 7.52% ± 13.59%. 4.5% of primary CTVs, 4.1% of secondary CTVs, and 26.8% tertiary CTVs didn't meet the V95 > 95% constraint on non-adapted IMPT plans. All adaptive plans were able to meet the coverage constraint.

CONCLUSION

sCTs can be a useful tool for accurate proton dose calculation. Adaptive IMPT resulted in better CTV coverage, OAR sparing and lower NTCP for some OARs as compared with non-adaptive IMPT.

Review and recommendations on deformable image registration uncertainties for radiotherapy applications

Deformable image registration (DIR) is a versatile tool used in many applications in radiotherapy (RT). DIR algorithms have been implemented in many commercial treatment planning systems providing accessible and easy-to-use solutions. However, the geometric uncertainty of DIR can be large and difficult to quantify, resulting in barriers to clinical practice. Currently, there is no agreement in the RT community on how to quantify these uncertainties and determine thresholds that distinguish a good DIR result from a poor one. This review summarises the current literature on sources of DIR uncertainties and their impact on RT applications. Recommendations are provided on how to handle these uncertainties for patient-specific use, commissioning, and research. Recommendations are also provided for developers and vendors to help users to understand DIR uncertainties and make the application of DIR in RT safer and more reliable.

Technical note: Evaluation of deep learning based synthetic CTs clinical readiness for dose and NTCP driven head and neck adaptive proton therapy

BACKGROUND

Adaptive proton therapy workflows rely on accurate imaging throughout the treatment course. Our centre currently utilizes weekly repeat CTs (rCTs) for treatment monitoring and plan adaptations. However, deep learning-based methods have recently shown to successfully correct CBCT images, which suffer from severe imaging artifacts, and generate high quality synthetic CT (sCT) images which enable CBCT-based proton dose calculations.

PURPOSE

To compare daily CBCT-based sCT images to planning CTs (pCT) and rCTs of head and neck (HN) cancer patients to investigate the dosimetric accuracy of CBCT-based sCTs in a scenario mimicking actual clinical practice.

METHODS

Data of 56 HN cancer patients, previously treated with proton therapy was used to generate 1.962 sCT images, using a previously developed and trained deep convolutional neural network. Clinical IMPT treatment plans were recalculated on the pCT, weekly rCTs and daily sCTs. The dosimetric accuracy of sCTs was compared to same day rCTs and the initial planning CT. As a reference, rCTs were also compared to pCTs. The dose difference between sCTs and rCTs/pCT was quantified by calculating the D98 difference for target volumes and Dmean difference for organs-at-risk. To investigate the clinical relevancy of possible dose differences, NTCP values were calculated for dysphagia and xerostomia.

RESULTS

For target volumes, only minor dose differences were found for sCT versus rCT and sCT versus pCT, with dose differences mostly within ±1.5%. Larger dose differences were observed in OARs, where a general shift towards positive differences was found, with the largest difference in the left parotid gland. Delta NTCP values for grade 2 dysphagia and xerostomia were within ±2.5% for 90% of the sCTs.

CONCLUSIONS

Target doses showed high similarity between rCTs and sCTs. Further investigations are required to identify the origin of the dose differences at OAR levels and its relevance in clinical decision making.

The Applications and Pitfalls of Cone-Beam Computed Tomography-Based Synthetic Computed Tomography for Adaptive Evaluation in Pencil-Beam Scanning Proton Therapy

PURPOSE

The study evaluates the efficacy of cone-beam computed tomography (CBCT)-based synthetic CTs (sCT) as a potential alternative to verification CT (vCT) for enhanced treatment monitoring and early adaptation in proton therapy.

METHODS

Seven common treatment sites were studied. Two sets of sCT per case were generated: direct-deformed (DD) sCT and image-correction (IC) sCT. The image qualities and dosimetric impact of the sCT were compared to the same-day vCT.

RESULTS

The sCT agreed with vCT in regions of homogeneous tissues such as the brain and breast; however, notable discrepancies were observed in the thorax and abdomen. The sCT outliers existed for DD sCT when there was an anatomy change and for IC sCT in low-density regions. The target coverage exhibited less than a 5% variance in most DD and IC sCT cases when compared to vCT. The Dmax of serial organ-at-risk (OAR) in sCT plans shows greater deviation from vCT than small-volume dose metrics (D0.1cc). The parallel OAR volumetric and mean doses remained consistent, with average deviations below 1.5%.

CONCLUSION

The use of sCT enables precise treatment and prompt early adaptation for proton therapy. The quality assurance of sCT is mandatory in the early stage of clinical implementation.

A systematic review of volumetric image guidance in proton therapy

In recent years, proton therapy centres have begun to shift from conventional 2D-kV imaging to volumetric imaging systems for image guided proton therapy (IGPT). This is likely due to the increased commercial interest and availability of volumetric imaging systems, as well as the shift from passively scattered proton therapy to intensity modulated proton therapy. Currently, there is no standard modality for volumetric IGPT, leading to variation between different proton therapy centres. This article reviews the reported clinical use of volumetric IGPT, as available in published literature, and summarises their utilisation and workflow where possible. In addition, novel volumetric imaging systems are also briefly summarised highlighting their potential benefits for IGPT and the challenges that need to be overcome before they can be used clinically.

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.

Image-based scatter correction for cone-beam CT using flip swin transformer U-shape network

BACKGROUND

Cone beam computed tomography (CBCT) plays an increasingly important role in image-guided radiation therapy. However, the image quality of CBCT is severely degraded by excessive scatter contamination, especially in the abdominal region, hindering its further applications in radiation therapy.

PURPOSE

To restore low-quality CBCT images contaminated by scatter signals, a scatter correction algorithm combining the advantages of convolutional neural networks (CNN) and Swin Transformer is proposed.

METHODS

In this paper a scatter correction model for CBCT image, the Flip Swin Transformer U-shape network (FSTUNet) model, is proposed. In this model, the advantages of CNN in texture detail and Swin Transformer in global correlation are used to accurately extract shallow and deep features, respectively. Instead of using the original Swin Transformer tandem structure, we build the Flip Swin Transformer Block to achieve a more powerful inter-window association extraction. The validity and clinical relevance of the method is demonstrated through extensive experiments on a Monte Carlo (MC) simulation dataset and frequency split dataset generated by a validated method, respectively.

RESULT

Experimental results on the MC simulated dataset show that the root mean square error of images corrected by the method is reduced from over 100 HU to about 7 HU. Both the structural similarity index measure (SSIM) and the universal quality index (UQI) are close to 1. Experimental results on the frequency split dataset demonstrate that the method not only corrects shading artifacts but also exhibits a high degree of structural consistency. In addition, comparison experiments show that FSTUNet outperforms UNet, Deep Residual Convolutional Neural Network (DRCNN), DSENet, Pix2pixGAN, and 3DUnet methods in both qualitative and quantitative metrics.

CONCLUSIONS

Accurately capturing the features at different levels is greatly beneficial for reconstructing high-quality scatter-free images. The proposed FSTUNet method is an effective solution to CBCT scatter correction and has the potential to improve the accuracy of CBCT image-guided radiation therapy.

A hybrid method of correcting CBCT for proton range estimation with deep learning and deformable image registration

Objective. This study aimed to develop a novel method for generating synthetic CT (sCT) from cone-beam CT (CBCT) of the abdomen/pelvis with bowel gas pockets to facilitate estimation of proton ranges.Approach. CBCT, the same-day repeat CT, and the planning CT (pCT) of 81 pediatric patients were used for training (n= 60), validation (n= 6), and testing (n= 15) of the method. The proposed method hybridizes unsupervised deep learning (CycleGAN) and deformable image registration (DIR) of the pCT to CBCT. The CycleGAN and DIR are respectively applied to generate the geometry-weighted (high spatial-frequency) and intensity-weighted (low spatial-frequency) components of the sCT, thereby each process deals with only the component weighted toward its strength. The resultant sCT is further improved in bowel gas regions and other tissues by iteratively feeding back the sCT to adjust incorrect DIR and by increasing the contribution of the deformed pCT in regions of accurate DIR.Main results. The hybrid sCT was more accurate than deformed pCT and CycleGAN-only sCT as indicated by the smaller mean absolute error in CT numbers (28.7 ± 7.1 HU versus 38.8 ± 19.9 HU/53.2 ± 5.5 HU;P≤ 0.012) and higher Dice similarity of the internal gas regions (0.722 ± 0.088 versus 0.180 ± 0.098/0.659 ± 0.129;P≤ 0.002). Accordingly, the hybrid method resulted in more accurate proton range for the beams intersecting gas pockets (11 fields in 6 patients) than the individual methods (the 90th percentile error in 80% distal fall-off, 1.8 ± 0.6 mm versus 6.5 ± 7.8 mm/3.7 ± 1.5 mm;P≤ 0.013). The gamma passing rates also showed a significant dosimetric advantage by the hybrid method (99.7 ± 0.8% versus 98.4 ± 3.1%/98.3 ± 1.8%;P≤ 0.007).Significance. The hybrid method significantly improved the accuracy of sCT and showed promises in CBCT-based proton range verification and adaptive replanning of abdominal/pelvic proton therapy even when gas pockets are present in the beam path.

CBCT-Based Dose Monitoring and Adaptive Planning Triggers in Head and Neck PBS Proton Therapy

PURPOSE

To investigate the feasibility of using cone-beam computed tomography (CBCT)-derived synthetic CTs to monitor the daily dose and trigger a plan review for adaptive proton therapy (APT) in head and neck cancer (HNC) patients.

METHODS

For 84 HNC patients treated with proton pencil-beam scanning (PBS), same-day CBCT and verification CT (vfCT) pairs were retrospectively collected. The ground truth CT (gtCT) was created by deforming the vfCT to the same-day CBCT, and it was then used as a dosimetric baseline and for establishing plan review trigger recommendations. Two different synthetic CT algorithms were tested; the corrected CBCT (corrCBCT) was created using an iterative image correction method and the virtual CT (virtCT) was created by deforming the planning CT to the CBCT, followed by a low-density masking process. Clinical treatment plans were recalculated on the image sets for evaluation.

RESULTS

Plan review trigger criteria for adaptive therapy were established after closely reviewing the cohort data. Compared to the vfCT, the corrCBCT and virtCT reliably produced dosimetric data more similar to the gtCT. The average discrepancy in D99 for high-risk clinical target volumes (CTV) was 1.1%, 0.7%, and 0.4% and for standard-risk CTVs was 1.8%, 0.5%, and 0.5% for the vfCT, corrCBCT, and virtCT, respectively.

CONCLUSION

Streamlined APT has been achieved with the proposed plan review criteria and CBCT-based synthetic CT workflow.

ScatterNet for projection-based 4D cone-beam computed tomography intensity correction of lung cancer patients

Background and purpose: In radiotherapy, dose calculations based on 4D cone beam CTs (4DCBCTs) require image intensity corrections. This retrospective study compared the dose calculation accuracy of a deep learning, projection-based scatter correction workflow (ScatterNet), to slower workflows: conventional 4D projection-based scatter correction (CBCTcor) and a deformable image registration (DIR)-based method (4DvCT). Materials and methods: For 26 lung cancer patients, planning CTs (pCTs), 4DCTs and CBCT projections were available. ScatterNet was trained with pairs of raw and corrected CBCT projections. Corrected projections from ScatterNet and the conventional workflow were reconstructed using MA-ROOSTER, yielding 4DCBCTSN and 4DCBCTcor. The 4DvCT was generated by 4DCT to 4DCBCT DIR, as part of the 4DCBCTcor workflow. Robust intensity modulated proton therapy treatment plans were created on free-breathing pCTs. 4DCBCTSN was compared to 4DCBCTcor and the 4DvCT in terms of image quality and dose calculation accuracy (dose-volume-histogram parameters and 3 % /3 mm gamma analysis). Results: 4DCBCTSN resulted in an average mean absolute error of 87 HU and 102 HU when compared to 4DCBCTcor and 4DvCT respectively. High agreement was observed in targets with median dose differences of 0.4 Gy (4DCBCTSN-4DCBCTcor) and 0.3 Gy (4DCBCTSN-4DvCT). The gamma analysis showed high average 3 % /3 mm pass rates of 96 % for both 4DCBCTSN vs. 4DCBCTcor and 4DCBCTSN vs. 4DvCT. Conclusions: Accurate 4D dose calculations are feasible for lung cancer patients using ScatterNet for 4DCBCT correction. Average scatter correction times could be reduced from 10 min (4DCBCTcor) to 3.9 s , showing the clinical suitability of the proposed deep learning-based method.

New technique and application of truncated CBCT processing in adaptive radiotherapy for breast cancer

OBJECTIVE

A generative adversarial network (TCBCTNet) was proposed to generate synthetic computed tomography (sCT) from truncated low-dose cone-beam computed tomography (CBCT) and planning CT (pCT). The sCT was applied to the dose calculation of radiotherapy for patients with breast cancer.

METHODS

The low-dose CBCT and pCT images of 80 female thoracic patients were used for training. The CBCT, pCT, and replanning CT (rCT) images of 20 thoracic patients and 20 patients with breast cancer were used for testing. All patients were fixed in the same posture with a vacuum pad. The CBCT images were scanned under the Fast Chest M20 protocol with a 50% reduction in projection frames compared with the standard Chest M20 protocol. Rigid registration was performed between pCT and CBCT, and deformation registration was performed between rCT and CBCT. In the training stage of the TCBCTNet, truncated CBCT images obtained from complete CBCT images by simulation were used. The input of the CBCT→CT generator was truncated CBCT and pCT, and TCBCTNet was applied to patients with breast cancer after training. The accuracy of the sCT was evaluated by anatomy and dosimetry and compared with the generative adversarial network with UNet and ResNet as the generators (named as UnetGAN, ResGAN).

RESULTS

The three models could improve the image quality of CBCT and reduce the scattering artifacts while preserving the anatomical geometry of CBCT. For the chest test set, TCBCTNet achieved the best mean absolute error (MAE, 21.18±3.76 HU), better than 23.06±3.90 HU in UnetGAN and 22.47±3.57 HU in ResGAN. When applied to patients with breast cancer, TCBCTNet performance decreased, and MAE was 25.34±6.09 HU. Compared with rCT, sCT by TCBCTNet showed consistent dose distribution and subtle absolute dose differences between the target and the organ at risk. The 3D gamma pass rates were 98.98%±0.64% and 99.69%±0.22% at 2 mm/2% and 3 mm/3%, respectively. Ablation experiments confirmed that pCT and content loss played important roles in TCBCTNet.

CONCLUSIONS

High-quality sCT images could be synthesized from truncated low-dose CBCT and pCT by using the proposed TCBCTNet model. In addition, sCT could be used to accurately calculate the dose distribution for patients with breast cancer.

CBCT-to-CT Translation Using Registration-Based Generative Adversarial Networks in Patients with Head and Neck Cancer

Recently, deep learning with generative adversarial networks (GANs) has been applied in multi-domain image-to-image translation. This study aims to improve the image quality of cone-beam computed tomography (CBCT) by generating synthetic CT (sCT) that maintains the patient’s anatomy as in CBCT, while having the image quality of CT. As CBCT and CT are acquired at different time points, it is challenging to obtain paired images with aligned anatomy for supervised training. To address this limitation, the study incorporated a registration network (RegNet) into GAN during training. RegNet can dynamically estimate the correct labels, allowing supervised learning with noisy labels. The study developed and evaluated the approach using imaging data from 146 patients with head and neck cancer. The results showed that GAN trained with RegNet performed better than those trained without RegNet. Specifically, in the UNIT model trained with RegNet, the mean absolute error (MAE) was reduced from 40.46 to 37.21, the root mean-square error (RMSE) was reduced from 119.45 to 108.86, the peak signal-to-noise ratio (PSNR) was increased from 28.67 to 29.55, and the structural similarity index (SSIM) was increased from 0.8630 to 0.8791. The sCT generated from the model had fewer artifacts and retained the anatomical information as in CBCT.

A correlated sampling-based Monte Carlo simulation for fast CBCT iterative scatter correction

BACKGROUND

In recent years, cone-beam computed tomography (CBCT) has played an important role in medical imaging. However, the applications of CBCT are limited due to the severe scatter contamination. Conventional Monte Carlo (MC) simulation can provide accurate scatter estimation for scatter correction, but the expensive computational cost has always been the bottleneck of MC method in clinical application.

PURPOSE

In this work, an MC simulation method combined with a variance reduction technique called correlated sampling is proposed for fast iterative scatter correction.

METHODS

Correlated sampling exploits correlation between similar simulation systems to reduce the variance of interest quantities. Specifically, conventional MC simulation is first performed on the scatter-contaminated CBCT to generate the initial scatter signal. In the subsequent correction iterations, scatter estimation is then updated by applying correlated MC sampling to the latest corrected CBCT images by reusing the random number sequences of the task-related photons in conventional MC. Afterward, the corrected projections obtained by subtracting the scatter estimation from raw projections are utilized for FDK reconstruction. These steps are repeated until an adequate scatter correction is obtained. The performance of the proposed framework is evaluated by the accuracy of the scatter estimation, the quality of corrected CBCT images and efficiency.

RESULTS

Overall, the difference in mean absolute percentage error between scatter estimation with and without correlated sampling is 0.25% for full-fan case and 0.34% for half-fan case, respectively. In simulation studies, scatter artifacts are substantially eliminated, where the mean absolute error value is reduced from 15 to 2 HU in full-fan case and from 53 to 13 HU in half-fan case. Scatter-to-primary ratio is reduced to 0.02 for full-fan and 0.04 for half-fan, respectively. In phantom study, the contrast-to-noise ratio (CNR) is increased by a factor of 1.63, and the contrast is increased by a factor of 1.77. As for clinical studies, the CNR is improved by 11% and 14% for half-fan and full-fan, respectively. The contrast after correction is increased by 19% for half-fan and 44% for full-fan. Furthermore, root mean square error is also effectively reduced, especially from 78 to 4 HU for full-fan. Experimental results demonstrate that the figure of merit is improved between 23 and 43 folds when using correlated sampling. The proposed method takes less than 25 s for the whole iterative scatter correction process.

CONCLUSIONS

The proposed correlated sampling-based MC simulation method can achieve fast and accurate scatter correction for CBCT, making it suitable for real-time clinical use.

Deep learning-based 4D-synthetic CTs from sparse-view CBCTs for dose calculations in adaptive proton therapy

BACKGROUND

Time-resolved 4D cone beam-computed tomography (4D-CBCT) allows a daily assessment of patient anatomy and respiratory motion. However, 4D-CBCTs suffer from imaging artifacts that affect the CT number accuracy and prevent accurate proton dose calculations. Deep learning can be used to correct CT numbers and generate synthetic CTs (sCTs) that can enable CBCT-based proton dose calculations.

PURPOSE

In this work, sparse view 4D-CBCTs were converted into 4D-sCT utilizing a deep convolutional neural network (DCNN). 4D-sCTs were evaluated in terms of image quality and dosimetric accuracy to determine if accurate proton dose calculations for adaptive proton therapy workflows of lung cancer patients are feasible.

METHODS

A dataset of 45 thoracic cancer patients was utilized to train and evaluate a DCNN to generate 4D-sCTs, based on sparse view 4D-CBCTs reconstructed from projections acquired with a 3D acquisition protocol. Mean absolute error (MAE) and mean error were used as metrics to evaluate the image quality of single phases and average 4D-sCTs against 4D-CTs acquired on the same day. The dosimetric accuracy was checked globally (gamma analysis) and locally for target volumes and organs-at-risk (OARs) (lung, heart, and esophagus). Furthermore, 4D-sCTs were also compared to 3D-sCTs. To evaluate CT number accuracy, proton radiography simulations in 4D-sCT and 4D-CTs were compared in terms of range errors. The clinical suitability of 4D-sCTs was demonstrated by performing a 4D dose reconstruction using patient specific treatment delivery log files and breathing signals.

RESULTS

4D-sCTs resulted in average MAEs of 48.1 ± 6.5 HU (single phase) and 37.7 ± 6.2 HU (average). The global dosimetric evaluation showed gamma pass ratios of 92.3% ± 3.2% (single phase) and 94.4% ± 2.1% (average). The clinical target volume showed high agreement in D₉₈ between 4D-CT and 4D-sCT, with differences below 2.4% for all patients. Larger dose differences were observed in mean doses of OARs (up to 8.4%). The comparison with 3D-sCTs showed no substantial image quality and dosimetric differences for the 4D-sCT average. Individual 4D-sCT phases showed slightly lower dosimetric accuracy. The range error evaluation revealed that lung tissues cause range errors about three times higher than the other tissues.

CONCLUSION

In this study, we have investigated the accuracy of deep learning-based 4D-sCTs for daily dose calculations in adaptive proton therapy. Despite image quality differences between 4D-sCTs and 3D-sCTs, comparable dosimetric accuracy was observed globally and locally. Further improvement of 3D and 4D lung sCTs could be achieved by increasing CT number accuracy in lung tissues.

Improving cone-beam CT quality using a cycle-residual connection with a dilated convolution-consistent generative adversarial network

Objective.Cone-Beam CT (CBCT) often results in severe image artifacts and inaccurate HU values, meaning poor quality CBCT images cannot be directly applied to dose calculation in radiotherapy. To overcome this, we propose a cycle-residual connection with a dilated convolution-consistent generative adversarial network (Cycle-RCDC-GAN). Approach. The cycle-consistent generative adversarial network (Cycle-GAN) was modified using a dilated convolution with different expansion rates to extract richer semantic features from input images. Thirty pelvic patients were used to investigate the effect of synthetic CT (sCT) from CBCT, and 55 head and neck patients were used to explore the generalizability of the model. Three generalizability experiments were performed and compared: the pelvis trained model was applied to the head and neck; the head and neck trained model was applied to the pelvis, and the two datasets were trained together. Main results. The mean absolute error (MAE), the root mean square error (RMSE), peak signal to noise ratio (PSNR), the structural similarity index (SSIM), and spatial nonuniformity (SNU) assessed the quality of the sCT generated from CBCT. Compared with CBCT images, the MAE improved from 28.81 to 18.48, RMSE from 85.66 to 69.50, SNU from 0.34 to 0.30, and PSNR from 31.61 to 33.07, while SSIM improved from 0.981 to 0.989. The sCT objective indicators of Cycle-RCDC-GAN were better than Cycle-GAN’s. The objective metrics for generalizability were also better than Cycle-GAN’s. Significance. Cycle-RCDC-GAN enhances CBCT image quality and has better generalizability than Cycle-GAN, which further promotes the application of CBCT in radiotherapy.

Dosimetric evaluation of cone-beam CT-based synthetic CTs in pediatric patients undergoing intensity-modulated proton therapy

PURPOSE

To evaluate dosimetric changes detected using synthetic computed tomography (sCT) derived from online cone-beam CTs (CBCT) in pediatric patients treated using intensity-modulated proton therapy (IMPT).

METHODS

Ten pediatric patients undergoing IMPT and aligned daily using proton gantry-mounted CBCT were identified for retrospective analysis with treated anatomical sites fully encompassed in the CBCT field of view. Dates were identified when the patient received both a CBCT and a quality assurance CT (qCT) for routine dosimetric evaluation. sCTs were generated based on a deformable registration between the initial plan CT (pCT) and CBCT. The clinical IMPT plans were re-computed on the same day qCT and sCT, and dosimetric changes due to tissue change or response from the initial plan were computed using each image. Linear regression analysis was performed to determine the correlation between dosimetric changes detected using the qCT and the sCT. Gamma analysis was also used to compare the dose distributions computed on the qCT and sCT.

RESULTS

The correlation coefficients (p-values) between qCTs and sCTs for changes detected in target coverage, overall maximum dose, and organ at risk dose were 0.97 (< .001), 0.84 (.002) and 0.91 (< .001), respectively. Mean ± SD gamma pass rates of the sCT-based dose compared to the qCT-based dose at 3%/3 mm, 3%/2 mm, and 2%/2 mm criteria were 96.5%±4.5%, 93.2%±6.3%, and 91.3%±7.8%, respectively. Pass rates tended to be lower for targets near lung.

CONCLUSION

While insufficient for re-planning, sCTs provide approximate dosimetry without administering additional imaging dose in pediatric patients undergoing IMPT. Dosimetric changes detected using sCTs are correlated with changes detected using clinically-standard qCTs; however, residual differences in dosimetry remain a limitation. Further improvements in sCT image quality may both improve online dosimetric evaluation and reduce imaging dose for pediatric patients by reducing the need for routine qCTs.

A two-step method to improve image quality of CBCT with phantom-based supervised and patient-based unsupervised learning strategies

Objective.In this study, we aimed to develop deep learning framework to improve cone-beam computed tomography (CBCT) image quality for adaptive radiation therapy (ART) applications.Approach.Paired CBCT and planning CT images of 2 pelvic phantoms and 91 patients (15 patients for testing) diagnosed with prostate cancer were included in this study. First, well-matched images of rigid phantoms were used to train a U-net, which is the supervised learning strategy to reduce serious artifacts. Second, the phantom-trained U-net generated intermediate CT images from the patient CBCT images. Finally, a cycle-consistent generative adversarial network (CycleGAN) was trained with intermediate CT images and deformed planning CT images, which is the unsupervised learning strategy to learn the style of the patient images for further improvement. When testing or applying the trained model on patient CBCT images, the intermediate CT images were generated from the original CBCT image by U-net, and then the synthetic CT images were generated by the generator of CycleGAN with intermediate CT images as input. The performance was compared with conventional methods (U-net/CycleGAN alone trained with patient images) on the test set.Results.The proposed two-step method effectively improved the CBCT image quality to the level of CT scans. It outperformed conventional methods for region-of-interest contouring and HU calibration, which are important to ART applications. Compared with the U-net alone, it maintained the structure of CBCT. Compared with CycleGAN alone, our method improved the accuracy of CT number and effectively reduced the artifacts, making it more helpful for identifying the clinical target volume.Significance.This novel two-step method improves CBCT image quality by combining phantom-based supervised and patient-based unsupervised learning strategies. It has immense potential to be integrated into the ART workflow to improve radiotherapy accuracy.

Empirical scatter correction (ESC): CBCT scatter artifact reduction without prior information

BACKGROUND

The image quality of cone-beam CT (CBCT) scans severely suffers from scattered radiation if no countermeasures are taken. Scatter artifacts may induce cupping and streak artifacts and lead to a reduced image contrast and wrong CT values of the reconstructed volumes. Established software-based approaches for a correction of scattered radiation typically rely on prior knowledge of the CT system, scan parameters, the scanned object, or all of the aforementioned.

PURPOSE

This study proposes a simple and effective post-processing software-based correction method of scatter artifacts in CBCT scans without specific prior knowledge.

METHODS

We propose the empirical scatter correction (ESC) which generates scatter-like basis images from each projection image by convolution operations. A linear combination of these basis images is subtracted from the original projection image. The logarithm is taken and an FDK reconstruction is performed. The coefficients needed for the linear combination are determined automatically by a downhill simplex algorithm such that the resulting reconstructed images show no scatter artifacts. We demonstrate the potential of ESC by correcting simulated volumes with Monte Carlo scatter artifacts, a head phantom scan performed on our table-top CBCT, and a pelvis scan from a Varian Edge CBCT scanner.

RESULTS

ESC is able to improve the image quality of CBCT scans which is shown on the basis of our simulations and on measured data. For a simulated head CT, the CT value difference to the scatter-free reference image was as low as -6 HU after using ESC whereas the uncorrected data deviated by more than -200 HU from the reference data. Simulations of thorax and abdomen CT scans show that although scatter artifacts are not fully removed, anatomical features which were hard to discover prior to the correction become clearly visible and better segmentable with ESC. Similar results are obtained in the phantom measurement where a comparison to a slit scan of our head phantom shows only small differences. The CT values in soft tissue are improved in this measurement, as well. In soft tissue areas with severe scatter artifacts the CT values agree well with those of the slit scan (difference to slit scan: 35 HU corrected, -289 HU uncorrected). Scatter artifacts in measured patient data can also be reduced using the proposed empirical scatter correction. The results are comparable to those achieved with designated correction algorithms installed on the Varian Edge CBCT system.

CONCLUSIONS

ESC allows to reduce artifacts caused by patient scatter solely based on the projection data. This article is protected by copyright. All rights reserved.

Dosimetric assessment of patient dose calculation on a deep learning-based synthesized computed tomography image for adaptive radiotherapy

Purpose

Dose computation using cone beam computed tomography (CBCT) images is inaccurate for the purpose of adaptive treatment planning. The main goal of this study is to assess the dosimetric accuracy of synthetic computed tomography (CT)‐based calculation for adaptive planning in the upper abdominal region. We hypothesized that deep learning‐based synthetically generated CT images will produce comparable results to a deformed CT (CTdef) in terms of dose calculation, while displaying a more accurate representation of the daily anatomy and therefore superior dosimetric accuracy.

Methods

We have implemented a cycle‐consistent generative adversarial networks (CycleGANs) architecture to synthesize CT images from the daily acquired CBCT image with minimal error. CBCT and CT images from 17 liver stereotactic body radiation therapy (SBRT) patients were used to train, test, and validate the algorithm.

Results

The synthetically generated images showed increased signal‐to‐noise ratio, contrast resolution, and reduced root mean square error, mean absolute error, noise, and artifact severity. Superior edge matching, sharpness, and preservation of anatomical structures from the CBCT images were observed for the synthetic images when compared to the CTdef registration method. Three verification plans (CBCT, CTdef, and synthetic) were created from the original treatment plan and dose volume histogram (DVH) statistics were calculated. The synthetic‐based calculation shows comparatively similar results to the CTdef‐based calculation with a maximum mean deviation of 1.5%.

Conclusions

Our findings show that CycleGANs can produce reliable synthetic images for the adaptive delivery framework. Dose calculations can be performed on synthetic images with minimal error. Additionally, enhanced image quality should translate into better daily alignment, increasing treatment delivery accuracy.

Please Place Your Seat in the Full Upright Position: A Technical Framework for Landing Upright Radiation Therapy in the 21st Century

Delivering radiotherapy to patients in an upright position can allow for increased patient comfort, reduction in normal tissue irradiation, or reduction of machine size and complexity. This paper gives an overview of the requirements for the delivery of contemporary arc and modulated radiation therapy to upright patients. We explore i) patient positioning and immobilization, ii) simulation imaging, iii) treatment planning and iv) online setup and image guidance. Treatment chairs have been designed to reproducibly position seated patients for treatment and can be augmented by several existing immobilisation systems or promising emerging technologies such as soft robotics. There are few solutions for acquiring CT images for upright patients, however, cone beam computed tomography (CBCT) scans of upright patients can be produced using the imaging capabilities of standard Linacs combined with an additional patient rotation device. While these images will require corrections to make them appropriate for treatment planning, several methods indicate the viability of this approach. Treatment planning is largely unchanged apart from translating gantry rotation to patient rotation, allowing for a fixed beam with a patient rotating relative to it. Rotation can be provided by a turntable during treatment delivery. Imaging the patient with the same machinery as used in treatment could be advantageous for online plan adaption. While the current focus is using clinical linacs in existing facilities, developments in this area could also extend to lower-cost and mobile linacs and heavy ion therapy.

Management of Motion and Anatomical Variations in Charged Particle Therapy: Past, Present, and Into the Future

The major aim of radiation therapy is to provide curative or palliative treatment to cancerous malignancies while minimizing damage to healthy tissues. Charged particle radiotherapy utilizing carbon ions or protons is uniquely suited for this task due to its ability to achieve highly conformal dose distributions around the tumor volume. For these treatment modalities, uncertainties in the localization of patient anatomy due to inter- and intra-fractional motion present a heightened risk of undesired dose delivery. A diverse range of mitigation strategies have been developed and clinically implemented in various disease sites to monitor and correct for patient motion, but much work remains. This review provides an overview of current clinical practices for inter and intra-fractional motion management in charged particle therapy, including motion control, current imaging and motion tracking modalities, as well as treatment planning and delivery techniques. We also cover progress to date on emerging technologies including particle-based radiography imaging, novel treatment delivery methods such as tumor tracking and FLASH, and artificial intelligence and discuss their potential impact towards improving or increasing the challenge of motion mitigation in charged particle therapy.

X-ray CT adaptation based on a 2D-3D deformable image registration framework using simulated in-room proton radiographies

The aim of this work is to investigate in-room proton radiographies to compensate realistic rigid and non-rigid transformations in clinical-like scenarios based on 2D-3D deformable image registration (DIR) framework towards future clinical implementation of adaptive radiation therapy (ART). Monte Carlo simulations of proton radiographies (pRads) based on clinical x-ray CT of a head and neck, and a brain tumor patients are simulated for two different detector configurations (i.e. integration-mode and list-mode detectors) including high and low proton statistics. A realistic deformation, derived from cone beam CT of the patient, is applied to the treatment planning CT. Rigid inaccuracies in patient positioning are also applied and the effect of small, medium and large fields of view (FOVs) is investigated. A stopping criterion, as desirable in realistic scenarios devoid of ground truth proton CT (pCT), is proposed and investigated. Results show that rigid and non-rigid transformations can be compensated based on a limited number of low dose pRads. The root mean square error with respect to the pCT shows that the 2D-3D DIR of the treatment planning CT based on 10 pRads from integration-mode data and 2 pRads from list-mode data is capable of achieving comparable accuracy (∼90% and >90%, respectively) to conventional 3D-3D DIR. The dice similarity coefficient over the segmented regions of interest also verifies the improvement in accuracy prior to and after 2D-3D DIR. No relevant changes in accuracy are found between high and low proton statistics except for 2 pRads from integration-mode data. The impact of FOV size is negligible. The convergence of the metric adopted for the stopping criterion indicates the optimal convergence of the 2D-3D DIR. This work represents a further step towards the potential implementation of ART in proton therapy. Further computational optimization is however required to enable extensive clinical validation.

Onboard cone-beam CT-based replan evaluation for head and neck proton therapy

Purpose

Quality assurance computed tomography (QACT) is the current clinical practice in proton therapy to evaluate the needs for replan. QACT could falsely indicate replan because of setup issues that would be solved on the treatment machine. Deforming the treatment planning CT (TPCT) to the pretreatment CBCT may eliminate this issue. We investigated the performance of replan evaluation based on deformed TPCT (TPCTdir) for proton head and neck (H&N) therapy.

Methods and materials

Twenty‐eight H&N datasets along with pretreatment CBCT and QACT were used to validate the method. The changes in body volume were analyzed between the no‐replan and replan groups. The dose on the TPCTdir, the deformed QACT (QACTdir), and the QACT were calculated by applying the clinical plans to these image sets. Dosimetric parameters’ changes, including ΔD95, ΔDmean, and ΔD1 for the clinical target volumes (CTVs) were calculated. Receiver operating characteristic curves for replan evaluation based on ΔD95 on QACT and TPCTdir were calculated, using ΔD95 on QACTdir as the reference. A threshold for replan based on ΔD95 on TPCTdir is proposed. The specificities for the proposed method were calculated.

Results

The changes in the body contour were 95.8 ± 83.8 cc versus 305.0 ± 235.0 cc (p < 0.01) for the no‐replan and replan groups, respectively. The ΔD95, ΔDmean, and ΔD1 are all comparable for all the evaluations. The differences between TPCTdir and QACTdir evaluations were 0.30% ± 0.86%, 0.00 ± 0.22 Gy, and −0.17 ± 0.61 Gy for CTV ΔD95, ΔDmean, and ΔD1, respectively. The corresponding differences between the QACT and QACTdir were 0.12% ± 1.1%, 0.02 ± 0.32 Gy, and −0.01 ± 0.71 Gy. CTV ΔD95 > 2.6% in TPCTdir was chosen as the threshold to trigger QACT/replan. The corresponding specificity was 94% and 98% for the clinical practice and the proposed method, respectively.

Conclusions

The replan evaluation based on TPCTdir provides better specificity than that based on the QACT.

A Comparison Study Between CNN-Based Deformed Planning CT and CycleGAN-Based Synthetic CT Methods for Improving iCBCT Image Quality

Purpose

The aim of this study is to compare two methods for improving the image quality of the Varian Halcyon cone-beam CT (iCBCT) system through the deformed planning CT (dpCT) based on the convolutional neural network (CNN) and the synthetic CT (sCT) generation based on the cycle-consistent generative adversarial network (CycleGAN).

Methods

A total of 190 paired pelvic CT and iCBCT image datasets were included in the study, out of which 150 were used for model training and the remaining 40 were used for model testing. For the registration network, we proposed a 3D multi-stage registration network (MSnet) to deform planning CT images to agree with iCBCT images, and the contours from CT images were propagated to the corresponding iCBCT images through a deformation matrix. The overlap between the deformed contours (dpCT) and the fixed contours (iCBCT) was calculated for purposes of evaluating the registration accuracy. For the sCT generation, we trained the 2D CycleGAN using the deformation-registered CT-iCBCT slicers and generated the sCT with corresponding iCBCT image data. Then, on sCT images, physicians re-delineated the contours that were compared with contours of manually delineated iCBCT images. The organs for contour comparison included the bladder, spinal cord, femoral head left, femoral head right, and bone marrow. The dice similarity coefficient (DSC) was used to evaluate the accuracy of registration and the accuracy of sCT generation.

Results

The DSC values of the registration and sCT generation were found to be 0.769 and 0.884 for the bladder (p < 0.05), 0.765 and 0.850 for the spinal cord (p < 0.05), 0.918 and 0.923 for the femoral head left (p > 0.05), 0.916 and 0.921 for the femoral head right (p > 0.05), and 0.878 and 0.916 for the bone marrow (p < 0.05), respectively. When the bladder volume difference in planning CT and iCBCT scans was more than double, the accuracy of sCT generation was significantly better than that of registration (DSC of bladder: 0.859 vs. 0.596, p < 0.05).

Conclusion

The registration and sCT generation could both improve the iCBCT image quality effectively, and the sCT generation could achieve higher accuracy when the difference in planning CT and iCBCT was large.

CT-on-Rails Versus In-Room CBCT for Online Daily Adaptive Proton Therapy of Head-and-Neck Cancers

PURPOSE

To compare the efficacy of CT-on-rails versus in-room CBCT for daily adaptive proton therapy.

METHODS

We analyzed a cohort of ten head-and-neck patients with daily CBCT and corresponding virtual CT images. The necessity of moving the patient after a CT scan is the most significant difference in the adaptation workflow, leading to an increased treatment execution uncertainty σ. It is a combination of the isocenter-matching σ_i and random patient movements induced by the couch motion σ_m. The former is assumed to never exceed 1 mm. For the latter, we studied three different scenarios with σ_m = 1, 2, and 3 mm. Accordingly, to mimic the adaptation workflow with CT-on-rails, we introduced random offsets after Monte-Carlo-based adaptation but before delivery of the adapted plan.

RESULTS

There were no significant differences in accumulated dose-volume histograms and dose distributions for σ_m = 1 and 2 mm. Offsets with σ_m = 3 mm resulted in underdosage to CTV and hot spots of considerable volume.

CONCLUSION

Since σ_m typically does not exceed 2 mm for in-room CT, there is no clinically significant dosimetric difference between the two modalities for online adaptive therapy of head-and-neck patients. Therefore, in-room CT-on-rails can be considered a good alternative to CBCT for adaptive proton therapy.

Adaptive proton therapy

Radiation therapy treatments are typically planned based on a single image set, assuming that the patient's anatomy and its position relative to the delivery system remains constant during the course of treatment. Similarly, the prescription dose assumes constant biological dose-response over the treatment course. However, variations can and do occur on multiple time scales. For treatment sites with significant intra-fractional motion, geometric changes happen over seconds or minutes, while biological considerations change over days or weeks. At an intermediate timescale, geometric changes occur between daily treatment fractions. Adaptive radiation therapy is applied to consider changes in patient anatomy during the course of fractionated treatment delivery. While traditionally adaptation has been done off-line with replanning based on new CT images, online treatment adaptation based on on-board imaging has gained momentum in recent years due to advanced imaging techniques combined with treatment delivery systems. Adaptation is particularly important in proton therapy where small changes in patient anatomy can lead to significant dose perturbations due to the dose conformality and finite range of proton beams. This review summarizes the current state-of-the-art of on-line adaptive proton therapy and identifies areas requiring further research.

Planning CT-guided robust and fast cone-beam CT scatter correction using a local filtration technique

PURPOSE

Cone-beam CT (CBCT) has been widely utilized in image-guided radiotherapy. Planning CT (pCT)-aided CBCT scatter correction could further enhance image quality and extend CBCT application to dose calculation and adaptive planning. Nevertheless, existing pCT-based approaches demand accurate registration between pCT and CBCT, leading to limited imaging performance and increased computational cost when large anatomical discrepancies exist. In this work, we proposed a robust and fast CBCT scatter correction method using local filtration technique and rigid registration between pCT and CBCT (LF-RR).

METHODS

First of all, the pCT was rigidly registered with CBCT, then forward projection was performed on registered pCT to create scatter-free projections. The raw scatter signals were obtained via subtracting the scatter-free projections from the measured CBCT projections. Based on frequency and intensity threshold criteria, reliable scatter signals were selected from the raw scatter signals, and further filtered for global scatter estimation via local filtration technique. Finally, corrected CBCT was reconstructed with the projections generated by subtracting the scatter estimation from the raw CBCT projections using FDK algorithm. The LF-RR method was evaluated via comparison with another pCT-based scatter correction method based on Median and Gaussian filters (MG method).

RESULTS

Proposed method was first validated on an anthropomorphic pelvis phantom, and showed satisfied performance on scatter removal even when anatomical mismatches were intentionally created on pCT. The quantitative analysis was further performed on four clinical CBCT images. Compared with the uncorrected CBCT, CBCT corrected by MG with rigid registration (MG-RR), MG with deformable registration (MG-DR), and LF-RR reduced the CT number error from 79 ± 35 to 25 ± 18 , 17 ± 13 and 7 ± 3 HU for adipose and from 115 ± 61 to 36 ± 22 , 30 ± 24 , 7 ± 3 HU for muscle, respectively. After correction, the spatial non-uniformity (SNU) of CBCT corrected with MG-RR, MG-DR and LF-RR was 51 ± 13 , 60 ± 21 , and 21 ± 9 HU for adipose, and 50 ± 22 , 57 ± 41 , and 25 ± 6 HU for muscle, respectively. Meanwhile, the contrast-to-noise ratio (CNR) between muscle and adipose was increased by a factor of 2.70, 2.89 and 2.56, respectively. Using the LF-RR method, the scatter correction of 656 projections can be finished within 10 s and the corrected volumetric images (200 slices) can be obtained within 2 min.

CONCLUSION

We developed a fast and robust pCT-based CBCT scatter correction method which exploits the local-filtration technique to promote the accuracy of scatter estimation and is resistant to pCT-to-CBCT registration uncertainties. Both phantom and patient studies showed the superiority of the proposed correction in imaging accuracy and computational efficiency, indicating promisingfuture clinical application.

Generating synthetic CT from low-dose cone-beam CT by using generative adversarial networks for adaptive radiotherapy

OBJECTIVE

To develop high-quality synthetic CT (sCT) generation method from low-dose cone-beam CT (CBCT) images by using attention-guided generative adversarial networks (AGGAN) and apply these images to dose calculations in radiotherapy.

METHODS

The CBCT/planning CT images of 170 patients undergoing thoracic radiotherapy were used for training and testing. The CBCT images were scanned under a fast protocol with 50% less clinical projection frames compared with standard chest M20 protocol. Training with aligned paired images was performed using conditional adversarial networks (so-called pix2pix), and training with unpaired images was carried out with cycle-consistent adversarial networks (cycleGAN) and AGGAN, through which sCT images were generated. The image quality and Hounsfield unit (HU) value of the sCT images generated by the three neural networks were compared. The treatment plan was designed on CT and copied to sCT images to calculated dose distribution.

RESULTS

The image quality of sCT images by all the three methods are significantly improved compared with original CBCT images. The AGGAN achieves the best image quality in the testing patients with the smallest mean absolute error (MAE, 43.5 ± 6.69), largest structural similarity (SSIM, 93.7 ± 3.88) and peak signal-to-noise ratio (PSNR, 29.5 ± 2.36). The sCT images generated by all the three methods showed superior dose calculation accuracy with higher gamma passing rates compared with original CBCT image. The AGGAN offered the highest gamma passing rates (91.4 ± 3.26) under the strictest criteria of 1 mm/1% compared with other methods. In the phantom study, the sCT images generated by AGGAN demonstrated the best image quality and the highest dose calculation accuracy.

CONCLUSIONS

High-quality sCT images were generated from low-dose thoracic CBCT images by using the proposed AGGAN through unpaired CBCT and CT images. The dose distribution could be calculated accurately based on sCT images in radiotherapy.

Image-based shading correction for narrow-FOV truncated pelvic CBCT with deep convolutional neural networks and transfer learning

Purpose: Cone beam computed tomography (CBCT) is a standard solution for in‐room image guidance for radiation therapy. It is used to evaluate and compensate for anatomopathological changes between the dose delivery plan and the fraction delivery day. CBCT is a fast and versatile solution, but it suffers from drawbacks like low contrast and requires proper calibration to derive density values. Although these limitations are even more prominent with in‐room customized CBCT systems, strategies based on deep learning have shown potential in improving image quality. As such, this article presents a method based on a convolutional neural network and a novel two‐step supervised training based on the transfer learning paradigm for shading correction in CBCT volumes with narrow field of view (FOV) acquired with an ad hoc in‐room system.

Methods: We designed a U‐Net convolutional neural network, trained on axial slices of corresponding CT/CBCT couples. To improve the generalization capability of the network, we exploited two‐stage learning using two distinct data sets. At first, the network weights were trained using synthetic CBCT scans generated from a public data set, and then only the deepest layers of the network were trained again with real‐world clinical data to fine‐tune the weights. Synthetic data were generated according to real data acquisition parameters. The network takes a single grayscale volume as input and outputs the same volume with corrected shading and improved HU values.

Results: Evaluation was carried out with a leave‐one‐out cross‐validation, computed on 18 unique CT/CBCT pairs from six different patients from a real‐world dataset. Comparing original CBCT to CT and improved CBCT to CT, we obtained an average improvement of 6 dB on peak signal‐to‐noise ratio (PSNR), +2% on structural similarity index measure (SSIM). The median interquartile range (IQR) Hounsfield unit (HU) difference between CBCT and CT improved from 161.37 (162.54) HU to 49.41 (66.70) HU. Region of interest (ROI)‐based HU difference was narrowed by 75% in the spongy bone (femoral head), 89% in the bladder, 85% for fat, and 83% for muscle. The improvement in contrast‐to‐noise ratio for these ROIs was about 67%.

Conclusions: We demonstrated that shading correction obtaining CT‐compatible data from narrow‐FOV CBCTs acquired with a customized in‐room system is possible. Moreover, the transfer learning approach proved particularly beneficial for such a shading correction approach.

Experimental validation of daily adaptive proton therapy

Anatomical changes during proton therapy require rapid treatment plan adaption to mitigate the associated dosimetric impact. This in turn requires a highly efficient workflow that minimizes the time between imaging and delivery. At the Paul Scherrer Institute, we have developed an online adaptive workflow, which is specifically designed for treatments in the skull-base/cranium, with the focus set on simplicity and minimizing changes to the conventional workflow. The dosimetric and timing performance of this daily adaptive proton therapy (DAPT) workflow has been experimentally investigated using an in-house developed DAPT software and specifically developed anthropomorphic phantom. After a standard treatment preparation, which includes the generation of a template plan, the treatment can then be adapted each day, based on daily imaging acquired on an in-room CT. The template structures are then rigidly propagated to this CT and the daily plan is fully re-optimized using the same field arrangement, DVH constraints and optimization settings of the template plan. After a dedicated plan QA, the daily plan is delivered. To minimize the time between imaging and delivery, clinically integrated software for efficient execution of all online adaption steps, as well as tools for comprehensive and automated QA checks, have been developed. Film measurements of an end-to-end validation of a multi-fraction DAPT treatment showed high agreement to the calculated doses. Gamma pass rates with a 3%/3 mm criteria were >92% when comparing the measured dose to the template plan. Additionally, a gamma pass rate >99% was found comparing measurements to the Monte Carlo dose of the daily plans reconstructed from the logfile, accumulated over the delivered fractions. With this, we experimentally demonstrate that the described adaptive workflow can be delivered accurately in a timescale similar to a standard delivery.

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 Supervised and Unsupervised Approaches for the Generation of Synthetic CT from Cone-Beam CT

Due to major artifacts and uncalibrated Hounsfield units (HU), cone-beam computed tomography (CBCT) cannot be used readily for diagnostics and therapy planning purposes. This study addresses image-to-image translation by convolutional neural networks (CNNs) to convert CBCT to CT-like scans, comparing supervised to unsupervised training techniques, exploiting a pelvic CT/CBCT publicly available dataset. Interestingly, quantitative results were in favor of supervised against unsupervised approach showing improvements in the HU accuracy (62% vs. 50%), structural similarity index (2.5% vs. 1.1%) and peak signal-to-noise ratio (15% vs. 8%). Qualitative results conversely showcased higher anatomical artifacts in the synthetic CBCT generated by the supervised techniques. This was motivated by the higher sensitivity of the supervised training technique to the pixel-wise correspondence contained in the loss function. The unsupervised technique does not require correspondence and mitigates this drawback as it combines adversarial, cycle consistency, and identity loss functions. Overall, two main impacts qualify the paper: (a) the feasibility of CNN to generate accurate synthetic CT from CBCT images, which is fast and easy to use compared to traditional techniques applied in clinics; (b) the proposal of guidelines to drive the selection of the better training technique, which can be shifted to more general image-to-image translation.

A Deep Unsupervised Learning Model for Artifact Correction of Pelvis Cone-Beam CT

Purpose

In recent years, cone-beam computed tomography (CBCT) is increasingly used in adaptive radiation therapy (ART). However, compared with planning computed tomography (PCT), CBCT image has much more noise and imaging artifacts. Therefore, it is necessary to improve the image quality and HU accuracy of CBCT. In this study, we developed an unsupervised deep learning network (CycleGAN) model to calibrate CBCT images for the pelvis to extend potential clinical applications in CBCT-guided ART.

Methods

To train CycleGAN to generate synthetic PCT (sPCT), we used CBCT and PCT images as inputs from 49 patients with unpaired data. Additional deformed PCT (dPCT) images attained as CBCT after deformable registration are utilized as the ground truth before evaluation. The trained uncorrected CBCT images are converted into sPCT images, and the obtained sPCT images have the characteristics of PCT images while keeping the anatomical structure of CBCT images unchanged. To demonstrate the effectiveness of the proposed CycleGAN, we use additional nine independent patients for testing.

Results

We compared the sPCT with dPCT images as the ground truth. The average mean absolute error (MAE) of the whole image on testing data decreased from 49.96 ± 7.21HU to 14.6 ± 2.39HU, the average MAE of fat and muscle ROIs decreased from 60.23 ± 7.3HU to 16.94 ± 7.5HU, and from 53.16 ± 9.1HU to 13.03 ± 2.63HU respectively.

Conclusion

We developed an unsupervised learning method to generate high-quality corrected CBCT images (sPCT). Through further evaluation and clinical implementation, it can replace CBCT in ART.

Range probing as a quality control tool for CBCT-based synthetic CTs: In vivo application for head and neck cancer patients

Purpose

Cone‐beam CT (CBCT)‐based synthetic CTs (sCT) produced with a deep convolutional neural network (DCNN) show high image quality, suggesting their potential usability in adaptive proton therapy workflows. However, the nature of such workflows involving DCNNs prevents the user from having direct control over their output. Therefore, quality control (QC) tools that monitor the sCTs and detect failures or outliers in the generated images are needed.

This work evaluates the potential of using a range‐probing (RP)‐based QC tool to verify sCTs generated by a DCNN. Such a RP QC tool experimentally assesses the CT number accuracy in sCTs.

Methods

A RP QC dataset consisting of repeat CTs (rCT), CBCTs, and RP acquisitions of seven head and neck cancer patients was retrospectively assessed. CBCT‐based sCTs were generated using a DCNN. The CT number accuracy in the sCTs was evaluated by computing relative range errors between measured RP fields and RP field simulations based on rCT and sCT images.

Results

Mean relative range errors showed agreement between measured and simulated RP fields, ranging from −1.2% to 1.5% in rCTs, and from −0.7% to 2.7% in sCTs.

Conclusions

The agreement between measured and simulated RP fields suggests the suitability of sCTs for proton dose calculations. This outcome brings sCTs generated by DCNNs closer toward clinical implementation within adaptive proton therapy treatment workflows. The proposed RP QC tool allows for CT number accuracy assessment in sCTs and can provide means of in vivo range verification.

Measurement-based range evaluation for quality assurance of CBCT-based dose calculations in adaptive proton therapy

PURPOSE

The implementation of volumetric in-room imaging for online adaptive radiotherapy makes extensive testing of this image data for treatment planning necessary. Especially for proton beams the higher sensitivity to stopping power properties of the tissue results in more stringent requirements. Current approaches mainly focus on recalculation of the plans on the new image data, lacking experimental verification, and ignoring the impact on the plan re-optimization process. The aim of this study was to use gel and film dosimetry coupled with a three-dimensional (3D) printed head phantom (based on the planning CT of the patient) for 3D range verification of intensity-corrected cone beam computed tomography (CBCT) image data for adaptive proton therapy.

METHODS

Single field uniform dose pencil beam scanning proton plans were optimized for three different patients on the patients' planning CT (planCT) and the patients' intensity-corrected CBCT (scCBCT) for the same target volume using the same optimization constraints. The CBCTs were corrected on projection level using the planCT as a prior. The dose optimized on planCT and recalculated on scCBCT was compared in terms of proton range differences (80% distal fall-off, recalculation). Moreover, the dose distribution resulting from recalculation of the scCBCT-optimized plan on the planCT and the original planCT dose distribution were compared (simulation). Finally, the two plans of each patient were irradiated on the corresponding patient-specific 3D printed head phantom using gel dosimetry inserts for one patient and film dosimetry for all three patients. Range differences were extracted from the measured dose distributions. The measured and the simulated range differences were corrected for range differences originating from the initial plans and evaluated.

RESULTS

The simulation approach showed high agreement with the standard recalculation approach. The median values of the range differences of these two methods agreed within 0.1 mm and the interquartile ranges (IQRs) within 0.3 mm for all three patients. The range differences of the film measurement were accurately matching with the simulation approach in the film plane. The median values of these range differences deviated less than 0.1 mm and the IQRs less than 0.4 mm. For the full 3D evaluation of the gel range differences, the median value and IQR matched those of the simulation approach within 0.7 and 0.5 mm, respectively. scCBCT- and planCT-based dose distributions were found to have a range agreement better than 3 mm (median and IQR) for all considered scenarios (recalculation, simulation, and measurement).

CONCLUSIONS

The results of this initial study indicate that an online adaptive proton workflow based on scatter-corrected CBCT image data for head irradiations is feasible. The novel presented measurement- and simulation-based method was shown to be equivalent to the standard literature recalculation approach. Additionally, it has the capability to catch effects of image differences on the treatment plan optimization. This makes the measurement-based approach particularly interesting for quality assurance of CBCT-based online adaptive proton therapy. The observed uncertainties could be kept within those of the registration and positioning. The proposed validation could also be applied for other alternative in-room images, e.g. for MR-based pseudoCTs.

Cone beam CT based validation of neural network generated synthetic CTs for radiotherapy in the head region

PURPOSE

In the past years, many different neural network-based conversion techniques for synthesizing computed tomographys (sCTs) from MR images have been published. While the model's performance can be checked during the training against the test set, test datasets can never represent the whole population. Conversion errors can still occur for special cases, for example, for unusual anatomical situations. Therefore, the performance of sCT conversion needs to be verified on a patient specific level, especially in the absence of a planning CT (pCT). In this study, the capability of cone-beam CTs (CBCTs) for the validation of sCTs generated by a neural network was investigated.

METHODS

41 patients with tumors in the head region were selected. 20 of them were used for model training and 10 for validation. Different implementations of CycleGAN (with/without identity and feature loss) were used to generate sCTs. The pixel (MAE, RMSE, PSNR) and geometric error (DICE, Sensitivity, Specificity) values were reported to identify the best model. VMAT plans were created for the remaining 11 patients on the pCTs. These plans were re-calculated on sCTs and CBCTs. An automatic density overriding method ( C B C T RS ) and a population-based dose calculation method ( C B C T Pop ) were employed for CBCT-based dose calculation. The dose distributions were analysed using 3D global gamma analysis, applying a threshold of 10% with respect to the prescribed dose. Differences in DVH metrics for the PTV and the organs-at-risk were compared among the dose distributions based on pCTs, sCTs, and CBCTs.

RESULTS

The best model was the CycleGAN without identity and feature matching loss. Including the identity loss led to a metric decrease of 10% for DICE and a metric increase of 20-60 HU for MAE. Using the 2%/2 mm gamma criterion and pCT as reference, the mean gamma pass rates were 99.0  ±  0.4% for sCTs. Mean gamma pass rate values comparing pCT and CBCT were 99.0  ±  0.8% and 99.1  ±  0.8% for the C B C T RS and C B C T Pop , respectively. The mean gamma pass rates comparing sCT and CBCT resulted in 98.4  ±  1.6% and 99.2  ±  0.6% for C B C T RS and C B C T Pop , respectively. The differences between the gamma-pass-rates of the sCT and two CBCT-based methods were not significant. The majority of deviations of the investigated DVH metrices between sCTs and CBCTs were within 2%.

CONCLUSION

The dosimetric results demonstrate good agreement between sCT, CBCT, and pCT based calculations. A properly applied CBCT conversion method can serve as a tool for quality assurance procedures in an MR only radiotherapy workflow for head patients. Dosimetric deviations of DVH metrics between sCT and CBCTs of larger than 2% should be followed up. A systematic shift of approximately 1% should be taken into account when using the C B C T RS approach in an MR only workflow.

Accuracy of dose calculation on iterative CBCT for head and neck radiotherapy

PURPOSE

To evaluate the feasibility of the use of iterative cone-beam computed tomography (CBCT) for dose calculation in the head and neck region.

METHODS

This study includes phantom and clinical studies. All acquired CBCT images were reconstructed with Feldkamp-Davis-Kress algorithm-based CBCT (FDK-CBCT) and iterative CBCT (iCBCT) algorithm. The Hounsfield unit (HU) consistency between the head and body phantoms was determined in both reconstruction techniques. Volumetric modulated arc therapy (VMAT) plans were generated for 16 head and neck patients on a planning CT scan, and the doses were recalculated on FDK-CBCT and iCBCT with Anisotropic Analytical Algorithm (AAA) and Acuros XB (AXB). As a comparison of the accuracy of dose calculations, the absolute dosimetric difference and 1%/1 mm gamma passing rate analysis were analyzed.

RESULTS

The difference in the mean HU values between the head and body phantoms was larger for FDK-CBCT (max value: 449.1 HU) than iCBCT (260.0 HU). The median dosimetric difference from the planning CT were <1.0% for both FDK-CBCT and iCBCT but smaller differences were found with iCBCT (planning target volume D_50%: 0.38% (0.15-0.59%) for FDK-CBCT, 0.28% (0.13-0.49%) for iCBCT, AAA; 0.14% (0.04-0.19%) for FDK-CBCT, 0.07% (0.02-0.20%) for iCBCT). The mean gamma passing rate was significantly better in iCBCT than FDK-CBCT (AAA: 98.7% for FDK-CBCT, 99.4% for iCBCT; AXB: 96.8% for FDK_CBCT, 97.5% for iCBCT).

CONCLUSION

The iCBCT-based dose calculation in VMAT for head and neck cancer was accurate compared to FDK-CBCT.

CBCT-Based Adaptive Assessment Workflow for Intensity Modulated Proton Therapy for Head and Neck Cancer

Purpose

Anatomical changes and patient setup uncertainties during intensity modulated proton therapy (IMPT) of head and neck (HN) cancers demand frequent evaluation of delivered dose. This work investigated a cone-beam computed tomography (CBCT) and deformable image registration based therapy workflow to demonstrate the feasibility of proton dose calculation on synthetic computed tomography (sCT) for adaptive IMPT treatment of HN cancer.

Materials and Methods

Twenty-one patients with HN cancer were enrolled in this study, a retrospective institutional review board protocol. They had previously been treated with volumetric modulated arc therapy and had daily iterative CBCT. For each patient, robust optimization (RO) IMPT plans were generated using ±3 mm patient setup and ±3% proton range uncertainties. The sCTs were created and the weekly delivered dose was recalculated using an adaptive dose accumulation workflow in which the planning computed tomography (CT) was deformably registered to CBCTs and Hounsfield units transferred from the planning CT. Accumulated doses from ±3 mm/±3% RO-IMPT plans were evaluated using clinical dose-volume constraints for targets (clinical target volume, or CTV) and organs at risk.

Results

Evaluation of weekly recalculated dose on sCTs showed that most of the patient plans maintained target dose coverage. The primary CTV remained covered by the V95 > 95% (95% of the volume receiving more than 95% of the prescription dose) worst-case scenario for 84.5% of the weekly fractions. The oral cavity accumulated mean dose remained lower than the worst-case scenario for all patients. Parotid accumulated mean dose remained within the uncertainty bands for 18 of the 21 patients, and all were kept lower than RO-IMPT worst-case scenario for 88.7% and 84.5% for left and right parotids, respectively.

Conclusion

This study demonstrated that RO-IMPT plans account for most setup and anatomical uncertainties, except for large weight-loss changes that need to be tracked throughout the treatment course. We showed that sCTs could be a powerful decision tool for adaptation of these cases in order to reduce workload when using repeat CTs.