Joost M Verburg

CT metal artifact reduction method correcting for beam hardening and missing projections

We present and validate a computed tomography (CT) metal artifact reduction method that is effective for a wide spectrum of clinical implant materials. Projections through low-Z implants such as titanium were corrected using a novel physics correction algorithm that reduces beam hardening errors. In the case of high-Z implants (dental fillings, gold, platinum), projections through the implant were considered missing and regularized iterative reconstruction was performed. Both algorithms were combined if multiple implant materials were present. For comparison, a conventional projection interpolation method was implemented. In a blinded and randomized evaluation, ten radiation oncologists ranked the quality of patient scans on which the different methods were applied. For scans that included low-Z implants, the proposed method was ranked as the best method in 90% of the reviews. It was ranked superior to the original reconstruction (p = 0.0008), conventional projection interpolation (p < 0.0001) and regularized limited data reconstruction (p = 0.0002). All reviewers ranked the method first for scans with high-Z implants, and better as compared to the original reconstruction (p < 0.0001) and projection interpolation (p = 0.004). We conclude that effective reduction of CT metal artifacts can be achieved by combining algorithms tailored to specific types of implant materials.

Energy- and time-resolved detection of prompt gamma-rays for proton range verification

In this work, we present experimental results of a novel prompt gamma-ray detector for proton beam range verification. The detection system features an actively shielded cerium-doped lanthanum(III) bromide scintillator, coupled to a digital data acquisition system. The acquisition was synchronized to the cyclotron radio frequency to separate the prompt gamma-ray signals from the later-arriving neutron-induced background. We designed the detector to provide a high energy resolution and an effective reduction of background events, enabling discrete proton-induced prompt gamma lines to be resolved. Measuring discrete prompt gamma lines has several benefits for range verification. As the discrete energies correspond to specific nuclear transitions, the magnitudes of the different gamma lines have unique correlations with the proton energy and can be directly related to nuclear reaction cross sections. The quantification of discrete gamma lines also enables elemental analysis of tissue in the beam path, providing a better prediction of prompt gamma-ray yields. We present the results of experiments in which a water phantom was irradiated with proton pencil-beams in a clinical proton therapy gantry. A slit collimator was used to collimate the prompt gamma-rays, and measurements were performed at 27 positions along the path of proton beams with ranges of 9, 16 and 23 g cm(-2) in water. The magnitudes of discrete gamma lines at 4.44, 5.2 and 6.13 MeV were quantified. The prompt gamma lines were found to be clearly resolved in dimensions of energy and time, and had a reproducible correlation with the proton depth-dose curve. We conclude that the measurement of discrete prompt gamma-rays for in vivo range verification of clinical proton beams is feasible, and plan to further study methods and detector designs for clinical use.

Proton range verification through prompt gamma-ray spectroscopy

We present an experimental study of a novel method to verify the range of proton therapy beams. Differential cross sections were measured for 15 prompt gamma-ray lines from proton-nuclear interactions with (12)C and (16)O at proton energies up to 150 MeV. These cross sections were used to model discrete prompt gamma-ray emissions along proton pencil-beams. By fitting detected prompt gamma-ray counts to these models, we simultaneously determined the beam range and the oxygen and carbon concentration of the irradiated matter. The performance of the method was assessed in two phantoms with different elemental concentrations, using a small scale prototype detector. Based on five pencil-beams with different ranges delivering 5 × 10(8) protons and without prior knowledge of the elemental composition at the measurement point, the absolute range was determined with a standard deviation of 1.0-1.4 mm. Relative range shifts at the same dose level were detected with a standard deviation of 0.3-0.5 mm. The determined oxygen and carbon concentrations also agreed well with the actual values. These results show that quantitative prompt gamma-ray measurements enable knowledge of nuclear reaction cross sections to be used for precise proton range verification in the presence of tissue with an unknown composition.

Simulation of prompt gamma-ray emission during proton radiotherapy

