Biomedical subjects

D Baltas

Publications and source records attributed to D Baltas.

At least 19 recordsLinked to original sources

PubMed2006

Polymer gel dosimetry for the TG-43 dosimetric characterization of a new 125I interstitial brachytherapy seed.

In this work, a polymer gel-magnetic resonance (MR) imaging method is employed for the dosimetric characterization of a new 125I low dose rate seed (IsoSeed model I25.S17). Two vials filled with PABIG gel were prepared in-house and one new seed as well as one commercially available 125I seed of similar dose rate and well-known dosimetric parameters (IsoSeed model I25.S06) were positioned in each vial. Both seeds in each vial were MR scanned simultaneously on days 11 and 26 after implantation. The data obtained from the known seed in each vial are used to calibrate the gel dose response which, for the prolonged irradiation duration necessitated by the investigated dose rates, depends on the overall irradiation time. Data for this study are presented according to the AAPM TG-43 dosimetric formalism. Polymer gel results concerning the new seed are compared to corresponding, published dosimetric results obtained, for the purpose of the new seed clinical implementation, by our group using the established methods of Monte Carlo (MC) simulation and thermo-luminescence dosimetry (TLD). Polymer gel dosimetry yields an average dose rate constant value of lambda = (0.921 +/- 0.031) cGy h(-1) U(-1) relative to (MC)lambda = (0.929 +/- 0.014) cGy h(-1) U(-1), (TLD)lambda = (0.951 +/- 0.044) cGy h(-1) U(-1) and the average value of Lambda = (0.940 +/- 0.051) cGy h(-1) U(-1) proposed for the clinical implementation of the new seed. Results for radial dose function, g(L)(r), and anisotropy function, F(r, theta), also agree with corresponding MC calculations within experimental uncertainties which are smaller for the polymer gel method compared to TLD. It is concluded that the proposed polymer gel-magnetic resonance imaging methodology could be used at least as a supplement to the established techniques for the dosimetric characterization of new low energy and low dose rate interstitial brachytherapy seeds.

Brachytherapy↗

PubMed2006

Comparison of radiation shielding requirements for HDR brachytherapy using 169Yb and 192Ir sources.

169Yb has received a renewed focus lately as an alternative to 192Ir sources for high dose rate (HDR) brachytherapy. Following the results of a recent work by our group which proved 169Yb to be a good candidate for HDR prostate brachytherapy, this work seeks to quantify the radiation shielding requirements for 169Yb HDR brachytherapy applications in comparison to the corresponding requirements for the current 192Ir HDR brachytherapy standard. Monte Carlo simulation (MC) is used to obtain 169Yb and 192Ir broad beam transmission data through lead and concrete. Results are fitted to an analytical equation which can be used to readily calculate the barrier thickness required to achieve a given dose rate reduction. Shielding requirements for a HDR brachytherapy treatment room facility are presented as a function of distance, occupancy, dose limit, and facility workload, using analytical calculations for both 169Yb and 192Ir HDR sources. The barrier thickness required for 169Yb is lower than that for 192Ir by a factor of 4-5 for lead and 1.5-2 for concrete. Regarding 169Yb HDR brachytherapy applications, the lead shielding requirements do not exceed 15 mm, even in highly conservative case scenarios. This allows for the construction of a lead door in most cases, thus avoiding the construction of a space consuming, specially designed maze. The effects of source structure, attenuation by the patient, and scatter conditions within an actual treatment room on the above-noted findings are also discussed using corresponding MC simulation results.

Brachytherapy↗

PubMed2006

Radiological evaluation of the metal-bone interface of a porous tantalum monoblock acetabular component.

Between January 1998 and December 1998, 82 consecutive patients (86 hips) underwent total hip arthroplasty using a trabecular metal monoblock acetabular component. All patients had a clinical and radiological follow-up evaluation at six, 12 and 24 weeks, 12 months, and then annually thereafter. On the initial post-operative radiograph 25 hips had a gap between the outer surface of the component and the acetabular host bed which ranged from 1 to 5 mm. All patients were followed up clinically and radiologically for a mean of 7.3 years (7 to 7.5). The 25 hips with the 1 to 5 mm gaps were studied for component migration at two years using the Einzel-Bild-Roentgen-Analyse (EBRA) digital measurement method. At 24 weeks all the post-operative gaps were filled with bone and no acetabular component had migrated. The radiographic outcome of all 86 components showed no radiolucent lines and no evidence of lysis. No acetabular implant was revised. There were no dislocations or other complications. The bridging of the interface gaps (up to 5 mm) by the trabecular metal monoblock acetabular component indicates the strong osteoconductive, and possibly osteoinductive, properties of trabecular metal.

