Valentina Giacometti

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.

Software platform for simulation of a prototype proton CT scanner

Purpose: Proton computed tomography (pCT) is a promising imaging technique to substitute or at least complement x‐ray CT for more accurate proton therapy treatment planning as it allows calculating directly proton relative stopping power from proton energy loss measurements. A proton CT scanner with a silicon‐based particle tracking system and a five‐stage scintillating energy detector has been completed. In parallel a modular software platform was developed to characterize the performance of the proposed pCT. Method: The modular pCT software platform consists of (1) a Geant4‐based simulation modeling the Loma Linda proton therapy beam line and the prototype proton CT scanner, (2) water equivalent path length (WEPL) calibration of the scintillating energy detector, and (3) image reconstruction algorithm for the reconstruction of the relative stopping power (RSP) of the scanned object. In this work, each component of the modular pCT software platform is described and validated with respect to experimental data and benchmarked against theoretical predictions. In particular, the RSP reconstruction was validated with both experimental scans, water column measurements, and theoretical calculations. Results: The results show that the pCT software platform accurately reproduces the performance of the existing prototype pCT scanner with a RSP agreement between experimental and simulated values to better than 1.5%. Conclusions: The validated platform is a versatile tool for clinical proton CT performance and application studies in a virtual setting. The platform is flexible and can be modified to simulate not yet existing versions of pCT scanners and higher proton energies than those currently clinically available.

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.

Results from a Prototype Proton-CT Head Scanner

Abstract We are exploring low-dose proton radiography and computed tomography (pCT) as techniques to improve the accuracy of proton treatment planning and to provide artifact-free images for verification and adaptive therapy at the time of treatment. Here we report on comprehensive beam test results with our prototype pCT head scanner. The detector system and data acquisition attain a sustained rate of more than a million protons individually measured per second, allowing a full CT scan to be completed in six minutes or less of beam time. In order to assess the performance of the scanner for proton radiography as well as computed tomography, we have performed numerous scans of phantoms at the Northwestern Medicine Chicago Proton Center including a custom phantom designed to assess the spatial resolution, a phantom to assess the measurement of relative stopping power, and a dosimetry phantom. Some images, performance, and dosimetry results from those phantom scans are presented together with a description of the instrument, the data acquisition system, and the calibration methods.

Development of a high resolution voxelised head phantom for medical physics applications

Computational anthropomorphic phantoms have become an important investigation tool for medical imaging and dosimetry for radiotherapy and radiation protection. The development of computational phantoms with realistic anatomical features contribute significantly to the development of novel methods in medical physics. For many applications, it is desirable that such computational phantoms have a real-world physical counterpart in order to verify the obtained results. In this work, we report the development of a voxelised phantom, the HIGH_RES_HEAD, modelling a paediatric head based on the commercial phantom 715-HN (CIRS). HIGH_RES_HEAD is unique for its anatomical details and high spatial resolution (0.18×0.18mm2 pixel size). The development of such a phantom was required to investigate the performance of a new proton computed tomography (pCT) system, in terms of detector technology and image reconstruction algorithms. The HIGH_RES_HEAD was used in an ad-hoc Geant4 simulation modelling the pCT system. The simulation application was previously validated with respect to experimental results. When compared to a standard spatial resolution voxelised phantom of the same paediatric head, it was shown that in pCT reconstruction studies, the use of the HIGH_RES_HEAD translates into a reduction from 2% to 0.7% of the average relative stopping power difference between experimental and simulated results thus improving the overall quality of the head phantom simulation. The HIGH_RES_HEAD can also be used for other medical physics applications such as treatment planning studies. A second version of the voxelised phantom was created that contains a prototypic base of skull tumour and surrounding organs at risk.

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.

Track reconstruction with the silicon strip tracker of the proton CT Phase 2 scanner

We present algorithms for the proton track reconstruction for the Phase II proton CT scanner designed and built by the pCT collaboration. The entrance and exit telescopes of the scanner consist of four planes of silicon sensors with horizontal and vertical strips for proton track reconstruction. The goal of the track reconstruction is to allow usage of submillimeter voxels for subsequent image reconstruction. A dedicated algorithm was developed that recovers protons that traversed the dead areas at sensor joints and the dead or masked strips. The track reconstruction was tested on the data collected with 200 MeV protons from the proton synchrotron of the Loma Linda University Medical Center. Test results show that the reconstruction accuracy is close to the geometric limit of the strip sensors and multiple Coulomb scattering of the protons in air, sensors and scanning object.

Dosimetric evaluation of proton CT using a prototype proton CT scanner

These instructions provide guidelines for preparing manuscripts for submission to the Conference Record (CR) of the 2016 IEEE Nuclear Science Symposium and Medical Imaging Conference. If you are using Microsoft Word to prepare your manuscript, you should use this document as a template. Define all symbols used in the abstract. Do not cite references in the abstract. Do not delete the blank line immediately above the abstract; it sets the footnote at the bottom of this column. Proton computed tomography has been suggested as an imaging technique alternative to x-ray CT for proton therapy treatment planning and image guidance. Dose, image quality, and range prediction accuracy are closely related and the calculation of the dose delivered during a proton CT scan is crucial for characterizing a proton CT scanner. Here we report on the dosimetric evaluation of proton CT scanner using a prototype built by the pCT collaboration between Loma Linda University, the University of California Santa Cruz, and Baylor University. The Catphan CTP554 16 cm acrylic dose phantom, representing a human head, was used to measure the dose to water during a typical proton CT scan at the Northwestern Medicine Chicago Proton Center in Warrenville, Illinois. A Farmer ionization chamber was installed in two locations (centre and periphery of the phantom) and the accumulated charge was measured with the 35040 Advanced Therapy Dosimeter (Fluke Biomedical). The proton CT scanner was exposed to a uniform beam profile of 200 MeV protons on the uniform scanning clinical proton beam line. At a rate of 1 M protons sec it takes 7 minutes in order collect about 400 M protons during a continuous 360 deg rotation. The proton fluence at the level of upstream inner tracker planes during the 7 minutes of exposure was estimated to be 1.4 M protons/cm2. The dose to water was 1.46 mGy in the peripheral location and 1.48 mGy in the central location, demonstrating that it is practically uniform across the phantom. According to these results, we can conclude that, based on an initial dosimetric characterization of a prototype pCT scanner conducted on a clinical proton beam line, proton CT is a promising modality for low-dose image guidance and adaptive proton therapy.

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.