Reconstruction for proton computed tomography: a practical approach

Proton Computed Tomography (CT) has important implication for both image-guided diagnosis and radiation therapy. For diagnosis, the fact that the patient dose committed by proton CT compares favorably with that delivered by traditional x-ray CT, for the same density resolution and contrast, may be exploited in dose-critical clinical settings. Proton CT is also the most appropriate imaging method to perform planning and verification of proton-based radiation treatment, since proton stopping power maps constructed by table-based transformation of x-ray CT images only render approximate stopping power estimates. In proton CT, sharp features become blurred by the phenomenon of multiple Coulomb scattering (MCS), resulting in a resolution of around 3 to 5 mm. Studies showed that the spatial resolution of proton radiography and CT can be improved to about 1-2 mm by tracking individual protons in coincidence as they enter and exit the imaged object This paper describes a new practical implementation of this approach. We first bin the captured protons into slots of similar tracks. Optionally, proton energy statistics can be collected within each bin to obtain further parameters for tissue characterization. The envelope of path uncertainty due to MCS can be be modeled as a banana-shaped curve. The 3D reconstruction, using either filter-backprojection or iterative algorithms, can be performed rapidly on graphics hardware, using a slice blurring technique to model the MCS uncertainty curve.