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.

PTCOG Gastrointestinal Subcommittee Lower Gastrointestinal Tract Malignancies Consensus Statement

Purpose

Radiotherapy delivery in the definitive management of lower gastrointestinal (LGI) tract malignancies is associated with substantial risk of acute and late gastrointestinal (GI), genitourinary, dermatologic, and hematologic toxicities. Advanced radiation therapy techniques such as proton beam therapy (PBT) offer optimal dosimetric sparing of critical organs at risk, achieving a more favorable therapeutic ratio compared with photon therapy.

Materials and Methods

The international Particle Therapy Cooperative Group GI Subcommittee conducted a systematic literature review, from which consensus recommendations were developed on the application of PBT for LGI malignancies.

Results

Eleven recommendations on clinical indications for which PBT should be considered are presented with supporting literature, and each recommendation was assessed for level of evidence and strength of recommendation. Detailed technical guidelines pertaining to simulation, treatment planning and delivery, and image guidance are also provided.

Conclusion

PBT may be of significant value in select patients with LGI malignancies. Additional clinical data are needed to further elucidate the potential benefits of PBT for patients with anal cancer and rectal cancer.

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.

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.

Monitoring of Interfractional Proton Range Verification and Dosimetric Impact Based on Daily CBCT for Pediatric Patients with Pelvic Tumors

(1) Background: Synthetic CT images of the pelvis were generated from daily CBCT images to monitor changes in water equivalent path length (WEPL) and determine the dosimetric impact of anatomy changes along the proton beam's path; (2) Methods: Ten pediatric patients with pelvic tumors treated using proton therapy with daily CBCT were included. The original planning CT was deformed to the same-day CBCT to generate synthetic CT images for WEPL comparison and dosimetric evaluation; (3) Results: WEPL changes of 20 proton fields at the distal edge of the CTV ranged from 0.1 to 12 mm with a median of 2.5 mm, and 75th percentile of 5.1 mm for (the original CT-rescanned CT) and ranged from 0.3 to 10.1 mm with a median of 2.45 mm and 75th percentile of 4.8 mm for (the original CT-synthetic CT). The dosimetric impact was due to proton range pullback or overshoot, which led to reduced coverage in CTV Dmin averaging 12.1% and 11.3% in the rescanned and synthetic CT verification plans, respectively; (4) Conclusions: The study demonstrated that synthetic CT generated by deforming the original planning CT to daily CBCT can be used to quantify proton range changes and predict adverse dosimetric scenarios without the need for excessive rescanned CT scans during large interfractional variations in adaptive proton therapy of pediatric pelvic tumors.

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.

Evaluating synthetic computed tomography images for adaptive radiotherapy decision making in head and neck cancer

Background and purpose

Adaptive radiotherapy (ART) decision-making benefits from dosimetric information to supplement image inspection when assessing the significance of anatomical changes. This study evaluated a dosimetry-based clinical decision workflow for ART utilizing deformable registration of the original planning computed tomography (CT) image to the daily Cone Beam CT (CBCT) to replace the need for a replan CT for dose estimation.

Materials and methods

We used 12 retrospective Head & Neck patient cases having a ground truth - a replan CT (rCT) in response to anatomical changes apparent in the daily CBCT - to evaluate the accuracy of dosimetric assessment conducted on synthetic CTs (sCT) generated by deforming the original planning CT Hounsfield Units to the daily CBCT anatomy.The original plan was applied to the sCT and dosimetric accuracy of the sCT was assessed by analyzing plan objectives for targets and organs-at-risk compared to calculations on the ground-truth rCT. Three commercial DIR algorithms were compared.

Results

For the best-performing algorithms, the majority of dose metrics calculated on the sCTs differed by less than 4 Gy (5.7% of 70 Gy prescription dose). An uncertainty of ±2.5 Gy (3.6% of 70 Gy prescription) is recommended as a conservative tolerance when evaluating dose metrics on sCTs for head and neck.

Conclusions

