Search PubMed⌕ Search

Biomedical subjects

D Baltas

Publications and source records attributed to D Baltas.

At least 19 recordsLinked to original sources

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↗

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↗

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↗

Dosimetry comparison of 192Ir sources.

192Ir sources besides being widely utilized in the field of conventional brachytherapy also find use in contemporary peripheral and coronal intravascular applications. In this study, the same Monte Carlo simulation code and input data were used to investigate differences between the dose rate distributions of the most commonly used 192Ir sources in the cm and mm distance range. Findings are discussed in view of differences in source and encapsulation dimensions as well as structural details. Results are presented in the AAPM TG-43 formalism, as generalized by AAPM TG-60, for five 192Ir HDR source designs as well as an LDR seed and an LDR wire source. Dose rate constants of the sources at r0 = 1 cm and r0 = 2 mm were found proportional to the corresponding geometry factors along the transverse source bisectors and an equation of the form lambda r0(cGyh(-1) U(-1)) = 1.12 x G(r0,90 degrees) provides results within clinical accuracy (less than 2%) for any 192Ir source. Radial dose functions do not depend significantly on source and encapsulation geometry and agree within 2% with that of a point 192Ir source. Anisotropy is of importance for accurate dosimetry at the cm distance range but it does not affect dose rate in the mm distance range significantly. At such short radial distances the source geometry factor defines the shape of isodose lines. Dose uniformity at given distances from the sources is strongly dependent on source dimensions as indicated by dose rate profiles in polar and Cartesian coordinates.

Anisotropy↗

A multiobjective gradient-based dose optimization algorithm for external beam conformal radiotherapy.

A multiobjective gradient-based algorithm has been developed for the purpose of dose distribution optimization in external beam conformal radiotherapy. This algorithm is based on the concept of gathering the values of all objectives into a single value. The weighting factors of the composite objective values are varied in different steps, allowing the reconstruction of the trade-off surfaces (three or more objectives) or curves (two objectives) which define the boundary between the feasible and non-feasible domain regions. The analysis of these curves allows the decision-maker to select the solution that best fits the clinical goals. In contrast to all the other algorithms, our method provides not a single solution but a sample of solutions representing all possible clinical importance factors (weights) for the objectives used. The application of this algorithm to two test cases shows that a correct selection for the importance factors to multiply the individual objectives in the global objective value is not trivial and that the location and shape of the boundary region between the feasible and non-feasible solution regions are case dependent. Provided that the individual objective functions are analytically differentiable and that the number of objectives is the range of two to three, the computation times are acceptable for clinical use. Furthermore, the optimization for a unique combination of importance factors within the aggregate objective function is performed in less than 1 min.

Algorithms↗

A new algorithm for autoreconstruction of catheters in computed tomography-based brachytherapy treatment planning.

This paper describes innovative software for automatic reconstruction, which we term autoreconstruction, of plastic and metallic brachytherapy catheters using computed tomography (CT) data. No such automatic facility has previously existed in any treatment planning software. The patient data consists of a set of post-implantation CT images with the catheters in situ in their final positions. This new software solves those difficulties which arise when the catheters are intersecting or when loop techniques are used. With the software algorithms, catheter reconstruction time is significantly reduced and accuracy is also improved when compared with that achieved using the classical manual method of CT-slice-by-CT-slice reconstruction.

Algorithms↗

Dosimetry of 192Ir wires for LDR interstitial brachytherapy following the AAPM TG-43 dosimetric formalism.

Implementation of the AAPM Task Group 43 dosimetric formalism for 192Ir wires used as interstitial sources in low dose-rate (LDR) brachytherapy applications is investigated. Geometry factors, dose-rate constant values, radial dose functions, and anisotropy functions, utilized in this formalism, were calculated for various lengths of all commercially available wire source designs by means of a well-established Monte Carlo simulation code and an improved modification of the Sievert integral method. Results are presented in the form of look up tables that allow interpolation for dose-rate calculations around all practically used wire lengths, with accuracy acceptable for clinical applications.

Anisotropy↗

Dosimetry close to an 192Ir HDR source using N-vinylpyrrolidone based polymer gels and magnetic resonance imaging.