The measurement of prompt gamma rays emitted from proton-induced nuclear reactions has been proposed as a method to verify in vivo the range of a clinical proton radiotherapy beam. A good understanding of the prompt gamma-ray emission during proton therapy is key to develop a clinically feasible technique, as it can facilitate accurate simulations and uncertainty analysis of gamma detector designs. Also, the gamma production cross-sections may be incorporated as prior knowledge in the reconstruction of the proton range from the measurements. In this work, we performed simulations of proton-induced nuclear reactions with the main elements of human tissue, carbon-12, oxygen-16 and nitrogen-14, using the nuclear reaction models of the GEANT4 and MCNP6 Monte Carlo codes and the dedicated nuclear reaction codes TALYS and EMPIRE. For each code, we made an effort to optimize the input parameters and model selection. The results of the models were compared to available experimental data of discrete gamma line cross-sections. Overall, the dedicated nuclear reaction codes reproduced the experimental data more consistently, while the Monte Carlo codes showed larger discrepancies for a number of gamma lines. The model differences lead to a variation of the total gamma production near the end of the proton range by a factor of about 2. These results indicate a need for additional theoretical and experimental study of proton-induced gamma emission in human tissue.

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.

Range verification of passively scattered proton beams based on prompt gamma time patterns

We propose a proton range verification technique for passive scattering proton therapy systems where spread out Bragg peak (SOBP) fields are produced with rotating range modulator wheels. The technique is based on the correlation of time patterns of the prompt gamma ray emission with the range of protons delivering the SOBP. The main feature of the technique is the ability to verify the proton range with a single point of measurement and a simple detector configuration. We performed four-dimensional (time-dependent) Monte Carlo simulations using TOPAS to show the validity and accuracy of the technique. First, we validated the hadronic models used in TOPAS by comparing simulations and prompt gamma spectrometry measurements published in the literature. Second, prompt gamma simulations for proton range verification were performed for the case of a water phantom and a prostate cancer patient. In the water phantom, the proton range was determined with 2 mm accuracy with a full ring detector configuration for a dose of ~2.5 cGy. For the prostate cancer patient, 4 mm accuracy on range determination was achieved for a dose of ~15 cGy. The results presented in this paper are encouraging in view of a potential clinical application of the technique.

Dosimetric accuracy of proton therapy for chordoma patients with titanium implants.

PURPOSE To investigate dosimetric errors in proton therapy treatment planning due to titanium implants, and to determine how these affect postoperative passively scattered proton therapy for chordoma patients with orthopedic hardware. METHODS The presence of titanium hardware near the tumor may affect the dosimetric accuracy of proton therapy. Artifacts in the computed tomography (CT) scan can cause errors in the proton stopping powers used for dose calculation, which are derived from CT numbers. Also, clinical dose calculation algorithms may not accurately simulate proton beam transport through the implants, which have very different properties as compared to human tissue. The authors first evaluated the impact of these two main issues. Dose errors introduced by metal artifacts were studied using phantoms with and without titanium inserts, and patient scans on which a metal artifact reduction method was applied. Pencil-beam dose calculations were compared to models of nuclear interactions in titanium and Monte Carlo simulations. Then, to assess the overall impact on treatment plans for chordoma, the authors compared the original clinical treatment plans to recalculated dose distributions employing both metal artifact reduction and Monte Carlo methods. RESULTS Dose recalculations of clinical proton fields showed that metal artifacts cause range errors up to 6 mm distal to regions affected by CT artifacts. Monte Carlo simulations revealed dose differences >10% in the high-dose area, and range differences up to 10 mm. Since these errors are mostly local in nature, the large number of fields limits the impact on target coverage in the chordoma treatment plans to a small decrease of dose homogeneity. CONCLUSIONS In the presence of titanium implants, CT metal artifacts and the approximations of pencil-beam dose calculations cause considerable errors in proton dose calculation. The spatial distribution of the errors however limits the overall impact on passively scattered proton therapy for chordoma.

The impact of low-Z and high-Z metal implants in IMRT: A Monte Carlo study of dose inaccuracies in commercial dose algorithms

