R. F. Hurley

Water-equivalent path length calibration of a prototype proton CT scanner

PURPOSE The authors present a calibration method for a prototype proton computed tomography (pCT) scanner. The accuracy of these measurements depends upon careful calibration of the energy detector used to measure the residual energy of the protons that passed through the object. METHODS A prototype pCT scanner with a cesium iodide (CsI(Tl)) crystal calorimeter was calibrated by measuring the calorimeter response for protons of 200 and 100 MeV initial energies undergoing degradation in polystyrene plates of known thickness and relative stopping power (RSP) with respect to water. Calibration curves for the two proton energies were obtained by fitting a second-degree polynomial to the water-equivalent path length versus calorimeter response data. Using the 100 MeV calibration curve, the RSP values for a variety of tissue-equivalent materials were measured and compared to values obtained from a standard depth-dose range shift measurement using a water-tank. A cylindrical water phantom was scanned with 200 MeV protons and its RSP distribution was reconstructed using the 200 MeV calibration. RESULTS It is shown that this calibration method produces measured RSP values of various tissue-equivalent materials that agree to within 0.5% of values obtained using an established water-tank method. The mean RSP value of the water phantom reconstruction was found to be 0.995 ± 0.006. CONCLUSIONS The method presented provides a simple and reliable procedure for calibration of a pCT scanner.

A Fast Experimental Scanner for Proton CT: Technical Performance and First Experience With Phantom Scans

We report on the design, fabrication, and first tests of a tomographic scanner developed for proton computed tomography (pCT) of head-sized objects. After extensive preclinical testing, pCT is intended to be employed in support of proton therapy treatment planning and pre-treatment verification in patients undergoing particle-beam therapy. The scanner consists of two silicon-strip telescopes that track individual protons before and after the phantom, and a novel multistage scintillation detector that measures a combination of the residual energy and range of the proton, from which we derive the water equivalent path length (WEPL) of the protons in the scanned object. The set of WEPL values and the associated paths of protons passing through the object over a 360 ° angular scan are processed by an iterative, parallelizable reconstruction algorithm that runs on modern GP-GPU hardware. In order to assess the performance of the scanner, we have performed tests with 200 MeV protons from the synchrotron of the Loma Linda University Medical Center and the IBA cyclotron of the Northwestern Medicine Chicago Proton Center. Our first objective was calibration of the instrument, including tracker channel maps and alignment as well as the WEPL calibration. Then we performed the first CT scans on a series of phantoms. The very high sustained rate of data acquisition, exceeding one million protons per second, allowed a full 360 ° scan to be completed in less than 10 minutes, and reconstruction of a CATPHAN 404 phantom verified accurate reconstruction of the proton relative stopping power in a variety of materials.

Novel scintillation detector design and performance for proton radiography and computed tomography

PURPOSE Proton computed tomography (pCT) will enable accurate prediction of proton and ion range in a patient while providing the benefit of lower radiation exposure than in x-ray CT. The accuracy of the range prediction is essential for treatment planning in proton or ion therapy and depends upon the detector used to evaluate the water-equivalent path length (WEPL) of a proton passing through the object. A novel approach is presented for an inexpensive WEPL detector for pCT and proton radiography. METHODS A novel multistage detector with an aperture of 10 × 37.5 cm was designed to optimize the accuracy of the WEPL measurements while simplifying detector construction and the performance requirements of its components. The design of the five-stage detector was optimized through simulations based on the geant4 detector simulation toolkit, and the fabricated prototype was calibrated in water-equivalent millimeters with 200 MeV protons in the research beam line of the clinical proton synchrotron at Loma Linda University Medical Center. A special polystyrene step phantom was designed and built to speed up and simplify the calibration procedure. The calibrated five-stage detector was tested in the 200 MeV proton beam as part of the pCT head scanner, using a water phantom and polystyrene slabs to verify the WEPL reconstruction accuracy. RESULTS The beam-test results demonstrated excellent performance of the new detector, in good agreement with the simulation results. The WEPL measurement accuracy is about 3.0 mm per proton in the 0-260 mm WEPL range required for a pCT head scan with a 200 MeV proton beam. CONCLUSIONS The new multistage design approach to WEPL measurements for proton CT and radiography has been prototyped and tested. The test results show that the design is competitive with much more expensive calorimeter and range-counter designs.

An evaluation of spatial resolution of a prototype proton CT scanner