In this work, the utilization of polymer gel-MRI dosimetry for measurements at distances relevant to clinical brachytherapy and intravascular applications [i.e., in the mm range, where steep three-dimensional (3-D) dose gradients exist] is investigated using N-vinylpyrrolidone-based gels. Transverse axis radial dose distributions, dose distributions parallel to the source axis, and 2-D dose distributions around the commonly used microSelectron 192Ir HDR source are measured for single source dwell position irradiations. Experimental results are found in good agreement with verified Monte Carlo calculations, even for distances less than 3 mm from the source. The effect of various MRI parameters, such as slice thickness, slice mispositioning, and in-plane resolution, on the accuracy of the method is also investigated. Possible limitations of the method are discussed, and its' overall potential in brachytherapy dosimetry is evaluated. Experimental 2-D dose distributions for an intravascular application following the Paris irradiation protocol are compared to corresponding commercial treatment planning system calculations. Results suggest that polymer gel-MRI dosimetry is capable of experimentally verifying dose distributions in relevant clinical intravascular applications.

Anisotropy↗

Monte Carlo dosimetry of the selectSeed 125I interstitial brachytherapy seed.

This work provides full dosimetric data for the new selectSeed 125I prostate seed source to be distributed by Nucletron B.V. The AAPM TG-43 dosimetric formalism and the new 1999 NIST air kerma strength calibration standard have been followed. Air kerma strength, dose rate constant, radial dose functions, anisotropy functions, and anisotropy factors were calculated using Monte Carlo simulation. Corresponding calculations were also performed for the commercially available 6711 seed source, which is of similar design, for reasons of comparison. The calculated dose rate constant of the selectSeed was 0.954+/-0.005 cGy h(-1) U(-1) compared to 0.953+/-0.005 cGy h(-1) U(-1) for the 6711 source design. The latter value for the 6711 source suggests that the correction factor proposed by NIST for conversion of dose rate constants to the new 1999 NIST calibration standard may be overestimated by 2-3%. Radial dose functions of the two sources were found in good agreement for radial distances up to 4 cm, the selectSeed being less penetrating at greater radial distances (approximately 4% at 10 cm). The selectSeed source presents similar anisotropy characteristics with the 6711 source design. For both source designs, a distance and polar angle dependent discontinuity of anisotropy function values was observed owing to the dose contribution of radioactivity distributed on the ends of the cylindrical source cores. Variation of dosimetric parameters with possible variation in radioactive silver halide coating thickness of the silver source core of the new source was also investigated.

Anisotropy↗

Correcting organ motion artifacts in x-ray CT systems based on tracking of motion phase by the spatial overlap correlator. II. Experimental study.

This paper presents the experimental part of an investigation on tracking and eliminating organ motion artifacts in x-ray CT cardiac applications with emphasis on imaging coronary calcification. The system methodology consists of a software implementation of the spatial overlap correlator (SSOC) concept in x-ray CT scanners to track the net amplitude and phase of organ motion during the CT data acquisition process. A coherent sinogram synthesis (CSS) method is then used to identify the repeated phases of a periodic organ motion from the information provided by the SSOC process and hence synthesize a new sinogram with no motion effects. Since the SSOC scheme is capable of tracking cardiac motion, it identifies also the projection points associated with minimum amplitude cardiac motion effects. These points are used to identify a 180 degrees plus the fan angle sinogram for image reconstruction. This leads to a retrospective gating (RG) scheme that is based on the output of the SSOC process. Performance comparison of the proposed methodology with the retrospective ECG gating using real data sets with phantoms and human patients provides a performance assessment of the merits of the proposed methods. Real results demonstrate that the new methodology eliminates the requirement for ECG gating. Moreover, the CSS and the new RG methods do not require breath holding and they can be implemented in x-ray CT scanners to image coronary calcification and the heart's ventricles.

Algorithms↗

Beta versus gamma dosimetry close to Ir-192 brachytherapy sources.

The relative importance of the dose rate component owing to the beta spectrum emitted by 192Ir brachytherapy sources at the short radial distances of interest in intravascular and endobronchial applications is investigated. Separate dosimetric calculations, using Monte Carlo simulations, were performed for the gamma and beta dose rate components of an 192Ir ideal point source as well as real 192Ir source designs used in clinical practice including wire and seed sources and both Nucletron and Varian, old and new, high dose rate (HDR) source designs. A significant dose rate enhancement due to the beta spectrum emitted by 192Ir, greater than 50% for radial distances r<2 mm, was observed for an ideal point source. For real source designs, however, the magnitude of this enhancement was found to depend strongly on the sources' geometric as well as compositional details of the active core and encapsulation. A detectable effect was found for the majority of the investigated sources at radial distances less than 1 mm, but overall findings suggest that the contribution of beta particles is not significant in 192Ir clinical intravascular applications that are currently carried out. However, since treatment of vessels with smaller diameters, in the future, may lead to the development of 192Ir sources and catheters of reduced diameters, the potential effect of the beta spectrum in terms of dose enhancement to tissues in close proximity to 192Ir sources should not be ignored.

