A comparison of direct reconstruction algorithms in proton computed tomography
- 6 March 2020
- journal article
- research article
- Published by IOP Publishing in Physics in Medicine & Biology
- Vol. 65 (10), 105010
- https://doi.org/10.1088/1361-6560/ab7d53
Abstract
Several direct algorithms have been proposed to take into account the non-linear path of protons in the reconstruction of a proton CT (pCT) image. This paper presents a comparison between five of them, in terms of spatial resolution and relative stopping power (RSP) accuracy. Our comparison includes (1) a distance-driven algorithm extending the filtered backprojection to non-linear trajectories (DD), (2) an algorithm reconstructing a pCT image from optimized projections (ML), (3) a backproject-then-filter approach using a 2D cone filter (BTF), (4) a differentiated backprojection algorithm based on the inversion of the Hilbert transform (DBP), and (5) an algorithm using a 2D directional ramp filter (DR). We have simulated a single tracking pCT set-up using Geant4 through GATE, with a proton source and two position, direction and energy detectors upstream and downstream from the object. Tracker uncertainties were added on the position and direction measurements. A Catphan 528 phantom and a spiral phantom were simulated to measure the spatial resolution and a Gammex 467 phantom was used for the RSP accuracy. Each proton's trajectory was estimated using a most likely path (MLP) formalism. The spatial resolution was evaluated using the frequency corresponding to an MTF of 10% of its peak value and the RSP accuracy using the mean values in the inserts of the Gammex phantom. In terms of spatial resolution, it was shown that, for ideal trackers, the DR and BTF methods offer a slightly better resolution since each proton is directly binned in the image grid according to its MLP. However, all methods but the ML show comparable resolution when using realistic trackers. Regarding the RSP, three algorithms (DR, DD and BTF) show a mean relative error inside the inserts about 0.1%. As the DR and BTF methods are more computationally expensive, the DD — which allows the same spatial resolution in realistic conditions and the same accuracy — and the DBP — which has a fairly good accuracy (< 0.2%) and allows reconstruction from truncated data — can be used for a reduced reconstruction time.Funding Information
- H2020 Marie Skłodowska-Curie Actions (753370)
- Fondation pour la Recherche Médicale (DIC20161236452)
This publication has 32 references indexed in Scilit:
- Towards task-based assessment of CT performance: System and object MTF across different reconstruction algorithmsMedical Physics, 2012
- Comprehensive analysis of proton range uncertainties related to patient stopping-power-ratio estimation using the stoichiometric calibrationPhysics in Medicine & Biology, 2012
- Range uncertainties in proton therapy and the role of Monte Carlo simulationsPhysics in Medicine & Biology, 2012
- GATE V6: a major enhancement of the GATE simulation platform enabling modelling of CT and radiotherapyPhysics in Medicine & Biology, 2011
- Total variation superiorization schemes in proton computed tomography image reconstructionMedical Physics, 2010
- A more accurate reconstruction system matrix for quantitative proton computed tomographyMedical Physics, 2009
- A maximum likelihood proton path formalism for application in proton computed tomographyMedical Physics, 2008
- Image reconstruction via the finite Hilbert transform of the derivative of the backprojectionMedical Physics, 2007
- Practical cone-beam algorithmJournal of the Optical Society of America A, 1984
- Algebraic Reconstruction Techniques (ART) for three-dimensional electron microscopy and X-ray photographyJournal of Theoretical Biology, 1970