PURPOSE To evaluate the spatial resolution of proton CT using both a prototype proton CT scanner and Monte Carlo simulations. METHODS A custom cylindrical edge phantom containing twelve tissue-equivalent inserts with four different compositions at varying radial displacements from the axis of rotation was developed for measuring the modulation transfer function (MTF) of a prototype proton CT scanner. Two scans of the phantom, centered on the axis of rotation, were obtained with a 200 MeV, low-intensity proton beam: one scan with steps of 4°, and one scan with the phantom continuously rotating. In addition, Monte Carlo simulations of the phantom scan were performed using scanners idealized to various degrees. The data were reconstructed using an iterative projection method with added total variation superiorization based on individual proton histories. Edge spread functions in the radial and azimuthal directions were obtained using the oversampling technique. These were then used to obtain the modulation transfer functions. The spatial resolution was defined by the 10% value of the modulation transfer function (MTF10%) in units of line pairs per centimeter (lp/cm). Data from the simulations were used to better understand the contributions of multiple Coulomb scattering in the phantom and the scanner hardware, as well as the effect of discretization of proton location. RESULTS The radial spatial resolution of the prototype proton CT scanner depends on the total path length, W, of the proton in the phantom, whereas the azimuthal spatial resolution depends both on W and the position, u-, at which the most-likely path uncertainty is evaluated along the path. For protons contributing to radial spatial resolution, W varies with the radial position of the edge, whereas for protons contributing to azimuthal spatial resolution, W is approximately constant. For a pixel size of 0.625 mm, the radial spatial resolution of the image reconstructed from the fully idealized simulation data ranged between 6.31 ± 0.36 lp/cm for W = 197 mm i.e., close to the center of the phantom, and 13.79 ± 0.36 lp/cm for W = 97 mm, near the periphery of the phantom. The azimuthal spatial resolution ranged from 6.99 ± 0.23 lp/cm at u- = 75 mm (near the center) to 11.20 ± 0.26 lp/cm at u- = 20 mm (near the periphery). Multiple Coulomb scattering limits the radial spatial resolution for path lengths greater than approximately 130 mm, and the azimuthal spatial resolution for positions of evaluation greater than approximately 40 mm for W = 199 mm. The radial spatial resolution of the image reconstructed from data from the 4° stepped experimental scan ranged from 5.11 ± 0.61 lp/cm for W = 197 mm to 8.58 ± 0.50 lp/cm for W = 97 mm. In the azimuthal direction, the spatial resolution ranged from 5.37 ± 0.40 lp/cm at u- = 75 mm to 7.27 ± 0.39 lp/cm at u- = 20 mm. The continuous scan achieved the same spatial resolution as that of the stepped scan. CONCLUSIONS Multiple Coulomb scattering in the phantom is the limiting physical factor of the achievable spatial resolution of proton CT; additional loss of spatial resolution in the prototype system is associated with scattering in the proton tracking system and inadequacies of the proton path estimate used in the iterative reconstruction algorithm. Improvement in spatial resolution may be achievable by improving the most likely path estimate by incorporating information about high and low density materials, and by minimizing multiple Coulomb scattering in the proton tracking system.

A novel detector for 2D ion detection in low-pressure gas and its applications

A detector array for single-ion registration was developed and studied. The detector array was operated with various working gases, including propane, air, and water vapor at low gas pressure (a few mbar) in a limited Geiger mode, providing single ion sensitivity. The patterned detector structure, comprising sub-millimeter hole diameter, millimeter pitch, ground and 2D readout electrodes, were manufactured using standard PCB technology. A glass cathode was utilized to allow the limited discharge mode of operation. Essentially, the new detector is a combination of a hole-type micro-pattern detector and a resistive plate counter operating at reverse polarity in a low pressure gas. It was demonstrated that the unique combined properties of these well established individual detectors in the low pressure gas environment allows for single ion registration through ion-impact induced discharge confined inside an individual cell (hole) of the detector array. Possible applications of the new ion detector in gas chromatography - mass spectrometry and track structure imaging for radiation protection and hadron therapy are briefly discussed.

Space carving and filtered back-projection as preconditioners for proton computed tomography reconstruction

This paper considers how to determine the boundary of an object by comparing two methods: space carving (SC) and filtered back-projection (FBP). Determination of the boundary is an important first step in proton CT. The boundary is used to set up the large sparse linear system of equations, which are then solved to determine the relative stopping power of each element (voxel) in the object. For instance, to find the path of the proton through the object, the entry and exit points on the boundary must first be found. The boundary also becomes important in the iterative reconstruction, as only voxels inside the object are updated, to reduce computational complexity and prevent external artifacts from forming. SC and FBP are compared on speed and boundary results for four cases: (1) noiseless simulated data, (2) noisy simulated data, (3) a real pediatric head phantom, and (4) a real rat. The usefulness and potential of both methods are discussed and future directions are outlined.

