|
Llosa, G., Barrio, J., Lacasta, C., Bisogni, M. G., Del Guerra, A., Marcatili, S., et al. (2010). Characterization of a PET detector head based on continuous LYSO crystals and monolithic, 64-pixel silicon photomultiplier matrices. Phys. Med. Biol., 55(23), 7299–7315.
Abstract: The characterization of a PET detector head based on continuous LYSO crystals and silicon photomultiplier (SiPM) arrays as photodetectors has been carried out for its use in the development of a small animal PET prototype. The detector heads are composed of a continuous crystal and a SiPM matrix with 64 pixels in a common substrate, fabricated specifically for this project. Three crystals of 12 mm x 12 mm x 5 mm size with different types of painting have been tested: white, black and black on the sides but white on the back of the crystal. The best energy resolution, obtained with the white crystal, is 16% FWHM. The detector response is linear up to 1275 keV. Tests with different position determination algorithms have been carried out with the three crystals. The spatial resolution obtained with the center of gravity algorithm is around 0.9 mm FWHM for the three crystals. As expected, the use of this algorithm results in the displacement of the reconstructed position toward the center of the crystal, more pronounced in the case of the white crystal. A maximum likelihood algorithm has been tested that can reconstruct correctly the interaction position of the photons also in the case of the white crystal.
|
|
|
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.
|
|
|
Blume, M., Navab, N., & Rafecas, M. (2012). Joint image and motion reconstruction for PET using a B-spline motion model. Phys. Med. Biol., 57(24), 22pp.
Abstract: We present a novel joint image and motion reconstruction method for PET. The method is based on gated data and reconstructs an image together with amotion function. The motion function can be used to transform the reconstructed image to any of the input gates. All available events (from all gates) are used in the reconstruction. The presented method uses a B-spline motion model, together with a novel motion regularization procedure that does not need a regularization parameter (which is usually extremely difficult to adjust). Several image and motion grid levels are used in order to reduce the reconstruction time. In a simulation study, the presented method is compared to a recently proposed joint reconstruction method. While the presented method provides comparable reconstruction quality, it is much easier to use since no regularization parameter has to be chosen. Furthermore, since the B-spline discretization of the motion function depends on fewer parameters than a displacement field, the presented method is considerably faster and consumes less memory than its counterpart. The method is also applied to clinical data, for which a novel purely data-driven gating approach is presented.
|
|
|
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.
|
|
|
Solevi, P. et al, Oliver, J. F., Gillam, J. E., & Rafecas, M. (2013). A Monte-Carlo based model of the AX-PET demonstrator and its experimental validation. Phys. Med. Biol., 58(16), 5495–5510.
Abstract: AX-PET is a novel PET detector based on axially oriented crystals and orthogonal wavelength shifter (WLS) strips, both individually read out by silicon photo-multipliers. Its design decouples sensitivity and spatial resolution, by reducing the parallax error due to the layered arrangement of the crystals. Additionally the granularity of AX-PET enhances the capability to track photons within the detector yielding a large fraction of inter-crystal scatter events. These events, if properly processed, can be included in the reconstruction stage further increasing the sensitivity. Its unique features require dedicated Monte-Carlo simulations, enabling the development of the device, interpreting data and allowing the development of reconstruction codes. At the same time the non-conventional design of AX-PET poses several challenges to the simulation and modeling tasks, mostly related to the light transport and distribution within the crystals and WLS strips, as well as the electronics readout. In this work we present a hybrid simulation tool based on an analytical model and a Monte-Carlo based description of the AX-PET demonstrator. It was extensively validated against experimental data, providing excellent agreement.
|
|
|
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.
|
|
|
Brzezinski, K., Oliver, J. F., Gillam, J., & Rafecas, M. (2014). Study of a high-resolution PET system using a Silicon detector probe. Phys. Med. Biol., 59(20), 6117–6140.
Abstract: A high-resolution silicon detector probe, in coincidence with a conventional PET scanner, is expected to provide images of higher quality than those achievable using the scanner alone. Spatial resolution should improve due to the finer pixelization of the probe detector, while increased sensitivity in the probe vicinity is expected to decrease noise. A PET-probe prototype is being developed utilizing this principle. The system includes a probe consisting of ten layers of silicon detectors, each a 80 x 52 array of 1 x 1 x 1 mm(3) pixels, to be operated in coincidence with a modern clinical PET scanner. Detailed simulation studies of this system have been performed to assess the effect of the additional probe information on the quality of the reconstructed images. A grid of point sources was simulated to study the contribution of the probe to the system resolution at different locations over the field of view (FOV). A resolution phantom was used to demonstrate the effect on image resolution for two probe positions. A homogeneous source distribution with hot and cold regions was used to demonstrate that the localized improvement in resolution does not come at the expense of the overall quality of the image. Since the improvement is constrained to an area close to the probe, breast imaging is proposed as a potential application for the novel geometry. In this sense, a simplified breast phantom, adjacent to heart and torso compartments, was simulated and the effect of the probe on lesion detectability, through measurements of the local contrast recovery coefficient-to-noise ratio (CNR), was observed. The list-mode ML-EM algorithm was used for image reconstruction in all cases. As expected, the point spread function of the PET-probe system was found to be non-isotropic and vary with position, offering improvement in specific regions. Increase in resolution, of factors of up to 2, was observed in the region close to the probe. Images of the resolution phantom showed visible improvement in resolution when including the probe in the simulations. The image quality study demonstrated that contrast and spill-over ratio in other areas of the FOV were not sacrificed for this enhancement. The CNR study performed on the breast phantom indicates increased lesion detectability provided by the probe.
|
|
|
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.
|
|
|
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
|
|