Synthetic CTs present a valuable addition to the adaptive radiotherapy workflow, and synthetic CT dose estimates can be effectively used in addition to the current practice of visually inspecting the overlay of the planning CT and CBCT to assess the significance of anatomical change.

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.

Evaluation of Dose Calculation Based on Cone-Beam CT Using Different Measuring Correction Methods for Head and Neck Cancer Patients

Purpose: To investigate and compare 2 cone-beam computed tomography (CBCT) correction methods for CBCT-based dose calculation. Materials and Methods: Routine CBCT image sets of 12 head and neck cancer patients who received volumetric modulated arc therapy (VMAT) treatment were retrospectively analyzed. The CBCT images obtained using an on-board imager (OBI) at the first treatment fraction were firstly deformable registered and padded with the kVCT images to provide enough anatomical information about the tissues for dose calculation. Then, 2 CBCT correction methods were developed and applied to correct CBCT Hounsfield unit (HU) values. One method (HD method) is based on protocol-specific CBCT HU to physical density (HD) curve, and the other method (HM method) is based on histogram matching (HM) of HU value. The corrected CBCT images (CBCT_HD and CBCT_HM for HD and HM methods) were imported into the original planning system for dose calculation based on the HD curve of kVCT (the planning CT). The dose computation result was analyzed and discussed to compare these 2 CBCT-correction methods. Results: Dosimetric parameters, such as the D_mean, D_max and D_5% of the target volume in CBCT plan doses, were higher than those in the kVCT plan doses; however, the deviations were less than 2%. The D_2%, in parallel organs such as the parotid glands, the deviations from the CBCT_HM plan dose were less than those of the CBCT_HD plan dose. The differences were statistically significant (P < .05). Meanwhile, the V₃₀ value based on the HM method was better than that based on the HD method in the oral cavity region (P = .016). In addition, we also compared the γ passing rates of kVCT plan doses with the 2 CBCT plan doses, and negligible differences were found. Conclusion: The HM method was more suitable for head and neck cancer patients than the HD one. Furthermore, with the CBCT_HM-based method, the dose calculation result better matches the kVCT-based dose calculation.

Evaluating Proton Dose and Associated Range Uncertainty Using Daily Cone-Beam CT

Purpose

This study aimed to quantitatively evaluate the range uncertainties that arise from daily cone-beam CT (CBCT) images for proton dose calculation compared to CT using a measurement-based technique.

Methods

For head and thorax phantoms, wedge-shaped intensity-modulated proton therapy (IMPT) treatment plans were created such that the gradient of the wedge intersected and was measured with a 2D ion chamber array. The measured 2D dose distributions were compared with 2D dose planes extracted from the dose distributions using the IMPT plan calculated on CT and CBCT. Treatment plans of a thymoma cancer patient treated with breath-hold (BH) IMPT were recalculated on 28 CBCTs and 9 CTs, and the resulting dose distributions were compared.

Results

The range uncertainties for the head phantom were determined to be 1.2% with CBCT, compared to 0.5% for CT, whereas the range uncertainties for the thorax phantom were 2.1% with CBCT, compared to 0.8% for CT. The doses calculated on CBCT and CT were similar with similar anatomy changes. For the thymoma patient, the primary source of anatomy change was the BH uncertainty, which could be up to 8 mm in the superior-inferior (SI) direction.

Conclusion

We developed a measurement-based range uncertainty evaluation method with high sensitivity and used it to validate the accuracy of CBCT-based range and dose calculation. Our study demonstrated that the CBCT-based dose calculation could be used for daily dose validation in selected proton patients.

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.

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.

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

Purpose

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

Methods

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

Results

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

Conclusion

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

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.

Deep learning based synthetic-CT generation in radiotherapy and PET: A review

