|
|
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.
|
|
|
|
Ballester, F., Granero, D., Perez-Calatayud, J., Venselaar, J. L. M., & Rivard, M. J. (2010). Study of encapsulated Tm-170 sources for their potential use in brachytherapy. Med. Phys., 37(4), 1629–1637.
Abstract: Methods: The authors have assumed a theoretical Tm-170 cylindrical source encapsulated with stainless steel and typical dimensions taken from the currently available HDR Ir-192 brachytherapy sources. The dose-rate distribution was calculated for this source using the GEANT4 Monte Carlo (MC) code considering both photon and electron Tm-170 spectra. The AAPM TG-43 U1 brachytherapy dosimetry parameters were derived. To study general properties of Tm-170 encapsulated sources, spherical sources encapsulated with stainless steel and platinum were also studied. Moreover, the influence of small variations in the active core and capsule dimensions on the dosimetric characteristics was assessed. Treatment times required for a Tm-170 source were compared to those for Ir-192 and Yb-169 for the same contained activity. Results: Due to the energetic beta spectrum and the large electron yield, the bremsstrahlung contribution to the dose was of the same order of magnitude as from the emitted gammas and characteristic x rays. Moreover, the electron spectrum contribution to the dose was significant up to 4 mm from the source center compared to the photon contribution. The dose-rate constant Lambda of the cylindrical source was 1.23 cGy h(-1) U-1. The behavior of the radial dose function showed promise for applications in brachytherapy. Due to the electron spectrum, the anisotropy was large for r < 6 mm. Variations in manufacturing tolerances did not significantly influence the final dosimetry data when expressed in cGy h(-1) U-1. For typical capsule dimensions, maximum reference dose rates of about 0.2, 10, and 2 Gy min(-1) would then be obtained for Tm-170, Ir-192, and Yb-169, respectively, resulting in treatment times greater than those for HDR Ir-192 brachytherapy. Conclusions: The dosimetric characteristics of source designs exploiting the low photon energy of Tm-170 were studied for potential application in HDR-brachytherapy. Dose-rate distributions were obtained for cylindrical and simplified spherical Tm-170 source designs (stainless steel and platinum capsule materials) using MC calculations. Despite the high activity of Tm-170, calculated treatment times were much longer than for Ir-192.
|
|
|
|
Rivard, M. J., Granero, D., Perez-Calatayud, J., & Ballester, F. (2010). Influence of photon energy spectra from brachytherapy sources on Monte Carlo simulations of kerma and dose rates in water and air. Med. Phys., 37(2), 869–876.
Abstract: Methods: For Ir-192, I-125, and Pd-103, the authors considered from two to five published spectra. Spherical sources approximating common brachytherapy sources were assessed. Kerma and dose results from GEANT4, MCNP5, and PENELOPE-2008 were compared for water and air. The dosimetric influence of Ir-192, I-125, and Pd-103 spectral choice was determined. Results: For the spectra considered, there were no statistically significant differences between kerma or dose results based on Monte Carlo code choice when using the same spectrum. Water-kerma differences of about 2%, 2%, and 0.7% were observed due to spectrum choice for Ir-192, I-125, and Pd-103, respectively (independent of radial distance), when accounting for photon yield per Bq. Similar differences were observed for air-kerma rate. However, their ratio (as used in the dose-rate constant) did not significantly change when the various photon spectra were selected because the differences compensated each other when dividing dose rate by air-kerma strength. Conclusions: Given the standardization of radionuclide data available from the National Nuclear Data Center (NNDC) and the rigorous infrastructure for performing and maintaining the data set evaluations, NNDC spectra are suggested for brachytherapy simulations in medical physics applications.
|
|
|
|
Granero, D., Vijande, J., Ballester, F., & Rivard, M. J. (2011). Dosimetry revisited for the HDR Ir-192 brachytherapy source model mHDR-v2. Med. Phys., 38(1), 487–494.
Abstract: Purpose: Recently, the manufacturer of the HDR Ir-192 mHDR-v2 brachytherapy source reported small design changes (referred to herein as mHDR-v2r) that are within the manufacturing tolerances but may alter the existing dosimetric data for this source. This study aimed to (1) check whether these changes affect the existing dosimetric data published for this source; (2) obtain new dosimetric data in close proximity to the source, including the contributions from 192Ir electrons and considering the absence of electronic equilibrium; and (3) obtain scatter dose components for collapsed cone treatment planning system implementation. Methods: Three different Monte Carlo (MC) radiation transport codes were used: MCNP5, PENELOPE2008, and GEANT4. The source was centrally positioned in a 40 cm radius water phantom. Absorbed dose and collision kerma were obtained using 0.1 mm (0.5 mm) thick voxels to provide high-resolution dosimetry near (far from) the source. Dose-rate distributions obtained with the three MC codes were compared. Results: Simulations of mHDR-v2 and mHDR-v2r designs performed with three radiation transport codes showed agreement typically within 0.2% for r >= 0.25 cm. Dosimetric contributions from source electrons were significant for r<0.25 cm. The dose-rate constant and radial dose function were similar to those from previous MC studies of the mHDR-v2 design. The 2D anisotropy function also coincided with that of the mHDR-v2 design for r >= 0.25 cm. Detailed results of dose distributions and scatter components are presented for the modified source design. Conclusions: Comparison of these results to prior MC studies showed agreement typically within 0.5% for r >= 0.25 cm. If dosimetric data for r<0.25 cm are not needed, dosimetric results from the prior MC studies will be adequate. c 2011 American Association of Physicists in Medicine.
|
|
|
|
Granero, D., Perez-Calatayud, J., Vijande, J., Ballester, F., & Rivard, M. J. (2014). Limitations of the TG-43 formalism for skin high-dose-rate brachytherapy dose calculations. Med. Phys., 41(2), 021703–8pp.
Abstract: Purpose: In skin high-dose-rate (HDR) brachytherapy, sources are located outside, in contact with, or implanted at some depth below the skin surface. Most treatment planning systems use the TG-43 formalism, which is based on single-source dose superposition within an infinite water medium without accounting for the true geometry in which conditions for scattered radiation are altered by the presence of air. The purpose of this study is to evaluate the dosimetric limitations of the TG-43 formalism in HDR skin brachytherapy and the potential clinical impact. Methods: Dose rate distributions of typical configurations used in skin brachytherapy were obtained: a 5 cm x 5 cm superficial mould; a source inside a catheter located at the skin surface with and without backscatter bolus; and a typical interstitial implant consisting of an HDR source in a catheter located at a depth of 0.5 cm. Commercially available HDR Co-60 and Ir-192 sources and a hypothetical Yb-169 source were considered. The Geant4Monte Carlo radiation transport code was used to estimate dose rate distributions for the configurations considered. These results were then compared to those obtained with the TG-43 dose calculation formalism. In particular, the influence of adding bolus material over the implant was studied. Results: For a 5 cm x 5 cm Ir-192 superficial mould and 0.5 cm prescription depth, dose differences in comparison to the TG-43 method were about -3%. When the source was positioned at the skin surface, dose differences were smaller than -1% for Co-60 and Ir-192, yet -3% for Yb-169. For the interstitial implant, dose differences at the skin surface were -7% for Co-60, -0.6% for Ir-192, and -2.5% for Yb-169. Conclusions: This study indicates the following: (i) for the superficial mould, no bolus is needed; (ii) when the source is in contact with the skin surface, no bolus is needed for either Co-60 and Ir-192. For lower energy radionuclides like Yb-169, bolus may be needed; and (iii) for the interstitial case, at least a 0.1 cm bolus is advised for Co-60 to avoid underdosing superficial target layers. For Ir-192 and Yb-169, no bolus is needed. For those cases where no bolus is needed, its use might be detrimental as the lack of radiation scatter may be beneficial to the patient, although the 2% tolerance for dose calculation accuracy recommended in the AAPM TG-56 report is not fulfilled.
|
|
|
|
Ballester, F., Tedgren, A. C., Granero, D., Haworth, A., Mourtada, F., Fonseca, G. P., et al. (2015). A generic high-dose rate Ir-192 brachytherapy source for evaluation of model-based dose calculations beyond the TG-43 formalism. Med. Phys., 42(6), 3048–3062.
Abstract: Purpose: In order to facilitate a smooth transition for brachytherapy dose calculations from the American Association of Physicists in Medicine (AAPM) Task Group No. 43 (TG-43) formalism to model-based dose calculation algorithms (MBDCAs), treatment planning systems (TPSs) using a MBDCA require a set of well-defined test case plans characterized by Monte Carlo (MC) methods. This also permits direct dose comparison to TG-43 reference data. Such test case plans should be made available for use in the software commissioning process performed by clinical end users. To this end, a hypothetical, generic high-dose rate (HDR) Ir-192 source and a virtual water phantom were designed, which can be imported into a TPS. Methods: A hypothetical, generic HDR Ir-192 source was designed based on commercially available sources as well as a virtual, cubic water phantom that can be imported into any TPS in DICOM format. The dose distribution of the generic Ir-192 source when placed at the center of the cubic phantom, and away from the center under altered scatter conditions, was evaluated using two commercial MBDCAs [Oncentra (R) Brachy with advanced collapsed-cone engine (ACE) and BrachyVision AcuRos (TM)]. Dose comparisons were performed using state-of-the-art MC codes for radiation transport, including ALGEBRA, BrachyDose, GEANT4, MCNP5, MCNP6, and pENELopE2008. The methodologies adhered to recommendations in the AAPM TG-229 report on high-energy brachytherapy source dosimetry. TG-43 dosimetry parameters, an along-away dose-rate table, and primary and scatter separated (PSS) data were obtained. The virtual water phantom of (201)(3) voxels (1 mm sides) was used to evaluate the calculated dose distributions. Two test case plans involving a single position of the generic HDR Ir-192 source in this phantom were prepared: (i) source centered in the phantom and (ii) source displaced 7 cm laterally from the center. Datasets were independently produced by different investigators. MC results were then compared against dose calculated using TG-43 and MBDCA methods. Results: TG-43 and PSS datasets were generated for the generic source, the PSS data for use with the ACE algorithm. The dose-rate constant values obtained from seven MC simulations, performed independently using different codes, were in excellent agreement, yielding an average of 1.1109 +/- 0.0004 cGy/(h U) (k = 1, Type A uncertainty). MC calculated dose-rate distributions for the two plans were also found to be in excellent agreement, with differences within type A uncertainties. Differences between commercial MBDCA and MC results were test, position, and calculation parameter dependent. On average, however, these differences were within 1% for ACUROS and 2% for ACE at clinically relevant distances. Conclusions: A hypothetical, generic HDR Ir-192 source was designed and implemented in two commercially available TPSs employing different MBDCAs. Reference dose distributions for this source were benchmarked and used for the evaluation of MBDCA calculations employing a virtual, cubic water phantom in the form of a CT DICOM image series. The implementation of a generic source of identical design in all TPSs using MBDCAs is an important step toward supporting univocal commissioning procedures and direct comparisons between TPSs.
|
|
|
|
Candela-Juan, C., Vijande, J., Garcia-Martinez, T., Niatsetski, Y., Nauta, G., Schuurman, J., et al. (2015). Comparison and uncertainty evaluation of different calibration protocols and ionization chambers for low-energy surface brachytherapy dosimetry. Med. Phys., 42(8), 4954–4964.
Abstract: Purpose: A surface electronic brachytherapy (EBT) device is in fact an x-ray source collimated with specific applicators. Low-energy (<100 kVp) x-ray beam dosimetry faces several challenges that need to be addressed. A number of calibration protocols have been published for x-ray beam dosimetry. The media in which measurements are performed are the fundamental difference between them. The aim of this study was to evaluate the surface dose rate of a low-energy x-ray source with small field applicators using different calibration standards and different small-volume ionization chambers, comparing the values and uncertainties of each methodology. Methods: The surface dose rate of the EBT unit Esteya (Elekta Brachytherapy, The Netherlands), a 69.5 kVp x-ray source with applicators of 10, 15, 20, 25, and 30 mm diameter, was evaluated using the AAPM TG-61 (based on air kerma) and International Atomic Energy Agency (IAEA) TRS-398 (based on absorbed dose to water) dosimetry protocols for low-energy photon beams. A plane parallel T34013 ionization chamber (PTW Freiburg, Germany) calibrated in terms of both absorbed dose to water and air kerma was used to compare the two dosimetry protocols. Another PTW chamber of the same model was used to evaluate the reproducibility between these chambers. Measurements were also performed with two different Exradin A20 (Standard Imaging, Inc., Middleton, WI) chambers calibrated in terms of air kerma. Results: Differences between surface dose rates measured in air and in water using the T34013 chamber range from 1.6% to 3.3%. No field size dependence has been observed. Differences are below 3.7% when measurements with the A20 and the T34013 chambers calibrated in air are compared. Estimated uncertainty (with coverage factor k = 1) for the T34013 chamber calibrated in water is 2.2%-2.4%, whereas it increases to 2.5% and 2.7% for the A20 and T34013 chambers calibrated in air, respectively. The output factors, measured with the PTW chambers, differ by less than 1.1% for any applicator size when compared to the output factors that were measured with the A20 chamber. Conclusions: Measurements using both dosimetric protocols are consistent, once the overall uncertainties are considered. There is also consistency between measurements performed with both chambers calibrated in air. Both the T34013 and A20 chambers have negligible stem effect. Any x-ray surface brachytherapy system, including Esteya, can be characterized using either one of these calibration protocols and ionization chambers. Having less correction factors, lower uncertainty, and based on measurements, performed in closer to clinical conditions, the TRS-398 protocol seems to be the preferred option.
|
|
|
|
Candela-Juan, C., Niatsetski, Y., van der Laarse, R., Granero, D., Ballester, F., Perez-Calatayud, J., et al. (2016). Design and characterization of a new high-dose-rate brachytherapy Valencia applicator for larger skin lesions. Med. Phys., 43(4), 1639–1648.
Abstract: Purpose: The aims of this study were (i) to design a new high-dose-rate (HDR) brachytherapy applicator for treating surface lesions with planning target volumes larger than 3 cm in diameter and up to 5 cm in size, using the microSelectron-HDR or Flexitron afterloader (Elekta Brachytherapy) with a Ir-192 source; (ii) to calculate by means of the Monte Carlo (MC) method the dose distribution for the new applicator when it is placed against a water phantom; and (iii) to validate experimentally the dose distributions in water. Methods: The PENELOPE2008 MC code was used to optimize dwell positions and dwell times. Next, the dose distribution in a water phantom and the leakage dose distribution around the applicator were calculated. Finally, MC data were validated experimentally for a 192Ir mHDR-v2 source by measuring (i) dose distributions with radiochromic EBT3 films (ISP); (ii) percentage depth-dose (PDD) curve with the parallel-plate ionization chamber Advanced Markus (PTW); and (iii) absolute dose rate with EBT3 films and the PinPoint T31016 (PTW) ionization chamber. Results: The new applicator is made of tungsten alloy (Densimet) and consists of a set of interchangeable collimators. Three catheters are used to allocate the source at prefixed dwell positions with preset weights to produce a homogenous dose distribution at the typical prescription depth of 3 mm in water. The same plan is used for all available collimators. PDD, absolute dose rate per unit of air kerma strength, and off-axis profiles in a cylindrical water phantom are reported. These data can be used for treatment planning. Leakage around the applicator was also scored. The dose distributions, PDD, and absolute dose rate calculated agree within experimental uncertainties with the doses measured: differences of MC data with chamber measurements are up to 0.8% and with radiochromic films are up to 3.5%. Conclusions: The new applicator and the dosimetric data provided here will be a valuable tool in clinical practice, making treatment of large skin lesions simpler, faster, and safer. Also the dose to surrounding healthy tissues is minimal.
|
|
|
|
Granero, D., Candela-Juan, C., Vijande, J., Ballester, F., Perez-Calatayud, J., Jacob, D., et al. (2016). Technical Note: Dosimetry of Leipzig and Valencia applicators without the plastic cap. Med. Phys., 43(5), 2087–4pp.
Abstract: Purpose: High dose rate (HDR) brachytherapy for treatment of small skin lesions using the Leipzig and Valencia applicators is a widely used technique. These applicators are equipped with an attachable plastic cap to be placed during fraction delivery to ensure electronic equilibrium and to prevent secondary electrons from reaching the skin surface. The purpose of this study is to report on the dosimetric impact of the cap being absent during HDR fraction delivery, which has not been explored previously in the literature. Methods: GEANT4 Monte Carlo simulations (version 10.0) have been performed for the Leipzig and Valencia applicators with and without the plastic cap. In order to validate the Monte Carlo simulations, experimental measurements using radiochromic films have been done. Results: Dose absorbed within 1 mm of the skin surface increases by a factor of 1500% for the Leipzig applicators and of 180% for the Valencia applicators. Deeper than 1 mm, the overdosage flattens up to a 10% increase. Conclusions: Differences of treating with or without the plastic cap are significant. Users must check always that the plastic cap is in place before any treatment in order to avoid overdosage of the skin. Prior to skin HDR fraction delivery, the timeout checklist should include verification of the cap placement. (C) 2016 American Association of Physicists in Medicine.
|
|
|
|
Albiol, A., Corbi, A., & Albiol, F. (2017). Automatic intensity windowing of mammographic images based on a perceptual metric. Med. Phys., 44(4), 1369–1378.
Abstract: Purpose: Initial auto-adjustment of the window level WL and width WW applied to mammographic images. The proposed intensity windowing (IW) method is based on the maximization of the mutual information (MI) between a perceptual decomposition of the original 12-bit sources and their screen displayed 8-bit version. Besides zoom, color inversion and panning operations, IW is the most commonly performed task in daily screening and has a direct impact on diagnosis and the time involved in the process. Methods: The authors present a human visual system and perception-based algorithm named GRAIL (Gabor-relying adjustment of image levels). GRAIL initially measures a mammogram's quality based on the MI between the original instance and its Gabor-filtered derivations. From this point on, the algorithm performs an automatic intensity windowing process that outputs the WL/WW that best displays each mammogram for screening. GRAIL starts with the default, high contrast, wide dynamic range 12-bit data, and then maximizes the graphical information presented in ordinary 8-bit displays. Tests have been carried out with several mammogram databases. They comprise correlations and an ANOVA analysis with the manual IW levels established by a group of radiologists. A complete MATLAB implementation of GRAIL is available at . Results: Auto-leveled images show superior quality both perceptually and objectively compared to their full intensity range and compared to the application of other common methods like global contrast stretching (GCS). The correlations between the human determined intensity values and the ones estimated by our method surpass that of GCS. The ANOVA analysis with the upper intensity thresholds also reveals a similar outcome. GRAIL has also proven to specially perform better with images that contain micro-calcifications and/or foreign X-ray-opaque elements and with healthy BI-RADS A-type mammograms. It can also speed up the initial screening time by a mean of 4.5 s per image. Conclusions: A novel methodology is introduced that enables a quality-driven balancing of the WL/WW of mammographic images. This correction seeks the representation that maximizes the amount of graphical information contained in each image. The presented technique can contribute to the diagnosis and the overall efficiency of the breast screening session by suggesting, at the beginning, an optimal and customized windowing setting for each mammogram.
|
|
