Proof of concept image artifact reduction by energy-modulated proton computed tomography (EMpCT)

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 ourenergy paintingmethod, 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.

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.

Relative stopping power resolution in time-of-flight proton CT

Objective. Proton computed tomography (CT) is similar to x-ray CT but relies on protons rather than photons to form an image. In its most common operation mode, the measured quantity is the amount of energy that a proton has lost while traversing the imaged object from which a relative stopping power map can be obtained via tomographic reconstruction. To this end, a calorimeter which measures the energy deposited by protons downstream of the scanned object has been studied or implemented as energy detector in several proton CT prototypes. An alternative method is to measure the proton’s residual velocity and thus its kinetic energy via the time of flight (TOF) between at least two sensor planes. In this work, we study the RSP resolution, seen as image noise, which can be expected from TOF proton CT systems.

Approach. We rely on physics models on the one hand and statistical models of the relevant uncertainties on the other to derive closed form expressions for the noise in projection images. The TOF measurement error scales with the distance between the TOF sensor planes and is reported as velocity error in ps/m. We use variance reconstruction to obtain noise maps of a water cylinder phantom given the scanner characteristics and additionally reconstruct noise maps for a calorimeter-based proton CT system as reference. We use Monte Carlo simulations to verify our model and to estimate the noise due to multiple Coulomb scattering inside the object. We also provide a comparison of TOF helium and proton CT.

Main results. We find that TOF proton CT with 30 ps m−1 velocity error reaches similar image noise as a calorimeter-based proton CT system with 1% energy error (1 sigma error). A TOF proton CT system with a 50 ps m−1 velocity error produces slightly less noise than a 2% calorimeter system. Noise in a reconstructed TOF proton CT image is spatially inhomogeneous with a marked increase towards the object periphery. Our modelled noise was consistent with Monte Carlo simulated images. TOF helium CT offers lower RSP noise at equal fluence, but is less advantageous at equal imaging dose.

Significance. This systematic study of image noise in TOF proton CT can serve as a guide for future developments of this alternative solution for estimating the residual energy of protons and helium ions after the scanned object.

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.

Comparative accuracy and resolution assessment of two prototype proton computed tomography scanners

BACKGROUND

Improving the accuracy of relative stopping power (RSP) in proton therapy may allow reducing range margins. Proton computed tomography (pCT) has been shown to provide state-of-the-art RSP accuracy estimation, and various scanner prototypes have recently been built. The different approaches used in scanner design are expected to impact spatial resolution and RSP accuracy.

PURPOSE

The goal of this study was to perform the first direct comparison, in terms of spatial resolution and RSP accuracy, of two pCT prototype scanners installed at the same facility and by using the same image reconstruction algorithm.

METHODS

A phantom containing cylindrical inserts of known RSP was scanned at the phase-II pCT prototype of the U.S. pCT collaboration and at the commercially oriented ProtonVDA scanner. Following distance-driven binning filtered backprojection reconstruction, the radial edge spread function of high-density inserts was used to estimate the spatial resolution. RSP accuracy was evaluated by the mean absolute percent error (MAPE) over the inserts. No direct imaging dose estimation was possible, which prevented a comparison of the two scanners in terms of RSP noise.

RESULTS

In terms of RSP accuracy, both scanners achieved the same MAPE of 0.72% when excluding the porous sinus insert from the evaluation. The ProtonVDA scanner reached a better overall MAPE when all inserts and the body of the phantom were accounted for (0.81%), compared to the phase-II scanner (1.14%). The spatial resolution with the phase-II scanner was found to be 0.61 lp/mm, while for the ProtonVDA scanner somewhat lower at 0.46 lp/mm.

CONCLUSIONS

The comparison between two prototype pCT scanners operated in the same clinical facility showed that they both fulfill the requirement of an RSP accuracy of about 1%. Their spatial resolution performance reflects the different design choices of either a scanner with full tracking capabilities (phase-II) or of a more compact tracker system which only provides the positions of protons but not their directions (ProtonVDA). This article is protected by copyright. All rights reserved.

Evaluation of the impact of a scanner prototype on proton CT and helium CT image quality and dose efficiency with Monte Carlo simulation