Beta Particles↗

Reconstruction and navigation system for intraoperative brachytherapy using the flab technique for colorectal tumor bed irradiation.

PURPOSE: To present the development of a new navigation and reconstruction system based on an electromagnetic free-hand tracker and on CT imaging for treatment planning of intraoperative high-dose-rate brachytherapy (IORT-HDRB) in the sacral region. Our aim is to improve accuracy and to enable individualized treatment planning and dose documentation to be performed for IORT-HDRB using a flab technique. METHODS AND MATERIALS: The material consists of an electromagnetic 3D tracker system, a PC workstation with Microsoft Windows NT 4.0 operating system, and a recognition program for continuous speech. In addition, we designed an external reference system constructed of titanium and Perspex, which is positioned in the pelvis, and a special digitizer pen for reconstruction of the flab geometry. The flab design incorporates a series of silicon 10-mm-diameter spherical pellets. Measurements were made with a pelvic phantom in order to study the accuracy of the system. The reconstruction results are stored and can be exported via network or floppy to our different treatment planning systems. RESULTS: Our results for the reconstruction of a flab with six catheters and a total of 100 spherical pellets give mean errors in the range (2.5 +/- 0.6) mm to (3.5 +/- 0.8) mm depending on the positions of the pelvic phantom and transmitter relative to the operation table. These errors are calculated by comparing the reconstruction results of our system with those using a CT-based reconstruction of the flab geometry. For the accuracy of the navigation system for the pelvic phantom, we obtained mean errors in the range (2.2 +/- 0.7) mm to (3. 1 +/- 1.0) mm. CONCLUSIONS: The new system we have developed enables navigation and reconstruction within the surgical environment with a clinically acceptable level of accuracy. It offers the possibility of individualized treatment planning and effective documentation of the 3D dose distribution in IORT-HDRB using a flab technique.

Brachytherapy↗

3D interstitial HDR brachytherapy combined with 3D external beam radiotherapy and androgen deprivation for prostate cancer. Preliminary results.

BACKGROUND: Evaluation of feasibility, tolerance and efficiency for a new 3D interstitial HDR brachytherapy technique combined with 3D external beam radiotherapy and androgen deprivation for prostate cancer. PATIENTS AND METHODS: Between January 1997 and August 1998 we treated 35 patients with stage cT1-3 N0 M0 prostate cancer. Thirty-two patients with a follow-up of 12 to 28 months (median: 18 months) were evaluated. After ultrasound-guided transrectal implantation of 4 non-parallel needles, CT based 3D brachytherapy treatment planning ("Offenbach system") was performed. All patients received 4 fractions brachytherapy using a fractional dose of 5 or 7 Gy. Time between each fraction was 14 days. After brachytherapy 3D external irradiation followed up to 39.6 or 45.0 Gy. All patients received androgen deprivation, starting 2 to 19 months before brachytherapy, ending 3 months after 3D external radiotherapy. RESULTS: Posttreatment PSA levels dropped to < 1.5 ng/ml in 29/32 patients (91%). In 25 patients PSA levels were < 0.5 ng/ml, in 4 patients 0.5 to 1.5 ng/ml. In 2 patients we noted biochemical relapse. Transrectal implantation was very well tolerated. Grade 3 acute urinary toxicity occurred in 1 patient. We noted no Grade 2 or higher acute gastrointestinal toxicity. One patient developed a Grade 3 late urinary toxicity. No patient showed late gastrointestinal side effects. All 140 dose-volume histograms for 3D HDR brachytherapy were analyzed. CONCLUSIONS: The new 3D HDR brachytherapy technique, combined with 3D external irradiation and androgen deprivation, is a feasible, so far well-tolerated and effective treatment in the short-time follow-up of median 18 months.

Adenocarcinoma↗

Navigation system for interstitial brachytherapy.

