|
|
Aguiar, P., Rafecas, M., Ortuño, J. E., Kontaxakis, G., Santos, A., Pavia, J., et al. (2010). Geometrical and Monte Carlo projectors in 3D PET reconstruction. Med. Phys., 37(11), 5691–5702.
Abstract: Purpose: In the present work, the authors compare geometrical and Monte Carlo projectors in detail. The geometrical projectors considered were the conventional geometrical Siddon ray-tracer (S-RT) and the orthogonal distance-based ray-tracer (OD-RT), based on computing the orthogonal distance from the center of image voxel to the line-of-response. A comparison of these geometrical projectors was performed using different point spread function (PSF) models. The Monte Carlo-based method under consideration involves an extensive model of the system response matrix based on Monte Carlo simulations and is computed off-line and stored on disk. Methods: Comparisons were performed using simulated and experimental data of the commercial small animal PET scanner rPET. Results: The results demonstrate that the orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions yield better images in terms of contrast and spatial resolution than those obtained after using the conventional method and the multiray-based S-RT. Furthermore, the Monte Carlo-based method yields slight improvements in terms of contrast and spatial resolution with respect to these geometrical projectors. Conclusions: The orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions represent satisfactory alternatives to factorizing the system matrix or to the conventional on-the-fly ray-tracing methods for list-mode reconstruction, where an extensive modeling based on Monte Carlo simulations is unfeasible.
|
|
|
|
Vijande, J., Granero, D., Perez-Calatayud, J., & Ballester, F. (2013). Monte Carlo dosimetric study of the medium dose rate CSM40 source. Appl. Radiat. Isot., 82, 283–288.
Abstract: The Cs-137 medium dose rate (MDR) CSM40 source model (Eckert & Ziegler BEBIG, Germany) is in clinical use but no dosimetric dataset has been published. This study aims to obtain dosimetric data for the CSM40 source for its use in clinical practice as required by the American Association of Physicists in Medicine (AAPM) and the European Society for Radiotherapy and Oncology (ESTRO). Penelope2008 and Geant4 Monte Carlo codes were used to characterize this source dosimetrically. It was located in an unbounded water phantom with composition and mass density as recommended by AAPM and ESTRO. Due to the low photon energies of Cs-137, absorbed dose was approximated by collisional kerma. Additional simulations were performed to obtain the air-kerma strength, sic. Mass-energy absorption coefficients in water and air were consistently derived and used to calculate collisional kerma. Results performed with both radiation transport codes showed agreement typically within 0.05%. Dose rate constant, radial dose function and anisotropy function are provided for the CSM40 and compared with published data for other commercially available Cs-137 sources. An uncertainty analysis has been performed. The data provided by this study can be used as input data and verification in the treatment planning systems. (C) 2013 Elsevier Ltd. All rights reserved.
|
|
|
|
Gimenez-Alventosa, V., Vijande, J., Ballester, F., & Perez-Calatayud, J. (2016). Transit dose comparisons for Co-60 and Ir-192 HDR sources. J. Radiol. Prot., 36(4), 858–864.
Abstract: The goal of this study is to evaluate the ambient dose due to the transit of high dose rate (HDR) Co-60 sources along a transfer tube as compared to Ir-192 ones in a realistic clinical scenario. This goal is accomplished by evaluating air-kerma differences with Monte Carlo calculations using PENELOPE2011. Scatter from both the afterloader and the patient was not taken into account. Two sources, mHDR-v2 and Flexisource Co-60, (Elekta Brachytherapy, Veenendaal, the Netherlands) have been considered. These sources were simulated within a standard transfer tube located in an infinite air phantom. The movement of the source was included by displacing their positions along the connecting tube from z = – 75 cm to z = + 75 cm and combining them. Since modern afterloaders like Flexitron (Elekta) or Saginova (BEBIG GmbH) are able to use equally 192Ir and 60Co sources, it was assumed that both sources are displaced with equal speed. Typical HDR source activity content values were provided by the manufacturer. 2D distributions were obtained with type-A uncertainties (k = 2) less than 0.01%. From those, the air-kerma ratio Co-60/Ir-192 was evaluated weighted by their corresponding typical activities. It was found that it varies slowly with distance (less than 10% variation at 75 cm) but strongly in time due to the shorter half-life of the 192Ir (73.83 d). The maximum ratio is located close to the tube. It reaches a value of 0.57 when the typical activity of the sources at the time when they were installed by the vendor was used. Such ratio increases up to 1.28 at the end of the recommended working life (90 d) of the Ir-192 source. Co-60/Ir-192 air-kerma ratios are almost constant (0.51-0.57) in the vicinity of the source-tube with recent installed sources. Nevertheless, air-kerma ratios increase rapidly (1.15-1.29) whenever the Ir-192 is approaching the end of its life. In case of a medical event requiring the medical staff to access the treatment room, these ratios indicate that the dosimetric impact on the medical team will be lower, with a few exceptions, in the case of Co-60-based HDR brachytherapy as compared to Ir-192-based one when typical air-kerma strength values are considered.
|
|
|
|
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.
|
|
|
|
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.
|
|
|
|
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.
|
|
|
|
Calatayud-Jordan, J., Candela-Juan, C., Palma, J. D., Pujades-Claumarchirant, M. C., Soriano, A., Gracia-Ochoa, M., et al. (2021). Influence of the simultaneous calibration of multiple ring dosimeters on the individual absorbed dose. J. Radiol. Prot., 41(2), 384–397.
Abstract: Ring dosimeters for personal dosimetry are calibrated in accredited laboratories following ISO 4037-3 guidelines. The simultaneous irradiation of multiple dosimeters would save time, but has to be carefully studied, since the scattering conditions could change and influence the absorbed dose in nearby dosimeters. Monte Carlo simulations using PENELOPE-2014 were performed to explore the need to increase the uncertainty of H-p (0.07) in the simultaneous irradiation of three and five DXT-RAD 707H-2 (Thermo Scientific) ring dosimeters with beam qualities: N-30, N-80 and N-300. Results show that the absorbed dose in each dosimeter is compatible with each of the others and with the reference simulation (a single dosimeter), with a coverage probability of 95% (k = 2). Comparison with experimental data yielded consistent results with the same coverage probability. Therefore, five ring dosimeters can be simultaneously irradiated with beam qualities ranging, at least, between N-30 and N-300 with a negligible impact on the uncertainty of H-p (0.07).
|
|
|
|
Gimenez-Alventosa, V., Gimenez, V., & Oliver, S. (2021). PenRed: An extensible and parallel Monte-Carlo framework for radiation transport based on PENELOPE. Comput. Phys. Commun., 267, 108065–12pp.
Abstract: Monte Carlo methods provide detailed and accurate results for radiation transport simulations. Unfortunately, the high computational cost of these methods limits its usage in real-time applications. Moreover, existing computer codes do not provide a methodology for adapting these kinds of simulations to specific problems without advanced knowledge of the corresponding code system, and this restricts their applicability. To help solve these current limitations, we present PenRed, a general-purpose, standalone, extensible and modular framework code based on PENELOPE for parallel Monte Carlo simulations of electron-photon transport through matter. It has been implemented in C++ programming language and takes advantage of modern object-oriented technologies. In addition, PenRed offers the capability to read and process DICOM images as well as to construct and simulate image-based voxelized geometries, so as to facilitate its usage in medical applications. Our framework has been successfully verified against the original PENELOPE Fortran code. Furthermore, the implemented parallelism has been tested showing a significant improvement in the simulation time without any loss in precision of results. Program summary Program title: PenRed: Parallel Engine for Radiation Energy Deposition. CPC Library link to program files: https://doi .org /10 .17632/rkw6tvtngy.1 Licensing provision: GNU Affero General Public License (AGPL). Programming language: C++ standard 2011. Nature of problem: Monte Carlo simulations usually require a huge amount of computation time to achieve low statistical uncertainties. In addition, many applications necessitate particular characteristics or the extraction of specific quantities from the simulation. However, most available Monte Carlo codes do not provide an efficient parallel and truly modular structure which allows users to easily customise their code to suit their needs without an in-depth knowledge of the code system. Solution method: PenRed is a fully parallel, modular and customizable framework for Monte Carlo simulations of the passage of radiation through matter. It is based on the PENELOPE [1] code system, from which inherits its unique physics models and tracking algorithms for charged particles. PenRed has been coded in C++ following an object-oriented programming paradigm restricted to the C++11 standard. Our engine implements parallelism via a double approach: on the one hand, by using standard C++ threads for shared memory, improving the access and usage of the memory, and, on the other hand, via the MPI standard for distributed memory infrastructures. Notice that both kinds of parallelism can be combined together in the same simulation. Moreover, both threads and MPI processes, can be balanced using the builtin load balance system (RUPER-LB [30]) to maximise the performance on heterogeneous infrastructures. In addition, PenRed provides a modular structure with methods designed to easily extend its functionality. Thus, users can create their own independent modules to adapt our engine to their needs without changing the original modules. Furthermore, user extensions will take advantage of the builtin parallelism without any extra effort or knowledge of parallel programming. Additional comments including restrictions and unusual features: PenRed has been compiled in linux systems withg++ of GCC versions 4.8.5, 7.3.1, 8.3.1 and 9; clang version 3.4.2 and intel C++ compiler (icc) version 19.0.5.281. Since it is a C++11-standard compliant code, PenRed should be able to compile with any compiler with C++11 support. In addition, if the code is compiled without MPI support, it does not require any non standard library. To enable MPI capabilities, the user needs to install whatever available MPI implementation, such as openMPI [24] or mpich [25], which can be found in the repositories of any linux distribution. Finally, to provide DICOM processing support, PenRed can be optionally compiled using the dicom toolkit (dcmtk) [32] library. Thus, PenRed has only two optional dependencies, an MPI implementation and the dcmtk library.
|
|
|
|
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).
|
|
