Validation of the proton range accuracy and optimization of CT calibration curves utilizing range probing

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.

A Monte-Carlo-based study of a single-2D-detector proton-radiography system

PURPOSE

To assess the feasibility of a proton radiography (pRG) system based on a single thin pixelated detector for water-equivalent path length (WEPL) and relative stopping power (RSP) measurements.

METHODS

A model of a pRG system consisting of a single pixelated detector measuring energy deposition and proton fluence was investigated in a Geant4-based Monte Carlo study. At the position directly after an object traversed by a broad proton beam, spatial 2D distributions are calculated of the energy deposition in, and the number of protons entering the detector. Their ratio relates to the 2D distribution of the average stopping power of protons in the detector. The system response is calibrated against the residual range in water of the protons to provide the 2D distribution of the WEPL of the object. The WEPL distribution is converted into the distribution of the RSP of the object. Simulations have been done, where the system has been tested on 13 samples of homogeneous materials of which the RSPs have been calculated and compared with RSPs determined from simulations of residual-range-in-water, which we refer to as reference RSPs.

RESULTS

For both human-tissue- and non-human-tissue-equivalent materials, the RSPs derived with the detector agree with the reference values within 1%.

CONCLUSION

The study shows that a pRG system based on one thin pixelated detection screen has the potential to provide RSP predictions with an accuracy of 1%.

Analysis of Performance and Failure Modes of the IROC Proton Liver Phantom

Purpose

To analyze trends in institutional performance and failure modes for the Imaging and Radiation Oncology Core's (IROC's) proton liver phantom.

Materials and Methods

Results of 66 phantom irradiations from 28 institutions between 2015 and 2020 were retrospectively analyzed. Univariate analysis and random forest models were used to associate irradiation conditions with phantom results. Phantom results included pass/fail classification, average thermoluminescent dosimeter (TLD) ratio of both targets, and percentage of pixels passing gamma of both targets. The following categories were evaluated in terms of how they predicted these outcomes: irradiation year, treatment planning system (TPS), TPS algorithm, treatment machine, number of irradiations, treatment technique, motion management technique, number of isocenters, and superior-inferior extent (in cm) of the 90% TPS isodose line for primary target 1 (PTV1) and primary target 2 (PTV2). In addition, failures were categorized by failure mode.

Results

Average pass rate was approximately 52% and average TLD ratio for both targets had slightly improved. As the treatment field increased to cover the target, the pass rate statistically significantly fell. Lower pass rates were observed for Mevion machines, scattered irradiation techniques, and gating and internal target volume (ITV) motion management techniques. Overall, the accuracy of the random forest modeling of the phantom results was approximately 73% ± 14%. The most important predictor was the superior-inferior extent for both targets and irradiation year. Three failure modes dominated the failures of the phantom: (1) systematic underdosing, (2) poor localization in the superior-inferior direction, and (3) range error. Only 44% of failures have similar failure modes between the 2 targets.

Conclusion

Improvement of the proton liver phantom has been observed; however, the pass rate remains the lowest among all IROC phantoms. Through various analysis techniques, range uncertainty, motion management, and underdosing are the main culprits of failures of the proton liver phantom. Clinically, careful consideration of the influences of liver proton therapy is needed to improve phantom performance and patient outcome.

An evaluation of the use of DirectSPR images for proton planning in the RayStation treatment planning software

