Design optimization of a pixel-based range telescope for proton computed tomography

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.

Exploration of differentiability in a proton computed tomography simulation framework

Objective.Gradient-based optimization using algorithmic derivatives can be a useful technique to improve engineering designs with respect to a computer-implemented objective function. Likewise, uncertainty quantification through computer simulations can be carried out by means of derivatives of the computer simulation. However, the effectiveness of these techniques depends on how 'well-linearizable' the software is. In this study, we assess how promising derivative information of a typical proton computed tomography (pCT) scan computer simulation is for the aforementioned applications.Approach.This study is mainly based on numerical experiments, in which we repeatedly evaluate three representative computational steps with perturbed input values. We support our observations with a review of the algorithmic steps and arithmetic operations performed by the software, using debugging techniques.Main results.The model-based iterative reconstruction (MBIR) subprocedure (at the end of the software pipeline) and the Monte Carlo (MC) simulation (at the beginning) were piecewise differentiable. However, the observed high density and magnitude of jumps was likely to preclude most meaningful uses of the derivatives. Jumps in the MBIR function arose from the discrete computation of the set of voxels intersected by a proton path, and could be reduced in magnitude by a 'fuzzy voxels' approach. The investigated jumps in the MC function arose from local changes in the control flow that affected the amount of consumed random numbers. The tracking algorithm solves an inherently non-differentiable problem.Significance.Besides the technical challenges of merely applying AD to existing software projects, the MC and MBIR codes must be adapted to compute smoother functions. For the MBIR code, we presented one possible approach for this while for the MC code, this will be subject to further research. For the tracking subprocedure, further research on surrogate models is necessary.

Uncertainty-aware spot rejection rate as quality metric for proton therapy using a digital tracking calorimeter

Objective.Proton therapy is highly sensitive to range uncertainties due to the nature of the dose deposition of charged particles. To ensure treatment quality, range verification methods can be used to verify that the individual spots in a pencil beam scanning treatment fraction match the treatment plan. This study introduces a novel metric for proton therapy quality control based on uncertainties in range verification of individual spots.Approach.We employ uncertainty-aware deep neural networks to predict the Bragg peak depth in an anthropomorphic phantom based on secondary charged particle detection in a silicon pixel telescope designed for proton computed tomography. The subsequently predicted Bragg peak positions, along with their uncertainties, are compared to the treatment plan, rejecting spots which are predicted to be outside the 95% confidence interval. The such-produced spot rejection rate presents a metric for the quality of the treatment fraction.Main results.The introduced spot rejection rate metric is shown to be well-defined for range predictors with well-calibrated uncertainties. Using this method, treatment errors in the form of lateral shifts can be detected down to 1 mm after around 1400 treated spots with spot intensities of 1 × 107protons. The range verification model used in this metric predicts the Bragg peak depth to a mean absolute error of 1.107 ± 0.015 mm.Significance.Uncertainty-aware machine learning has potential applications in proton therapy quality control. This work presents the foundation for future developments in this area.

Feasibility study of a proton CT system based on 4D-tracking and residual energy determination via time-of-flight

Objective.For dose calculations in ion beam therapy, it is vital to accurately determine the relative stopping power (RSP) distribution within the treatment volume. A suitable imaging modality to achieve the required RSP accuracy is proton computed tomography (pCT), which usually uses a tracking system and a separate residual energy (or range) detector to directly measure the RSP distribution. This work investigates the potential of a novel pCT system based on a single detector technology, namely low gain avalanche detectors (LGADs). LGADs are fast 4D-tracking detectors, which can be used to simultaneously measure the particle position and time with precise timing and spatial resolution. In contrast to standard pCT systems, the residual energy is determined via a time-of-flight (TOF) measurement between different 4D-tracking stations.Approach.To show the potential of using 4D-tracking for proton imaging, we studied and optimized the design parameters for a realistic TOF-pCT system using Monte Carlo simulations. We calculated the RSP accuracy and RSP resolution inside the inserts of the CTP404 phantom and compared the results to a simulation of an ideal pCT system.Main results.After introducing a dedicated calibration procedure for the TOF calorimeter, RSP accuracies less than 0.6% could be achieved. We also identified the design parameters with the strongest impact on the RSP resolution and proposed a strategy to further improve the image quality.Significance.This comprehensive study of the most important design aspects for a novel TOF-pCT system could help guide future hardware developments and, once implemented, improve the quality of treatment planning in ion beam therapy.

