Relative Stopping Power Research Articles

Development of a Tumor Growth-responsive Robust Planning Method for Head and Neck Cancer in Carbon-ion Radiotherapy.

The current study aimed to evaluate a treatment planning method that is robust against tumor growth and to assess its effectiveness in particle therapy for head and neck cancer. The proposed method optimizes dose distribution by replacing the relative stopping power ratio (rSPR) of the clinical target volume (CTV) cavity region with a tumor-equivalent rSPR (Condition 1). The optimized initial treatment plan template was then recalculated using in-room CT images acquired in the same treatment position, and the doses to the tumor and organs at risk were compared with those in the initial treatment plan. We evaluated this method in 10 patients with head and neck cancer treated with carbon ion radiotherapy. To evaluate the effectiveness of the proposed method, we compared it to the initial treatment plan without the replacement (Condition 2). CTV V95% reduction relative to that of the initial treatment plan at the end of treatment was 1.3%±2.9% and 2.6%±3.7% for Condition 1 and Condition 2, respectively, with Condition 1 (Replacement condition) providing better CTV coverage. Subgroup analysis showed a higher change in target coverage in the mucosal melanoma than in the adenoid cystic carcinoma, suggesting that this was influenced by the rate of tumor growth. The proposed treatment planning method, which adjusts for tumor growth by modifying the rSPR in the cavity region, improved robustness in carbon ion radiotherapy for head and neck cancer. Condition 1 (with replacement) achieved better CTV coverage than Condition 2 (without replacement), particularly in fast-growing tumors like mucosal melanoma. This method ensures consistent dose delivery to tumors while maintaining safe doses to organs at risk, offering potential for improved treatment outcomes.

First test beam of the DMAPS-based proton tracker at the pMBRT facility at the Curie Institute

Objective.Proton radiotherapy's efficacy relies on an accurate relative stopping power (RSP) map of the patient to optimise the treatment plan and minimize uncertainties. Currently, a conversion of a Hounsfield Units map obtained by a common x-ray computed tomography (CT) is used to compute the RSP. This conversion is one of the main limiting factors for proton radiotherapy. To bypass this conversion a direct RSP map could be obtained by performing a proton CT (pCT). The focal point of this article is to present a proof of concept of the potential of fast pixel technologies for proton tracking at clinical facilities.Approach.A two-layer tracker based on the TJ-Monopix1, a depleted monolithic active pixel sensor (DMAPS) chip initially designed for the ATLAS, was tested at the proton minibeam radiotherapy beamline at the Curie Institute. The chips were subjected to 100 MeV protons passing through the single slit collimator (0.4×20mm2) with fluxes up to1.3×107p s-1 cm-2. The performance of the proton tracker was evaluated with GEANT4 simulations.Main results.The beam profile and dispersion in air were characterized by an opening of 0.031 mm cm-1, and aσx=0.172mm at the position of the slit. The results of the proton tracking show how the TJ-Monopix1 chip can effectively track protons in a clinical environment, achieving a tracking purity close to 70%, and a position resolution below 0.5 mm; confirming the chip's ability to handle high proton fluxes with a competitive performance.Significance.This work suggests that DMAPS technologies can be a cost-effective alternative for proton imaging. Additionally, the study identifies areas where further optimization of chip design is required to fully leverage these technologies for clinical ion imaging applications.

Proton dose calculation with LSTM networks in presence of a magnetic field