An important source of uncertainty in proton therapy treatment planning is the assignment of stopping-power ratio (SPR) from CT data. A commercial product is now available that creates an SPR map directly from dual-energy CT (DECT). This paper investigates the use of this new product in proton treatment planning and compares the results to the current method of assigning SPR based on a single-energy CT (SECT). Two tissue surrogate phantoms were CT scanned using both techniques. The SPRs derived from single-energy CT and by DirectSPR™ were compared to measured values. SECT-based values agreed with measurements within 4% except for low density lung and high density bone, which differed by 13% and 8%, respectively. DirectSPR™ values were within 2% of measured values for all tissues studied. Both methods were also applied to scanned containers of three types of animal tissue, and the expected range of protons of two different energies was calculated in the treatment planning system and compared to the range measured using a multi-layer ion chamber. The average difference between range measurements and calculations based on SPR maps from dual- and single-energy CT, respectively, was 0.1 mm (0.07%) versus 2.2 mm (1.5%). Finally, a phantom was created using a layer of various tissue surrogate plugs on top of a 2D ion chamber array. Dose measurements on this array were compared to predictions using both single- and dual-energy CTs and SPR maps. While standard gamma pass rates for predictions based on DECT-derived SPR maps were slightly higher than those based on single-energy CT, the differences were generally modest for this measurement setup. This study showed that SPR maps created by the commercial product from dual-energy CT can successfully be used in RayStation to generate proton dose distributions and that these predictions agree well with measurements.

Considerations for Upright Particle Therapy Patient Positioning and Associated Image Guidance

Particle therapy is a rapidly growing field in cancer therapy. Worldwide, over 100 centers are in operation, and more are currently in construction phase. The interest in particle therapy is founded in the superior target dose conformity and healthy tissue sparing achievable through the particles’ inverse depth dose profile. This physical advantage is, however, opposed by increased complexity and cost of particle therapy facilities. Particle therapy, especially with heavier ions, requires large and costly equipment to accelerate the particles to the desired treatment energy and steer the beam to the patient. A significant portion of the cost for a treatment facility is attributed to the gantry, used to enable different beam angles around the patient for optimal healthy tissue sparing. Instead of a gantry, a rotating chair positioning system paired with a fixed horizontal beam line presents a suitable cost-efficient alternative. Chair systems have been used already at the advent of particle therapy, but were soon dismissed due to increased setup uncertainty associated with the upright position stemming from the lack of dedicated image guidance systems. Recently, treatment chairs gained renewed interest due to the improvement in beam delivery, commercial availability of vertical patient CT imaging and improved image guidance systems to mitigate the problem of anatomical motion in seated treatments. In this review, economical and clinical reasons for an upright patient positioning system are discussed. Existing designs targeted for particle therapy are reviewed, and conclusions are drawn on the design and construction of chair systems and associated image guidance. Finally, the different aspects from literature are channeled into recommendations for potential upright treatment layouts, both for retrofitting and new facilities.

Proton Therapy in the Management of Pancreatic Cancer

Simple Summary

Radiation treatment plays a pivotal a role in the management of pancreatic cancer, mainly in the neoadjuvant setting for borderline resectable tumors and in the definitive setting for unresectable localized disease. Most of the studies on pancreatic cancer use photon-based radiation therapy modalities. However, there is a growing interest in the application of protons therapy for gastrointestinal cancers. This review summarizes the literature on the use of proton therapy in the management of pancreatic cancer.

Abstract

Radiation therapy plays a central role in the treatment of pancreatic cancer. While generally shown to be feasible, proton irradiation, particularly when an ablative dose is planned, remains a challenge, especially due to tumor motion and the proximity to organs at risk, like the stomach, duodenum, and bowel. Clinically, standard doses of proton radiation treatment have not been shown to be statistically different from photon radiation treatment in terms of oncologic outcomes and toxicity rates as per non-randomized comparative studies. Fractionation schedules and concurrent chemotherapy combinations are yet to be optimized for proton therapy and are the subject of ongoing trials.

A comparison of proton stopping power measured with proton CT and x-ray CT in fresh postmortem porcine structures

PURPOSE

Currently, calculations of proton range in proton therapy patients are based on a conversion of CT Hounsfield units of patient tissues into proton relative stopping power. Uncertainties in this conversion necessitate larger proximal and distal planned target volume margins. Proton CT can potentially reduce these uncertainties by directly measuring proton stopping power. We aim to demonstrate proton CT imaging with complex porcine samples, to analyze in detail three-dimensional regions of interest, and to compare proton stopping powers directly measured by proton CT to those determined from x-ray CT scans.