The role of Monte Carlo simulation in understanding the performance of proton computed tomography

Proton computed tomography (pCT) is a promising tomographic imaging modality allowing direct reconstruction of proton relative stopping power (RSP) required for proton therapy dose calculation. In this review article, we aim at highlighting the role of Monte Carlo (MC) simulation in pCT studies. After describing the requirements for performing proton computed tomography and the various pCT scanners actively used in recent research projects, we present an overview of available MC simulation platforms. The use of MC simulations in the scope of investigations of image reconstruction, and for the evaluation of optimal RSP accuracy, precision and spatial resolution omitting detector effects is then described. In the final sections of the review article, we present specific applications of realistic MC simulations of an existing pCT scanner prototype, which we describe in detail.

The accuracy of helium ion CT based particle therapy range prediction: an experimental study comparing different particle and x-ray CT modalities

This work provides a quantitative assessment of helium ion CT (HeCT) for particle therapy treatment planning. For the first time, HeCT based range prediction accuracy in a heterogeneous tissue phantom is presented and compared to single-energy x-ray CT (SECT), dual-energy x-ray CT (DECT) and proton CT (pCT). HeCT and pCT scans were acquired using the US pCT collaboration prototype particle CT scanner at the Heidelberg Ion-Beam Therapy Center. SECT and DECT scans were done with a Siemens Somatom Definition Flash and converted to RSP. A Catphan CTP404 module was used to study the RSP accuracy of HeCT. A custom phantom of 20 cm diameter containing several tissue equivalent plastic cubes was used to assess the spatial resolution of HeCT and compare it to DECT. A clinically realistic heterogeneous tissue phantom was constructed using cranial slices from a pig head placed inside a cylindrical phantom (ø150 mm). A proton beam (84.67 mm range) depth-dose measurement was acquired using a stack of GafchromicTM EBT-XD films in a central dosimetry insert in the phantom. CT scans of the phantom were acquired with each modality, and proton depth-dose estimates were simulated based on the reconstructions. The RSP accuracy of HeCT for the plastic phantom was found to be 0.3 ± 0.1%. The spatial resolution for HeCT of the cube phantom was 5.9 ± 0.4 lp cm-1for central, and 7.6 ± 0.8 lp cm-1for peripheral cubes, comparable to DECT spatial resolution (7.7 ± 0.3 lp cm-1and 7.4 ± 0.2 lp cm-1, respectively). For the pig head, HeCT, SECT, DECT and pCT predicted range accuracy was 0.25%, -1.40%, -0.45% and 0.39%, respectively. In this study, HeCT acquired with a prototype system showed potential for particle therapy treatment planning, offering RSP accuracy, spatial resolution, and range prediction accuracy comparable to that achieved with a commercial DECT scanner. Still, technical improvements of HeCT are needed to enable clinical implementation.

Investigating particle track topology for range telescopes in particle radiography using convolutional neural networks

BACKGROUND

Proton computed tomography (pCT) and radiography (pRad) are proposed modalities for improved treatment plan accuracy and in situ treatment validation in proton therapy. The pCT system of the Bergen pCT collaboration is able to handle very high particle intensities by means of track reconstruction. However, incorrectly reconstructed and secondary tracks degrade the image quality. We have investigated whether a convolutional neural network (CNN)-based filter is able to improve the image quality.

MATERIAL AND METHODS