Objective.To present a long short-term memory (LSTM) network-based dose calculation method for magnetic resonance (MR)-guided proton therapy.Approach.35 planning computed tomography (CT) images of prostate cancer patients were collected for Monte Carlo (MC) dose calculation under a perpendicular 1.5 T magnetic field. Proton pencil beams (PB) at three energies (150, 175, and 200 MeV) were simulated (7560 PBs at each energy). A 3D relative stopping power cuboid covering the extent of the PB dose was extracted and given as input to the LSTM model, yielding a 3D predicted PB dose. Three single-energy (SE) LSTM models were trained separately on the corresponding 150/175/200 MeV datasets and a multi-energy (ME) LSTM model with an energy embedding layer was trained on either the combined dataset with three energies or a continuous energy (CE) dataset with 1 MeV steps ranging from 125 to 200 MeV. For each model, training and validation involved 25 patients and 10 patients were for testing. Two single field uniform dose prostate treatment plans were optimized and recalculated with MC and the CE model.Results.Test results of all PBs from the three SE models showed a mean gamma passing rate (2%/2 mm, 10% dose cutoff) above 99.9% with an average center-of-mass (COM) discrepancy below 0.4 mm between predicted and simulated trajectories. The ME model showed a mean gamma passing rate exceeding 99.8% and a COM discrepancy of less than 0.5 mm at the three energies. Treatment plan recalculation by the CE model yielded gamma passing rates of 99.6% and 97.9%. The inference time of the models was 9-10 ms per PB.Significance.LSTM models for proton dose calculation in a magnetic field were developed and showed promising accuracy and efficiency for prostate cancer patients.

Development of a proton CT imaging system using scintillator-based range detection.

The accuracy of proton therapy and preclinical proton irradiation experiments is susceptible to proton range uncertainties, which partly stem from the inaccurate conversion between CT numbers and relative stopping power (RSP). Proton computed tomography (PCT) can reduce these uncertainties by directly acquiring RSP maps. This study aims to develop a novel PCT imaging system based on scintillator-based proton range detection for accurate RSP reconstruction. The proposed PCT system consists of a pencil-beam brass collimator with a 1mm aperture, an object stage capable of translation and 360° rotation, a plastic scintillator for dose-to-light conversion, and a complementary metal oxide semiconductor (CMOS) camera for light distribution acquisition. A calibration procedure based on Monte Carlo (MC) simulation was implemented to convert the obtained light ranges into water equivalent ranges. The water equivalent path lengths (WEPLs) of the imaged object were determined by calculating the differences in proton ranges obtained with and without the object in the beam path. To validate the WEPL calculation, measurements of WEPLs for eight tissue-equivalent inserts were conducted. PCT imaging was performed on a custom-designed phantom and a mouse, utilizing both 60 and 360 projections. The filtered back projection (FBP) algorithm was employed to reconstruct the RSP from WEPLs. Image quality was assessed based on the reconstructed RSP maps and compared to reference and simulation-based reconstructions. The differences between the calibrated and reference ranges of 110-150 MeV proton beams were within 0.18mm. The WEPLs of eight tissue-equivalent inserts were measured with accuracies better than 1%. Phantom experiments exhibited good agreement with reference and simulation-based reconstructions, demonstrating average RSP errors of 1.26%, 1.38%, and 0.38% for images reconstructed with 60 projections, 60 projections after penalized weighted least-squares algorithm denoising, and 360 projections, respectively. Mouse experiments provided clear observations of mouse contours and major tissue types. MC simulation estimated an imaging dose of 3.44 cGy for decent RSP reconstruction. The proposed PCT imaging system enables RSP map acquisition with high accuracy and has the potential to improve dose calculation accuracy in proton therapy and preclinical proton irradiation experiments.

Characterisation of 3D-printable thermoplastics to be used as tissue-equivalent materials in photon and proton beam radiotherapy end-to-end quality assurance devices

Objective. To investigate the potential of 3D-printable thermoplastics as tissue-equivalent materials to be used in multimodal radiotherapy end-to-end quality assurance (QA) devices. Approach. Six thermoplastics were investigated: Polylactic Acid (PLA), Acrylonitrile Butadiene Styrene (ABS), Polyethylene Terephthalate Glycol (PETG), Polymethyl Methacrylate (PMMA), High Impact Polystyrene (HIPS) and StoneFil. Measurements of mass density ( ρ ), Relative Electron Density (RED), in a nominal 6 MV photon beam, and Relative Stopping Power (RSP), in a 210 MeV proton pencil-beam, were performed. Average Hounsfield Units (HU) were derived from CTs acquired with two independent scanners. The calibration curves of both scanners were used to predict average ρ, RED and RSP values and compared against the experimental data. Finally, measured data of ρ, RED and RSP was compared against theoretical values estimated for the thermoplastic materials and biological tissues. Main results. Overall, good ρ and RSP CT predictions were made; only PMMA and PETG showed differences >5%. The differences between experimental and CT predicted RED values were also <5% for PLA, ABS, PETG and PMMA; for HIPS and StoneFil higher differences were found (6.94% and 9.42/15.34%, respectively). Small HU variations were obtained in the CTs for all materials indicating good uniform density distribution in the samples production. ABS, PLA, PETG and PMMA showed potential equivalency for a variety of soft tissues (adipose tissue, skeletal muscle, brain and lung tissues, differences within 0.19%–8.35% for all properties). StoneFil was the closest substitute to bone, but differences were >10%. Theoretical calculations of all properties agreed with experimental values within 5% difference for most thermoplastics. Significance. Several 3D-printed thermoplastics were promising tissue-equivalent materials to be used in devices for end-to-end multimodal radiotherapy QA and may not require corrections in treatment planning systems’ dose calculations. Theoretical calculations showed promise in identifying thermoplastics matching target biological tissues before experiments are performed.