Acetabulum↗

PubMed2005

Monte Carlo and thermoluminescence dosimetry of the new IsoSeed model I25.S17 125I interstitial brachytherapy seed.

Monte Carlo simulation and experimental thermoluminescence dosimetry were utilized for the dosimetric characterization of the new IsoSeed model I25.S17 125I interstitial brachytherapy seed. The new seed design is similar to that of the selectSeed and 6711 seeds, with the exception of its molybdenum marker. Full dosimetric data are presented following the recommendations in the Update of the AAPM Task Group 43 report (TG-43U1). A difference of 3.3% was found between Monte Carlo dose rate constant results calculated by air kerma strengths from simulations using a point detector and a detector resembling the solid angle subtended to the seed by the Wide Angle Free Air Chamber (WAFAC) in the primary standard calibration geometry. Following the TG-43U1 recommendations, an average value of lambdaMC = (0.929 +/- 0.014) cGy h(-1) U(-1) was adopted for the new seed. This value was then averaged with the measured value of lambdaEXP = (0.951 +/- 0.044) cGy h(-1) U(-1) to yield the proposed dose rate constant for the new seed that is equal to lambda = (0.940 +/- 0.051) cGy h(-1) U(-1). The Monte Carlo calculated radial dose function and two-dimensional (2-D) anisotropy function results for the new seed were found in agreement with experimental results to within statistical uncertainty of repeated measurements. Monte Carlo simulations were also performed for 125I seeds of similar geometry and dimensions for the purpose of comparison. The new seed presents dosimetric characteristics that are very similar to that of the selectSeed. In comparison to the most extensively studied Amersham 6711 seed, the new one presents similar dosimetric characteristics with a slightly reduced dose rate constant (1.5%).

Anisotropy↗

PubMed2005

A dosimetric comparison of 169Yb versus 192Ir for HDR prostate brachytherapy.

For the purpose of evaluating the use of 169Yb for prostate High Dose Rate brachytherapy (HDR), a hypothetical 169Yb source is assumed with the exact same design of the new microSelectron source replacing the 192Ir active core by pure 169Yb metal. Monte Carlo simulation is employed for the full dosimetric characterization of both sources and results are compared following the AAPM TG-43 dosimetric formalism. Monte Carlo calculated dosimetry results are incorporated in a commercially available treatment planning system (SWIFT), which features an inverse treatment planning option based on a multiobjective dose optimization engine. The quality of prostate HDR brachytherapy using the real 192Ir and hypothetical 169Yb source is compared in a comprehensive analysis of different prostate implants in terms of the multiobjective dose optimization solutions as well as treatment quality indices such as Dose Volume Histograms (DVH) and the Conformal Index (COIN). Given that scattering overcompensates for absorption in intermediate photon energies and distances in the range of interest to prostate HDR brachytherapy, 169Yb proves at least equivalent to 192Ir irrespective of prostate volume. This has to be evaluated in view of the shielding requirements for the 169Yb energies that are minimal relative to that for 192Ir.

Biophysical Phenomena↗

PubMed2004

The effect of patient inhomogeneities in oesophageal 192Ir HDR brachytherapy: a Monte Carlo and analytical dosimetry study.

