Home | << 1 2 3 4 >> |
![]() |
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.
|
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.
|
Ros Garcia, A., Barrio, J., Etxebeste, A., Garcia-Lopez, J., Jimenez-Ramos, M. C., Lacasta, C., et al. (2020). MACACO II test-beam with high energy photons. Phys. Med. Biol., 65(24), 245027–12pp.
Abstract: The IRIS group at IFIC Valencia is developing a three-layer Compton camera for treatment monitoring in proton therapy. The system is composed of three detector planes, each made of a LaBr3<i monolithic crystal coupled to a SiPM array. Having obtained successful results with the first prototype (MACACO) that demonstrated the feasibility of the proposed technology, a second prototype (MACACO II) with improved performance has been developed, and is the subject of this work. The new system has an enhanced detector energy resolution which translates into a higher spatial resolution of the telescope. The image reconstruction method has also been improved with an accurate model of the sensitivity matrix. The device has been tested with high energy photons at the National Accelerator Centre (CNA, Seville). The tests involved a proton beam of 18 MeV impinging on a graphite target, to produce 4.4 MeV photons. Data were taken at different system positions of the telescope with the first detector at 65 and 160 mm from the target, and at different beam intensities. The measurements allowed successful reconstruction of the photon emission distribution at two target positions separated by 5 mm in different telescope configurations. This result was obtained both with data recorded in the first and second telescope planes (two interaction events) and, for the first time in beam experiments, with data recorded in the three planes (three interaction events).
Keywords: Compton imaging; Compton camera; proton therapy; LaBr3; test-beam; image reconstruction
|
Valdes-Cortez, C., Ballester, F., Vijande, J., Gimenez, V., Gimenez-Alventosa, V., Perez-Calatayud, J., et al. (2020). Depth-dose measurement corrections for the surface electronic brachytherapy beams of an Esteya(R) unit: a Monte Carlo study. Phys. Med. Biol., 65(24), 245026–12pp.
Abstract: Three different correction factors for measurements with the parallel-plate ionization chamber PTW T34013 on the Esteya electronic brachytherapy unit have been investigated. This chamber type is recommended by AAPM TG-253 for depth-dose measurements in the 69.5 kV x-ray beam generated by the Esteya unit. Monte Carlo simulations using the PENELOPE-2018 system were performed to determine the absorbed dose deposited in water and in the chamber sensitive volume at different depths with a Type A uncertainty smaller than 0.1%. Chamber-to-chamber differences have been explored performing measurements using three different chambers. The range of conical applicators available, from 10 to 30 mm in diameter, has been explored. Using a depth-independent global chamber perturbation correction factor without a shift of the effective point of measurement yielded differences between the absorbed dose to water and the corrected absorbed dose in the sensitive volume of the chamber of up to 1% and 0.6% for the 10 mm and 30 mm applicators, respectively. Calculations using a depth-dependent perturbation factor, including or excluding a shift of the effective point of measurement, resulted in depth-dose differences of about +/- 0.5% or less for both applicators. The smallest depth-dose differences were obtained when a shift of the effective point of measurement was implemented, being displaced 0.4 mm towards the center of the sensitive volume of the chamber. The correction factors were obtained with combined uncertainties of 0.4% (k = 2). Uncertainties due to chamber-to-chamber differences are found to be lower than 2%. The results emphasize the relevance of carrying out detailed Monte Carlo studies for each electronic brachytherapy device and ionization chamber used for its dosimetry.
Keywords: electronic brachytherapy; eBT; dosimetry; ionization chamber; Monte Carlo
|
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.
|
Oliver, J. F., & Rafecas, M. (2010). Improving the singles rate method for modeling accidental coincidences in high-resolution PET. Phys. Med. Biol., 55(22), 6951–6971.
Abstract: Random coincidences ('randoms') are one of the main sources of image degradation in PET imaging. In order to correct for this effect, an accurate method to estimate the contribution of random events is necessary. This aspect becomes especially relevant for high-resolution PET scanners where the highest image quality is sought and accurate quantitative analysis is undertaken. One common approach to estimate randoms is the so-called singles rate method (SR) widely used because of its good statistical properties. SR is based on the measurement of the singles rate in each detector element. However, recent studies suggest that SR systematically overestimates the correct random rate. This overestimation can be particularly marked for low energy thresholds, below 250 keV used in some applications and could entail a significant image degradation. In this work, we investigate the performance of SR as a function of the activity, geometry of the source and energy acceptance window used. We also investigate the performance of an alternative method, which we call 'singles trues' (ST) that improves SR by properly modeling the presence of true coincidences in the sample. Nevertheless, in any real data acquisition the knowledge of which singles are members of a true coincidence is lost. Therefore, we propose an iterative method, STi, that provides an estimation based on ST but which only requires the knowledge of measurable quantities: prompts and singles. Due to inter-crystal scatter, for wide energy windows ST only partially corrects SR overestimations. While SR deviations are in the range 86-300% (depending on the source geometry), the ST deviations are systematically smaller and contained in the range 4-60%. STi fails to reproduce the ST results, although for not too high activities the deviation with respect to ST is only a few percent. For conventional energy windows, i.e. those without inter-crystal scatter, the ST method corrects the SR overestimations, and deviations from the true random rate are of the order of 1% or less. In addition, in the case of conventional energy window STi results reproduce ST results and therefore the former can be used to obtain the true random rate.
|
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).
Keywords: proton therapy; GATE; Monte Carlo simulations; J-PET; medical imaging
|
Hueso-Gonzalez, F., Berthold, J., Wohlfahrt, P., Bortfeld, T., Khamfongkhruea, C., Tattenberg, S., et al. (2024). Inter-center comparison of proton range verification prototypes with an anthropomorphic head phantom. Phys. Med. Biol., 69(22), 225010–14pp.
Abstract: Objective. To compare in reproducible and equalized conditions the performance of two independent proton range verification systems based on prompt gamma-ray detectors from two different proton therapy centers. Approach. An anthropomorphic head phantom with calibrated stopping power, serving as ground truth, was irradiated with comparable treatment plans, spot positions and energies in both facilities. Clinical beam current, tumor contour and dose were used. The absolute range measurement was compared to the expected value according to the ground truth. The statistical precision was assessed by repeating each measurement ten times. Sensitivity to relative range shifts was evaluated by introducing 2 mm and 5 mm plastic slabs on half of the field. Main results. The resulting absolute range accuracy was within 2.4 mm in all cases. Relative range shifts were detected with deviations lower than 14%. Significance. The performance of both systems was deemed worthy of clinical application for the detection of range deviations. This study represents the first comparison of independent prompt gamma-ray-based proton range verification systems under equalized conditions with realistic treatment fields and beam currents.
|
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.
Keywords: PET; high-resolution imaging; Si detectors; PET insert
|
Gimenez-Alventosa, V., Gimenez, V., Ballester, F., Vijande, J., & Andreo, P. (2020). Monte Carlo calculation of beam quality correction factors for PTW cylindrical ionization chambers in photon beams. Phys. Med. Biol., 65(20), 205005–11pp.
Abstract: The beam quality correction factork(Q)for megavoltage photon beams has been calculated for eight PTW (Freiburg, Germany) ionization chambers (Farmer chambers PTW30010, PTW30011, PTW30012, and PTW30013, Semiflex 3D chambers PTW31021, PTW31010, and PTW31013, and the PinPoint 3D chamber PTW31016). Simulations performed on the widely used NE-2571 ionization chamber have been used to benchmark the results. The Monte Carlo code PENELOPE/penEasy was used to calculate the absorbed dose to a point in water and the absorbed dose to the active air volume of the chambers for photon beams in the range 4 to 24 MV. Of the nine ionization chambers analysed, only five are included in the current version of the International Code of Practice for dosimetry based on standards of absorbed dose to water (IAEA TRS 398). The values reported in this work agree with those in the literature within the uncertainty estimates and are to be included in the average values of the data obtained by different working groups for the forthcoming update of TRS 398.
Keywords: TRS 398; Monte Carlo; dosimetry; ionization chambers; MV photon beams
|