A direct comparison of multi-energy x-ray and proton CT for imaging and relative stopping power estimation of plastic and ex-vivo phantoms

Objective.Proton therapy administers a highly conformal dose to the tumour region, necessitating accurate prediction of the patient's 3D map of proton relative stopping power (RSP) compared to water. This remains challenging due to inaccuracies inherent in single-energy computed tomography (SECT) calibration. Recent advancements in spectral x-ray CT (xCT) and proton CT (pCT) have shown improved RSP estimation compared to traditional SECT methods. This study aims to provide the first comparison of the imaging and RSP estimation performance among dual-energy CT (DECT) and photon-counting CT (PCCT) scanners, and a pCT system prototype.Approach.Two phantoms were scanned with the three systems for their performance characterisation: a plastic phantom, filled with water and containing four plastic inserts and a wood insert, and a heterogeneous biological phantom, containing a formalin-stabilised bovine specimen. RSP maps were generated by converting CT numbers to RSP using a calibration based on low- and high-energy xCT images, while pCT utilised a distance-driven filtered back projection algorithm for RSP reconstruction. Spatial resolution, noise, and RSP accuracy were compared across the resulting images.Main results.All three systems exhibited similar spatial resolution of around 0.54 lp/mm for the plastic phantom. The PCCT images were less noisy than the DECT images at the same dose level. The lowest mean absolute percentage error (MAPE) of RSP,(0.28±0.07)%, was obtained with the pCT system, compared to MAPE values of(0.51±0.08)%and(0.80±0.08)%for the DECT- and PCCT-based methods, respectively. For the biological phantom, the xCT-based methods resulted in higher RSP values in most of the voxels compared to pCT.Significance.The pCT system yielded the most accurate estimation of RSP values for the plastic materials, and was thus used to benchmark the xCT calibration performance on the biological phantom. This study underlined the potential benefits and constraints of utilising such a novelex-vivophantom for inter-centre surveys in future.

MRI-only based material mass density and relative stopping power estimation via deep learning for proton therapy: a preliminary study

Magnetic Resonance Imaging (MRI) is increasingly being used in treatment planning due to its superior soft tissue contrast, which is useful for tumor and soft tissue delineation compared to computed tomography (CT). However, MRI cannot directly provide mass density or relative stopping power (RSP) maps, which are required for calculating proton radiotherapy doses. Therefore, the integration of artificial intelligence (AI) into MRI-based treatment planning to estimate mass density and RSP directly from MRI has generated significant interest. A deep learning (DL) based framework was developed to establish a voxel-wise correlation between MR images and mass density as well as RSP. To facilitate the study, five tissue substitute phantoms were created, representing different tissues such as skin, muscle, adipose tissue, 45% hydroxyapatite (HA), and spongiosa bone. The composition of these phantoms was based on information from ICRP reports. Additionally, two animal tissue phantoms, simulating pig brain and liver, were prepared for DL training purposes. The phantom study involved the development of two DL models. The first model utilized clinical T1 and T2 MRI scans as input, while the second model incorporated zero echo time (ZTE) MRI scans. In the patient application study, two more DL models were trained: one using T1 and T2 MRI scans as input, and another model incorporating synthetic dual-energy computed tomography (sDECT) images to provide accurate bone tissue information. The DECT empirical model was used as a reference to evaluate the proposed models in both phantom and patient application studies. The DECT empirical model was selected as the reference for evaluating the proposed models in both phantom and patient application studies. In the phantom study, the DL model based on T1, and T2 MRI scans demonstrated higher accuracy in estimating mass density and RSP for skin, muscle, adipose tissue, brain, and liver. The mean absolute percentage errors (MAPE) were 0.42%, 0.14%, 0.19%, 0.78%, and 0.26% for mass density, and 0.30%, 0.11%, 0.16%, 0.61%, and 0.23% for RSP, respectively. The DL model incorporating ZTE MRI further improved the accuracy of mass density and RSP estimation for 45% HA and spongiosa bone, with MAPE values of 0.23% and 0.09% for mass density, and 0.19% and 0.07% for RSP, respectively. These results demonstrate the feasibility of using an MRI-only approach combined with DL methods for mass density and RSP estimation in proton therapy treatment planning. By employing this approach, it is possible to obtain the necessary information for proton radiotherapy directly from MRI scans, eliminating the need for additional imaging modalities.