The effect of patient inhomogeneities surrounding the oesophagus on the dosimetry planning of an upper thoracic oesophageal 192Ir HDR brachytherapy treatment is studied. The MCNPX Monte Carlo code is used for dosimetry in a patient-equivalent phantom geometry and results are compared in terms of isodose contours as well as dose volume histograms with corresponding calculations by a contemporary treatment planning system software featuring a full TG-43 dose calculation algorithm (PLATO BPS version 14.2.4). It is found that the presence of patient inhomogeneities does not alter the delivery of the planned dose distribution to the planning treatment volume. Regarding the organs at risk, the common practice of current treatment planning systems (TPSs) to consider the patient geometry as a homogeneous water medium leads to a dose overestimation of up to 13% to the spinal cord and an underestimation of up to 15% to the sternum bone. These findings which correspond to the dose region of about 5-10% of the prescribed dose could only be of significance when brachytherapy is used as a boost to external beam therapy. Additionally, an analytical dosimetry model, which is efficient in calculating dose in mathematical phantoms containing inhomogeneity shells of materials of radiobiological interest, is utilized for dosimetry in the patient-equivalent inhomogeneous phantom geometry. Analytical calculations in this work are in good agreement with corresponding Monte Carlo results within the bone inhomogeneities of spinal cord and sternum bone but, like treatment planning system calculations, the model fails to predict the dose distribution in the proximal lung surface as well as within the lungs just as the TPS does, due to its inherent limitation in treating lateral scatter and backscatter radiation.

Algorithms↗

PubMed2004

Evaluation of a TG-43 compliant analytical dosimetry model in clinical 192Ir HDR brachytherapy treatment planning and assessment of the significance of source position and catheter reconstruction uncertainties.

A simple, time efficient, analytical model incorporating heterogeneities and body dimensions around a point 192Ir source is generalized for accurate dosimetry around commercially available 192Ir brachytherapy sources. The generalized model was verified in dosimetry of a clinical 192Ir high dose rate prostate monotherapy application, involving 16 catheters and 83 source dwell positions, through comparison with corresponding treatment planning system data. The computational time efficiency and accuracy of the proposed model allowed the assessment of the impact that uncertainties in source dwell positions and catheter reconstruction may have on dose distributions, and how these could potentially affect the clinical outcome. Results revealed that a 0.1 cm catheter reconstruction uncertainty and a 0.15 cm source position uncertainty along the catheter lead to a dose uncertainty of less than 2% for doses lower than 200% of the prescribed dose, reaching up to 5% for points lying in close proximity to the catheters. These uncertainties were found to have no impact (less than 1%) on dose volume histogram results of both the planning target volume and the urethra. A catheter reconstruction uncertainty as high as 0.2 cm results in a dose uncertainty greater than 2%, reaching up to 9%, only for points inside the 150% contour. However, even in this case, the impact on dose volume histogram calculations is less than 3%.

Brachytherapy↗

PubMed2004

A monte carlo dosimetry study of vaginal 192Ir brachytherapy applications with a shielded cylindrical applicator set.

