Experimental validation of a multi-energy x-ray adapted scatter separation method.

Both in radiography and computed tomography (CT), recently emerged energy-resolved x-ray photon counting detectors enable the identification and quantification of individual materials comprising the inspected object. However, the approaches used for these operations require highly accurate x-ray images. The accuracy of the images is severely compromised by the presence of scattered radiation, which leads to a loss of spatial contrast and, more importantly, a bias in radiographic material imaging and artefacts in CT. The aim of the present study was to experimentally evaluate a recently introduced partial attenuation spectral scatter separation approach (PASSSA) adapted for multi-energy imaging. For this purpose, a prototype x-ray system was used. Several radiographic acquisitions of an anthropomorphic thorax phantom were performed. Reference primary images were obtained via the beam-stop (BS) approach. The attenuation images acquired from PASSSA-corrected data showed a substantial increase in local contrast and internal structure contour visibility when compared to uncorrected images. A substantial reduction of scatter induced bias was also achieved. Quantitatively, the developed method proved to be in relatively good agreement with the BS data. The application of the proposed scatter correction technique lowered the initial normalized root-mean-square error (NRMSE) of 45% between the uncorrected total and the reference primary spectral images by a factor of 9, thus reducing it to around 5%.

A novel scatter separation method for multi-energy x-ray imaging.

X-ray imaging coupled with recently emerged energy-resolved photon counting detectors provides the ability to differentiate material components and to estimate their respective thicknesses. However, such techniques require highly accurate images. The presence of scattered radiation leads to a loss of spatial contrast and, more importantly, a bias in radiographic material imaging and artefacts in computed tomography (CT). The aim of the present study was to introduce and evaluate a partial attenuation spectral scatter separation approach (PASSSA) adapted for multi-energy imaging. This evaluation was carried out with the aid of numerical simulations provided by an internal simulation tool, Sindbad-SFFD. A simplified numerical thorax phantom placed in a CT geometry was used. The attenuation images and CT slices obtained from corrected data showed a remarkable increase in local contrast and internal structure detectability when compared to uncorrected images. Scatter induced bias was also substantially decreased. In terms of quantitative performance, the developed approach proved to be quite accurate as well. The average normalized root-mean-square error between the uncorrected projections and the reference primary projections was around 23%. The application of PASSSA reduced this error to around 5%. Finally, in terms of voxel value accuracy, an increase by a factor >10 was observed for most inspected volumes-of-interest, when comparing the corrected and uncorrected total volumes.

An empirical model for calculation of the collimator contamination dose in therapeutic proton beams.

Collimators are used as lateral beam shaping devices in proton therapy with passive scattering beam lines. The dose contamination due to collimator scattering can be as high as 10% of the maximum dose and influences calculation of the output factor or monitor units (MU). To date, commercial treatment planning systems generally use a zero-thickness collimator approximation ignoring edge scattering in the aperture collimator and few analytical models have been proposed to take scattering effects into account, mainly limited to the inner collimator face component. The aim of this study was to characterize and model aperture contamination by means of a fast and accurate analytical model. The entrance face collimator scatter distribution was modeled as a 3D secondary dose source. Predicted dose contaminations were compared to measurements and Monte Carlo simulations. Measurements were performed on two different proton beam lines (a fixed horizontal beam line and a gantry beam line) with divergent apertures and for several field sizes and energies. Discrepancies between analytical algorithm dose prediction and measurements were decreased from 10% to 2% using the proposed model. Gamma-index (2%/1 mm) was respected for more than 90% of pixels. The proposed analytical algorithm increases the accuracy of analytical dose calculations with reasonable computation times.

Monte Carlo simulation of prompt γ-ray emission in proton therapy using a specific track length estimator.

