A more accurate reconstruction system matrix for quantitative proton computed tomography
- 9 September 2009
- journal article
- research article
- Published by Wiley in Medical Physics
- Vol. 36 (10), 4511-4518
- https://doi.org/10.1118/1.3218759
Abstract
An accurate system matrix is required for quantitative protonCT (pCT) image reconstruction with iterative projection algorithms. The system matrix is composed of chord lengths of individual proton path intersections with reconstruction pixels. In previous work, reconstructions were performed assuming constant intersection chord lengths, which led to systematic errors of the reconstructedprotonstopping powers. The purpose of the present work was to introduce a computationally efficient variable intersection chord length in order to improve the accuracy of the system matrix. An analytical expression that takes into account the discrete stepping nature of the pCT most likely path (MLP) reconstruction procedure was created to describe an angle-dependent effective mean chord length function. A pCT dataset was simulated with GEANT4 using a parallel beam of 200 MeV protons intersecting a computerized head phantom consisting of tissue-equivalent materials with known relative stopping power. The phantom stopping powers were reconstructed with the constant chord length, exact chord length, and effective mean chord length approaches, in combination with the algebraic reconstruction technique. Relative stopping power errors were calculated for each anatomical phantom region and compared for the various methods. It was found that the error of approximately 10% in the mean reconstructedstopping power value for a given anatomical region, resulting from a system matrix with a constant chord length, could be reduced to less than 0.5% with either the effective mean chord length or exact chord length approaches. Reconstructions with the effective mean chord length were found to be approximately 20% faster than reconstructions with an exact chord length. The effective mean chord length method provides the possibility for more accurate, computationally efficient quantitative pCT reconstructions.S. N. Penfold, A. B. Rosenfeld, R. W. Schulte and K. E. SchuberKeywords
Funding Information
- Cancer Institute NSW (08/RSA/1-02)
This publication has 16 references indexed in Scilit:
- A maximum likelihood proton path formalism for application in proton computed tomographyMedical Physics, 2008
- A pencil beam algorithm for proton dose calculationsPhysics in Medicine & Biology, 1996
- The calibration of CT Hounsfield units for radiotherapy treatment planningPhysics in Medicine & Biology, 1996
- The relation between x-ray CT numbers and charged particle stopping powers and its significance for radiotherapy treatment planningPhysics in Medicine & Biology, 1983
- Thomas Bayes's Bayesian InferenceJournal of the Royal Statistical Society. Series A (General), 1982
- Treatment planning for heavy ion radiotherapyInternational Journal of Radiation Oncology*Biology*Physics, 1979
- Three methods for reconstructing objects from x rays: A comparative studyComputer Graphics and Image Processing, 1973
- Algebraic Reconstruction Techniques (ART) for three-dimensional electron microscopy and X-ray photographyJournal of Theoretical Biology, 1970
- Multiple Scattering with Energy LossPhysical Review B, 1948
- Zur Theorie des Durchgangs schneller Korpuskularstrahlen durch MaterieAnnalen der Physik, 1930