METHODS

We have used a prototype proton imaging system with single proton tracking to acquire proton radiography and proton CT images of a sample of porcine pectoral girdle and ribs, and a pig's head. We also acquired close in time x-ray CT scans of the same samples and compared proton stopping power measurements from the two modalities. In the case of the pig's head, we obtained x-ray CT scans from two different scanners and compared results from high-dose and low-dose settings.

RESULTS

Comparing our reconstructed proton CT images with images derived from x-ray CT scans, we find agreement within 1% to 2% for soft tissues and discrepancies of up to 6% for compact bone. We also observed large discrepancies, up to 40%, for cavitated regions with mixed content of air, soft tissue, and bone, such as sinus cavities or tympanic bullae.

CONCLUSIONS

Our images and findings from a clinically realistic proton CT scanner demonstrate the potential for proton CT to be used for low-dose treatment planning with reduced margins.

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.

Can a ToF-PET photon attenuation reconstruction test stopping-power estimations in proton therapy? A phantom study

Objective. The aim of the phantom study was to validate and to improve the computed tomography (CT) images used for the dose computation in proton therapy. It was tested, if the joint reconstruction of activity and attenuation images of time-of-flight PET (ToF-PET) scans could improve the estimation of the proton stopping-power.Approach. The attenuation images, i.e. CT images with 511 keV gamma-rays (γCTs), were jointly reconstructed with activity maps from ToF-PET scans. Theβ⁺activity was produced with FDG and in a separate experiment with proton-induced radioactivation. The phantoms contained slabs of tissue substitutes. The use of theγCTs for the prediction of the beam stopping in proton therapy was based on a linear relationship between theγ-ray attenuation, the electron density, and the stopping-power of fast protons.Main results. The FDG based experiment showed sufficient linearity to detect a bias of bony tissue in the heuristic look-up table, which maps between x-ray CT images and proton stopping-power.γCTs can be used for dose computation, if the electron density of one type of tissue is provided as a scaling factor. A possible limitation is imposed by the spatial resolution, which is inferior by a factor of 2.5 compared to the one of the x-ray CT.γCTs can also be derived from off-line, ToF-PET scans subsequent to the application of a proton field with a hypofractionated dose level.Significance. γCTs are a viable tool to support the estimation of proton stopping with radiotracer-based ToF-PET data from diagnosis or staging. This could be of higher potential relevance in MRI-guided proton therapy.γCTs could form an alternative approach to make use of in-beam or off-line PET scans of proton-inducedβ⁺activity with possible clinical limitations due to the low number of coincidence counts.

Experimental assessment of inter-centre variation in stopping-power and range prediction in particle therapy

PURPOSE

Experimental assessment of inter-centre variation and absolute accuracy of stopping-power-ratio (SPR) prediction within 17 particle therapy centres of the European Particle Therapy Network.

MATERIAL AND METHODS

A head and body phantom with seventeen tissue-equivalent materials were scanned consecutively at the participating centres using their individual clinical CT scan protocol and translated into SPR with their in-house CT-number-to-SPR conversion. Inter-centre variation and absolute accuracy in SPR prediction were quantified for three tissue groups: lung, soft tissues and bones. The integral effect on range prediction for typical clinical beams traversing different tissues was determined for representative beam paths for the treatment of primary brain tumours as well as lung and prostate cancer.

RESULTS

An inter-centre variation in SPR prediction (2σ) of 8.7%, 6.3% and 1.5% relative to water was determined for bone, lung and soft-tissue surrogates in the head setup, respectively. Slightly smaller variations were observed in the body phantom (6.2%, 3.1%, 1.3%). This translated into inter-centre variation of integral range prediction (2σ) of 2.9%, 2.6% and 1.3% for typical beam paths of prostate-, lung- and primary brain-tumour treatments, respectively. The absolute error in range exceeded 2% in every fourth participating centre. The consideration of beam hardening and the execution of an independent HLUT validation had a positive effect, on average.