A Monte Carlo (MC) variance reduction technique is developed for prompt-γ emitters calculations in proton therapy. Prompt-γ emitted through nuclear fragmentation reactions and exiting the patient during proton therapy could play an important role to help monitoring the treatment. However, the estimation of the number and the energy of emitted prompt-γ per primary proton with MC simulations is a slow process. In order to estimate the local distribution of prompt-γ emission in a volume of interest for a given proton beam of the treatment plan, a MC variance reduction technique based on a specific track length estimator (TLE) has been developed. First an elemental database of prompt-γ emission spectra is established in the clinical energy range of incident protons for all elements in the composition of human tissues. This database of the prompt-γ spectra is built offline with high statistics. Regarding the implementation of the prompt-γ TLE MC tally, each proton deposits along its track the expectation of the prompt-γ spectra from the database according to the proton kinetic energy and the local material composition. A detailed statistical study shows that the relative efficiency mainly depends on the geometrical distribution of the track length. Benchmarking of the proposed prompt-γ TLE MC technique with respect to an analogous MC technique is carried out. A large relative efficiency gain is reported, ca. 10(5).

Technical Note: Experimental carbon ion range verification in inhomogeneous phantoms using prompt gammas.

PURPOSE: The purpose of this study was to experimentally assess the possibility to monitor carbon ion range variations--due to tumor shift and/or elongation or shrinking--using prompt-gamma (PG) emission with inhomogeneous phantoms. Such a study is related to the development of PG monitoring techniques to be used in a carbon ion therapy context. METHODS: A 95 MeV/u carbon ion beam was used to irradiate phantoms with a variable density along the ion path to mimic the presence of bone and lung in homogeneous humanlike tissue. PG profiles were obtained after a longitudinal scan of the phantoms. A setup comprising a narrow single-slit collimator and two detectors placed at 90° with respect to the beam axis was used. The time of flight technique was applied to allow the selection between PG and background events. RESULTS: Using the positions at 50% entrance and 50% falloff of the PG profiles, a quantity called prompt-gamma profile length (PGPL) is defined. It is possible to observe shifts in the PGPL when there are absolute ion range shifts as small as 1-2 mm. Quantitatively, for an ion range shift of -1.33 ± 0.46 mm (insertion of a Teflon slab), a PGPL difference of -1.93 ± 0.58 mm and -1.84 ± 1.27 mm is obtained using a BaF2 and a NaI(Tl) detector, respectively. In turn, when an ion range shift of 4.59 ± 0.42 mm (insertion of a lung-equivalent material slab) is considered, the difference is of 4.10 ± 0.54 and 4.39 ± 0.80 mm for the same detectors. CONCLUSIONS: Herein, experimental evidence of the usefulness of employing PG to monitor carbon ion range using inhomogeneous phantoms is presented. Considering the homogeneous phantom as reference, the results show that the information provided by the PG emission allows for detecting ion range shifts as small as 1-2 mm. When considering the expected PG emission from an energy slice in a carbon ion therapy scenario, the experimental setup would allow to retrieve the same PGPL as the high statistics of the full experimental dataset in 58% of the times. However, this success rate increases to 93% when using a better optimized setup by means of Monte Carlo simulations.

Absolute prompt-gamma yield measurements for ion beam therapy monitoring.

Prompt-gamma emission detection is a promising technique for hadrontherapy monitoring purposes. In this regard, obtaining prompt-gamma yields that can be used to develop monitoring systems based on this principle is of utmost importance since any camera design must cope with the available signal. Herein, a comprehensive study of the data from ten single-slit experiments is presented, five consisting in the irradiation of either PMMA or water targets with lower and higher energy carbon ions, and another five experiments using PMMA targets and proton beams. Analysis techniques such as background subtraction methods, geometrical normalization, and systematic uncertainty estimation were applied to the data in order to obtain absolute prompt-gamma yields in units of prompt-gamma counts per incident ion, unit of field of view, and unit of solid angle. At the entrance of a PMMA target, where the contribution of secondary nuclear reactions is negligible, prompt-gamma counts per incident ion, per millimetre and per steradian equal to (124 ± 0.7stat ± 30sys) × 10(-6) for 95 MeV u(-1) carbon ions, (79 ± 2stat ± 23sys) × 10(-6) for 310 MeV u(-1) carbon ions, and (16 ± 0.07stat ± 1sys) × 10(-6) for 160 MeV protons were found for prompt gammas with energies higher than 1 MeV. This shows a factor 5 between the yields of two different ions species with the same range in water (160 MeV protons and 310 MeV u(-1) carbon ions). The target composition was also found to influence the prompt-gamma yield since, for 300/310 MeV u(-1) carbon ions, a 42% greater yield ((112 ± 1stat ± 22sys) × 10(-6) counts ion(-1) mm(-1) sr(-1)) was obtained with a water target compared to a PMMA one.