Characterization of a novel proton-CT scanner based on Silicon and LaBr3\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_3$$\\end{document}(Ce) detectors

Treatment planning systems at proton-therapy centres entirely use X-ray computed tomography (CT) as primary imaging technique to infer the proton treatment doses to tumour and healthy tissues. However, proton stopping powers in the body, as derived from X-ray images, suffer from important proton-range uncertainties. In order to reduce this uncertainty in range, one could use proton-CT images instead. The main goal of this work is to test the capabilities of a newly-developed proton-CT scanner, based on the use of a set of tracking detectors and a high energy resolution scintillator for the residual energy of the protons. Different custom-made phantoms were positioned at the field of view of the scanner and were irradiated with protons at the CCB proton-therapy center in Krakow. We measured with the phantoms at different angles and produced sinograms that were used to obtain reconstructed images by Filtered Back-Projection. The obtained images were used to determine the capabilities of our scanner in terms of spatial resolution and proton Relative Stopping Power (RSP) mapping and validate its use as proton-CT scanner. The results show that the scanner can produce medium-high quality images, with spatial resolution better than 2 mm in radiography, below 3 mm in tomography and resolving power in the RSP comparable to other state-of-the-art pCT scanners.

1202: Comparison of performance of spectral x-ray CT and proton CT for relative stopping power estimation

1202: Comparison of performance of spectral x-ray CT and proton CT for relative stopping power estimation

Impact of iodinated oil in proton therapy on relative stopping power of liver post-cTACE

Objective. Conventional transarterial chemoembolization (cTACE) is a common treatment for hepatocellular carcinoma (HCC), often with unsatisfactory local controls. Combining cTACE with radiotherapy shows a promise for unresectable large HCC, with proton therapy preserving healthy liver tissue. However, the proton therapy benefits are subject to the accuracy of tissue relative stopping power (RSP) prediction. The RSP values are typically derived from computed tomography (CT) images using stoichiometric calibration. Lipiodol deposition significantly increases CT numbers in liver regions of post-cTACE. Hence, it is necessary to evaluate the accuracy of RSP in liver regions of post-cTACE. Approach. Liver, water, and iodinated oil samples were prepared. Some liver samples contained iodinated oil. The water equivalent path length (WEPL) of sample was measured through the pullbacks of spread-out Bragg peak (SOBP) depth-dose profiles scanned in a water tank with and without sample in the beam path. Measured RSP values were compared to estimated RSP values derived from the CT number based on the stoichiometric calibration method. Main results. The measured RSP of water was 0.991, confirming measurement system calibration. After removing the RSP contribution from container walls, the pure iodinated oil and liver samples had RSP values of 1.12 and 1.06, while the liver samples mixed with varying oil volumes (5 ml, 10 ml, 15 ml) showed RSP values of 1.05, 1.05 and 1.06. Using the stoichiometric calibration method, pure iodinated oil and liver samples had RSP values of 2.79 and 1.06. Liver samples mixed with iodinated oil (5 ml, 10 ml, 15 ml) had calculated RSP values of 1.21, 1.34, and 1.46. The RSP discrepancy reached 149.1% for pure iodinated oil. Significance. Iodinated oil notably raises CT numbers in liver tissue. However, there is almost no effect on its RSP value. Proton treatment of post-cTACE HCC patients can therefore be overshooting if no proper measures are taken against this specific effect.

