Feasibility of the J-PET to monitor the range of therapeutic proton beams.

PURPOSE: The aim of this work is to investigate the feasibility of the Jagiellonian Positron Emission Tomography (J-PET) scanner for intra-treatment proton beam range monitoring. METHODS: The Monte Carlo simulation studies with GATE and PET image reconstruction with CASToR were performed in order to compare six J-PET scanner geometries. We simulated proton irradiation of a PMMA phantom with a Single Pencil Beam (SPB) and Spread-Out Bragg Peak (SOBP) of various ranges. The sensitivity and precision of each scanner were calculated, and considering the setup's cost-effectiveness, we indicated potentially optimal geometries for the J-PET scanner prototype dedicated to the proton beam range assessment. RESULTS: The investigations indicate that the double-layer cylindrical and triple-layer double-head configurations are the most promising for clinical application. We found that the scanner sensitivity is of the order of 10-5 coincidences per primary proton, while the precision of the range assessment for both SPB and SOBP irradiation plans was found below 1 mm. Among the scanners with the same number of detector modules, the best results are found for the triple-layer dual-head geometry. The results indicate that the double-layer cylindrical and triple-layer double-head configurations are the most promising for the clinical application, CONCLUSIONS:: We performed simulation studies demonstrating that the feasibility of the J-PET detector for PET-based proton beam therapy range monitoring is possible with reasonable sensitivity and precision enabling its pre-clinical tests in the clinical proton therapy environment. Considering the sensitivity, precision and cost-effectiveness, the double-layer cylindrical and triple-layer dual-head J-PET geometry configurations seem promising for future clinical application.

Detection of range shifts in proton beam therapy using the J-PET scanner: a patient simulation study.

OBJECTIVE: The Jagiellonian PET (J-PET) technology, based on plastic scintillators, has been proposed as a cost effective tool for detecting range deviations during proton therapy. This study investigates the feasibility of using J-PET for range monitoring by means of a detailed Monte Carlo simulation study of 95 patients who underwent proton therapy at the Cyclotron Centre Bronowice (CCB) in Krakow, Poland. Approach: Discrepancies between prescribed and delivered treatments were artificially introduced in the simulations by means of shifts in patient positioning and in the Hounsfield unit to the relative proton stopping power calibration curve. A dual-layer, cylindrical J-PET geometry was simulated in an in-room monitoring scenario and a triple-layer, dual-head geometry in an in-beam protocol. The distribution of range shifts in reconstructed PET activity was visualised in the beam's eye view. Linear prediction models were constructed from all patients in the cohort, using the mean shift in reconstructed PET activity as a predictor of the mean proton range deviation. Main results: Maps of deviations in the range of reconstructed PET distributions showed agreement with those of deviations in dose range in most patients. The linear prediction model showed a good fit, with coefficient of determination r^2 = 0.84 (in-room) and 0.75 (in-beam). Residual standard error was below 1 mm: 0.33 mm (in-room) and 0.23 mm (in-beam). Significance: The precision of the proposed prediction models shows the sensitivity of the proposed J-PET scanners to shifts in proton range for a wide range of clinical treatment plans. Furthermore, it motivates the use of such models as a tool for predicting proton range deviations and opens up new prospects for investigations into the use of intra-treatment PET images for predicting clinical metrics that aid in the assessment of the quality of delivered treatment. .

GPU accelerated Monte Carlo scoring of positron emitting isotopes produced during proton therapy for PET verification.

Objective.Verification of delivered proton therapy treatments is essential for reaping the many benefits of the modality, with the most widely proposedin vivoverification technique being the imaging of positron emitting isotopes generated in the patient during treatment using positron emission tomography (PET). The purpose of this work is to reduce the computational resources and time required for simulation of patient activation during proton therapy using the GPU accelerated Monte Carlo code FRED, and to validate the predicted activity against the widely used Monte Carlo code GATE.Approach.We implement a continuous scoring approach for the production of positron emitting isotopes within FRED version 5.59.9. We simulate treatment plans delivered to 95 head and neck patients at Centrum Cyklotronowe Bronowice using this GPU implementation, and verify the accuracy using the Monte Carlo toolkit GATE version 9.0.Main results.We report an average reduction in computational time by a factor of 50 when using a local system with 2 GPUs as opposed to a large compute cluster utilising between 200 to 700 CPU threads, enabling simulation of patient activity within an average of 2.9 min as opposed to 146 min. All simulated plans are in good agreement across the two Monte Carlo codes. The two codes agree within a maximum of 0.95σon a voxel-by-voxel basis for the prediction of 7 different isotopes across 472 simulated fields delivered to 95 patients, with the average deviation over all fields being 6.4 × 10-3σ.Significance.The implementation of activation calculations in the GPU accelerated Monte Carlo code FRED provides fast and reliable simulation of patient activation following proton therapy, allowing for research and development of clinical applications of range verification for this treatment modality using PET to proceed at a rapid pace.

ProTheRaMon - a GATE simulation framework for proton therapy range monitoring using PET imaging.