A track length estimator method for dose calculations in low-energy X-ray irradiations: implementation, properties and performance.

The track length estimator (TLE) method, an "on-the-fly" fluence tally in Monte Carlo (MC) simulations, recently implemented in GATE 6.2, is known as a powerful tool to accelerate dose calculations in the domain of low-energy X-ray irradiations using the kerma approximation. Overall efficiency gains of the TLE with respect to analogous MC were reported in the literature for regions of interest in various applications (photon beam radiation therapy, X-ray imaging). The behaviour of the TLE method in terms of statistical properties, dose deposition patterns, and computational efficiency compared to analogous MC simulations was investigated. The statistical properties of the dose deposition were first assessed. Derivations of the variance reduction factor of TLE versus analogous MC were carried out, starting from the expression of the dose estimate variance in the TLE and analogous MC schemes. Two test cases were chosen to benchmark the TLE performance in comparison with analogous MC: (i) a small animal irradiation under stereotactic synchrotron radiation therapy conditions and (ii) the irradiation of a human pelvis during a cone beam computed tomography acquisition. Dose distribution patterns and efficiency gain maps were analysed. The efficiency gain exhibits strong variations within a given irradiation case, depending on the geometrical (voxel size, ballistics) and physical (material and beam properties) parameters on the voxel scale. Typical values lie between 10 and 10(3), with lower levels in dense regions (bone) outside the irradiated channels (scattered dose only), and higher levels in soft tissues directly exposed to the beams.

Split exponential track length estimator for Monte-Carlo simulations of small-animal radiation therapy.

We propose the split exponential track length estimator (seTLE), a new kerma-based method combining the exponential variant of the TLE and a splitting strategy to speed up Monte Carlo (MC) dose computation for low energy photon beams. The splitting strategy is applied to both the primary and the secondary emitted photons, triggered by either the MC events generator for primaries or the photon interactions generator for secondaries. Split photons are replaced by virtual particles for fast dose calculation using the exponential TLE. Virtual particles are propagated by ray-tracing in voxelized volumes and by conventional MC navigation elsewhere. Hence, the contribution of volumes such as collimators, treatment couch and holding devices can be taken into account in the dose calculation.We evaluated and analysed the seTLE method for two realistic small animal radiotherapy treatment plans. The effect of the kerma approximation, i.e. the complete deactivation of electron transport, was investigated. The efficiency of seTLE against splitting multiplicities was also studied. A benchmark with analog MC and TLE was carried out in terms of dose convergence and efficiency.The results showed that the deactivation of electrons impacts the dose at the water/bone interface in high dose regions. The maximum and mean dose differences normalized to the dose at the isocenter were, respectively of 14% and 2% . Optimal splitting multiplicities were found to be around 300. In all situations, discrepancies in integral dose were below 0.5% and 99.8% of the voxels fulfilled a 1%/0.3 mm gamma index criterion. Efficiency gains of seTLE varied from 3.2 × 10(5) to 7.7 × 10(5) compared to analog MC and from 13 to 15 compared to conventional TLE.In conclusion, seTLE provides results similar to the TLE while increasing the efficiency by a factor between 13 and 15, which makes it particularly well-suited to typical small animal radiation therapy applications.

Design optimisation of a TOF-based collimated camera prototype for online hadrontherapy monitoring.