Technical note: Radiological clinical equivalence for phantom materials in carbon ion therapy.

As carbon ion radiotherapy increases in use, there are limited phantom materials for heterogeneous or anthropomorphic phantom measurements. This work characterized the radiological clinical equivalence of several phantom materials in a therapeutic carbon ion beam. Eight materials were tested for radiological material-equivalence in a carbon ion beam. The materials were computed tomography (CT)-scanned to obtain Hounsfield unit (HU) values, then irradiated in a monoenergetic carbon ion beam to determine relative linear stopping power (RLSP). The corresponding HU and RLSP for each phantom material were compared to clinical carbon ion calibration curves. For absorbed dose comparison, ion chamber measurements were made in the center of a carbon ion spread-out Bragg peak (SOBP) in water and in the phantom material, evaluating whether the material perturbed the absorbed dose measurement beyond what was predicted by the HU-RLSP relationship. Polyethylene, solid water (Gammex and Sun Nuclear), Blue Water (Standard Imaging), and Techtron HPV had measured RLSP values that agreed within ±4.2% of RLSP values predicted by the clinical calibration curve. Measured RLSP for acrylic was 7.2% different from predicted. The agreement for balsa wood and cork varied between samples. Ion chamber measurements in the phantom materials were within 0.1% of ion chamber measurements in water for most materials (solid water, Blue Water, polyethylene, and acrylic), and within 1.9% for the rest of the materials (balsa wood, cork, and Techtron HPV). Several phantom materials (Blue Water, polyethylene, solid water [Gammex and Sun Nuclear], and Techtron HPV) are suitable for heterogeneous phantom measurements for carbon ion therapy. Low density materials should be carefully characterized due to inconsistencies between samples.

First experimental time-of-flight-based proton radiography using low gain avalanche diodes

Objective. Ion computed tomography (iCT) is an imaging modality for the direct determination of the relative stopping power (RSP) distribution within a patient’s body. Usually, this is done by estimating the path and energy loss of ions traversing the scanned volume utilising a tracking system and a separate residual energy detector. This study, on the other hand, introduces the first experimental study of a novel iCT approach based on time-of-flight (TOF) measurements, the so-called Sandwich TOF-iCT concept, which in contrast to any other iCT systems, does not require a residual energy detector for the RSP determination. Approach. A small Sandwich TOF-iCT demonstrator was built based on low gain avalanche diodes (LGADs), which are 4D-tracking detectors that allow to simultaneously measure the particle position and time-of-arrival with a precision better than 100 μm and 100 ps, respectively. Using this demonstrator, the material and energy-dependent TOF was measured for several homogeneous PMMA slabs in order to calibrate the acquired TOF against the corresponding water equivalent thickness (WET). With this calibration, two proton radiographs (pRads) of a small aluminium stair phantom were recorded at MedAustron using 83 MeV and 100.4 MeV protons. Main results. Due to the simplified WET calibration models used in this very first experimental study of this novel approach, the difference between the measured and theoretical WET ranged between 37.09% and 51.12%. Nevertheless, the first TOF-based pRad was successfully recorded showing that LGADs are suitable detector candidates for Sandwich TOF-iCT. Significance. While the system parameters and WET estimation algorithms require further optimization, this work was an important first step to realize Sandwich TOF-iCT. Due to its compact and cost-efficient design, Sandwich TOF-iCT has the potential to make iCT more feasible and attractive for clinical application, which, eventually, could enhance the treatment planning quality.

Study of energy loss, range, and stopping time for proton in germanium and copper materials

Abstract In this study, the energy loss, range, and stopping time were calculated using the Bethe theory for protons in germanium and copper materials over an energy range of 0.01–1,000 MeV. The mathematical equations were programmed and implemented using MATLAB2013 to obtain the results for the relative mass stopping power and range, illustrated in tables and graphs. The computed results from the Bethe and Ziegler equations were compared with the experimental results from Pstar. The results showed a difference between the experimental and theoretical results at the beginning of the range used in the Bethe equation because it is a quantized equation that is not suitable for low energies. Hence, a correction term was proposed for the Bethe equation to correct the results at low energies, which produced negative values. On the other hand, the Ziegler equation showed good agreement along the range used as it is a semi-empirical equation. The range and stopping time of the incident particle increased with increasing incident particle energy in the mentioned materials.

