|
|
Robert, C., Dedes, G., Battistoni, G., Bohlen, T. T., Buvat, I., Cerutti, F., et al. (2013). Distributions of secondary particles in proton and carbon-ion therapy: a comparison between GATE/Geant4 and FLUKA Monte Carlo codes. Phys. Med. Biol., 58(9), 2879–2899.
Abstract: Monte Carlo simulations play a crucial role for in-vivo treatment monitoring based on PET and prompt gamma imaging in proton and carbon-ion therapies. The accuracy of the nuclear fragmentation models implemented in these codes might affect the quality of the treatment verification. In this paper, we investigate the nuclear models implemented in GATE/Geant4 and FLUKA by comparing the angular and energy distributions of secondary particles exiting a homogeneous target of PMMA. Comparison results were restricted to fragmentation of O-16 and C-12. Despite the very simple target and set-up, substantial discrepancies were observed between the two codes. For instance, the number of high energy (>1 MeV) prompt gammas exiting the target was about twice as large with GATE/Geant4 than with FLUKA both for proton and carbon ion beams. Such differences were not observed for the predicted annihilation photon production yields, for which ratios of 1.09 and 1.20 were obtained between GATE and FLUKA for the proton beam and the carbon ion beam, respectively. For neutrons and protons, discrepancies from 14% (exiting protons-carbon ion beam) to 57% (exiting neutrons-proton beam) have been identified in production yields as well as in the energy spectra for neutrons.
|
|
|
|
Cabello, J., Etxebeste, A., Llosa, G., & Ziegler, S. I. (2015). Simulation study of PET detector limitations using continuous crystals. Phys. Med. Biol., 60(9), 3673–3694.
Abstract: Continuous crystals can potentially obtain better intrinsic detector spatial resolution compared to pixelated crystals, additionally providing depth of interaction (DoI) information from the light distribution. To achieve high performance sophisticated interaction position estimation algorithms are required. There are a number of algorithms in the literature applied to different crystal dimensions and different photodetectors. However, the different crystal properties and photodetector array geometries have an impact on the algorithm performance. In this work we analysed, through Monte Carlo simulations, different combinations of realistic crystals and photodetector parameters to better understand their influence on the interaction position estimation accuracy, with special emphasis on the DoI. We used an interaction position estimation based on an analytical model for the present work. Different photodetector granulation schemes were investigated. The impact of the number of crystal faces readout by photodetectors was studied by simulating scenarios with one and two photodetectors. In addition, crystals with different levels of reflection and aspect ratios (AR) were analysed. Results showed that the impact of photodetector granularity is mainly shown near the edges and specially in the corners of the crystal. The resulting intrinsic spatial resolution near the centre with a 12 x 12 x 10 mm(3) LYSO crystal was 0.7-0.9 mm, while the average spatial resolution calculated on the entire crystal was 0.77 +/- 0.18 mm for all the simulated geometries with one and two photodetectors. Having front and back photodetectors reduced the DoI bias (Euclidean distance between estimated DoI and real DoI) and improved the transversal resolution near the corners. In scenarios with one photodetector, small AR resulted in DoI inaccuracies for absorbed events at the entrance of the crystal. These inaccuracies were slightly reduced either by increasing the AR or reducing the amount of reflected light, and highly mitigated using two photodetectors. Using one photodetector, we obtained a piecewise DoI error model with a DoI resolution of 0.4-0.9 mm for a 1.2 AR crystal, and we observed that including a second photodetector or reducing the amount of reflections reduced the DoI bias but did not significantly improve the DoI resolution. Translating the piecewise DoI error model obtained in this study to image reconstruction we obtained a spatial resolution variability of 0.39 mm using 85% of the FoV, compared to 2.59 mm and 1.87 mm without DoI correction or with a dual layer system, respectively.
|
|
|
|
Torres-Espallardo, I., Diblen, F., Rohling, H., Solevi, P., Gillam, J., Watts, D., et al. (2015). Evaluation of resistive-plate-chamber-based TOF-PET applied to in-beam particle therapy monitoring. Phys. Med. Biol., 60(9), N187–N208.
Abstract: Particle therapy is a highly conformal radiotherapy technique which reduces the dose deposited to the surrounding normal tissues. In order to fully exploit its advantages, treatment monitoring is necessary to minimize uncertainties related to the dose delivery. Up to now, the only clinically feasible technique for the monitoring of therapeutic irradiation with particle beams is Positron Emission Tomography (PET). In this work we have compared a Resistive Plate Chamber (RPC)-based PET scanner with a scintillation-crystal-based PET scanner for this application. In general, the main advantages of the RPC-PET system are its excellent timing resolution, low cost, and the possibility of building large area systems. We simulated a partial-ring scannerbeam monitoring, which has an intrinsically low positron yield compared to diagnostic PET. In addition, for in-beam PET there is a further data loss due to the partial ring configuration. In order to improve the performance of the RPC-based scanner, an improved version of the RPC detector (modifying the thickness of the gas and glass layers), providing a larger sensitivity, has been simulated and compared with an axially extended version of the crystal-based device. The improved version of the RPC shows better performance than the prototype, but the extended version of the crystal-based PET outperforms all other options. based on an RPC prototype under construction within the Fondazione per Adroterapia Oncologica (TERA). For comparison with the crystal-based PET scanner we have chosen the geometry of a commercially available PET scanner, the Philips Gemini TF. The coincidence time resolution used in the simulations takes into account the current achievable values as well as expected improvements of both technologies. Several scenarios (including patient data) have been simulated to evaluate the performance of different scanners. Initial results have shown that the low sensitivity of the RPC hampers its application to hadron
|
|
|
|
Cabello, J., & Rafecas, M. (2012). Comparison of basis functions for 3D PET reconstruction using a Monte Carlo system matrix. Phys. Med. Biol., 57(7), 1759–1777.
Abstract: In emission tomography, iterative statistical methods are accepted as the reconstruction algorithms that achieve the best image quality. The accuracy of these methods relies partly on the quality of the system response matrix (SRM) that characterizes the scanner. The more physical phenomena included in the SRM, the higher the SRM quality, and therefore higher image quality is obtained from the reconstruction process. High-resolution small animal scanners contain as many as 10(3)-10(4) small crystal pairs, while the field of view (FOV) is divided into hundreds of thousands of small voxels. These two characteristics have a significant impact on the number of elements to be calculated in the SRM. Monte Carlo (MC) methods have gained popularity as a way of calculating the SRM, due to the increased accuracy achievable, at the cost of introducing some statistical noise and long simulation times. In the work presented here the SRM is calculated using MC methods exploiting the cylindrical symmetries of the scanner, significantly reducing the simulation time necessary to calculate a high statistical quality SRM and the storage space necessary. The use of cylindrical symmetries makes polar voxels a convenient basis function. Alternatively, spherically symmetric basis functions result in improved noise properties compared to cubic and polar basis functions. The quality of reconstructed images using polar voxels, spherically symmetric basis functions on a polar grid, cubic voxels and post-reconstruction filtered polar and cubic voxels is compared from a noise and spatial resolution perspective. This study demonstrates that polar voxels perform as well as cubic voxels, reducing the simulation time necessary to calculate the SRM and the disk space necessary to store it. Results showed that spherically symmetric functions outperform polar and cubic basis functions in terms of noise properties, at the cost of slightly degraded spatial resolution, larger SRM file size and longer reconstruction times. However, we demonstrate that post-reconstruction smoothing, usually applied in emission imaging to reduce the level of noise, can produce a spatial resolution degradation of similar to 50%, while spherically symmetric basis functions produce a degradation of only similar to 6%, compared to polar and cubic voxels, at the same noise level. Therefore, the image quality trade-off obtained with blobs is higher than that obtained with cubic or polar voxels.
|
|
|
|
Gillam, J. E., Solevi, P., Oliver, J. F., & Rafecas, M. (2013). Simulated one-pass list-mode: an approach to on-the-fly system matrix calculation. Phys. Med. Biol., 58(7), 2377–2394.
Abstract: In the development of prototype systems for positron emission tomography a valid and robust image reconstruction algorithm is required. However, prototypes often employ novel detector and system geometries which may change rapidly under optimization. In addition, developing systems generally produce highly granular, or possibly continuous detection domains which require some level of on-the-fly calculation for retention of measurement precision. In this investigation a new method of on-the-fly system matrix calculation is proposed that provides advantages in application to such list-mode systems in terms of flexibility in system modeling. The new method is easily adaptable to complicated system geometries and available computational resources. Detection uncertainty models are used as random number generators to produce ensembles of possible photon trajectories at image reconstruction time for each datum in the measurement list. However, the result of this approach is that the system matrix elements change at each iteration in a non-repetitive manner. The resulting algorithm is considered the simulation of a one-pass list (SOPL) which is generated and the list traversed during image reconstruction. SOPL alters the system matrix in use at each iteration and so behavior within the maximum likelihood-expectation maximization algorithm was investigated. A two-pixel system and a small two dimensional imaging model are used to illustrate the process and quantify aspects of the algorithm. The two-dimensional imaging system showed that, while incurring a penalty in image resolution, in comparison to a non-random equal-computation counterpart, SOPL provides much enhanced noise properties. In addition, enhancement in system matrix quality is straightforward (by increasing the number of samples in the ensemble) so that the resolution penalty can be recovered when desired while retaining improvement in noise properties. Finally the approach is tested and validated against a standard (highly accurate) system matrix using experimental data from a prototype system-the AX-PET.
|
|
|
|
Piersanti, L., Bellini, F., Bini, F., Collamati, F., De Lucia, E., Durante, M., et al. (2014). Measurement of charged particle yields from PMMA irradiated by a 220 MeV/u C-12 beam. Phys. Med. Biol., 59(7), 1857–1872.
Abstract: The radiation used in hadrontherapy treatments interacts with the patient body producing secondary particles, either neutral or charged, that can be used for dose and Bragg peak monitoring and to provide a fast feedback on the treatment plans. Recent results obtained from the authors on simplified setups (mono-energetic primary beams interacting with homogeneous tissue like target) have already indicated the correlation that exists between the flux of these secondaries coming from the target (e.g. protons and photons) and the position of the primary beam Bragg peak. In this paper, the measurements of charged particle fluxes produced by the interaction of a 220 MeV/u carbon ion beam at GSI, Darmstadt, with a polymethyl methacrylate target are reported. The emission region of protons (p), deuterons (d) and tritons (t) has been characterized using a drift chamber while the particle time-of-flight, used to compute the kinetic energy spectra, was measured with a LYSO scintillator.The energy released in the LYSO crystal was used for particle identification purposes. The measurements were repeated with the setup at 60 degrees and 90 degrees with respect to the primary beam direction. The accuracy on the fragments emission profile reconstruction and its relationship with the Bragg peak position have been studied. Based on the acquired experimental evidence, a method to monitor the dose profile and the position of the Bragg peak inside the target is proposed.
|
|
|
|
Cabello, J., & Wells, K. (2010). The spatial resolution of silicon-based electron detectors in beta-autoradiography. Phys. Med. Biol., 55(6), 1677–1699.
Abstract: Thin tissue autoradiography is an imaging modality where ex-vivo tissue sections are placed in direct contact with autoradiographic film. These tissue sections contain a radiolabelled ligand bound to a specific biomolecule under study. This radioligand emits beta- or beta+ particles ionizing silver halide crystals in the film. High spatial resolution autoradiograms are obtained using low energy radioisotopes, such as H-3 where an intrinsic 0.1-1 μm spatial resolution can be achieved. Several digital alternatives have been presented over the past few years to replace conventional film but their spatial resolution has yet to equal film, although silicon-based imaging technologies have demonstrated higher sensitivity compared to conventional film. It will be shown in this work how pixel size is a critical parameter for achieving high spatial resolution for low energy uncollimated beta imaging. In this work we also examine the confounding factors impeding silicon-based technologies with respect to spatial resolution. The study considers charge diffusion in silicon and detector noise, and this is applied to a range of radioisotopes typically used in autoradiography. Finally an optimal detector geometry to obtain the best possible spatial resolution for a specific technology and a specific radioisotope is suggested.
|
|
|
|
Ortega, P. G., Torres-Espallardo, I., Cerutti, F., Ferrari, A., Gillam, J. E., Lacasta, C., et al. (2015). Noise evaluation of Compton camera imaging for proton therapy. Phys. Med. Biol., 60(5), 1845–1863.
Abstract: Compton Cameras emerged as an alternative for real-time dose monitoring techniques for Particle Therapy (PT), based on the detection of prompt-gammas. As a consequence of the Compton scattering process, the gamma origin point can be restricted onto the surface of a cone (Compton cone). Through image reconstruction techniques, the distribution of the gamma emitters can be estimated, using cone-surfaces backprojections of the Compton cones through the image space, along with more sophisticated statistical methods to improve the image quality. To calculate the Compton cone required for image reconstruction, either two interactions, the last being photoelectric absorption, or three scatter interactions are needed. Because of the high energy of the photons in PT the first option might not be adequate, as the photon is not absorbed in general. However, the second option is less efficient. That is the reason to resort to spectral reconstructions, where the incoming. energy is considered as a variable in the reconstruction inverse problem. Jointly with prompt gamma, secondary neutrons and scattered photons, not strongly correlated with the dose map, can also reach the imaging detector and produce false events. These events deteriorate the image quality. Also, high intensity beams can produce particle accumulation in the camera, which lead to an increase of random coincidences, meaning events which gather measurements from different incoming particles. The noise scenario is expected to be different if double or triple events are used, and consequently, the reconstructed images can be affected differently by spurious data. The aim of the present work is to study the effect of false events in the reconstructed image, evaluating their impact in the determination of the beam particle ranges. A simulation study that includes misidentified events (neutrons and random coincidences) in the final image of a Compton Telescope for PT monitoring is presented. The complete chain of detection, from the beam particle entering a phantom to the event classification, is simulated using FLUKA. The range determination is later estimated from the reconstructed image obtained from a two and three-event algorithm based on Maximum Likelihood Expectation Maximization. The neutron background and random coincidences due to a therapeutic-like time structure are analyzed for mono-energetic proton beams. The time structure of the beam is included in the simulations, which will affect the rate of particles entering the detector.
|
|
|
|
Etxebeste, A., Dauvergne, D., Fontana, M., Letang, J. M., Llosa, G., Muñoz, E., et al. (2020). CCMod: a GATE module for Compton camera imaging simulation. Phys. Med. Biol., 65(5), 055004–17pp.
Abstract: Compton cameras are gamma-ray imaging systems which have been proposed for a wide variety of applications such as medical imaging, nuclear decommissioning or homeland security. In the design and optimization of such a system Monte Carlo simulations play an essential role. In this work, we propose a generic module to perform Monte Carlo simulations and analyses of Compton Camera imaging which is included in the open-source GATE/Geant4 platform. Several digitization stages have been implemented within the module to mimic the performance of the most commonly employed detectors (e.g. monolithic blocks, pixelated scintillator crystals, strip detectors...). Time coincidence sorter and sequence coincidence reconstruction are also available in order to aim at providing modules to facilitate the comparison and reproduction of the data taken with different prototypes. All processing steps may be performed during the simulation (on-the-fly mode) or as a post-process of the output files (offline mode). The predictions of the module have been compared with experimental data in terms of energy spectra, angular resolution, efficiency and back-projection image reconstruction. Consistent results within a 3-sigma interval were obtained for the energy spectra except for low energies where small differences arise. The angular resolution measure for incident photons of 1275 keV was also in good agreement between both data sets with a value close to 13 degrees. Moreover, with the aim of demonstrating the versatility of such a tool the performance of two different Compton camera designs was evaluated and compared.
|
|
|
|
Etxebeste, A., Barrio, J., Bernabeu, J., Lacasta, C., Llosa, G., Muñoz, E., et al. (2019). Study of sensitivity and resolution for full ring PET prototypes based on continuous crystals and analytical modeling of the light distribution. Phys. Med. Biol., 64(3), 035015–17pp.
Abstract: Sensitivity and spatial resolution are the main parameters to maximize in the performance of a PET scanner. For this purpose, detectors consisting of a combination of continuous crystals optically coupled to segmented photodetectors have been employed. With the use of continuous crystals the sensitivity is increased with respect to the pixelated crystals. In addition, spatial resolution is no longer limited to the crystal size. The main drawback is the difficulty in determining the interaction position. In this work, we present the characterization of the performance of a full ring based on cuboid continuous crystals coupled to SiPMs. To this end, we have employed the simulations developed in a previous work for our experimental detector head. Sensitivity could be further enhanced by using tapered crystals. This enhancement is obtained by increasing the solid angle coverage, reducing the wedge-shaped gaps between contiguous detectors. The performance of the scanners based on both crystal geometries was characterized following NEMA NU 4-2008 standardized protocol in order to compare them. An average sensitivity gain over the entire axial field of view of 13.63% has been obtained with tapered geometry while similar performance of the spatial resolution has been proven with both scanners. The activity at which NECR and true peak occur is smaller and the peak value is greater for tapered crystals than for cuboid crystals. Moreover, a higher degree of homogeneity was obtained in the sensitivity map due to the tighter packing of the crystals, which reduces the gaps and results in a better recovery of homogeneous regions than for the cuboid configuration. Some of the results obtained, such as spatial resolution, depend on the interaction position estimation and may vary if other method is employed.
|
|