A durable recommendation for brachytherapy treatment planning systems to account for the effect of tissue, applicator and shielding material heterogeneities exists. As different proposed approaches have not been integrated in clinical treatment planning routine yet, currently utilized systems disregard or, most commonly, do not fully account for the aforementioned effects. Therefore, it is of interest to evaluate the efficacy of current treatment planning in clinical applications susceptible to errors due to heterogeneities. In this work the effect of the internal structure as well as the shielding used with a commercially available cylindrical shielded applicator set (Nucletron part # 084.320) for vaginal and rectum treatments is studied using three-dimensional Monte Carlo simulation for a clinical treatment plan involving seven source dwell positions of the classic microSelectron HDR 192Ir source. Results are compared to calculations of a treatment planning system (Plato BPS v.14.2.7), which assumes homogeneous water medium and applies a constant, multiplicative transmission factor only at points lying in the shadow of the shield. It is found that the internal structure of the applicator (which includes stainless steel, air and plastic materials) with no shield loaded does not affect the dose distribution relative to homogeneous water. In the unshielded side of the applicator with a 90 degrees, 180 degrees, or 270 degrees tungsten alloy shield loaded, an overestimation of treatment planning system calculations relative to Monte Carlo results was observed which is both shield and position dependent. While significant (up to 15%) at increased distances, which are not of major clinical importance, this overestimation does not affect dose prescription distances by more than 3%. The inverse effect of approx. 3% dose increase at dose prescription distances is observed for stainless steel shields. Regarding the shielded side of the applicator, it is shown that the default treatment planning system transmission factors for tungsten alloy result in a consistent dose over-estimation thus constituting a safe approach given the nature of associated clinical applications. Stainless steel is shown to be an ineffective shielding material with transmission factors reaching up to 0.68 at increased distances irrespective of shield geometry.

Body Burden↗

PubMed2003

In vivo thermoluminescence dosimetry dose verification of transperineal 192Ir high-dose-rate brachytherapy using CT-based planning for the treatment of prostate cancer.

PURPOSE: To evaluate the potential of in vivo thermoluminescence dosimetry to estimate the accuracy of dose delivery in conformal high-dose-rate brachytherapy of prostate cancer. METHODS AND MATERIALS: A total of 50 LiF, TLD-100 cylindrical rods were calibrated in the dose range of interest and used as a batch for all fractions. Fourteen dosimeters for every treatment fraction were loaded in a plastic 4F catheter that was fixed in either one of the 6F needles implanted for treatment purposes or in an extra needle implanted after consulting with the patient. The 6F needles were placed either close to the urethra or in the vicinity of the median posterior wall of the prostate. Initial results are presented for 18 treatment fractions in 5 patients and compared to corresponding data calculated using the commercial treatment planning system used for the planning of the treatments based on CT images acquired postimplantation. RESULTS: The maximum observed mean difference between planned and delivered dose within a single treatment fraction was 8.57% +/- 2.61% (root mean square [RMS] errors from 4.03% to 9.73%). Corresponding values obtained after averaging results over all fractions of a patient were 6.88% +/- 4.93% (RMS errors from 4.82% to 7.32%). Experimental results of each fraction corresponding to the same patient point were found to agree within experimental uncertainties. CONCLUSIONS: Experimental results indicate that the proposed method is feasible for dose verification purposes and suggest that dose delivery in transperineal high-dose-rate brachytherapy after CT-based planning can be of acceptable accuracy.

Brachytherapy↗

PubMed2003

An analytical dosimetry model as a step towards accounting for inhomogeneities and bounded geometries in 192Ir brachytherapy treatment planning.

A simple analytical dose rate calculation model based on primary and scatter separation that treats 192Ir as a monoenergetic source by use of appropriate attenuation and mass energy absorption coefficients is documented for accurate dosimetry in water. This model is then generalized and tested for use in any homogeneous tissue material of radiobiological interest using scatter to primary ratios calculated in water with a material density scaling to account for the difference in the scattering properties of these materials and water. The potential of the analytical model for predicting the effect of the interference of an inhomogeneity is then evaluated by comparison with corresponding Monte Carlo calculations. It is found that regardless of the inhomogeneity dimensions and position relative to the source, the model is capable of increased accuracy (better than 2%) in calculating the primary dose rate at any point not only for low-Z tissue materials but also for high-Z shielding materials where a severe hardening of the primary photons occurs. Overall, for low-Z tissue inhomogeneities the proposed model succeeds in correcting dosimetry results towards the right direction compared to commercial treatment planning systems that currently ignore the effect of phantom dimensions and inhomogeneity interference. Regarding high-Z shielding materials the proposed model accurately predicts the dose reduction just beyond the inhomogeneity (for example it predicts a dose reduction of 47% just behind a tungsten alloy cylinder of 1 cm diameter and 2 mm thickness placed at 1.4 cm away from an 192Ir source, in agreement with corresponding results in the literature) but does not account for the increasing contribution of the laterally scattered photons with increasing distance from the bounded inhomogeneity.

Algorithms↗

PubMed2003

Global convergence analysis of fast multiobjective gradient-based dose optimization algorithms for high-dose-rate brachytherapy.

We consider the problem of the global convergence of gradient-based optimization algorithms for interstitial high-dose-rate (HDR) brachytherapy dose optimization using variance-based objectives. Possible local minima could lead to only sub-optimal solutions. We perform a configuration space analysis using a representative set of the entire non-dominated solution space. A set of three prostate implants is used in this study. We compare the results obtained by conjugate gradient algorithms, two variable metric algorithms and fast-simulated annealing. For the variable metric algorithm BFGS from numerical recipes, large fluctuations are observed. The limited memory L-BFGS algorithm and the conjugate gradient algorithm FRPR are globally convergent. Local minima or degenerate states are not observed. We study the possibility of obtaining a representative set of non-dominated solutions using optimal solution rearrangement and a warm start mechanism. For the surface and volume dose variance and their derivatives, a method is proposed which significantly reduces the number of required operations. The optimization time, ignoring a preprocessing step, is independent of the number of sampling points in the planning target volume. Multiobjective dose optimization in HDR brachytherapy using L-BFGS and a new modified computation method for the objectives and derivatives has been accelerated, depending on the number of sampling points, by a factor in the range 10-100.

Algorithms↗

PubMed2003

A hybrid evolutionary algorithm for multi-objective anatomy-based dose optimization in high-dose-rate brachytherapy.

Multiple objectives must be considered in anatomy-based dose optimization for high-dose-rate brachytherapy and a large number of parameters must be optimized to satisfy often competing objectives. For objectives expressed solely in terms of dose variances, deterministic gradient-based algorithms can be applied and a weighted sum approach is able to produce a representative set of non-dominated solutions. As the number of objectives increases, or non-convex objectives are used, local minima can be present and deterministic or stochastic algorithms such as simulated annealing either cannot be used or are not efficient. In this case we employ a modified hybrid version of the multi-objective optimization algorithm NSGA-II. This, in combination with the deterministic optimization algorithm, produces a representative sample of the Pareto set. This algorithm can be used with any kind of objectives, including non-convex, and does not require artificial importance factors. A representation of the trade-off surface can be obtained with more than 1000 non-dominated solutions in 2-5 min. An analysis of the solutions provides information on the possibilities available using these objectives. Simple decision making tools allow the selection of a solution that provides a best fit for the clinical goals. We show an example with a prostate implant and compare results obtained by variance and dose-volume histogram (DVH) based objectives.

Algorithms↗

PubMed2003

Monte Carlo dosimetry of a new 192Ir pulsed dose rate brachytherapy source.

A new microSelectron pulsed dose rate source has been designed, containing two active pellets instead of one inactive and one active pellet contained in the old design, to facilitate the incorporation of higher activity up to 74 GBq (2 Ci). In this work, Monte Carlo simulation is used to derive full dosimetric data following the AAPM TG-43 formalism, as well as the dose rate per unit air kerma strength data in Cartesian, "away and along" coordinates for both source designs. The calculated dose rate constant of the new PDR source design was found equal to lambda=(1.121 +/- 0.006) cGy h(-1) U(-1) compared to lambda = (1.124 +/- 0.006) cGy h(-1) U(-1) for the old design. Radial dose functions of the two sources calculated using the point source approximated geometry factors were found in close agreement (within 1%) except for radial distances under 2 mm. At polar angles close to the longitudinal source axis at the sources' distal end, the new design presents increased anisotropy (up to 10%) compared to the old one due to its longer active core. At polar angles close to the longitudinal source axis at the sources' drive wire end however, the old design presents increased anisotropy (up to 18%) due to attenuation of emitted photons through the inactive Ir pellet. These differences, also present in "away and along" dose rate results, necessitate the replacement of treatment planning input data for the new microSelectron pulsed dose rate source.

Anisotropy↗

PubMed2003

Are dose calculations during dose optimization in brachytherapy necessary?

We consider the problem of anatomy based dose optimization in brachytherapy. A calculation method for some objective functions and their derivatives is proposed which significantly reduces the number of required operations. The optimization in some cases, ignoring a preprocessing step, is independent of the number of sampling points. The idea is that some of the objectives and their derivatives used for dose optimization do not require the explicit calculation of dose values. Dose optimization with the new modified computation method for the objectives and derivatives is, depending on the number of sampling points, up to 100 times faster than the conventional method with dose calculation.

Algorithms↗

PubMed2003

3D dose verification in 192Ir HDR prostate monotherapy using polymer gels and MRI.

VIPAR polymer gels and 3D MRI techniques were evaluated for their ability to provide experimental verification of 3D dose distributions in a simulation of a 192Ir prostate monotherapy clinical application. A real clinical treatment plan was utilized, generated by post-irradiation, CT based calculations derived from Plato BPS and Swift treatment planning systems. The simulated treatment plan involved the use of 10 catheters and 39 source positions within a glass vessel of appropriate dimensions, homogeneously filled with the VIPAR gel. 3D high resolution MR scanning of the gel produced T2 relaxation time maps, from which 3D dose distributions were derived via an appropriate calibration procedure. Results were compared to corresponding dose distributions obtained from the Plato and Swift treatment planning systems. Quantitative comparison, on a point by point basis, was based on user adopted acceptance criteria of 5% dose-difference and 3 mm distance-to-agreement. Significant deviations between experimental and calculated dose distributions were found for doses lower than 50% due to the reduced dose resolution of the method in the low dose, low dose gradient region. Measurement errors were observed at 1.0-1.5 mm around each catheter due to MR imaging susceptibility artifacts. For most remaining points the acceptance criteria were fulfilled. Systematic offsets of the order of 1-2 mm, observed between measured and corresponding calculated isocontours at specific segments, are attributed to the 1 mm uncertainty in catheter reconstruction and 1 mm uncertainty in the alignment of the MR and CT imaging planes.

Algorithms↗

PubMed2002

Multiobjective anatomy-based dose optimization for HDR-brachytherapy with constraint free deterministic algorithms.

In high dose rate (HDR) brachytherapy, conventional dose optimization algorithms consider multiple objectives in the form of an aggregate function that transforms the multiobjective problem into a single-objective problem. As a result, there is a loss of information on the available alternative possible solutions. This method assumes that the treatment planner exactly understands the correlation between competing objectives and knows the physical constraints. This knowledge is provided by the Pareto trade-off set obtained by single-objective optimization algorithms with a repeated optimization with different importance vectors. A mapping technique avoids non-feasible solutions with negative dwell weights and allows the use of constraint free gradient-based deterministic algorithms. We compare various such algorithms and methods which could improve their performance. This finally allows us to generate a large number of solutions in a few minutes. We use objectives expressed in terms of dose variances obtained from a few hundred sampling points in the planning target volume (PTV) and in organs at risk (OAR). We compare two- to four-dimensional Pareto fronts obtained with the deterministic algorithms and with a fast-simulated annealing algorithm. For PTV-based objectives, due to the convex objective functions, the obtained solutions are global optimal. If OARs are included, then the solutions found are also global optimal, although local minima may be present as suggested.

Algorithms↗

PubMed2002

Brachytherapy dose-volume histogram computations using optimized stratified sampling methods.

A stratified sampling method for the efficient repeated computation of dose-volume histograms (DVHs) in brachytherapy is presented as used for anatomy based brachytherapy optimization methods. The aim of the method is to reduce the number of sampling points required for the calculation of DVHs for the body and the PTV. From the DVHs are derived the quantities such as Conformity Index COIN and COIN integrals. This is achieved by using partial uniform distributed sampling points with a density in each region obtained from a survey of the gradients or the variance of the dose distribution in these regions. The shape of the sampling regions is adapted to the patient anatomy and the shape and size of the implant. For the application of this method a single preprocessing step is necessary which requires only a few seconds. Ten clinical implants were used to study the appropriate number of sampling points, given a required accuracy for quantities such as cumulative DVHs, COIN indices and COIN integrals. We found that DVHs of very large tissue volumes surrounding the PTV, and also COIN distributions, can be obtained using a factor of 5-10 times smaller the number of sampling points in comparison with uniform distributed points.

Brachytherapy↗

PubMed2002

Thermoluminescent dosimetry of the selectseed 125I interstitial brachytherapy seed.

This work presents experimental dosimetry results for the new selectSeed 125I prostate seed design for use with the seedSelectron afterloading device, in accordance with the AAPM advisory that all new low energy interstitial brachytherapy seeds should undergo one Monte Carlo (MC) and at least one experimental dosimetry characterization. TLD dosimetry was performed using 120 cylindrical LiF TLD type-100 rods calibrated using a 6 MV photon beam. They were irradiated in solid water phantoms for the experimental determination of the seed dose rate constant, radial dose functions and anisotropy functions. MC simulations were performed for the determination of the TLDs relative energy response that was found position independent and equal to 1.40+/-0.03, and for the calculation of the ratio of dose in liquid water to dose in solid water that was found to be well described by Dliquidwater/Dsolidwater= 1.013*r+0.030 presenting only a minor dependence on polar angle. The selectSeed dose rate constant in liquid water was found equal to 0.938+/-0.065 cGy h(-1) U(-1), which agrees within experimental uncertainties with corresponding MC results of lambdaselect Seed=0.954+/-0.005 cGy h(-1) U(-1). The experimental radial dose and anisotropy function results were also found in good agreement with corresponding MC calculations.

Anisotropy↗