Hadrontherapy is an innovative radiation therapy modality for which one of the main key advantages is the target conformality allowed by the physical properties of ion species. However, in order to maximise the exploitation of its potentialities, online monitoring is required in order to assert the treatment quality, namely monitoring devices relying on the detection of secondary radiations. Herein is presented a method based on Monte Carlo simulations to optimise a multi-slit collimated camera employing time-of-flight selection of prompt-gamma rays to be used in a clinical scenario. In addition, an analytical tool is developed based on the Monte Carlo data to predict the expected precision for a given geometrical configuration. Such a method follows the clinical workflow requirements to simultaneously have a solution that is relatively accurate and fast. Two different camera designs are proposed, considering different endpoints based on the trade-off between camera detection efficiency and spatial resolution to be used in a proton therapy treatment with active dose delivery and assuming a homogeneous target.

Assessment and improvements of Geant4 hadronic models in the context of prompt-gamma hadrontherapy monitoring.

Monte Carlo simulations are nowadays essential tools for a wide range of research topics in the field of radiotherapy. They also play an important role in the effort to develop a real-time monitoring system for quality assurance in proton and carbon ion therapy, by means of prompt-gamma detection. The internal theoretical nuclear models of Monte Carlo simulation toolkits are of decisive importance for the accurate description of neutral or charged particle emission, produced by nuclear interactions between beam particles and target nuclei. We assess the performance of Geant4 nuclear models in the context of prompt-gamma emission, comparing them with experimental data from proton and carbon ion beams. As has been shown in the past and further indicated in our study, the prompt-gamma yields are consistently overestimated by Geant4 by a factor of about 100% to 200% over an energy range from 80 to 310 MeV/u for the case of (12)C, and to a lesser extent for 160 MeV protons. Furthermore, we focus on the quantum molecular dynamics (QMD) modeling of ion-ion collisions, in order to optimize its description of light nuclei, which are abundant in the human body and mainly anticipated in hadrontherapy applications. The optimization has been performed by benchmarking QMD free parameters with well established nuclear properties. In addition, we study the effect of this optimization on charged particle emission. With the usage of the proposed parameter values, discrepancies reduce to less than 70%, with the highest values being attributed to the nucleon-ion induced prompt-gammas. This conclusion, also confirmed by the disagreement we observe in the case of proton beams, indicates the need for further investigation on nuclear models which describe proton and neutron induced nuclear reactions.

Real-time proton beam range monitoring by means of prompt-gamma detection with a collimated camera.

Prompt-gamma profile was measured at WPE-Essen using 160 MeV protons impinging a movable PMMA target. A single collimated detector was used with time-of-flight (TOF) to reduce the background due to neutrons. The target entrance rise and the Bragg peak falloff retrieval precision was determined as a function of incident proton number by a fitting procedure using independent data sets. Assuming improved sensitivity of this camera design by using a greater number of detectors, retrieval precisions of 1 to 2 mm (rms) are expected for a clinical pencil beam. TOF improves the contrast-to-noise ratio and the performance of the method significantly.

Machine learning-based patient specific prompt-gamma dose monitoring in proton therapy.

Online dose monitoring in proton therapy is currently being investigated with prompt-gamma (PG) devices. PG emission was shown to be correlated with dose deposition. This relationship is mostly unknown under real conditions. We propose a machine learning approach based on simulations to create optimized treatment-specific classifiers that detect discrepancies between planned and delivered dose. Simulations were performed with the Monte-Carlo platform Gate/Geant4 for a spot-scanning proton therapy treatment and a PG camera prototype currently under investigation. The method first builds a learning set of perturbed situations corresponding to a range of patient translation. This set is then used to train a combined classifier using distal falloff and registered correlation measures. Classifier performances were evaluated using receiver operating characteristic curves and maximum associated specificity and sensitivity. A leave-one-out study showed that it is possible to detect discrepancies of 5 mm with specificity and sensitivity of 85% whereas using only distal falloff decreases the sensitivity down to 77% on the same data set. The proposed method could help to evaluate performance and to optimize the design of PG monitoring devices. It is generic: other learning sets of deviations, other measures and other types of classifiers could be studied to potentially reach better performance. At the moment, the main limitation lies in the computation time needed to perform the simulations.

Interaction vertex imaging (IVI) for carbon ion therapy monitoring: a feasibility study.