Objective.The use of ion computed tomography (CT) promises to yield improved relative stopping power (RSP) estimation as input to particle therapy treatment planning. Recently, proton CT (pCT) has been shown to yield RSP accuracy on par with state-of-the-art x-ray dual energy CT. There are however concerns that the lower spatial resolution of pCT compared to x-ray CT may limit its potential, which has spurred interest in the use of helium ion CT (HeCT). The goal of this study was to investigate image quality of pCT and HeCT in terms of noise, spatial resolution, RSP accuracy and imaging dose using a detailed Monte Carlo (MC) model of an existing ion CT prototype.Approach.Three phantoms were used in simulated pCT and HeCT scans allowing estimation of noise, spatial resolution and the scoring of dose. An additional phantom was used to evaluate RSP accuracy. The imaging dose required to achieve the same image noise in a water and a head phantom was estimated at both native spatial resolution, and in a scenario where the HeCT spatial resolution was reduced and matched to that of pCT using Hann windowing of the reconstruction filter. A variance reconstruction formalism was adapted to account for Hann windowing.Main results.We confirmed that the scanner prototype would produce higher spatial resolution for HeCT than pCT by a factor 1.8 (0.86 lp mm^-1versus 0.48 lp mm^-1at the center of a 20 cm water phantom). At native resolution, HeCT required a factor 2.9 more dose than pCT to achieve the same noise, while at matched resolution, HeCT required only 38% of the pCT dose. Finally, RSP mean absolute percent error (MAPE) was found to be 0.59% for pCT and 0.67% for HeCT.Significance.This work compared the imaging performance of pCT and HeCT when using an existing scanner prototype, with the spatial resolution advantage of HeCT coming at the cost of increased dose. When matching spatial resolution via Hann windowing, HeCT had a substantial dose advantage. Both modalities provided state-of-the-art RSP MAPE. HeCT might therefore help reduce the dose exposure of patients with comparable image noise to pCT, enhanced spatial resolution and acceptable RSP accuracy at the same time.

Energy-adaptive calculation of the most likely path in proton CT

This note addresses an issue faced by every proton computed tomography (CT) reconstruction software: the modelling and the parametrisation of the multiple Coulomb scattering power for the estimation of the most likely path (MLP) of each proton. The conventional approach uses a polynomial model parameterised as a function of depth for a given initial beam energy. This makes it cumbersome to implement a software that works for proton CT data acquired with an arbitrary beam energy or with energy modulation during acquisition. We propose a simple way to parametrise the scattering power based on the measured proton CT list-mode data only and derive a compact expression for the MLP based on a conventional MLP model. Our MLP does not require any parameter. The method assumes the imaged object to be homogeneous, as most conventional MLPs, but requires no information about the material as opposed to most conventional MLP expressions which often assume water to infer energy loss. Instead, our MLP automatically adapts itself to the energy-loss which actually occurred in the object and which is one of the measurements required for proton CT reconstruction. We validate our MLP method numerically and find excellent agreement with conventional MLP methods.

An empirical artifact correction for proton computed tomography

PURPOSE

To reduce image artifacts of proton computed tomography (pCT) from a preclinical scanner, for imaging of the relative stopping power (RSP) needed for particle therapy treatment planning using a simple empirical artifact correction method.

METHODS

We adapted and employed a correction method previously used for beam-hardening correction in x-ray CT which makes use of a single scan of a custom-built homogeneous phantom with known RSP. Exploiting the linearity of the filtered backprojection operation, a function was found which corrects water-equivalent path lengths (RSP line integrals) in experimental scans using a prototype pCT scanner. The correction function was applied to projection values of subsequent scans of a homogeneous water phantom, a sensitometric phantom with various inserts and an anthropomorphic head phantom. Data were acquired at two different incident proton energies to test the robustness of the method.

RESULTS

Inaccuracies in the detection process caused an offset and known ring artifacts in the water phantom which were considerably reduced using the proposed method. The mean absolute percentage error (MAPE) of mean RSP values of all inserts of the sensitometric phantom and the water phantom was reduced from 0.87% to 0.44% and from 0.86% to 0.48% for the two incident energies respectively. In the head phantom a clear reduction of artifacts was observed.

CONCLUSIONS

Image artifacts of experimental pCT scans with a prototype scanner could substantially be reduced both in homogeneous, heterogeneous and anthropomorphic phantoms. RSP accuracy was also improved.