Recently,deep learning (DL)-based methods for the generation of synthetic computed tomography (sCT) have received significant research attention as an alternative to classical ones. We present here a systematic review of these methods by grouping them into three categories, according to their clinical applications: (i) to replace computed tomography in magnetic resonance (MR) based treatment planning, (ii) facilitate cone-beam computed tomography based image-guided adaptive radiotherapy, and (iii) derive attenuation maps for the correction of positron emission tomography. Appropriate database searching was performed on journal articles published between January 2014 and December 2020. The DL methods' key characteristics were extracted from each eligible study, and a comprehensive comparison among network architectures and metrics was reported. A detailed review of each category was given, highlighting essential contributions, identifying specific challenges, and summarizing the achievements. Lastly, the statistics of all the cited works from various aspects were analyzed, revealing the popularity and future trends and the potential of DL-based sCT generation. The current status of DL-based sCT generation was evaluated, assessing the clinical readiness of the presented methods.

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.

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.

Comparison of weekly and daily online adaptation for head and neck intensity-modulated proton therapy

The high conformality of intensity-modulated proton therapy (IMPT) dose distributions causes treatment plans to be sensitive to geometrical changes during the course of a fractionated treatment. This can be addressed using adaptive proton therapy (APT). One important question in APT is the frequency of adaptations performed during a fractionated treatment, which is related to the question whether plan adaptation has to be done online or offline. The purpose of this work is to investigate the impact of weekly and daily online IMPT plan adaptation on the treatment quality for head and neck patients. A cohort of ten head and neck patients with daily acquired cone-beam CT (CBCT) images was evaluated retrospectively. Dose tracking of the IMPT treatment was performed for three scenarios: base plan with no adaptation (BP), weekly online adaptation (OA_W), and daily online adaptation (OA_D). Both adaptation schemes used an in-house developed online APT workflow, performing Monte Carlo (MC) dose calculations on scatter-corrected CBCTs. IMPT plan adaptation was achieved by only tuning the weights of a subset of beamlets, based on deformable image registration from the planning CT to each CBCT. Although OA_Dmitigated random delivery errors more effectively than OA_Won a fraction per fraction basis, both OA_Wand OA_Dachieved the clinical goals for all ten patients, while BP failed for six cases. In the high-risk CTV, accumulated values of D_98%ranged between 97.15% and 99.73% of the prescription dose for OA_D, with a median of 98.07%. For OA_W, values between 95.02% and 99.26% were obtained, with a median of 97.61% of the prescription dose. Otherwise, the dose to most organs at risk was similar for all three scenarios. Globally, our results suggest that OA_Wcould be used as an alternative approach to OA_Dfor most patients in order to reduce the clinical workload.

Assessment of daily dose accumulation for robustly optimized intensity modulated proton therapy treatment of prostate cancer

PURPOSE

To implement a daily CBCT based dose accumulation technique in order to assess ideal robust optimization (RO) parameters for IMPT treatment of prostate cancer.

METHODS

Ten prostate cancer patients previously treated with VMAT and having daily CBCT were included. First, RO-IMPT plans were created with ± 3 mm and ± 5 mm patient setup and ± 3% proton range uncertainties, respectively. Second, the planning CT (pCT) was deformably registered to the CBCT to create a synthetic CT (sCT). Both daily and weekly sampling strategies were employed to determine optimal dose accumulation frequency. Doses were recalculated on sCTs for both ± 3 mm/±3% and ± 5 mm/±3% uncertainties and were accumulated back to the pCT. Accumulated doses generated from ± 3 mm/±3% and ± 5 mm/±3% RO-IMPT plans were evaluated using the clinical dose volume constraints for CTV, bladder, and rectum.

RESULTS

Daily accumulated dose based on both ± 3mm/±3% and ±5 mm/±3% uncertainties for RO-IMPT plans resulted in satisfactory CTV coverage (RO-IMPT_3mm/3% CTV_V95 = 99.01 ± 0.87% vs. RO-IMPT_5mm/3% CTV_V95 = 99.81 ± 0.2%, P = 0.002). However, the accumulated dose based on ± 3 mm/3% RO-IMPT plans consistently provided greater OAR sparing than ±5 mm/±3% RO-IMPT plans (RO-IMPT_3mm/3% rectum_V65Gy = 2.93 ± 2.39% vs. RO-IMPT_5mm/3% rectum_V65Gy = 4.38 ± 3%, P < 0.01; RO-IMPT_3mm/3% bladder_V65Gy = 5.2 ± 7.12% vs. RO-IMPT_5mm/3% bladder_V65Gy = 7.12 ± 9.59%, P < 0.01). The gamma analysis showed high dosimetric agreement between weekly and daily accumulated dose distributions.