CONCLUSION

The large inter-centre variations in SPR and range prediction justify the currently clinically used margins accounting for range uncertainty, which are of the same magnitude as the inter-centre variation. This study underlines the necessity of higher standardisation in CT-number-to-SPR conversion.

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.

First-In-Human Validation of CT-Based Proton Range Prediction Using Prompt Gamma Imaging in Prostate Cancer Treatments

PURPOSE

Uncertainty in computed tomography (CT)-based range prediction substantially impairs the accuracy of proton therapy. Direct determination of the stopping-power ratio (SPR) from dual-energy CT (DECT) has been proposed (DirectSPR), and initial validation studies in phantoms and biological tissues have proven a high accuracy. However, a thorough validation of range prediction in patients has not yet been achieved by any means. Here, we present the first systematic validation of CT-based proton range prediction in patients using prompt gamma imaging (PGI).

METHODS AND MATERIALS

A PGI slit camera system with improved positioning accuracy, using a floor-based docking station, was used. Its overall uncertainty for range prediction validation was determined experimentally with both x-ray and beam measurements. The accuracy of range prediction in patients was determined from clinical PGI measurements during hypofractionated treatment of 5 patients with prostate cancer - in total 30 fractions with in-room control-CTs. For each pencil-beam-scanning spot, the range shift was obtained by comparing the PGI measurement to a control-CT-based PGI simulation. Three different SPR prediction approaches were applied in simulations: a standard CT-number-to-SPR conversion (Hounsfield look-up table [HLUT]), an adapted HLUT (DECT optimized), and DirectSPR. The spot-wise weighted mean range shift from all spots served as a measure for the accuracy of the respective range prediction approach.

RESULTS

A mean range prediction accuracy of 0.0% ± 0.5%, 0.3% ± 0.4%, and 1.8% ± 0.4% was obtained for DirectSPR, adapted HLUT, and standard HLUT, respectively. The overall validation uncertainty of the second-generation PGI slit camera is about 1 mm (2σ) for all approaches, which is smaller than the range prediction uncertainty for deep-seated tumors.

CONCLUSIONS

For the first time, range prediction accuracy was assessed in clinical routine using PGI range verification in prostate cancer treatments. Both DECT-derived range prediction approaches agree well with the measured proton range from PGI verification, whereas the standard HLUT approach differs relevantly. These results endorse the recent reduction of clinical safety margins in DirectSPR-based treatment planning in our institution.

Artificial intelligence supported single detector multi-energy proton radiography system

Proton radiography imaging was proposed as a promising technique to evaluate internal anatomical changes, to enable pre-treatment patient alignment, and most importantly, to optimize the patient specific CT number to stopping-power ratio conversion. The clinical implementation rate of proton radiography systems is still limited due to their complex bulky design, together with the persistent problem of (in)elastic nuclear interactions and multiple Coulomb scattering (i.e. range mixing). In this work, a compact multi-energy proton radiography system was proposed in combination with an artificial intelligence network architecture (ProtonDSE) to remove the persistent problem of proton scatter in proton radiography. A realistic Monte Carlo model of the Proteus^®One accelerator was built at 200 and 220 MeV to isolate the scattered proton signal in 236 proton radiographies of 80 digital anthropomorphic phantoms. ProtonDSE was trained to predict the proton scatter distribution at two beam energies in a 60%/25%/15% scheme for training, testing, and validation. A calibration procedure was proposed to derive the water equivalent thickness image based on the detector dose response relationship at both beam energies. ProtonDSE network performance was evaluated with quantitative metrics that showed an overall mean absolute percentage error below 1.4% ± 0.4% in our test dataset. For one example patient, detector dose to WET conversions were performed based on the total dose (ITotal), the primary proton dose (IPrimary), and the ProtonDSE corrected detector dose (ICorrected). The determined WET accuracy was compared with respect to the reference WET by idealistic raytracing in a manually delineated region-of-interest inside the brain. The error was determined 4.3% ± 4.1% forWET(ITotal),2.2% ± 1.4% forWET(IPrimary),and 2.5% ± 2.0% forWET(ICorrected).

