| Home | << 1 2 3 4 >> |
|
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
|
|
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
|
|
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
|
|
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.
|
|
Brzezinski, K. et al. (2023). Detection of range shifts in proton beam therapy using the J-PET scanner: a patient simulation study. Phys. Med. Biol., 68(14), 145016–17pp.
Abstract: Objective. The Jagiellonian positron emission tomography (J-PET) technology, based on plastic scintillators, has been proposed as a cost effective tool for detecting range deviations during proton therapy. This study investigates the feasibility of using J-PET for range monitoring by means of a detailed Monte Carlo simulation study of 95 patients who underwent proton therapy at the Cyclotron Centre Bronowice (CCB) in Krakow, Poland. Approach. Discrepancies between prescribed and delivered treatments were artificially introduced in the simulations by means of shifts in patient positioning and in the Hounsfield unit to the relative proton stopping power calibration curve. A dual-layer, cylindrical J-PET geometry was simulated in an in-room monitoring scenario and a triple-layer, dual-head geometry in an in-beam protocol. The distribution of range shifts in reconstructed PET activity was visualized in the beam's eye view. Linear prediction models were constructed from all patients in the cohort, using the mean shift in reconstructed PET activity as a predictor of the mean proton range deviation. Main results. Maps of deviations in the range of reconstructed PET distributions showed agreement with those of deviations in dose range in most patients. The linear prediction model showed a good fit, with coefficient of determination r (2) = 0.84 (in-room) and 0.75 (in-beam). Residual standard error was below 1 mm: 0.33 mm (in-room) and 0.23 mm (in-beam). Significance. The precision of the proposed prediction models shows the sensitivity of the proposed J-PET scanners to shifts in proton range for a wide range of clinical treatment plans. Furthermore, it motivates the use of such models as a tool for predicting proton range deviations and opens up new prospects for investigations into the use of intra-treatment PET images for predicting clinical metrics that aid in the assessment of the quality of delivered treatment.
|
|
Roser, J., Muñoz, E., Barrientos, L., Barrio, J., Bernabeu, J., Borja-Lloret, M., et al. (2020). Image reconstruction for a multi-layer Compton telescope: an analytical model for three interaction events. Phys. Med. Biol., 65(14), 145005–17pp.
Abstract: Compton Cameras are electronically collimated photon imagers suitable for sub-MeV to few MeV gamma-ray detection. Such features are desirable to enablein vivorange verification in hadron therapy, through the detection of secondary Prompt Gammas. A major concern with this technique is the poor image quality obtained when the incoming gamma-ray energy is unknown. Compton Cameras with more than two detector planes (multi-layer Compton Cameras) have been proposed as a solution, given that these devices incorporate more signal sequences of interactions to the conventional two interaction events. In particular, three interaction events convey more spectral information as they allow inferring directly the incident gamma-ray energy. A three-layer Compton Telescope based on continuous Lanthanum (III) Bromide crystals coupled to Silicon Photomultipliers is being developed at the IRIS group of IFIC-Valencia. In a previous work we proposed a spectral reconstruction algorithm for two interaction events based on an analytical model for the formation of the signal. To fully exploit the capabilities of our prototype, we present here an extension of the model for three interaction events. Analytical expressions of the sensitivity and the System Matrix are derived and validated against Monte Carlo simulations. Implemented in a List Mode Maximum Likelihood Expectation Maximization algorithm, the proposed model allows us to obtain four-dimensional (energy and position) images by using exclusively three interaction events. We are able to recover the correct spectrum and spatial distribution of gamma-ray sources when ideal data are employed. However, the uncertainties associated to experimental measurements result in a degradation when real data from complex structures are employed. Incorrect estimation of the incident gamma-ray interaction positions, and missing deposited energy associated with escaping secondaries, have been identified as the causes of such degradation by means of a detailed Monte Carlo study. As expected, our current experimental resolution and efficiency to three interaction events prevents us from correctly recovering complex structures of radioactive sources. However, given the better spectral information conveyed by three interaction events, we expect an improvement of the image quality of conventional Compton imaging when including such events. In this regard, future development includes the incorporation of the model assessed in this work to the two interaction events model in order to allow using simultaneously two and three interaction events in the image reconstruction.
|
|
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.
|
|
Muñoz, E., Barrio, J., Bernabeu, J., Etxebeste, A., Lacasta, C., Llosa, G., et al. (2018). Study and comparison of different sensitivity models for a two-plane Compton camera. Phys. Med. Biol., 63(13), 135004–19pp.
Abstract: Given the strong variations in the sensitivity of Compton cameras for the detection of events originating from different points in the field of view (FoV), sensitivity correction is often necessary in Compton image reconstruction. Several approaches for the calculation of the sensitivity matrix have been proposed in the literature. While most of these models are easily implemented and can be useful in many cases, they usually assume high angular coverage over the scattered photon, which is not the case for our prototype. In this work, we have derived an analytical model that allows us to calculate a detailed sensitivity matrix, which has been compared to other sensitivity models in the literature. Specifically, the proposed model describes the probability of measuring a useful event in a two-plane Compton camera, including the most relevant physical processes involved. The model has been used to obtain an expression for the system and sensitivity matrices for iterative image reconstruction. These matrices have been validated taking Monte Carlo simulations as a reference. In order to study the impact of the sensitivity, images reconstructed with our sensitivity model and with other models have been compared. Images have been reconstructed from several simulated sources, including point-like sources and extended distributions of activity, and also from experimental data measured with Na-22 sources. Results show that our sensitivity model is the best suited for our prototype. Although other models in the literature perform successfully in many scenarios, they are not applicable in all the geometrical configurations of interest for our system. In general, our model allows to effectively recover the intensity of point-like sources at different positions in the FoV and to reconstruct regions of homogeneous activity with minimal variance. Moreover, it can be employed for all Compton camera configurations, including those with low angular coverage over the scatterer.
Keywords: Compton camera imaging; MLEM; Monte Carlo simulations; image quality
|
|
Gimenez-Alventosa, V., Gimenez, V., Ballester, F., Vijande, J., & Andreo, P. (2018). Correction factors for ionization chamber measurements with the 'Valencia' and 'large field Valencia' brachytherapy applicators. Phys. Med. Biol., 63(12), 125004–10pp.
Abstract: Treatment of small skin lesions using HDR brachytherapy applicators is a widely used technique. The shielded applicators currently available in clinical practice are based on a tungsten-alloy cup that collimates the source-emitted radiation into a small region, hence protecting nearby tissues. The goal of this manuscript is to evaluate the correction factors required for dose measurements with a plane-parallel ionization chamber typically used in clinical brachytherapy for the 'Valencia' and 'large field Valencia' shielded applicators. Monte Carlo simulations have been performed using the PENELOPE-2014 system to determine the absorbed dose deposited in a water phantom and in the chamber active volume with a Type A uncertainty of the order of 0.1%. The average energies of the photon spectra arriving at the surface of the water phantom differ by approximately 10%, being 384 keV for the 'Valencia' and 343 keV for the 'large field Valencia'. The ionization chamber correction factors have been obtained for both applicators using three methods, their values depending on the applicator being considered. Using a depth-independent global chamber perturbation correction factor and no shift of the effective point of measurement yields depth-dose differences of up to 1% for the 'Valencia' applicator. Calculations using a depth-dependent global perturbation factor, or a shift of the effective point of measurement combined with a constant partial perturbation factor, result in differences of about 0.1% for both applicators. The results emphasize the relevance of carrying out detailed Monte Carlo studies for each shielded brachytherapy applicator and ionization chamber.
|