OPTIma: simplifying calorimetry for proton computed tomography in high proton flux environments

Objective. Proton computed tomography (pCT) offers a potential route to reducing range uncertainties for proton therapy treatment planning, however the current trend towards high current spot scanning treatment systems leads to high proton fluxes which are challenging for existing systems. Here we demonstrate a novel approach to energy reconstruction, referred to as ‘de-averaging’, which allows individual proton energies to be recovered using only a measurement of their integrated energy without the need for spatial information from the calorimeter. Approach. The method is evaluated in the context of the Optimising Proton Therapy through Imaging (OPTIma) system which uses a simple, relatively inexpensive, scintillator-based calorimeter that reports only the integrated energy deposited by all protons within a cyclotron period, alongside a silicon strip based tracking system capable of reconstructing individual protons in a high flux environment. GEANT4 simulations have been performed to examine the performance of such a system at a modern commercial cyclotron facility using a σ ≈ 10 mm beam for currents in the range 10–50 pA at the nozzle. Main results. Apart from low-density lung tissue, a discrepancy of less than 1% on the Relative Stopping Power is found for all other considered tissues when embedded within a 150 mm spherical Perspex phantom in the 10–30 pA current range, and for some tissues even up to 50 pA. Significance. By removing the need for the calorimeter system to provide spatial information, it is hoped that the de-averaging approach can facilitate clinically relevant, cost effective and less complex calorimeter systems for performing high current pCTs.

Energy painting: helium-beam radiography with thin detectors and multiple beam energies

Objective. Compact ion imaging systems based on thin detectors are a promising prospect for the clinical environment since they are easily integrated into the clinical workflow. Their measurement principle is based on energy deposition instead of the conventionally measured residual energy or range. Therefore, thin detectors are limited in the water-equivalent thickness range they can image with high precision. This article presents our energy painting method, which has been developed to render high precision imaging with thin detectors feasible even for objects with larger, clinically relevant water-equivalent thickness (WET) ranges. Approach. A detection system exclusively based on pixelated silicon Timepix detectors was used at the Heidelberg ion-beam therapy center to track single helium ions and measure their energy deposition behind the imaged object. Calibration curves were established for five initial beam energies to relate the measured energy deposition to WET. They were evaluated regarding their accuracy, precision and temporal stability. Furthermore, a 60 mm × 12 mm region of a wedge phantom was imaged quantitatively exploiting the calibrated energies and five different mono-energetic images. These mono-energetic images were combined in a pixel-by-pixel manner by averaging the WET-data weighted according to their single-ion WET precision (SIWP) and the number of contributing ions. Main result. A quantitative helium-beam radiograph of the wedge phantom with an average SIWP of 1.82(5) % over the entire WET interval from 150 mm to 220 mm was obtained. Compared to the previously used methodology, the SIWP improved by a factor of 2.49 ± 0.16. The relative stopping power value of the wedge derived from the energy-painted image matches the result from range pullback measurements with a relative deviation of only 0.4 %. Significance. The proposed method overcomes the insufficient precision for wide WET ranges when employing detection systems with thin detectors. Applying this method is an important prerequisite for imaging of patients. Hence, it advances detection systems based on energy deposition measurements towards clinical implementation.

A denoising method based on deep learning for proton radiograph using energy resolved dose function