Technical Note: CT calibration for proton treatment planning by cross-calibration with proton CT data

PURPOSE

This study explores the possibility of a new method for x-ray computed tomography (CT) calibration by means of cross-calibration with proton CT (pCT) data. The proposed method aims at a more accurate conversion of CT Hounsfield Units (HU) into proton stopping power ratio (SPR) relative to water to be used in proton-therapy treatment planning.

METHODS

X-ray CT scan was acquired on a synthetic anthropomorphic phantom, composed of different tissue equivalent materials (TEMs). A pCT apparatus was instead adopted to obtain a reference three-dimensional distribution of the phantom's SPR values. After rigid registration, the x-ray CT was artificially blurred to the same resolution of pCT. Then a scatter plot showing voxel-by-voxel SPR values as a function of HU was employed to link the two measurements and thus obtaining a cross-calibrated x-ray CT calibration curve. The cross-calibration was tested at treatment planning system and then compared with a conventional calibration based on exactly the same TEMs constituting the anthropomorphic phantom.

RESULTS

Cross-calibration provided an accurate SPR mapping, better than by conventional TEMs calibration. The dose distribution of single beams optimized on the reference SPR map was recomputed on cross-calibrated CT, showing, with respect to conventional calibration, minor deviation at the dose fall-off (lower than 1%).

CONCLUSIONS

The presented data demonstrated that, by means of reference pCT data, a heterogeneous phantom can be used for CT calibration, paving the way to the use of biological samples, with their accurate description of patients' tissues. This overcomes the limitations of conventional CT calibration requiring homogenous samples, only available by synthetic TEMs, which fail in accurately mimicking the properties of biological tissues. Once a heterogeneous biological sample is provided with its corresponding reference SPR maps, a cross-calibration procedure could be adopted by other PT centers, even when not equipped with a pCT system.

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

PURPOSE

The capability of proton therapy to provide highly conformal dose distributions is impaired by range uncertainties. The aim of this work is to apply range probing (RP), a form of a proton radiography-based quality control (QC) procedure for range accuracy assessment in head and neck cancer (HNC) patients in a clinical setting.

METHODS AND MATERIALS

This study included seven HNC patients. RP acquisition was performed using a multi-layer ionization chamber (MLIC). Per patient, two RP frames were acquired within the first two weeks of treatment, on days when a repeated CT scan was obtained. Per RP frame, integral depth dose (IDD) curves of 81 spots around the treatment isocenter were acquired. Range errors are determined as a discrepancy between calculated IDDs in the treatment planning system and measured residual ranges by the MLIC. Range errors are presented relative to the water equivalent path length of individual proton spots. In addition to reporting results for complete measurement frames, an analysis, excluding range error contributions due to anatomical changes, is presented.

RESULTS

Discrepancies between measured and calculated ranges are smaller when performing RP calculations on the day-specific patient anatomy rather than the planning CT. The patient-specific range evaluation shows an agreement between calculated and measured ranges for spots in anatomically consistent areas within 3% (1.5 standard deviation).

CONCLUSIONS

The results of an RP-based QC procedure implemented in the clinical practice for HNC patients have been demonstrated. The agreement of measured and simulated proton ranges confirms the 3% uncertainty margin for robust optimization. Anatomical variations show a predominant effect on range accuracy, motivating efforts towards the implementation of adaptive radiotherapy.

openPR - A computational tool for CT conversion assessment with proton radiography

PURPOSE

