Positron emission tomography is one of the most mature techniques for monitoring the particles range in hadron therapy, aiming to reduce treatment uncertainties and therefore the extent of safety ...margins in the treatment plan. In-beam PET monitoring has been already performed using inter-spill and post-irradiation data, i.e. while the particle beam is off or paused. The full beam acquisition procedure is commonly discarded because the particle spills abruptly increase the random coincidence rates and therefore the image noise. This is because random coincidences cannot be separated by annihilation photons originating from radioactive decays and cannot be corrected with standard random coincidence techniques due to the time correlation of the beam-induced background with the ion beam microstructure. The aim of this paper is to provide a new method to recover in-spill data to improve the images obtained with full-beam PET acquisitions. This is done by estimating the temporal microstructure of the beam and thus selecting input PET events that are less likely to be random ones. The PET detector we used was the one developed within the INSIDE project and tested at the CNAO synchrotron-based facility. The data were taken on a PMMA phantom irradiated with 72 MeV proton pencil beams. The obtained results confirm the possibility of improving the acquired PET data without any external signal coming from the synchrotron or ad hoc detectors.
During particle therapy irradiation, positron emitters with half-lives ranging from 2 to 20 min are generated from nuclear processes. The half-lives are such that it is possible either to detect the ...positron signal in the treatment room using an in-beam positron emission tomography (PET) system, right after the irradiation, or to quickly transfer the patient to a close PET/CT scanner. Since the activity distribution is spatially correlated with the dose, it is possible to use PET imaging as an indirect method to assure the quality of the dose delivery. In this work, we present a new dedicated PET system able to operate in-beam. The PET apparatus consists in two 10 cm × 10 cm detector heads. Each detector is composed of four scintillating matrices of 23 × 23 LYSO crystals. The crystal size is 1.9 mm × 1.9 mm × 16 mm. Each scintillation matrix is read out independently with a modularized acquisition system. The distance between the two opposing detector heads was set to 20 cm. The system has very low dead time per detector area and a 3 ns coincidence window, which is capable to sustain high single count rates and to keep the random counts relatively low. This allows a new full-beam monitoring modality that includes data acquisition also while the beam is on. The PET system was tested during the irradiation at the CATANA (INFN, Catania, Italy) cyclotron-based proton therapy facility. Four acquisitions with different doses and dose rates were analysed. In all cases the random to total coincidences ratio was equal or less than 25%. For each measurement we estimated the accuracy and precision of the activity range on a set of voxel lines within an irradiated PMMA phantom. Results show that the inclusion of data acquired during the irradiation, referred to as beam-on data, improves both the precision and accuracy of the range measurement with respect to data acquired only after irradiation. Beam-on data alone are enough to give precisions better than 1 mm when at least 5 Gy are delivered.
Treatment quality assessment is a crucial feature for both present and next-generation ion therapy facilities. Several approaches are being explored, based on prompt radiation emission or on PET ...signals by Formula: see text-decaying isotopes generated by beam interactions with the body. In-beam PET monitoring at synchrotron-based ion therapy facilities has already been performed, either based on inter-spill data only, to avoid the influence of the prompt radiation, or including both in-spill and inter-spill data. However, the PET images either suffer of poor statistics (inter-spill) or are more influenced by the background induced by prompt radiation (in-spill). Both those problems are expected to worsen for accelerators with improved duty cycle where the inter-spill interval is reduced to shorten the treatment time. With the aim of assessing the detector performance and developing techniques for background reduction, a test of an in-beam PET detector prototype was performed at the CNAO synchrotron-based ion therapy facility in full-beam acquisition modality. Data taken with proton beams impinging on PMMA phantoms showed the system acquisition capability and the resulting activity distribution, separately reconstructed for the in-spill and the inter-spill data. The coincidence time resolution for in-spill and inter-spill data shows a good agreement, with a slight deterioration during the spill. The data selection technique allows the identification and rejection of most of the background originated during the beam delivery. The activity range difference between two different proton beam energies (68 and 72 MeV) was measured and found to be in sub-millimeter agreement with the expected result. However, a slightly longer (2 mm) absolute profile length is obtained for in-spill data when compared to inter-spill data.
In this study a procedure for range verification in proton therapy by means of a planar in-beam PET system is presented. The procedure consists of two steps: the measurement of the β+-activity ...induced in the irradiated body by the proton beam and the comparison of these distributions with simulations. The experimental data taking was performed at the CNAO center in Pavia, Italy, irradiating plastic phantoms. For two different cases we demonstrate how a real-time feedback of the delivered treatment plan can be obtained with in-beam PET imaging.
The high dose conformity and healthy tissue sparing achievable in Particle Therapy when using C ions calls for safety factors in treatment planning, to prevent the tumor under-dosage related to the ...possible occurrence of inter-fractional morphological changes during a treatment. This limitation could be overcome by a range monitor, still missing in clinical routine, capable of providing on-line feedback. The Dose Profiler (DP) is a detector developed within the INnovative Solution for In-beam Dosimetry in hadronthErapy (INSIDE) collaboration for the monitoring of carbon ion treatments at the CNAO facility (Centro Nazionale di Adroterapia Oncologica) exploiting the detection of charged secondary fragments that escape from the patient. The DP capability to detect inter-fractional changes is demonstrated by comparing the obtained fragment emission maps in different fractions of the treatments enrolled in the first ever clinical trial of such a monitoring system, performed at CNAO. The case of a CNAO patient that underwent a significant morphological change is presented in detail, focusing on the implications that can be drawn for the achievable inter-fractional monitoring DP sensitivity in real clinical conditions. The results have been cross-checked against a simulation study.
Positron emission tomography (PET) is one of the most mature techniques for monitoring in ion beam therapy. PET allows to reconstruct the β+ activity generated in the patient by the nuclear ...interaction of the ions. Taking advantage of the spatial correlation between positron emitters created along the ions path and the dose distribution, it is possible to perform a quality control of the treatment. Usually, to reconstruct the activity generated within the irradiated volume, standard 3D PET reconstruction techniques are implemented. In this work, we explore a new reconstruction method (Straight Forward Reconstruction) particularly useful for reconstructing activity distribution generated by mono-energetic pencil beams. The method was validated by measuring the produced activity distribution with the DoPET system. Irradiations performed with mono-energetic pencil beams on phantoms mimicking human tissues were used for this study. Both reconstruction methods reach an accuracy in the reconstruction of the activity distribution width of the order of 1.5 mm for 2⋅108 primaries.
•PET is one of the most mature techniques for monitoring in ion beam therapy.•A new reconstruction method (Straight Forward Reconstruction) for PET is presented.•Activity distribution and accuracy by pencil beams irradiation are presented.
In this paper, we investigate the use of a positron emission tomography (PET) system to monitor the proton therapy. The monitoring procedure is based on the comparison between the <inline-formula> ...<tex-math notation="LaTeX">{ \beta }+ </tex-math></inline-formula> activity generated in the irradiated volume during the treatment, with the <inline-formula> <tex-math notation="LaTeX">{ \beta }+ </tex-math></inline-formula> activity distribution obtained with Monte Carlo (MC) simulation. The dedicated PET system is a dual head detection system; each head is composed of nine scintillating LYSO crystal matrices read out independently with a custom modularized acquisition system. Our experimental data were acquired at the Cyclotron Centre Bronowice, Institute Nuclear Physics in Kraków, Poland, and were simulated with the FLUKA MC code. Homogeneous and heterogeneous plastic phantoms were irradiated with monoenergetic 130 MeV protons. The capabilities of our PET system to distinguish different irradiated materials were investigated, and the proton pencil-beams were used as probes. Our focus was to analyze the activity width and the total activity event number in several cases. Irradiations were performed using either single pencil-beams one at a time, or two pencil-beams during the same data taking. The comparison of 1-D activity profile for experimental data and MC simulation were always in good agreement showing that, the treatment quality assessment in proton therapy can be based on <inline-formula> <tex-math notation="LaTeX">{ \beta }+ </tex-math></inline-formula> activity measurements.
Display omitted
•Range monitoring is necessary to fully exploit hadrontherapy advantages.•The Monte Carlo (MC) simulation tool for the INSIDE in-beam PET scanner is described.•The Monte Carlo tool ...provides images directly comparable with the experimental ones.•The comparison with in vivo experimental data validates the tool reliability.•The simulation information is already reliable after the first half of the treatment.
Hadrontherapy is a method for treating cancer with very targeted dose distributions and enhanced radiobiological effects. To fully exploit these advantages, in vivo range monitoring systems are required. These devices measure, preferably during the treatment, the secondary radiation generated by the beam-tissue interactions. However, since correlation of the secondary radiation distribution with the dose is not straightforward, Monte Carlo (MC) simulations are very important for treatment quality assessment.
The INSIDE project constructed an in-beam PET scanner to detect signals generated by the positron-emitting isotopes resulting from projectile-target fragmentation. In addition, a FLUKA-based simulation tool was developed to predict the corresponding reference PET images using a detailed scanner model.
The INSIDE in-beam PET was used to monitor two consecutive proton treatment sessions on a patient at the Italian Center for Oncological Hadrontherapy (CNAO). The reconstructed PET images were updated every 10 s providing a near real-time quality assessment. By half-way through the treatment, the statistics of the measured PET images were already significant enough to be compared with the simulations with average differences in the activity range less than 2.5 mm along the beam direction. Without taking into account any preferential direction, differences within 1 mm were found.
In this paper, the INSIDE MC simulation tool is described and the results of the first in vivo agreement evaluation are reported. These results have justified a clinical trial, in which the MC simulation tool will be used on a daily basis to study the compliance tolerances between the measured and simulated PET images.
Purpose:
A fully three-dimensional (3D) massively parallelizable list-mode ordered-subsets expectation-maximization (LM-OSEM) reconstruction algorithm has been developed for high-resolution PET ...cameras. System response probabilities are calculated online from a set of parameters derived from Monte Carlo simulations. The shape of a system response for a given line of response (LOR) has been shown to be asymmetrical around the LOR. This work has been focused on the development of efficient region-search techniques to sample the system response probabilities, which are suitable for asymmetric kernel models, including elliptical Gaussian models that allow for high accuracy and high parallelization efficiency. The novel region-search scheme using variable kernel models is applied in the proposed PET reconstruction algorithm.
Methods:
A novel region-search technique has been used to sample the probability density function in correspondence with a small dynamic subset of the field of view that constitutes the region of response (ROR). The ROR is identified around the LOR by searching for any voxel within a dynamically calculated contour. The contour condition is currently defined as a fixed threshold over the posterior probability, and arbitrary kernel models can be applied using a numerical approach. The processing of the LORs is distributed in batches among the available computing devices, then, individual LORs are processed within different processing units. In this way, both multicore and multiple many-core processing units can be efficiently exploited. Tests have been conducted with probability models that take into account the noncolinearity, positron range, and crystal penetration effects, that produced tubes of response with varying elliptical sections whose axes were a function of the crystal's thickness and angle of incidence of the given LOR. The algorithm treats the probability model as a 3D scalar field defined within a reference system aligned with the ideal LOR.
Results:
This new technique provides superior image quality in terms of signal-to-noise ratio as compared with the histogram-mode method based on precomputed system matrices available for a commercial small animal scanner. Reconstruction times can be kept low with the use of multicore, many-core architectures, including multiple graphic processing units.
Conclusions
: A highly parallelizable LM reconstruction method has been proposed based on Monte Carlo simulations and new parallelization techniques aimed at improving the reconstruction speed and the image signal-to-noise of a given OSEM algorithm. The method has been validated using simulated and real phantoms. A special advantage of the new method is the possibility of defining dynamically the cut-off threshold over the calculated probabilities thus allowing for a direct control on the trade-off between speed and quality during the reconstruction.