CONCLUSIONS

This study demonstrated that for RO-IMPT optimization, ±3mm/±3% uncertainty is sufficient to create plans that meet desired CTV coverage while achieving superior sparing to OARs when compared with ± 5 mm/±3% uncertainty. Furthermore, weekly dose accumulation can accurately estimate the overall dose delivered to prostate cancer patients.

Anthropomorphic lung phantom based validation of in-room proton therapy 4D-CBCT image correction for dose calculation

PURPOSE

Ventilation-induced tumour motion remains a challenge for the accuracy of proton therapy treatments in lung patients. We investigated the feasibility of using a 4D virtual CT (4D-vCT) approach based on deformable image registration (DIR) and motion-aware 4D CBCT reconstruction (MA-ROOSTER) to enable accurate daily proton dose calculation using a gantry-mounted CBCT scanner tailored to proton therapy.

METHODS

Ventilation correlated data of 10 breathing phases were acquired from a porcine ex-vivo functional lung phantom using CT and CBCT. 4D-vCTs were generated by (1) DIR of the mid-position 4D-CT to the mid-position 4D-CBCT (reconstructed with the MA-ROOSTER) using a diffeomorphic Morphons algorithm and (2) subsequent propagation of the obtained mid-position vCT to the individual 4D-CBCT phases. Proton therapy treatment planning was performed to evaluate dose calculation accuracy of the 4D-vCTs. A robust treatment plan delivering a nominal dose of 60Gy was generated on the average intensity image of the 4D-CT for an approximated internal target volume (ITV). Dose distributions were then recalculated on individual phases of the 4D-CT and the 4D-vCT based on the optimized plan. Dose accumulation was performed for 4D-vCT and 4D-CT using DIR of each phase to the mid position, which was chosen as reference. Dose based on the 4D-vCT was then evaluated against the dose calculated on 4D-CT both, phase-by-phase as well as accumulated, by comparing dose volume histogram (DVH) values (D_mean, D_2%, D_98%, D_95%) for the ITV, and by a 3D-gamma index analysis (global, 3%/3mm, 5Gy, 20Gy and 30Gy dose thresholds).

RESULTS

Good agreement was found between the 4D-CT and 4D-vCT-based ITV-DVH curves. The relative differences ((CT-vCT)/CT) between accumulated values of ITV D_mean, D_2%, D_95% and D_98% for the 4D-CT and 4D-vCT-based dose distributions were -0.2%, 0.0%, -0.1% and -0.1%, respectively. Phase specific values varied between -0.5% and 0.2%, -0.2% and 0.5%, -3.5% and 1.5%, and -5.7% and 2.3%. The relative difference of accumulated D_mean over the lungs was 2.3% and D_mean for the phases varied between -5.4% and 5.8%. The gamma pass-rates with 5Gy, 20Gy and 30Gy thresholds for the accumulated doses were 96.7%, 99.6% and 99.9%, respectively. Phase-by-phase comparison yielded pass-rates between 86% and 97%, 88% and 98%, and 94% and 100%.

CONCLUSIONS

Feasibility of the suggested 4D-vCT workflow using proton therapy specific imaging equipment was shown. Results indicate the potential of the method to be applied for daily 4D proton dose estimation.

Comparison of CBCT based synthetic CT methods suitable for proton dose calculations in adaptive proton therapy

