Technical Note: First report on an in vivo range probing quality control procedure for scanned proton beam therapy in head and neck cancer patients

Imaging lung tumor motion using integrated-mode proton radiography-A phantom study towards tumor tracking in proton radiotherapy

BACKGROUND

Motion of lung tumors during radiotherapy leads to decreased accuracy of the delivered dose distribution. This is especially true for proton radiotherapy due to the finite range of the proton beam. Methods for mitigating motion rely on knowing the position of the tumor during treatment.

PURPOSE

Proton radiography uses the treatment beam, at an energy high enough to traverse the patient, to produce a radiograph. This work shows the first results of using an integrated-mode proton radiography system to track the position of moving objects in an experimental phantom study; demonstrating the potential of using this method for measuring tumor motion.

METHODS

Proton radiographs of an anthropomorphic lung phantom, with a motor-driven tumor insert, were acquired approximately every 1 s, using tumor inserts of 10, 20, and 30 mm undergoing a known periodic motion. The proton radiography system used a monolithic scintillator block and digital cameras to capture the residual range of each pencil beam passing through the phantom. These ranges were then used to produce a water equivalent thickness map of the phantom. The centroid of the tumor insert in the radiographs was used to determine its position. This measured position was then compared to the known motion of the phantom to determine the accuracy.

RESULTS

Submillimeter accuracy on the measurement of the tumor insert was achieved when using a 30 mm tumor insert with a period of 24 s and was found to be improved for decreasing motion amplitudes with a mean absolute error (MAE) of 1.0, 0.9, and 0.7 mm for 20, 15, and 10 mm respectively. Using smaller tumor inserts reduced the accuracy with a MAE of 1.8 and 1.9 mm for a 20 and 10 mm insert respectively undergoing a periodic motion with an amplitude of 20 mm and a period of 24 s. Using a shorter period resulted in significant motion artifacts reducing the accuracy to a MAE of 2.2 mm for a 12 s period and 3.1 mm for a 6 s period for the 30 mm insert with an amplitude of 20 mm.

CONCLUSIONS

This work demonstrates that the position of a lung tumor insert in a realistic anthropomorphic phantom can be measured with high accuracy using proton radiographs. Results show that the accuracy of the position measurement is the highest for slower tumor motions due to a reduction in motion artifacts. This indicates that the primary obstacle to accurate measurement is the speed of the radiograph acquisition. Although the slower tumor motions used in this study are not clinically realistic, this work demonstrates the potential for using proton radiography for measuring tumor motion with an increased scanning speed that results in a decreased acquisition time.

A probabilistic evaluation of the Dutch robustness and model-based selection protocols for Head-and-Neck IMPT: A multi-institutional study

BACKGROUND AND PURPOSE

In the Netherlands, 2 protocols have been standardized for PT among the 3 proton centers: a robustness evaluation (RE) to ensure adequate CTV dose and a model-based selection (MBS) approach for IMPT patient-selection. This multi-institutional study investigates (i) inter-patient and inter-center variation of target dose from the RE protocol and (ii) the robustness of the MBS protocol against treatment errors for a cohort of head-and-neck cancer (HNC) patients treated in the 3 Dutch proton centers.

MATERIALS AND METHODS

Clinical treatment plans of 100 HNC patients were evaluated. Polynomial Chaos Expansion (PCE) was used to perform a comprehensive robustness evaluation per plan, enabling the probabilistic evaluation of 100,000 complete fractionated treatments. PCE allowed to derive scenario distributions of clinically relevant dosimetric parameters to assess CTV dose (D99.8%/D0.2%, based on a prior photon plan calibration) and tumour control probabilities (TCP) as well as the evaluation of the dose to OARs and normal tissue complication probabilities (NTCP) per center.

RESULTS