Proton imaging can be seen as a powerful technique for online monitoring of ion range during carbon ion therapy irradiations. Indeed, a large number of secondary protons are created during nuclear reactions, and many of these protons are likely to escape from the patient even for deep-seated tumors, carrying accurate information on the reaction vertex position. Two detection techniques have been considered: (i) double-proton detection by means of two forward-located trackers and (ii) single-proton detection in coincidence with the incoming carbon ion detected by means of a beam hodoscope. Geant4 simulations, validated by proton yield measurements performed at GANIL and GSI, show that ion-range monitoring is accessible on a pencil-beam basis with the single-proton imaging technique. Millimetric precision on the Bragg peak position is expected in the ideal case of homogeneous targets. The uncertainties in more realistic conditions should be investigated, in particular the influence of tissue heterogeneity in the very last part of the ion path (about 20 mm).

GATE as a GEANT4-based Monte Carlo platform for the evaluation of proton pencil beam scanning treatment plans.

Active scanning delivery systems take full advantage of ion beams to best conform to the tumor and to spare surrounding healthy tissues; however, it is also a challenging technique for quality assurance. In this perspective, we upgraded the GATE/GEANT4 Monte Carlo platform in order to recalculate the treatment planning system (TPS) dose distributions for active scanning systems. A method that allows evaluating the TPS dose distributions with the GATE Monte Carlo platform has been developed and applied to the XiO TPS (Elekta), for the IBA proton pencil beam scanning (PBS) system. First, we evaluated the specificities of each dose engine. A dose-conversion scheme that allows one to convert dose to medium into dose to water was implemented within GATE. Specific test cases in homogeneous and heterogeneous configurations allowed for the estimation of the differences between the beam models implemented in XiO and GATE. Finally, dose distributions of a prostate treatment plan were compared. In homogeneous media, a satisfactory agreement was generally obtained between XiO and GATE. The maximum stopping power difference of 3% occurred in a human tissue of 0.9 g cm(-3) density and led to a significant range shift. Comparisons in heterogeneous configurations pointed out the limits of the TPS dose calculation accuracy and the superiority of Monte Carlo simulations. The necessity of computing dose to water in our Monte Carlo code for comparisons with TPSs is also presented. Finally, the new capabilities of the platform are applied to a prostate treatment plan and dose differences between both dose engines are analyzed in detail. This work presents a generic method to compare TPS dose distributions with the GATE Monte Carlo platform. It is noteworthy that GATE is also a convenient tool for imaging applications, therefore opening new research possibilities for the PBS modality.

A Monte Carlo pencil beam scanning model for proton treatment plan simulation using GATE/GEANT4.

This work proposes a generic method for modeling scanned ion beam delivery systems, without simulation of the treatment nozzle and based exclusively on beam data library (BDL) measurements required for treatment planning systems (TPS). To this aim, new tools dedicated to treatment plan simulation were implemented in the Gate Monte Carlo platform. The method was applied to a dedicated nozzle from IBA for proton pencil beam scanning delivery. Optical and energy parameters of the system were modeled using a set of proton depth-dose profiles and spot sizes measured at 27 therapeutic energies. For further validation of the beam model, specific 2D and 3D plans were produced and then measured with appropriate dosimetric tools. Dose contributions from secondary particles produced by nuclear interactions were also investigated using field size factor experiments. Pristine Bragg peaks were reproduced with 0.7 mm range and 0.2 mm spot size accuracy. A 32 cm range spread-out Bragg peak with 10 cm modulation was reproduced with 0.8 mm range accuracy and a maximum point-to-point dose difference of less than 2%. A 2D test pattern consisting of a combination of homogeneous and high-gradient dose regions passed a 2%/2 mm gamma index comparison for 97% of the points. In conclusion, the generic modeling method proposed for scanned ion beam delivery systems was applicable to an IBA proton therapy system. The key advantage of the method is that it only requires BDL measurements of the system. The validation tests performed so far demonstrated that the beam model achieves clinical performance, paving the way for further studies toward TPS benchmarking. The method involves new sources that are available in the new Gate release V6.1 and could be further applied to other particle therapy systems delivering protons or other types of ions like carbon.