PURPOSE: To develop a computed tomography (CT) based electromagnetic navigation system for interstitial brachytherapy. This is especially designed for situations when needles have to be positioned adjacent to or within critical anatomical structures. In such instances interactive 3D visualisation of the needle positions is essential. METHODS AND MATERIALS: The material consisted of a Polhemus electromagnetic 3D digitizer, a Pentium 200 MHz laptop and a voice recognition for continuous speech. In addition, we developed an external reference system constructed of Perspex which could be positioned above the tumour region and attached to the patient using a non-invasive fixation method. A specially designed needle holder and patient bed were also developed. Measurements were made on a series of phantoms in order to study the efficacy and accuracy of the navigation system. RESULTS: The mean navigation accuracy of positioning the 20.0 cm length metallic needles within the phantoms was in the range 2.0-4.1 mm with a maximum of 5.4 mm. This is an improvement on the accuracy of a CT-guided technique which was in the range 6.1-11.3 mm with a maximum of 19.4 mm. The mean reconstruction accuracy of the implant geometry was 3.2 mm within a non-ferromagnetic environment. We found that although the needles were metallic this did not have a significant influence. We also found for our experimental setups that the CT table and operation table non-ferromagnetic parts had no significant influence on the navigation accuracy. CONCLUSIONS: This navigation system will be a very useful clinical tool for interstitial brachytherapy applications, particularly when critical structures have to be avoided. It also should provide a significant improvement on our existing technique.

Brachytherapy↗

Optimized bounding boxes for three-dimensional treatment planning in brachytherapy.

It is sometimes necessary to determine the optimal value for a direction dependent quantity. Using a search technique based on Powell's quadratic convergent method such an optimal direction can be approximated. The necessary geometric transformations in n-dimensional space are introduced. As an example we consider the approximation of the minimum bounding box of a set of three-dimensional points. Minimum bounding boxes can significantly improve accuracy and efficiency of the calculations in modern brachytherapy treatment planning of the volumes of objects or the dose distribution inside an object. A covariance matrix based approximation method for the minimum bounding box is compared with the results of the search method. The benefits of the use of optimal oriented bounding boxes in brachytherapy treatment planning systems are demonstrated and discussed.

Algorithms↗

Dose-volume histograms computation comparisons using conventional methods and optimized fast Fourier transforms algorithms for brachytherapy.

In anatomy based optimization procedures for large volume implants the calculation of dose-volume histograms (DVH) accounts for the major part of the time involved and can be as long as a few hours. This time is proportional to the number of seeds or source dwell positions required for the implant. A procedure for the calculation of brachytherapy seed dose distribution calculation employing fast Fourier transforms (FFT) and the convolution theorem has been described by others and was supposed to significantly improve the speed of the dose distribution computation. Using new significantly improved FFT algorithms and various other optimization techniques we have compared the calculated differential and integral DVHs in high dose rate (HDR) brachytherapy with a single stepping source using actual clinical implants. This is so that we could assess the efficiency and accuracy of the FFT method with that of conventional methods. Our results showed that the FFT based method of calculating DVHs in brachytherapy is comparable in speed with conventional dose calculation methods, but only for implants with more than 287 sources. It is therefore of limited practical use even for large implants. This result is in direct opposition to the claim by other authors.

Algorithms↗

Autoactivation of source dwell positions for HDR brachytherapy treatment planning.

The most accurate classical dose optimization algorithms in HDR brachytherapy strongly depend on an appropriate selection of source dwell positions which fulfill user-defined geometrical boundary conditions which are relative to patient anatomy. Most anatomical situations, such as for prostate and head and neck tumors, are complex and can require geometries with 5-15 catheters with 48 possible dwell positions per catheter depending on the tumor volume. The manual selection of dwell positions using visual checks by trial and error is very time consuming. This can only be improved by the use of a technique which automatically recognizes and selects the optimum dwell positions for each catheter. We have developed an algorithm, termed an autoactivation algorithm, which improves implant planning by providing a facility for the necessary automatic recognition of HDR source dwell positions.

Algorithms↗

Generation of uniformly distributed dose points for anatomy-based three-dimensional dose optimization methods in brachytherapy.

We have studied the accuracy of statistical parameters of dose distributions in brachytherapy using actual clinical implants. These include the mean, minimum and maximum dose values and the variance of the dose distribution inside the PTV (planning target volume), and on the surface of the PTV. These properties have been studied as a function of the number of uniformly distributed sampling points. These parameters, or the variants of these parameters, are used directly or indirectly in optimization procedures or for a description of the dose distribution. The accurate determination of these parameters depends on the sampling point distribution from which they have been obtained. Some optimization methods ignore catheters and critical structures surrounded by the PTV or alternatively consider as surface dose points only those on the contour lines of the PTV. D(min) and D(max) are extreme dose values which are either on the PTV surface or within the PTV. They must be avoided for specification and optimization purposes in brachytherapy. Using D(mean) and the variance of D which we have shown to be stable parameters, achieves a more reliable description of the dose distribution on the PTV surface and within the PTV volume than does D(min) and D(max). Generation of dose points on the real surface of the PTV is obligatory and the consideration of catheter volumes results in a realistic description of anatomical dose distributions.

Algorithms↗