For the CTV70.00, doses from the RE protocol were consistent with the clinical plan evaluation metrics used in the 3 centers. For the CTV54.25, D99.8% were consistent with the clinical plan evaluation metrics at center 1 and 2 while, for center 3, a reduction of 1 GyRBE was found on average. This difference did not impact modelled TCP at center 3. Differences between expected and nominal NTCP were below 0.3 percentage point for most patients.

CONCLUSION

The standardization of the RE and MBS protocol lead to comparable results in terms of TCP and the NTCPs. Still, significant inter-patient and inter-center variation in dosimetric parameters remained due to clinical practice differences at each institution. The MBS approach is a robust protocol to qualify patients for PT.

Opportunities and challenges of upright patient positioning in radiotherapy

Objective.Upright positioning has seen a surge in interest as a means to reduce radiotherapy (RT) cost, improve patient comfort, and, in selected cases, benefit treatment quality. In particle therapy (PT) in particular, eliminating the need for a gantry can present massive cost and facility footprint reduction. This review discusses the opportunities of upright RT in perspective of the open challenges.Approach.The clinical, technical, and workflow challenges that come with the upright posture have been extracted from an extensive literature review, and the current state of the art was collected in a synergistic perspective from photon and particle therapy. Considerations on future developments and opportunities are provided.Main results.Modern image guidance is paramount to upright RT, but it is not clear which modalities are essential to acquire in upright posture. Using upright MRI or upright CT, anatomical differences between upright/recumbent postures have been observed for nearly all body sites. Patient alignment similar to recumbent positioning was achieved in small patient/volunteer cohorts with prototype upright positioning systems. Possible clinical advantages, such as reduced breathing motion in upright position, have been reported, but limited cohort sizes prevent resilient conclusions on the treatment impact. Redesign of RT equipment for upright positioning, such as immobilization accessories for various body regions, is necessary, where several innovations were recently presented. Few clinical studies in upright PT have already reported promising outcomes for head&neck patients.Significance.With more evidence for benefits of upright RT emerging, several centers worldwide, particularly in PT, are installing upright positioning devices or have commenced upright treatment. Still, many challenges and open questions remain to be addressed to embed upright positioning firmly in the modern RT landscape. Guidelines, professionals trained in upright patient positioning, and large-scale clinical studies are required to bring upright RT to fruition.

Integrated-mode proton radiography with 2D lateral projections

Integrated-mode proton radiography leading to water equivalent thickness (WET) maps is an avenue of interest for motion management, patient positioning, andin vivorange verification. Radiographs can be obtained using a pencil beam scanning setup with a large 3D monolithic scintillator coupled with optical cameras. Established reconstruction methods either (1) involve a camera at the distal end of the scintillator, or (2) use a lateral view camera as a range telescope. Both approaches lead to limited image quality. The purpose of this work is to propose a third, novel reconstruction framework that exploits the 2D information provided by two lateral view cameras, to improve image quality achievable using lateral views. The three methods are first compared in a simulated Geant4 Monte Carlo framework using an extended cardiac torso (XCAT) phantom and a slanted edge. The proposed method with 2D lateral views is also compared with the range telescope approach using experimental data acquired with a plastic volumetric scintillator. Scanned phantoms include a Las Vegas (contrast), 9 tissue-substitute inserts (WET accuracy), and a paediatric head phantom. Resolution increases from 0.24 (distal) to 0.33 lp mm-1(proposed method) on the simulated slanted edge phantom, and the mean absolute error on WET maps of the XCAT phantom is reduced from 3.4 to 2.7 mm with the same methods. Experimental data from the proposed 2D lateral views indicate a 36% increase in contrast relative to the range telescope method. High WET accuracy is obtained, with a mean absolute error of 0.4 mm over 9 inserts. Results are presented for various pencil beam spacing ranging from 2 to 6 mm. This work illustrates that high quality proton radiographs can be obtained with clinical beam settings and the proposed reconstruction framework with 2D lateral views, with potential applications in adaptive proton therapy.

Technical note: development of a simulation framework, enabling the investigation of locally tuned single energy proton radiography