The CNN was trained by simulation and reconstruction of tens of millions of proton and helium tracks. The CNN filter was then compared to simple energy loss threshold methods using the Area Under the Receiver Operating Characteristics curve (AUROC), and by comparing the image quality and Water Equivalent Path Length (WEPL) error of proton and helium radiographs filtered with the same methods.

RESULTS

The CNN method led to a considerable improvement of the AUROC, from 74.3% to 97.5% with protons and from 94.2% to 99.5% with helium. The CNN filtering reduced the WEPL error in the helium radiograph from 1.03 mm to 0.93 mm while no improvement was seen in the CNN filtered pRads.

CONCLUSION

The CNN improved the filtering of proton and helium tracks. Only in the helium radiograph did this lead to improved image quality.

Helium radiography with a digital tracking calorimeter-a Monte Carlo study for secondary track rejection

Radiation therapy using protons and heavier ions is a fast-growing therapeutic option for cancer patients. A clinical system for particle imaging in particle therapy would enable online patient position verification, estimation of the dose deposition through range monitoring and a reduction of uncertainties in the calculation of the relative stopping power of the patient. Several prototype imaging modalities offer radiography and computed tomography using protons and heavy ions. A Digital Tracking Calorimeter (DTC), currently under development, has been proposed as one such detector. In the DTC 43 longitudinal layers of laterally stacked ALPIDE CMOS monolithic active pixel sensor chips are able to reconstruct a large number of simultaneously recorded proton tracks. In this study, we explored the capability of the DTC for helium imaging which offers favorable spatial resolution over proton imaging. Helium ions exhibit a larger cross section for inelastic nuclear interactions, increasing the number of produced secondaries in the imaged object and in the detector itself. To that end, a filtering process able to remove a large fraction of the secondaries was identified, and the track reconstruction process was adapted for helium ions. By filtering on the energy loss along the tracks, on the incoming angle and on the particle ranges, 97.5% of the secondaries were removed. After passing through 16 cm water, 50.0% of the primary helium ions survived; after the proposed filtering 42.4% of the primaries remained; finally after subsequent image reconstruction 31% of the primaries remained. Helium track reconstruction leads to more track matching errors compared to protons due to the increased available focus strength of the helium beam. In a head phantom radiograph, the Water Equivalent Path Length error envelope was 1.0 mm for helium and 1.1 mm for protons. This accuracy is expected to be sufficient for helium imaging for pre-treatment verification purposes.

A Geant4 simulation for three-dimensional proton imaging of microscopic samples

Proton imaging can be carried out on microscopic samples by focusing the beam to a diameter ranging from a few micrometers down to a few tens of nanometers, depending on the required beam intensity and spatial resolution. Three-dimensional (3D) imaging by tomography is obtained from proton transmission (STIM: Scanning Transmission Ion Microscopy) and/or X-ray emission (PIXE: Particle Induced X-ray Emission). In these experiments, the samples are dehydrated for under vacuum analysis. In situ quantification of nanoparticles has been carried out at CENBG in the frame of nanotoxicology studies, on cells and small organisms used as biological models, especially on Caenorhabditis elegans (C. elegans) nematodes. Tomography experiments reveal the distribution of mass density and chemical content (in g.cm^-3) within the analyzed volume. These density values are obtained using an inversion algorithm. To investigate the effect of this data reduction process, we defined different numerical phantoms, including a (dehydrated) C. elegans phantom whose geometry and density were derived from experimental data. A Monte Carlo simulation based on the Geant4 toolkit was developed. Using different simulation and reconstruction conditions, we compared the resulting tomographic images to the initial numerical reference phantom. A study of the relative error between the reconstructed and the reference images lead to the result that 20 protons per shot can be considered as an optimal number for 3D STIM imaging. Preliminary results for PIXE tomography are also presented, showing the interest of such numerical phantoms to produce reference data for future studies on X-ray signal attenuation in thick samples.