|
|
Roser, J., Barrientos, L., Bernabeu, J., Borja-Lloret, M., Muñoz, E., Ros, A., et al. (2022). Joint image reconstruction algorithm in Compton cameras. Phys. Med. Biol., 67(15), 155009–15pp.
Abstract: Objective. To demonstrate the benefits of using an joint image reconstruction algorithm based on the List Mode Maximum Likelihood Expectation Maximization that combines events measured in different channels of information of a Compton camera. Approach. Both simulations and experimental data are employed to show the algorithm performance. Main results. The obtained joint images present improved image quality and yield better estimates of displacements of high-energy gamma-ray emitting sources. The algorithm also provides images that are more stable than any individual channel against the noisy convergence that characterizes Maximum Likelihood based algorithms. Significance. The joint reconstruction algorithm can improve the quality and robustness of Compton camera images. It also has high versatility, as it can be easily adapted to any Compton camera geometry. It is thus expected to represent an important step in the optimization of Compton camera imaging.
|
|
|
|
Borys, D. et al, & Brzezinski, K. (2022). ProTheRaMon-a GATE simulation framework for proton therapy range monitoring using PET imaging. Phys. Med. Biol., 67(22), 224002–15pp.
Abstract: 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 (Borys et al 2022).
|
|
|
|
Borja-Lloret, M., Barrientos, L., Bernabeu, J., Lacasta, C., Muñoz, E., Ros, A., et al. (2023). Influence of the background in Compton camera images for proton therapy treatment monitoring. Phys. Med. Biol., 68(14), 144001–16pp.
Abstract: Objective. Background events are one of the most relevant contributions to image degradation in Compton camera imaging for hadron therapy treatment monitoring. A study of the background and its contribution to image degradation is important to define future strategies to reduce the background in the system. Approach. In this simulation study, the percentage of different kinds of events and their contribution to the reconstructed image in a two-layer Compton camera have been evaluated. To this end, GATE v8.2 simulations of a proton beam impinging on a PMMA phantom have been carried out, for different proton beam energies and at different beam intensities. Main results. For a simulated Compton camera made of Lanthanum (III) Bromide monolithic crystals, coincidences caused by neutrons arriving from the phantom are the most common type of background produced by secondary radiations in the Compton camera, causing between 13% and 33% of the detected coincidences, depending on the beam energy. Results also show that random coincidences are a significant cause of image degradation at high beam intensities, and their influence in the reconstructed images is studied for values of the time coincidence windows from 500 ps to 100 ns. Significance. Results indicate the timing capabilities required to retrieve the fall-off position with good precision. Still, the noise observed in the image when no randoms are considered make us consider further background rejection methods.
|
|
|
|
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.
|
|
