Purpose To demonstrate a proton‐imaging system based on well‐established fast scintillator technology to achieve high performance with low cost and complexity, with the potential of a straightforward translation into clinical use. Methods The system tracks individual protons through one (X, Y) scintillating fiber tracker plane upstream and downstream of the object and into a 13‐cm ‐thick scintillating block residual energy detector. The fibers in the tracker planes are multiplexed into silicon photomultipliers (SiPMs) to reduce the number of electronics channels. The light signal from the residual energy detector is collected by 16 photomultiplier tubes (PMTs). Only four signals from the PMTs are output from each event, which allows for fast signal readout. A robust calibration method of the PMT signal to residual energy has been developed to obtain accurate proton images. The development of patient‐specific scan patterns using multiple input energies allows for an image to be produced with minimal excess dose delivered to the patient. Results The calibration of signals in the energy detector produces accurate residual range measurements limited by intrinsic range straggling. We measured the water‐equivalent thickness (WET) of a block of solid water (physical thickness of 6.10 mm) with a proton radiograph. The mean WET from all pixels in the block was 6.13 cm (SD 0.02 cm). The use of patient‐specific scan patterns using multiple input energies enables imaging with a compact range detector. Conclusions We have developed a prototype clinical proton radiography system for pretreatment imaging in proton radiation therapy. We have optimized the system for use with pencil beam scanning systems and have achieved a reduction of size and complexity compared to previous designs.

Verification of patient-specific proton stopping powers obtained in the patient's treatment position can be used to reduce the distal and proximal margins needed in particle beam planning. Proton radiography can be used as a pretreatment instrument to verify integrated stopping power consistency with the treatment planning CT. Although a proton radiograph is a pixel by pixel representation of integrated stopping powers, the image may also be of high enough quality and contrast to be used for patient alignment. This investigation quantifies the accuracy and image quality of a prototype proton radiography system on a clinical proton delivery system. Methods: We have developed a clinical prototype proton radiography system designed for integration into efficient clinical workflows. We tested the images obtained by this system for water-equivalent thickness (WET) accuracy, image noise, and spatial resolution. We evaluated the WET accuracy by comparing the average WET and rms error in several regions of interest (ROI) on a proton radiograph of a custom peg phantom. We measured the spatial resolution on a CATPHAN Line Pair phantom and a custom edge phantom by measuring the 10% value of the modulation transfer function (MTF). In addition, we tested the ability to detect proton range errors due to anatomical changes in a patient with a customized CIRS pediatric head phantom and inserts of varying WET placed in the posterior fossae of the brain. We took proton radiographs of the phantom with each insert in place and created difference maps between the resulting images. Integrated proton range was measured from an ROI in the difference maps. Results: We measured the WET accuracy of the proton radiographic images to be AE0.2 mm (0.33%) from known values. The spatial resolution of the images was 0.6 lp/mm on the line pair phantom and 1.13 lp/mm on the edge phantom. We were able to detect anatomical changes producing changes in WET as low as 0.6 mm.

Proton CT (pCT) is a promising new imaging technique that can reconstruct relative stopping power (RSP) more accurately than x-ray CT in each cubic millimeter voxel of the patient. This, in turn, will result in better proton range accuracy and, therefore, smaller planned tumor volumes (PTV). The hardware description and some reconstructed images have previously been reported. In a series of two contributions, we focus on presenting the software algorithms that convert pCT detector data to the final reconstructed pCT images for application in proton treatment planning. There were several options on how to accomplish this, and we will describe our solutions at each stage of the data processing chain. In the first paper of this series, we present the data acquisition with the pCT tracking and energy-range detectors and how the data are preprocessed, including the conversion to the well-formatted track information from tracking data and water-equivalent path length from the data of a calibrated multi-stage energy-range detector. These preprocessed data are then used for the initial image formation with an FDK cone-beam CT algorithm. The output of data acquisition, preprocessing, and FDK reconstruction is presented along with illustrative imaging results for two phantoms, including a pediatric head phantom. The second paper in this series will demonstrate the use of iterative solvers in conjunction with the superiorization methodology to further improve the images resulting from the upfront FDK image reconstruction and the implementation of these algorithms on a hybrid CPU/GPU computer cluster. INDEX TERMS Proton computed tomography, data acquisition, preprocessing, initial image formation I. INTRODUCTION T HE interest in the technological development of proton computed tomography, called pCT herein, has increased in recent years due to its potential to reduce the range uncertainty problem in proton therapy [1], [2]. The range of proton beams in a given patient is associated with substantial uncertainties, including the conversion of Hounsfield units (HU) to relative stopping power (RSP) with respect to water, daily variation in the patient setup as well as in the distribution and composition of tissues. These uncertainties