OBJECTIVE: This paper reports on the implementation and shows examples of the use of the ProTheRaMon framework for simulating the delivery of proton therapy treatment plans and range monitoring using positron emission tomography (PET). ProTheRaMon offers complete processing of proton therapy treatment plans, patient CT geometries, and intra-treatment PET imaging, taking into account therapy and imaging coordinate systems and activity decay during the PET imaging protocol specific to a given proton therapy facility. We present the ProTheRaMon framework and illustrate its potential use case and data processing steps for a patient treated at the Cyclotron Centre Bronowice (CCB) proton therapy center in Krakow, Poland. APPROACH: The ProTheRaMon framework is based on GATE Monte Carlo software, the CASToR reconstruction package and in-house developed Python and bash scripts. The framework consists of five separated simulation and data processing steps, that can be further optimized according to the user's needs and specific settings of a given proton therapy facility and PET scanner design. MAIN RESULTS: ProTheRaMon is presented using example data from a patient treated at CCB and the J-PET scanner to demonstrate the application of the framework for proton therapy range monitoring. The output of each simulation and data processing stage is described and visualized. SIGNIFICANCE: We demonstrate that the ProTheRaMon simulation platform is a high-performance tool, capable of running on a computational cluster and suitable for multi-parameter studies, with databases consisting of large number of patients, as well as different PET scanner geometries and settings for range monitoring in a clinical environment. Due to its modular structure, the ProTheRaMon framework can be adjusted for different proton therapy centers and/or different PET detector geometries. It is available to the community via github.

Photogrammetry-Based Volume Measurement Framework for the Particle Density Estimation of LECA.

This paper presents a photogrammetry-based volume measurement framework for the particle density estimation of Lightweight expanded clay aggregate (LECA). The results are compared with computed tomography (CT) and Archimedes' method measurements. All of the steps required in order to apply the proposed approach are explained. Next, we discuss how the interpretation of open pores affects the results of volume measurements. We propose to process the shapes obtained from different methods by applying an Ambient Occlusion algorithm with the same threshold, t = 0.175. The difference between the CT and SfM methods is less than 0.006 g/cm3, proving that the photogrammetry-based approach is accurate enough. The Archimedes' method significantly overestimates the density of the particles. Nevertheless, its accuracy is acceptable for most engineering purposes. Additionally, we evaluate the accuracy of shape reconstruction (in terms of the Hausdorff distance). For 95% of the grain's surface, the maximum error is between 0.073 mm and 0.129 mm (depending on the grain shape). The presented approach is helpful for measuring the particle density of porous aggregates. The proposed methodology can be utilized in order to estimate intergranular porosity, which is valuable information for the calibration of DEM models.

Evaluation of a high resolution silicon PET insert module.

Conventional PET systems can be augmented with additional detectors placed in close proximity of the region of interest. We developed a high resolution PET insert module to evaluate the added benefit of such a combination. The insert module consists of two back-to-back 1 mm thick silicon sensors, each segmented into 1040 1 mm2 pads arranged in a 40 by 26 array. A set of 16 VATAGP7.1 ASICs and a custom assembled data acquisition board were used to read out the signal from the insert module. Data were acquired in slice (2D) geometry with a Jaszczak phantom (rod diameters of 1.2-4.8 mm) filled with 18F-FDG and the images were reconstructed with ML-EM method. Both data with full and limited angular coverage from the insert module were considered and three types of coincidence events were combined. The ratio of high-resolution data that substantially improves quality of the reconstructed image for the region near the surface of the insert module was estimated to be about 4%. Results from our previous studies suggest that such ratio could be achieved at a moderate technological expense by using an equivalent of two insert modules (an effective sensor thickness of 4 mm).

A high-resolution PET demonstrator using a silicon "magnifying glass".

To assist ongoing investigations of the limits of the tradeoff between spatial resolution and noise in PET imaging, several PET instruments based on silicon-pad detectors have been developed. The latest is a segment of a dual-ring device to demonstrate that excellent reconstructed image resolution can be achieved with a scanner that uses high-resolution detectors placed close to the object of interest or surrounding a small field-of-view in combination with detectors having modest resolution at larger radius. The outer ring of our demonstrator comprises conventional BGO block detectors scavenged from a clinical PET scanner and located at a 500mm radius around a 50mm diameter field-of-view. The inner detector-in contrast to the high-Z scintillator typically used in PET-is based on silicon-pad detectors located at 70mm nominal radius. Each silicon detector has 512 1.4mm x 1.4mm x 1mm detector elements in a 16 x 32 array and is read out using VATA GP7 ASICs (Gamma Medica-Ideas, Northridge, CA). Even though virtually all interactions of 511 keV annihilation photons in silicon are Compton-scatter, both high spatial resolution and reasonable sensitivity appears possible. The system has demonstrated resolution of ~0.7mm FWHM with Na-22 for coincidences having the highest intrinsic resolution (silicon-silicon) and 5-6mm FWHM for the lowest resolution BGO-BGO coincidences. Spatial resolution for images reconstructed from the mixed silicon-BGO coincidences is ~1.5mm FWHM demonstrating the "magnifying-glass" concept.