Real-time monitoring of the Bragg-peak position in ion therapy by means of single photon detection.

For real-time monitoring of the longitudinal position of the Bragg-peak during an ion therapy treatment, a novel non-invasive technique has been recently proposed that exploits the detection of prompt gamma-rays issued from nuclear fragmentation. Two series of experiments have been performed at the GANIL and GSI facilities with 95 and 305 MeV/u (12)C(6+) ion beams stopped in PMMA and water phantoms. In both experiments, a clear correlation was obtained between the carbon ion range and the prompt photon profile. Additionally, an extensive study has been performed to investigate whether a prompt neutron component may be correlated with the carbon ion range. No such correlation was found. The present paper demonstrates that a collimated set-up can be used to detect single photons by means of time-of-flight measurements, at those high energies typical for ion therapy. Moreover, the applicability of the technique both at cyclotron and at synchrotron facilities is shown. It is concluded that the detected photon count rates provide sufficiently high statistics to allow real-time control of the longitudinal position of the Bragg-peak under clinical conditions.

Simulation of dose deposition in stereotactic synchrotron radiation therapy: a fast approach combining Monte Carlo and deterministic algorithms.

A hybrid approach, combining deterministic and Monte Carlo (MC) calculations, is proposed to compute the distribution of dose deposited during stereotactic synchrotron radiation therapy treatment. The proposed approach divides the computation into two parts: (i) the dose deposited by primary radiation (coming directly from the incident x-ray beam) is calculated in a deterministic way using ray casting techniques and energy-absorption coefficient tables and (ii) the dose deposited by secondary radiation (Rayleigh and Compton scattering, fluorescence) is computed using a hybrid algorithm combining MC and deterministic calculations. In the MC part, a small number of particle histories are simulated. Every time a scattering or fluorescence event takes place, a splitting mechanism is applied, so that multiple secondary photons are generated with a reduced weight. The secondary events are further processed in a deterministic way, using ray casting techniques. The whole simulation, carried out within the framework of the Monte Carlo code Geant4, is shown to converge towards the same results as the full MC simulation. The speed of convergence is found to depend notably on the splitting multiplicity, which can easily be optimized. To assess the performance of the proposed algorithm, we compare it to state-of-the-art MC simulations, accelerated by the track length estimator technique (TLE), considering a clinically realistic test case. It is found that the hybrid approach is significantly faster than the MC/TLE method. The gain in speed in a test case was about 25 for a constant precision. Therefore, this method appears to be suitable for treatment planning applications.

[Ascaris lumbricoides--a cause of acute abdomen].

Ascaris lumbricoides infestation is the most common helminthic disease in many developing countries and may cause severe surgical complications, especially in children. We present a 3-year-old Ethiopian immigrant brought to hospital directly from the airport because of signs of acute abdomen. It was found to be caused by volvulus and necrosis of a loop of bowel impacted with worms. The necrotic bowel loop was resected and a temporary ileostomy was formed; recovery was uneventful. This disease and its severe complications are rare in Israel and the western world. Awareness of its occurrence is important in view of the present world-wide large-scale immigration from developing countries.

Subtotal gastrectomy in infancy for perforating necrotizing gastritis.

Necrotizing gastritis is a rare pathology causing a high rate of morbidity and mortality in the infants. A 4-week-old baby, with right hypoplastic kidney and ectopic ureter, was admitted because of profound septic shock and "coffee ground" vomitus. Aggressive treatment was started and hemodynamic stabilization was achieved. On the fourth admission day, ascites was noted and on the eighth day in a plain abdominal x-ray, free air was shown. At urgent explorative laparotomy, double posterior gastric wall perforations with extensive gastric wall necrosis were found, which required subtotal gastrectomy. The etiology and pathophysiology of this rare process are discussed with an emphasis on the difficulty in diagnosis of posterior gastric perforation into the lesser sac.

Persistent müllerian duct remnant in a newborn.

Müllerian duct remnant in a newborn is a rare finding. The presented case illustrates well the ultrasonic and radiologic features of the entity.