Objective. Proton radiograph has been broadly applied in proton radiotherapy which is affected by scattered protons which result in the lower spatial resolution of proton radiographs than that of x-ray images. Traditional image denoising method may lead to the change of water equivalent path length (WEPL) resulting in the lower WEPL measurement accuracy. In this study, we proposed a new denoising method of proton radiographs based on energy resolved dose function curves. Approach. Firstly, the corresponding relationship between the distortion of WEPL characteristic curve, and energy and proportion of scattered protons was established. Then, to improve the accuracy of proton radiographs, deep learning technique was used to remove scattered protons and correct deviated WEPL values. Experiments on a calibration phantom to prove the effectiveness and feasibility of this method were performed. In addition, an anthropomorphic head phantom was selected to demonstrate the clinical relevance of this technology and the denoising effect was analyzed. Main results. The curves of WEPL profiles of proton radiographs became smoother and deviated WEPL values were corrected. For the calibration phantom proton radiograph, the average absolute error of WEPL values decreased from 2.23 to 1.72, the mean percentage difference of all materials of relative stopping power decreased from 1.24 to 0.39, and the average relative WEPL corrected due to the denoising process was 1.06%. In addition, WEPL values correcting were also observed on the proton radiograph for anthropomorphic head phantom due to this denoising process. Significance. The experiments showed that this new method was effective for proton radiograph denoising and had greater advantages than end-to-end image denoising methods, laying the foundation for the implementation of precise proton radiotherapy.

Single energy CT-based mass density and relative stopping power estimation for proton therapy using deep learning method

BackgroundThe number of patients undergoing proton therapy has increased in recent years. Current treatment planning systems (TPS) calculate dose maps using three-dimensional (3D) maps of relative stopping power (RSP) and mass density. The patient-specific maps of RSP and mass density were obtained by translating the CT number (HU) acquired using single-energy computed tomography (SECT) with appropriate conversions and coefficients. The proton dose calculation uncertainty of this approach is 2.5%-3.5% plus 1 mm margin. SECT is the major clinical modality for proton therapy treatment planning. It would be intriguing to enhance proton dose calculation accuracy using a deep learning (DL) approach centered on SECT.ObjectivesThe purpose of this work is to develop a deep learning method to generate mass density and relative stopping power (RSP) maps based on clinical single-energy CT (SECT) data for proton dose calculation in proton therapy treatment.MethodsArtificial neural networks (ANN), fully convolutional neural networks (FCNN), and residual neural networks (ResNet) were used to learn the correlation between voxel-specific mass density, RSP, and SECT CT number (HU). A stoichiometric calibration method based on SECT data and an empirical model based on dual-energy CT (DECT) images were chosen as reference models to evaluate the performance of deep learning neural networks. SECT images of a CIRS 062M electron density phantom were used as the training dataset for deep learning models. CIRS anthropomorphic M701 and M702 phantoms were used to test the performance of deep learning models.ResultsFor M701, the mean absolute percentage errors (MAPE) of the mass density map by FCNN are 0.39%, 0.92%, 0.68%, 0.92%, and 1.57% on the brain, spinal cord, soft tissue, bone, and lung, respectively, whereas with the SECT stoichiometric method, they are 0.99%, 2.34%, 1.87%, 2.90%, and 12.96%. For RSP maps, the MAPE of FCNN on M701 are 0.85%, 2.32%, 0.75%, 1.22%, and 1.25%, whereas with the SECT reference model, they are 0.95%, 2.61%, 2.08%, 7.74%, and 8.62%. ConclusionThe results show that deep learning neural networks have the potential to generate accurate voxel-specific material property information, which can be used to improve the accuracy of proton dose calculation.Advances in knowledgeDeep learning-based frameworks are proposed to estimate material mass density and RSP from SECT with improved accuracy compared with conventional methods.

A systematic approach for calibrating a Monte Carlo code to a treatment planning system for obtaining dose, LET, variable proton RBE and out-of-field dose

Objective. While integration of variable relative biological effectiveness (RBE) has not reached full clinical implementation, the importance of having the ability to recalculate proton treatment plans in a flexible, dedicated Monte Carlo (MC) code cannot be understated . Here we provide a step-wise method for calibrating dose from a MC code to a treatment planning system (TPS), to obtain required parameters for calculating linear energy transfer (LET), variable RBE and in general enabling clinical realistic research studies beyond the capabilities of a TPS. Approach. Initially, Pristine Bragg peaks (PBP) were calculated in both the Eclipse TPS and the FLUKA MC code. A rearranged Bortfeld energy-range relation was applied to the initial energy of the beam to fine-tune the range of the MC code at 80% dose level distal to the PBP. The energy spread was adapted by dividing the TPS range by the MC range for dose level 80%–20% distal to the PBP. Density and relative proton stopping power were adjusted by comparing the TPS and MC for different Hounsfield units. To find the relationship of dose per primary particle from the MC to dose per monitor unit in the TPS, integration was applied to the area of the Bragg curve. The calibration was validated for spread-out Bragg peaks (SOBP) in water and patient treatment plans. Following the validation, variable RBE were calculated using established models. Main results. The PBPs ranges were within ± mm threshold, and a maximum of 5.5% difference for the SOBPs was observed. The patient validation showed excellent dose agreement between the TPS and MC, with the greatest differences for the lung tumor patient. Significance. A procedure for calibrating a MC code to a TPS was developed and validated. The procedure enables MC-based calculation of dose, LET, variable RBE, advanced (secondary) particle tracking and more from treatment plans.

