Proton radiography and proton computed tomography based on time-resolved dose measurements

We present a proof of principle study of proton radiography and proton computed tomography (pCT) based on time-resolved dose measurements. We used a prototype, two-dimensional, diode-array detector capable of fast dose rate measurements, to acquire proton radiographic images expressed directly in water equivalent path length (WEPL). The technique is based on the time dependence of the dose distribution delivered by a proton beam traversing a range modulator wheel in passive scattering proton therapy systems. The dose rate produced in the medium by such a system is periodic and has a unique pattern in time at each point along the beam path and thus encodes the WEPL. By measuring the time dose pattern at the point of interest, the WEPL to this point can be decoded. If one measures the time–dose patterns at points on a plane behind the patient for a beam with sufficient energy to penetrate the patient, the obtained 2D distribution of the WEPL forms an image. The technique requires only a 2D dosimeter array and it uses only the clinical beam for a fraction of second with negligible dose to patient. We first evaluated the accuracy of the technique in determining the WEPL for static phantoms aiming at beam range verification of the brain fields of medulloblastoma patients. Accurate beam ranges for these fields can significantly reduce the dose to the cranial skin of the patient and thus the risk of permanent alopecia. Second, we investigated the potential features of the technique for real-time imaging of a moving phantom. Real-time tumor tracking by proton radiography could provide more accurate validations of tumor motion models due to the more sensitive dependence of proton beam on tissue density compared to x-rays. Our radiographic technique is rapid (~100 ms) and simultaneous over the whole field, it can image mobile tumors without the problem of interplay effect inherently challenging for methods based on pencil beams. Third, we present the reconstructed pCT images of a cylindrical phantom containing inserts of different materials. As for all conventional pCT systems, the method illustrated in this work produces tomographic images that are potentially more accurate than x-ray CT in providing maps of proton relative stopping power (RSP) in the patient without the need for converting x-ray Hounsfield units to proton RSP. All phantom tests produced reasonable results, given the currently limited spatial and time resolution of the prototype detector. The dose required to produce one radiographic image, with the current settings, is ~0.7 cGy. Finally, we discuss a series of techniques to improve the resolution and accuracy of radiographic and tomographic images for the future development of a full-scale detector.

[1]  K. Yokoyama,et al.  Evaluation of Spatial Resolution for Heavy Ion CT System Based on the Measurement of Residual Range Distribution With HIMAC , 2009, IEEE Transactions on Nuclear Science.

[2]  R.W. Schulte,et al.  Proton computed tomography: Update on current status , 2007, 2007 IEEE Nuclear Science Symposium Conference Record.

[3]  Yurii Nesterov,et al.  Smooth minimization of non-smooth functions , 2005, Math. Program..

[4]  Simon Rit,et al.  Filtered backprojection proton CT reconstruction along most likely paths. , 2013, Medical physics.

[5]  Paola Coan,et al.  X-ray phase-contrast imaging: from pre-clinical applications towards clinics , 2013, Physics in medicine and biology.

[6]  Kenneth M. Hanson,et al.  COMPUTED TOMOGRAPHY USING PROTON ENERGY LOSS , 1981, Physics in Medicine and Biology.

[7]  H Paganetti,et al.  Adaptation of GEANT4 to Monte Carlo dose calculations based on CT data. , 2004, Medical physics.

[8]  Barbara Kaser-Hotz,et al.  First proton radiography of an animal patient. , 2004, Medical physics.

[9]  Tianfang Li,et al.  Density resolution of proton computed tomography. , 2005, Medical physics.

[10]  U. Schneider,et al.  Technical note: spatial resolution of proton tomography: impact of air gap between patient and detector. , 2012, Medical physics.

[11]  Cinzia Talamonti,et al.  Proton radiography for clinical applications , 2010 .

[12]  Hsiao-Ming Lu,et al.  Effect of tissue heterogeneity on an in vivo range verification technique for proton therapy , 2012, Physics in medicine and biology.

[13]  Joao Seco,et al.  Sensitivity study of proton radiography and comparison with kV and MV x-ray imaging using GEANT4 Monte Carlo simulations , 2011, Physics in medicine and biology.

[14]  Jan-Jakob Sonke,et al.  Adaptive radiotherapy for lung cancer. , 2010, Seminars in radiation oncology.

[15]  E Pedroni,et al.  Multiple Coulomb scattering and spatial resolution in proton radiography. , 1994, Medical physics.

[16]  On the use of a proton path probability map for proton computed tomography reconstructiona). , 2010, Medical physics.

[17]  A Ghebremedhin,et al.  Water-equivalent path length calibration of a prototype proton CT scanner. , 2012, Medical physics.

[18]  Hsiao-Ming Lu,et al.  A potential method for in vivo range verification in proton therapy treatment , 2008, Physics in medicine and biology.

[19]  Hsiao-Ming Lu,et al.  Optimization of current modulation function for proton spread-out Bragg peak fields. , 2006, Medical physics.

[20]  Cinzia Talamonti,et al.  The PRIMA collaboration: Preliminary results in FBP reconstruction of pCT data , 2013 .

[21]  Emmanuel J. Candès,et al.  NESTA: A Fast and Accurate First-Order Method for Sparse Recovery , 2009, SIAM J. Imaging Sci..

[22]  H. Paganetti.  Range uncertainties in proton therapy and the role of Monte Carlo simulations , 2012, Physics in medicine and biology.

[23]  H. Sadrozinski,et al.  Design of a proton computed tomography system for applications in proton radiation therapy , 2004, 2003 IEEE Nuclear Science Symposium. Conference Record (IEEE Cat. No.03CH37515).

[24]  Oliver Jäkel,et al.  Quantitative carbon ion beam radiography and tomography with a flat-panel detector , 2012, Physics in medicine and biology.

[25]  Tatsuaki Kanai,et al.  Application of heavy-ion CT , 2006, Physics in medicine and biology.

[26]  Marco Riboldi,et al.  Real-time tumour tracking in particle therapy: technological developments and future perspectives. , 2012, The Lancet. Oncology.

[27]  D. Prieels,et al.  Water equivalent path length measurement in proton radiotherapy using time resolved diode dosimetry. , 2011, Medical physics.

[28]  Hyungjoon Ryu,et al.  Density and spatial resolutions of proton radiography using a range modulation technique , 2008, Physics in medicine and biology.

[29]  Joao Seco,et al.  Proof of principle study of the use of a CMOS active pixel sensor for proton radiography. , 2011, Medical physics.

[30]  Wojtek Hajdas,et al.  Construction, test and operation of a proton range radiography system , 2011 .

[31]  J. Verburg,et al.  CT metal artifact reduction method correcting for beam hardening and missing projections , 2012, Physics in medicine and biology.

[32]  S N Penfold,et al.  A maximum likelihood proton path formalism for application in proton computed tomography. , 2008, Medical physics.

[33]  I Rinaldi,et al.  Experimental characterization of a prototype detector system for carbon ion radiography and tomography , 2013, Physics in medicine and biology.

[34]  A M Koehler,et al.  Proton Beam Radiography in Tumor Detection , 1973, Science.

[35]  P Zygmanski,et al.  The measurement of proton stopping power using proton-cone-beam computed tomography. , 2000, Physics in medicine and biology.