In-room imaging is a prerequisite for adaptive proton therapy. The use of onboard cone-beam computed tomography (CBCT) imaging, which is routinely acquired for patient position verification, can enable daily dose reconstructions and plan adaptation decisions. Image quality deficiencies though, hamper dose calculation accuracy and make corrections of CBCTs a necessity. This study compared three methods to correct CBCTs and create synthetic CTs that are suitable for proton dose calculations. CBCTs, planning CTs and repeated CTs (rCT) from 33 H&N cancer patients were used to compare a deep convolutional neural network (DCNN), deformable image registration (DIR) and an analytical image-based correction method (AIC) for synthetic CT (sCT) generation. Image quality of sCTs was evaluated by comparison with a same-day rCT, using mean absolute error (MAE), mean error (ME), Dice similarity coefficient (DSC), structural non-uniformity (SNU) and signal/contrast-to-noise ratios (SNR/CNR) as metrics. Dosimetric accuracy was investigated in an intracranial setting by performing gamma analysis and calculating range shifts. Neural network-based sCTs resulted in the lowest MAE and ME (37/2 HU) and the highest DSC (0.96). While DIR and AIC generated images with a MAE of 44/77 HU, a ME of -8/1 HU and a DSC of 0.94/0.90. Gamma and range shift analysis showed almost no dosimetric difference between DCNN and DIR based sCTs. The lower image quality of AIC based sCTs affected dosimetric accuracy and resulted in lower pass ratios and higher range shifts. Patient-specific differences highlighted the advantages and disadvantages of each method. For the set of patients, the DCNN created synthetic CTs with the highest image quality. Accurate proton dose calculations were achieved by both DCNN and DIR based sCTs. The AIC method resulted in lower image quality and dose calculation accuracy was reduced compared to the other methods.

Proton therapy for head and neck squamous cell carcinomas: A review of the physical and clinical challenges

The quality of radiation therapy has been shown to significantly influence the outcomes for head and neck squamous cell carcinoma (HNSCC) patients. The results of dosimetric studies suggest that intensity-modulated proton therapy (IMPT) could be of added value for HNSCC by being more effective than intensity-modulated (photon) radiation therapy (IMRT) for reducing side effects of radiation therapy. However, the physical properties of protons make IMPT more sensitive than photons to planning uncertainties. This could potentially have a negative effect on the quality of IMPT planning and delivery. For this review, the three French proton therapy centers collaborated to evaluate the differences between IMRT and IMPT. The review explored the effects of these uncertainties and their management for developing a robust and optimized IMPT treatment delivery plan to achieve clinical outcomes that are superior to those for IMRT. We also provide practical suggestions for the management of HNSCC carcinoma with IMPT. Because metallic dental implants can increase range uncertainties (3-10%), patient preparation for IMPT may require more systematic removal of in-field alien material than is done for IMRT. Multi-energy CT may be an alternative to calculate more accurately the dose distribution. The practical aspects that we describe are essential to guarantee optimal quality in radiation therapy in both model-based and randomized clinical trials.

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

Purpose

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

Methods

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

Results

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

Conclusion

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

Beam angle optimization using angular dependency of range variation assessed via water equivalent path length (WEPL) calculation for head and neck proton therapy

PURPOSE

To investigate angular sensitivity of proton range variation due to anatomic change in patients and patient setup error via water equivalent path length (WEPL) calculations.

METHODS

Proton range was estimated by calculating WEPL to the distal edge of target volume using planning CT (pCT) and weekly scatter-corrected cone-beam CT (CBCT) images of 11 head and neck patients. Range variation was estimated as the difference between the distal WEPLs calculated on pCT and scatter-corrected CBCT (cCBCT). This WEPL analysis was performed every five degrees ipsilaterally to the target. Statistics of the distal WEPL difference were calculated over the distal area to compare between different beam angles. Physician-defined contours were used for the WEPL calculation on both pCT and cCBCT, not considering local deformation of target volume. It was also tested if a couch kick (10°) can mitigate the range variation due to anatomic change and patient setup error.

RESULTS

For most of the patients considered, median, 75% quantile, and 95% quantile of the distal WEPL difference were largest for posterior oblique angles, indicating a higher chance of overdosing normal tissues at distal edge with these angles. Using a couch kick resulted in decrease in the WEPL difference for some posterior oblique angles.

CONCLUSIONS

It was demonstrated that the WEPL change has angular dependency for the cohort of head and neck cancer patients. Selecting beam configuration robust to anatomic change in patient and patient setup error may improve the treatment outcome of head and neck proton therapy.