PURPOSE The aim of the study was to evaluate the dosimetric impact of low-Z and high-Z metallic implants on IMRT plans. METHODS Computed tomography (CT) scans of three patients were analyzed to study effects due to the presence of Titanium (low-Z), Platinum and Gold (high-Z) inserts. To eliminate artifacts in CT images, a sinogram-based metal artifact reduction algorithm was applied. IMRT dose calculations were performed on both the uncorrected and corrected images using a commercial planning system (convolution/superposition algorithm) and an in-house Monte Carlo platform. Dose differences between uncorrected and corrected datasets were computed and analyzed using gamma index (Pγ<1) and setting 2 mm and 2% as distance to agreement and dose difference criteria, respectively. Beam specific depth dose profiles across the metal were also examined. RESULTS Dose discrepancies between corrected and uncorrected datasets were not significant for low-Z material. High-Z materials caused under-dosage of 20%-25% in the region surrounding the metal and over dosage of 10%-15% downstream of the hardware. Gamma index test yielded Pγ<1>99% for all low-Z cases; while for high-Z cases it returned 91% < Pγ<1< 99%. Analysis of the depth dose curve of a single beam for low-Z cases revealed that, although the dose attenuation is altered inside the metal, it does not differ downstream of the insert. However, for high-Z metal implants the dose is increased up to 10%-12% around the insert. In addition, Monte Carlo method was more sensitive to the presence of metal inserts than superposition/convolution algorithm. CONCLUSIONS The reduction in terms of dose of metal artifacts in CT images is relevant for high-Z implants. In this case, dose distribution should be calculated using Monte Carlo algorithms, given their superior accuracy in dose modeling in and around the metal. In addition, the knowledge of the composition of metal inserts improves the accuracy of the Monte Carlo dose calculation significantly.

Range verification of passively scattered proton beams using prompt gamma-ray detection.

We performed an experimental study to verify the range of passively scattered proton beams by detecting prompt gamma-rays emitted from proton-nuclear interactions. A method is proposed using a single scintillation detector positioned near the distal end of the irradiated target. Lead shielding was used to attenuate gamma-rays emitted along most of the entrance path of the beam. By synchronizing the prompt gamma-ray detector to the rotation of the range modulation wheel, the relation between the gamma emission from the distal part of the target and the range of the incident proton beam was determined. In experiments with a water phantom and an anthropomorphic head phantom, this relation was found to be sensitive to range shifts that were introduced. The wide opening angle of the detector enabled a sufficient signal-to-background ratio to be achieved in the presence of neutron-induced background from the scattering and collimating devices. Uniform range shifts were detected with a standard deviation of 0.1 mm to 0.2 mm at a dose level of 30 cGy to 50 cGy (RBE). The detectable magnitude of a range shift limited to a part of the treatment field area was approximately proportional to the ratio between the field area and the area affected by the range shift. We conclude that it is feasible to detect changes in the range of passively scattered proton beams using a relatively simple prompt gamma-ray detection system. The method can be employed for in vivo verification of the consistency of the delivered range in fractionated treatments.

Metal Artifact Reduction in CT: Where Are We After Four Decades?

Methods to overcome metal artifacts in computed tomography (CT) images have been researched and developed for nearly 40 years. When X-rays pass through a metal object, depending on its size and density, different physical effects will negatively affect the measurements, most notably beam hardening, scatter, noise, and the non-linear partial volume effect. These phenomena severely degrade image quality and hinder the diagnostic power and treatment outcomes in many clinical applications. In this paper, we first review the fundamental causes of metal artifacts, categorize metal object types, and present recent trends in the CT metal artifact reduction (MAR) literature. To improve image quality and recover information about underlying structures, many methods and correction algorithms have been proposed and tested. We comprehensively review and categorize these methods into six different classes of MAR: metal implant optimization, improvements to the data acquisition process, data correction based on physics models, modifications to the reconstruction algorithm (projection completion and iterative reconstruction), and image-based post-processing. The primary goals of this paper are to identify the strengths and limitations of individual MAR methods and overall classes, and establish a relationship between types of metal objects and the classes that most effectively overcome their artifacts. The main challenges for the field of MAR continue to be cases with large, dense metal implants, as well as cases with multiple metal objects in the field of view. Severe photon starvation is difficult to compensate for with only software corrections. Hence, the future of MAR seems to be headed toward a combined approach of improving the acquisition process with dual-energy CT, higher energy X-rays, or photon-counting detectors, along with advanced reconstruction approaches. Additional outlooks are addressed, including the need for a standardized evaluation system to compare MAR methods.