Range uncertainties remain a limitation for the confined dose distribution that proton therapy can offer. The uncertainty stems from the ambiguity when translating CT Hounsfield Units (HU) into proton stopping powers. Proton Radiography (PR) can be used to verify the proton range. Specifically, PR can be used as a quality-control tool for CBCT-based synthetic CTs. An essential part of the work illustrating the potential of PR has been conducted using multi-layer ionization chamber (MLIC) detectors and mono-energetic PR. Due to the dimensions of commercially available MLICs, clinical adoption is cumbersome. Here, we present a simulation framework exploring locally-tuned single energy (LTSE) proton radiography and corresponding potential compact PR detector designs. Based on a planning CT data set, the presented framework models the water equivalent thickness. Subsequently, it analyses the proton energies required to pass through the geometry within a defined ROI. In the final step, an LTSE PR is simulated using the MCsquare Monte Carlo code. In an anatomical head phantom, we illustrate that LTSE PR allows for a significantly shorter longitudinal dimension of MLICs. We compared PR simulations for two exemplary 30 × 30 mm2proton fields passing the phantom at a 90° angle at an anterior and a posterior location in an iso-centric setup. The longitudinal distance over which all spots per field range out is significantly reduced for LTSE PR compared to mono-energetic PR. In addition, we illustrate the difference in shape of integral depth dose (IDD) when using constrained PR energies. Finally, we demonstrate the accordance of simulated and experimentally acquired IDDs for an LTSE PR acquisition. As the next steps, the framework will be used to investigate the sensitivity of LTSE PR to various sources of errors. Furthermore, we will use the framework to systematically explore the dimensions of an optimized MLIC design for daily clinical use.

PTV-based VMAT vs. robust IMPT for head-and-neck cancer: A probabilistic uncertainty analysis of clinical plan evaluation with the Dutch model-based selection

BACKGROUND AND PURPOSE

In the Netherlands, head-and-neck cancer (HNC) patients are referred for proton therapy (PT) through model-based selection (MBS). However, treatment errors may compromise adequate CTV dose. Our aims are: (i) to derive probabilistic plan evaluation metrics on the CTV consistent with clinical metrics; (ii) to evaluate plan consistency between photon (VMAT) and proton (IMPT) planning in terms of CTV dose iso-effectiveness and (iii) to assess the robustness of the OAR doses and of the risk toxicities involved in the MBS.

MATERIALS AND METHODS

Sixty HNC plans (30 IMPT/30 VMAT) were included. A robustness evaluation with 100,000 treatment scenarios per plan was performed using Polynomial Chaos Expansion (PCE). PCE was applied to determine scenario distributions of clinically relevant dosimetric parameters, which were compared between the 2 modalities. Finally, PCE-based probabilistic dose parameters were derived and compared to clinical PTV-based photon and voxel-wise proton evaluation metrics.

RESULTS

Probabilistic dose to near-minimum volume v = 99.8% for the CTV correlated best with clinical PTV-D98% and VWmin-D98%,CTV doses for VMAT and IMPT respectively. IMPT showed slightly higher nominal CTV doses, with an average increase of 0.8 GyRBE in the median of the D99.8%,CTV distribution. Most patients qualified for IMPT through the dysphagia grade II model, for which an average NTCP gain of 10.5 percentages points (%-point) was found. For all complications, uncertainties resulted in moderate NTCP spreads lower than 3 p.p. on average for both modalities.

CONCLUSION

Despite the differences between photon and proton planning, the comparison between PTV-based VMAT and robust IMPT is consistent. Treatment errors had a moderate impact on NTCPs, showing that the nominal plans are a good estimator to qualify patients for PT.

Technical note: Flat panel proton radiography with a patient specific imaging field for accurate WEPL assessment

BACKGROUND