One of the main sources of uncertainty in proton therapy is the conversion of the Hounsfield Units of the planning CT to (relative) proton stopping powers. Proton radiography provides range error maps but these can be affected by other sources of errors as well as the CT conversion (e.g., residual misalignment). To better understand and quantify range uncertainty, it is desirable to measure the individual contributions and particularly those associated to the CT conversion.

METHODS

A workflow is proposed to carry out an assessment of the CT conversion solely on the basis of proton radiographs of real tissues measured with a multilayer ionization chamber (MLIC). The workflow consists of a series of four stages: (a) CT and proton radiography acquisitions, (b) CT and proton radiography registration in postprocessing, (c) sample-specific validation of the semi-empirical model both used in the registration and to estimate the water equivalent path length (WEPL), and (d) WEPL error estimation. The workflow was applied to a pig head as part of the validation of the CT calibration of the proton therapy center PARTICLE at UZ Leuven, Belgium.

RESULTS

The CT conversion-related uncertainty computed based on the well-established safety margin rule of 1.2 mm + 2.4% were overestimated by 71% on the pig head. However, the range uncertainty was very much underestimated where cavities were encountered by the protons. Excluding areas with cavities, the overestimation of the uncertainty was 500%. A correlation was found between these localized errors and HUs between -1000 and -950, suggesting that the underestimation was not a consequence of an inaccurate conversion but was probably rather due to the resolution of the CT leading to material mixing at interfaces. To reduce these errors, the CT calibration curve was adapted by increasing the HU interval corresponding to the air up to -950.

CONCLUSION

The application of the workflow as part of the validation of the CT conversion to RSPs showed an overall overestimation of the expected uncertainty. Moreover, the largest WEPL errors were found to be related to the presence of cavities which nevertheless are associated with low WEPL values. This suggests that the use of this workflow on patients or in a generalized study on different types of animal tissues could shed sufficient light on how the contributions to the CT conversion-related uncertainty add up to potentially reduce up to several millimeters the uncertainty estimations taken into account in treatment planning. All the algorithms required to perform the workflow were implemented in the computational tool named openPR which is part of openREGGUI, an open-source image processing platform for adaptive proton therapy.

Head and neck IMPT probabilistic dose accumulation: Feasibility of a 2 mm setup uncertainty setting

OBJECTIVE

To establish optimal robust optimization uncertainty settings for clinical head and neck cancer (HNC) patients undergoing 3D image-guided pencil beam scanning (PBS) proton therapy.

METHODS

We analyzed ten consecutive HNC patients treated with 70 and 54.25 Gy_RBE to the primary and prophylactic clinical target volumes (CTV) respectively using intensity-modulated proton therapy (IMPT). Clinical plans were generated using robust optimization with 5 mm/3% setup/range uncertainties (RayStation v6.1). Additional plans were created for 4, 3, 2 and 1 mm setup and 3% range uncertainty and for 3 mm setup and 3%, 2% and 1% range uncertainty. Systematic and random error distributions were determined for setup and range uncertainties based on our quality assurance program. From these, 25 treatment scenarios were sampled for each plan, each consisting of a systematic setup and range error and daily random setup errors. Fraction doses were calculated on the weekly verification CT closest to the date of treatment as this was considered representative of the daily patient anatomy.

RESULTS

Plans with a 2 mm/3% setup/range uncertainty setting adequately covered the primary and prophylactic CTV (V₉₅ ≥ 99% in 98.8% and 90.8% of the treatment scenarios respectively). The average organ-at-risk dose decreased with 1.1 Gy_RBE/mm setup uncertainty reduction and 0.5 Gy_RBE/1% range uncertainty reduction. Normal tissue complication probabilities decreased by 2.0%/mm setup uncertainty reduction and by 0.9%/1% range uncertainty reduction.

CONCLUSION

The results of this study indicate that margin reduction below 3 mm/3% is possible but requires a larger cohort to substantiate clinical introduction.