CBCT correction using a cycle-consistent generative adversarial network and unpaired training to enable photon and proton dose calculation

In presence of inter-fractional anatomical changes, clinical benefits are anticipated from image-guided adaptive radiotherapy. Nowadays, cone-beam CT (CBCT) imaging is mostly utilized during pre-treatment imaging for position verification. Due to various artifacts, image quality is typically not sufficient for photon or proton dose calculation, thus demanding accurate CBCT correction, as potentially provided by deep learning techniques. This work aimed at investigating the feasibility of utilizing a cycle-consistent generative adversarial network (cycleGAN) for prostate CBCT correction using unpaired training. Thirty-three patients were included. The network was trained to translate uncorrected, original CBCT images (CBCT_org) into planning CT equivalent images (CBCT_cycleGAN). HU accuracy was determined by comparison to a previously validated CBCT correction technique (CBCT_cor). Dosimetric accuracy was inferred for volumetric-modulated arc photon therapy (VMAT) and opposing single-field uniform dose (OSFUD) proton plans, optimized on CBCT_cor and recalculated on CBCT_cycleGAN. Single-sided SFUD proton plans were utilized to assess proton range accuracy. The mean HU error of CBCT_cycleGAN with respect to CBCT_cor decreased from 24 HU for CBCT_org to  -6 HU. Dose calculation accuracy was high for VMAT, with average pass-rates of 100%/89% for a 2%/1% dose difference criterion. For proton OSFUD plans, the average pass-rate for a 2% dose difference criterion was 80%. Using a (2%, 2 mm) gamma criterion, the pass-rate was 96%. 93% of all analyzed SFUD profiles had a range agreement better than 3 mm. CBCT correction time was reduced from 6-10 min for CBCT_cor to 10 s for CBCT_cycleGAN. Our study demonstrated the feasibility of utilizing a cycleGAN for CBCT correction, achieving high dose calculation accuracy for VMAT. For proton therapy, further improvements may be required. Due to unpaired training, the approach does not rely on anatomically consistent training data or potentially inaccurate deformable image registration. The substantial speed-up for CBCT correction renders the method particularly interesting for adaptive radiotherapy.

[Proton therapy for head and neck squamous cell carcinomas: From physics to clinic]

Intensity-modulated radiation therapy (IMRT) is presently the recommended technique for the treatment of locally advanced head and neck carcinomas. Proton therapy would allow to reduce the volume of irradiated normal tissue and, thus, to decrease the risk of late dysphagia, xerostomia, dysgeusia and hypothyroidism. An exhaustive research was performed with the search engine PubMed by focusing on the papers about the physical difficulties that slow down use of proton therapy for head and neck carcinomas. Range uncertainties in proton therapy (±3 %) paradoxically limit the use of the steep dose gradient in distality. Calibration uncertainties can be important in the treatment of head and neck cancer in the presence of materials of uncertain stoichiometric composition (such as with metal implants, dental filling, etc.) and complex heterogeneities. Dental management for example may be different with IMRT or proton therapy. Some uncertainties can be somewhat minimized at the time of optimization. Inter- and intrafractional variations and uncertainties in Hounsfield units/stopping power can be integrated in a robust optimization process. Additional changes in patient's anatomy (tumour shrinkage, changes in skin folds in the beam patch, large weight loss or gain) require rescanning. Dosimetric and small clinical studies comparing photon and proton therapy have well shown the interest of proton therapy for head and neck cancers. Intensity-modulated proton therapy is a promising treatment as it can reduce the substantial toxicity burden of patients with head and neck squamous cell carcinoma compared to IMRT. Robust optimization will allow to perform an optimal treatment and to use proton therapy in current clinical practice.

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

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

Current state and future applications of radiological image guidance for particle therapy