Carbon ion radiography with a composite ionization chamber detector

Range uncertainty in carbon ion therapy can diminish treatment efficacy because it may cause deviation from the planned dose distribution. The precise and accurate determination of relative stopping power (RSP) maps of carbon ions in the patient is a direct solution to this problem. To obtain RSP maps in patients undergoing carbon ion radiography, our team developed a preliminary prototype of a composite ionization chamber detection (CICD) system. The CICD prototype employs synchronously gated integral electronics with the ability to measure the depth-to-dose curve and the beam profile simultaneously. Carbon ion radiography experiments were performed on hemispherical, sloped, and stepped phantoms using the Heavy Ion Medical Machine (HIMM) beam. The beam energy was 190.19 MeV/μ and the beam spot full width at half maximum (FWHM) was 7.42 mm. The radiographic image of the sloped phantom, the thickness prediction accuracy of each pixel (2 mm) is 88.25%, its absolute mean error (AME) is 1.07 mm, and the maximum absolute deviation (MAD) is 2.64 mm. The prediction accuracy of the CICD prototype is mainly affected by electronic noise, with a noise-to-signal ratio (NSR) of about 14.36 dB. Carbon ion radiography simulations were performed in this study using Geant4 software to eliminate the effect of the electronic noise. The thickness prediction accuracy is 98.54%, 98.62%, and 99.07% per pixel for hemispherical, sloped and stepped phantoms, respectively, with AME of 0.09 mm, 0.27 mm, and 0.48 mm. Carbon ion radiography utilizing the CICD prototype scheme has the ability to refine the accuracy and resolution of radiographic images, consequently establishing a scientific foundation for diminishing the effects of range uncertainty and fully exploiting the advantages of precision particle therapy.

A stopping criterion for iterative proton CT image reconstruction based on correlated noise properties.

Whereas filtered back projection algorithms for voxel-based CT image reconstruction have noise properties defined by the filter, iterative algorithms must stop at some point in their convergence and do not necessarily produce consistent noise properties for images with different degrees of heterogeneity. A least-squares iterative algorithm for proton CT (pCT) image reconstruction converges toward a unique solution for relative stopping power (RSP) that optimally fits the protons. We present a stopping criterion that delivers solutions with the property that correlations of RSP noise between voxels are relatively low. This provides a method to produce pCT images with consistent noise properties useful for proton therapy treatment planning, which relies on summing RSP along lines of voxels. Consistent noise properties will also be useful for future studies of image quality using metrics such as contrast to noise ratio, and to compare RSP noise and dose of pCT with other modalities such as dual-energyCT. With simulated and real images with varying heterogeneity from a prototype clinical proton imaging system, we calculate average RSP correlations between voxel pairs in uniform regions-of-interest versus distance between voxels. We define a parameter r, the remaining distance to the unique solution relative to estimated RSP noise, and our stopping criterion is based on r falling below a chosenvalue. We find large correlations between voxels for larger values of r, and anticorrelations for smaller values. For r in the range of 0.5-1, voxels are relatively uncorrelated, and compared to smaller values of r have lower noise with only slight loss of spatialresolution. Iterative algorithms not using a specific metric or rationale for stopping iterations may produce images with an unknown and arbitrary level of convergence or smoothing. We resolve this issue by stopping iterations of a least-squares iterative algorithm when r reaches the range of 0.5-1. This defines a pCT image reconstruction method with consistent statistical properties optimal for clinical use, including for treatment planning with pCTimages.