Results from a pre-clinical head scanner for proton CT

We report on the first beam test results with our pre-clinical (Phase-II) head scanner developed for proton computed tomography (pCT). After extensive preclinical testing, pCT will be employed in support of proton therapy treatment planning and pre-treatment verification in patients undergoing treatment with particle beam therapy. The Phase-II pCT system consists of two silicon-strip telescopes that track individual protons before and after the phantom or patient, and a novel multistage scintillation detector that measures a combination of the residual energy and range of the proton, from which we derive the water equivalent path length (WEPL) of the protons in the scanned object. The set of WEPL values and associated paths of protons passing through the object over a 360° angular scan is processed by an iterative, parallelizable reconstruction algorithm that runs on modern GP-GPU hardware. In order to assess the performance of the scanner, we have performed beam tests with 200 MeV protons from the synchrotron of the Loma Linda University Medical Center. The first objective was the calibration of the instrument, including tracker channel maps and alignment as well as the WEPL calibration. Then we performed the first CT scans on a series of phantoms. The very high sustained rate of data acquisition, exceeding one million protons per second, allowed a full 360° scan to be completed in less than 10 minutes, and reconstruction of a CATPHAN 404 phantom verified accurate reconstruction of the proton relative stopping power in a variety of materials.

Monte Carlo simulations for the development a clinical proton CT scanner

In support of developing the next phase of a proton computed tomography (pCT) scanner with features making it applicable to clinical situations, much insight can be gained through Monte Carlo simulation using Geant4. Careful simulation of energy/range detectors, as well as silicon strip detectors (SSDs), has offered insights into the physical limitations placed on a pCT scanner. Simulation also offers the opportunity to evaluate different detector design schemes and regimes for reconstructing CT images using protons.

Characterisation of a track structure imaging detector

The spatial distribution of radiation-induced ionisations in sub-cellular structures plays an important role in the initial formation of radiation damage to biological tissues. Using the nanodosimetry approach, physical characteristics of the track structure can be measured and correlated to DNA damage. In this work, a novel nanodosimeter is presented, which detects positive ions produced by radiation interacting with a gas-sensitive volume in order to obtain a high resolution image of the radiation track structure. The characterisation of the detector prototype was performed and different configurations of the device were tested by varying the detector cathode material and the working gas. Preliminary results show that the ionisation cluster size distribution can be obtained with this approach. Further work is planned to improve the detector efficiency in order to register the complete three-dimensional track structure of ionising radiation.

SU-F-I-54: Spatial Resolution Studies in Proton CT Using a Phase-II Prototype Head Scanner

PURPOSE To characterize the modulation transfer function (MTF) of the pre-clinical (phase II) head scanner developed for proton computed tomography (pCT) by the pCT collaboration. To evaluate the spatial resolution achievable by this system. METHODS Our phase II proton CT scanner prototype consists of two silicon telescopes that track individual protons upstream and downstream from a phantom, and a 5-stage scintillation detector that measures a combination of the residual energy and range of the proton. Residual energy is converted to water equivalent path length (WEPL) of the protons in the scanned object. The set of WEPL values and associated paths of protons passing through the object over a 360° angular scan is processed by an iterative parallelizable reconstruction algorithm that runs on GP-GPU hardware. A custom edge phantom composed of water-equivalent polymer and tissue-equivalent material inserts was constructed. The phantom was first simulated in Geant4 and then built to perform experimental beam tests with 200 MeV protons at the Northwestern Medicine Chicago Proton Center. The oversampling method was used to construct radial and azimuthal edge spread functions and modulation transfer functions. The spatial resolution was defined by the 10% point of the modulation transfer function in units of lp/cm. RESULTS The spatial resolution of the image was found to be strongly correlated with the radial position of the insert but independent of the relative stopping power of the insert. The spatial resolution varies between roughly 4 and 6 lp/cm in both the the radial and azimuthal directions depending on the radial displacement of the edge. CONCLUSION The amount of image degradation due to our detector system is small compared with the effects of multiple Coulomb scattering, pixelation of the image and the reconstruction algorithm. Improvements in reconstruction will be made in order to achieve the theoretical limits of spatial resolution.