Proton computed tomography (pCT) has high accuracy and dose efficiency in producing spatial maps of the relative stopping power (RSP) required for treatment planning in proton therapy. With fluence-modulated pCT (FMpCT), prescribed noise distributions can be achieved, which allows to decrease imaging dose by employing object-specific dynamically modulated fluence during the acquisition. For FMpCT acquisitions we divide the image into region-of-interest (ROI) and non-ROI volumes. In proton therapy, the ROI volume would encompass all treatment beams. An optimization algorithm then calculates dynamically modulated fluence that achieves low prescribed noise inside the ROI and high prescribed noise elsewhere. It also produces a planned noise distribution, which is the expected noise map for that fluence, as calculated with a Monte Carlo simulation. The optimized fluence can be achieved by acquiring pCT images with grids of intensity modulated pencil beams. In this work, we interfaced the control system of a clinical proton beam line to deliver the optimized fluence. Using three phantoms we acquired images with uniform fluence, with a constant noise prescription, and with an FMpCT task. Image noise distributions as well as fluence maps were compared to the corresponding planned distributions as well as to the prescription. Furthermore, we propose a correction method that removes image artifacts stemming from the acquisition with pencil beams having a spatially varying energy distribution that is not seen in clinical operation. RSP accuracy of FMpCT scans was compared to uniform scans and was found to be comparable to standard pCT scans. While we identified technical improvements for future experimental acquisitions, in particular related to an unexpected pencil beam size reduction and a misalignment of the fluence pattern, agreement with the planned noise was satisfactory and we conclude that FMpCT optimized for specific image noise prescriptions is experimentally feasible.

Purpose: Currently, calculations of proton range in proton therapy patients are based on a conversion of CT Hounsfield units of patient tissues into proton relative stopping power. Uncertainties in this conversion necessitate larger proximal and distal planned target volume margins. Proton CT can potentially reduce these uncertainties by directly measuring proton stopping power. We aim to demonstrate proton CT imaging with complex porcine samples, to analyze in detail three-dimensional regions of interest, and to compare proton stopping powers directly measured by proton CT to those determined from x-ray CT scans. Methods: We have used a prototype proton imaging system with single proton tracking to acquire proton radiography and proton CT images of a sample of porcine pectoral girdle and ribs, and a pig's head. We also acquired close in time x-ray CT scans of the same samples and compared proton stopping power measurements from the two modalities. In the case of the pig's head, we obtained x-ray CT scans from two different scanners and compared results from high-dose and low-dose settings. Results: Comparing our reconstructed proton CT images with images derived from x-ray CT scans, we find agreement within 1% to 2% for soft tissues and discrepancies of up to 6% for compact bone. We also observed large discrepancies, up to 40%, for cavitated regions with mixed content of air, soft tissue, and bone, such as sinus cavities or tympanic bullae. Conclusions:Our images and findings from a clinically realistic proton CT scanner demonstrate the potential for proton CT to be used for low-dose treatment planning with reduced margins. K E Y W O R D Sdigitally reconstructed radiograph, iterative algorithm, proton computed tomography, proton imaging, proton radiography, relative stopping power 7998

This work provides a quantitative assessment of helium ion CT (HeCT) for particle therapy treatment planning. For the first time, HeCT based range prediction accuracy in a heterogeneous tissue phantom is presented and compared to single-energy x-ray CT (SECT), dual-energy x-ray CT (DECT) and proton CT (pCT). HeCT and pCT scans were acquired using the US pCT collaboration prototype particle CT scanner at the Heidelberg Ion-Beam Therapy Center. SECT and DECT scans were done with a Siemens Somatom Definition Flash and converted to RSP. A Catphan CTP404 module was used to study the RSP accuracy of HeCT. A custom phantom of 20 cm diameter containing several tissue equivalent plastic cubes was used to assess the spatial resolution of HeCT and compare it to DECT. A clinically realistic heterogeneous tissue phantom was constructed using cranial slices from a pig head placed inside a cylindrical phantom (ø150 mm). A proton beam (84.67 mm range) depth-dose measurement was acquired using a stack of GafchromicTM EBT-XD films in a central dosimetry insert in the phantom. CT scans of the phantom were acquired with each modality, and proton depth-dose estimates were simulated based on the reconstructions. The RSP accuracy of HeCT for the plastic phantom was found to be 0.3 ± 0.1%. The spatial resolution for HeCT of the cube phantom was 5.9 ± 0.4 lp cm−1 for central, and 7.6 ± 0.8 lp cm−1 for peripheral cubes, comparable to DECT spatial resolution (7.7 ± 0.3 lp cm−1 and 7.4 ± 0.2 lp cm−1, respectively). For the pig head, HeCT, SECT, DECT and pCT predicted range accuracy was 0.25%, −1.40%, −0.45% and 0.39%, respectively. In this study, HeCT acquired with a prototype system showed potential for particle therapy treatment planning, offering RSP accuracy, spatial resolution, and range prediction accuracy comparable to that achieved with a commercial DECT scanner. Still, technical improvements of HeCT are needed to enable clinical implementation.

No abstract

No abstract