In this review paper, we first give a short overview of radiological image guidance in photon radiotherapy, placing emphasis on the fact that linac based radiotherapy has outpaced particle therapy in the adoption of volumetric image guidance. While cone beam computed tomography (CBCT) has been an established technique in linac treatment rooms for almost two decades, the widespread adoption of volumetric image guidance in particle therapy, whether by means of CBCT or in-room CT imaging, is recent. This lag may be attributable to the bespoke nature and lower number of particle therapy installations, as well as the differences in geometry between those installations and linac treatment rooms. In addition, for particle therapy the so called shift invariance of the dose distribution rarely applies. An overview of the different volumetric image guidance solutions found at modern particle therapy facilities is provided, covering gantry, nozzle, C-arm, and couch-mounted CBCT as well different in-room CT configurations. A summary of the use of in-room volumetric imaging data beyond anatomy-based positioning is also presented as well as the necessary corrections to CBCT images for accurate water equivalent thickness calculation. Finally, the use of non-ionizing imaging modalities is discussed.

An automated planning strategy for near real-time adaptive proton therapy in prostate cancer

Proton therapy plans are very sensitive to anatomical changes such as density changes along the pencil-beam paths and changes in organ shape and location. Previously, we developed a restoration method which compensates for density changes along the pencil-beam paths but which is unable to adapt for anatomical changes. This study's purpose is to develop and evaluate an automated method for adaptation of IMPT plans in near real-time to the anatomy of the day. We developed an automated treatment plan adaptation method using (1) a restoration of spot positions (Bragg peaks) by adapting the energies to the new water equivalent path lengths; and (2) a spot addition to fully cover the target of the day, followed by a fast reference point method optimization of the spot weights resulting in a Pareto optimal plan for the daily anatomy. The method was developed and evaluated using 8-10 repeat CT scans of 11 prostate cancer patients, prescribing 55 Gy(RBE) (seminal vesicles and lymph nodes) with a boost to 74 Gy(RBE) (prostate). Applying the automated adaptation method resulted in a clinically acceptable target coverage (V _95% [Formula: see text] 98% and V _107% [Formula: see text] 2%) for 96% of the scans after a single iteration of adding 2500 spots. The other scans obtained target coverages with V _95% [Formula: see text] 98% and 2  <  V _107% [Formula: see text] 5%. When using two spot-addition iterations, all scans obtained clinically acceptable results. Compared to the restoration method the adaptation lowered the mean dose to rectum and bladder with median values of 6.2 Gy(RBE) and 4.7 Gy(RBE) respectively. The largest improvements were obtained for V _45Gy(RBE) for both rectum and bladder, with median differences of 10.3%-point and 10.8%-point respectively, and maximum differences up to 22%-point. The two adaptation steps took on average 7.3 s and 1.7 min respectively. No user interaction was needed, making this fast and fully automated method a first step towards online adaptive proton therapy.

Accuracy of radiotherapy dose calculations based on cone-beam CT: comparison of deformable registration and image correction based methods

Radiotherapy dose calculations based on cone-beam CT (CBCT) images can be inaccurate due to unreliable Hounsfield units (HU) in the CBCT. Deformable image registration of planning CT images to CBCT, and direct correction of CBCT image values are two methods proposed to allow heterogeneity corrected dose calculations based on CBCT. In this paper we compare the accuracy and robustness of these two approaches. CBCT images for 44 patients were used including pelvis, lung and head & neck sites. CBCT HU were corrected using a 'shading correction' algorithm and via deformable registration of planning CT to CBCT using either Elastix or Niftyreg. Radiotherapy dose distributions were re-calculated with heterogeneity correction based on the corrected CBCT and several relevant dose metrics for target and OAR volumes were calculated. Accuracy of CBCT based dose metrics was determined using an 'override ratio' method where the ratio of the dose metric to that calculated on a bulk-density assigned version of the same image is assumed to be constant for each patient, allowing comparison to the patient's planning CT as a gold standard. Similar performance is achieved by shading corrected CBCT and both deformable registration algorithms, with mean and standard deviation of dose metric error less than 1% for all sites studied. For lung images, use of deformed CT leads to slightly larger standard deviation of dose metric error than shading corrected CBCT with more dose metric errors greater than 2% observed (7% versus 1%).

Management of organ motion in scanned ion beam therapy