Proton radiography (PR) uses highly energetic proton beams to create images where energy loss is the main contrast mechanism. Water-equivalent path length (WEPL) measurements using flat panel PR (FP-PR) have potential for in vivo range verification. However, an accurate WEPL measurement via FP-PR requires irradiation with multiple energy layers, imposing high imaging doses.

PURPOSE

A FP-PR method is proposed for accurate WEPL determination based on a patient-specific imaging field with a reduced number of energies (n) to minimize imaging dose.

METHODS

Patient-specific FP-PRs were simulated and measured for a head and neck (HN) phantom. An energy selection algorithm estimated spot-wise the lowest energy required to cross the anatomy (Emin) using a water-equivalent thickness map. Starting from Emin, n was restricted to certain values (n = 26, 24, 22, …, 2 for simulations, n = 10 for measurements), resulting in patient-specific FP-PRs. A reference FP-PR with a complete set of energies was compared against patient-specific FP-PRs covering the whole anatomy via mean absolute WEPL differences (MAD), to evaluate the impact of the developed algorithm. WEPL accuracy of patient-specific FP-PRs was assessed using mean relative WEPL errors (MRE) with respect to measured multi-layer ionization chamber PRs (MLIC-PR) in the base of skull, brain, and neck regions.

RESULTS

MADs ranged from 2.1 mm (n = 26) to 21.0 mm (n = 2) for simulated FP-PRs, and 7.2 mm for measured FP-PRs (n = 10). WEPL differences below 1 mm were observed across the whole anatomy, except at the phantom surfaces. Measured patient-specific FP-PRs showed good agreement against MLIC-PRs, with MREs of 1.3 ± 2.0%, -0.1 ± 1.0%, and -0.1 ± 0.4% in the three regions of the phantom.

CONCLUSION

A method to obtain accurate WEPL measurements using FP-PR with a reduced number of energies selected for the individual patient anatomy was established in silico and validated experimentally. Patient-specific FP-PRs could provide means of in vivo range verification.

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.

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.

Optimizing calibration settings for accurate water equivalent path length assessment using flat panel proton radiography

OBJECTIVE

Proton range uncertainties can compromise the effectiveness of proton therapy treatments. Water equivalent path length (WEPL) assessment by flat panel detector proton radiography (FP-PR) can provide means of range uncertainty detection. Since WEPL accuracy intrinsically relies on the FP-PR calibration parameters, the purpose of this study is to establish an optimal calibration procedure that ensures high accuracy of WEPL measurements. To that end, several calibration settings were investigated.

APPROACH

FP-PR calibration datasets were obtained simulating PR fields with different proton energies, directed towards water-equivalent material slabs of increasing thickness. The parameters investigated were the spacing between energy layers (ΔE) and the increment in thickness of the water-equivalent material slabs (ΔX) used for calibration. 30 calibrations were simulated, as a result of combining ΔE = 9, 7, 5, 3, 1 MeV and ΔX = 10, 8, 5, 3, 2, 1 mm. FP-PRs through a CIRS electron density phantom were simulated, and WEPL images corresponding to each calibration were obtained. Ground truth WEPL values were provided by range probing multi-layer ionization chamber simulations on each insert of the phantom. Relative WEPL errors between FP-PR simulations and ground truth were calculated for each insert. Mean relative WEPL errors and standard deviations across all inserts were computed for WEPL images obtained with each calibration.

MAIN RESULTS

Large mean and standard deviations were found in WEPL images obtained with large ΔEvalues (ΔE = 9 or 7 MeV), for any ΔX. WEPL images obtained with ΔE ≤ 5 MeV and ΔX ≤ 5 mm resulted in a WEPL accuracy with mean values within ±0.5% and standard deviations around 1%.

SIGNIFICANCE

An optimal FP calibration in the framework of this study was established, characterized by 3 MeV ≤ ΔE ≤ 5 MeV and 2 mm ≤ ΔX ≤ 5 mm. Within these boundaries, highly accurate WEPL acquisitions using FP-PR are feasible and practical, holding the potential to assist future online range verification quality control procedures.

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.