Scanned ion beam therapy has special demands for treatment of intra-fractionally moving tumors such as lesions in lung or liver. Interplay effects between beam and organ motion can in those settings lead to under-dosage of the target volume. Dedicated treatment techniques such as gating or abdominal compression are required. In addition 4D treatment planning should be used to determine strategies for patient specific treatment planning such as an increased beam focus or the use of internal target volumes incorporating range changes.Several work packages of the Clinical Research Units 214 and 214/2 funded by the German Research Council investigated the management of organ motion in scanned ion beam therapy. A focus was laid on 4D treatment planning using TRiP4D and the development of motion mitigation strategies including their quality assurance. This review focuses on the activity in the second funding period covering adaptive treatment planning strategies, 4D treatment plan optimization, and the application of motion management in pre-clinical research on radiation therapy of cardiac arrhythmias.

Water equivalent path length calculations using scatter-corrected head and neck CBCT images to evaluate patients for adaptive proton therapy

Proton therapy has dosimetric advantages due to the well-defined range of the proton beam over photon radiotherapy. When the proton beams, however, are delivered to the patient in fractionated radiation treatment, the treatment outcome is affected by delivery uncertainties such as anatomic change in the patient and daily patient setup error. This study aims at establishing a method to evaluate the dosimetric impact of the anatomic change and patient setup error during head and neck proton therapy. Range variations due to the delivery uncertainties were assessed by calculating water equivalent path length (WEPL) to the distal edge of tumor volume using planning CT and weekly treatment cone-beam CT (CBCT) images. Specifically, mean difference and root mean squared deviation (RMSD) of the distal WEPLs were calculated as the weekly range variations. To accurately calculate the distal WEPLs, an existing CBCT scatter correction algorithm was used. An automatic rigid registration was used to align the planning CT and treatment CBCT images, simulating a six degree-of-freedom couch correction at treatments. The authors conclude that the dosimetric impact of the anatomic change and patient setup error was reasonably captured in the differences of the distal WEPL variation with a range calculation uncertainty of 2%. The proposed method to calculate the distal WEPL using the scatter-corrected CBCT images can be an essential tool to decide the necessity of re-planning in adaptive proton therapy.

Feasibility of automated proton therapy plan adaptation for head and neck tumors using cone beam CT images

Background

Intensity modulated proton therapy (IMPT) of head and neck (H&N) tumors may benefit from plan adaptation to correct for the dose perturbations caused by weight loss and tumor volume changes observed in these patients. As cone beam CT (CBCT) is increasingly considered in proton therapy, it may be possible to use available CBCT images following intensity correction for plan adaptation. This is the first study exploring IMPT plan adaptation on CBCT images corrected and delineated by deformable image registration of the planning CT (pCT) to the CBCT, yielding a virtual CT (vCT).

Methods

A Morphons algorithm was used to deform the pCTs and corresponding delineations of 9 H&N cancer patients to a weekly CBCT acquired within ±3 days of a control replanning CT scan (rpCT). The IMPT treatment plans were adapted using the vCT and the adapted and original plans were recalculated on the rpCT for dose/volume parameter evaluation of the impact of adaptation.

Results

On the rpCT, the adapted plans were equivalent to the original plans in terms of target volumes D₉₅ and V₉₅, but showed a significant reduction of D₂ in these volumes. OAR doses were mostly equivalent or reduced. In particular, the adapted plans did not reduce parotid gland D_mean, but the dose to the optical system. For three patients the spinal cord or brain stem received higher, though well below tolerance, maximum dose. Subsequent tightening of the treatment planning constraints for these OARs on new vCT-adapted plans did not degrade target coverage and yielded pCT equivalent plans on the vCT.

Conclusions

An offline automated procedure to generate an adapted IMPT plan on CBCT images was developed and investigated. When evaluating the adapted plan on a control rpCT we observed reduced D₂ in target volumes as major improvement. OAR sparing was only partially improved by the procedure. Despite potential limitations in the accuracy of the vCT approach, an improved quality of the adapted plans could be achieved.

Electronic supplementary material

The online version of this article (doi:10.1186/s13014-016-0641-7) contains supplementary material, which is available to authorized users.