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.

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.

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.

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.

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.

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.

CT imaging based digitally reconstructed radiographs and their application in brachytherapy.

The aim of our study was to develop an algorithm to simulate the digitally reconstructed radiograph (DRR) calculation process for different beam qualities (photon energies) in the range 50 keV to 12 MeV. This was achieved using volumetric anatomical data for the patient obtained from three-dimensional diagnostic CT images. These DRR images can be used in three-dimensional treatment planning for external beam radiotherapy as well as for brachytherapy in the same way as conventional radiographic films. The advantages of using such DRRs in modern 3D brachytherapy treatment planning are shown. A number of tools are described, illustrating that the application of DRRs in brachytherapy is very useful.

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.

Catheter autoreconstruction in computed tomography based brachytherapy treatment planning.

The aim of this study is to develop an automatic reconstruction of brachytherapy catheters using CT (computed tomography) data. Previously no such automatic facility has existed in any treatment planning software. To achieve this facility we have developed tools for the automatic reconstruction (which we term autoreconstruction) of plastic and metallic catheters. These algorithms overcome a number of difficulties which arise when a large number of catheters are present. These include situations with intersecting catheters and with loop techniques. The time required for the catheter reconstruction process using our autoreconstruction method is significantly reduced. The accuracy of our autoreconstruction is at least as high as the classical manual slice-by-slice method.

Comparison of calibration procedures for 192Ir high-dose-rate brachytherapy sources.

PURPOSE: To compare the efficacy of different calibration procedures for 192Ir high-dose-rate (HDR) brachytherapy sources and to determine their suitability in clinical practice. In addition the manufacturer's calibration is compared with our experimental measurements so that the accuracy of the source strength on the manufacturer certificate which is supplied with each new 192Ir source can be accessed. METHODS AND MATERIALS: We compared three types of calibration system: well-type chambers (HDR-1000 and SDS), cylindrical phantom, and plate phantom. The total number of measurements we obtained was 365. The number of sources used for the calibration procedure comparison was 20 and the number used for comparison with the manufacturer's calibration was 46. This study was made during the period 1989-1997. Also, Physikalisch-Technische Bundesanstalt (PTB) calibrated one of our sources using their PTB protocol so that the results could be compared with our own. RESULTS: The sensitivity of each system on scattering from the room walls was studied. It was found that different minimum lateral distances from the walls were required for the different systems tested: 15 cm and 25 cm for the well-type chambers, 75 cm for the cylindrical phantom, and 13 cm for the plate phantom. The minimum thickness required to reach phantom scattering saturation for the plate phantom setup is 24 cm. The influence of the applicator material used in the calibration setup was found to be 1.7% for the stainless steel dosimetry applicator compared to the plastic 5F applicator. The accuracy of source positioning within the applicator can lead to dosimetric errors of +/-1.2% for the radial distance of 8.0 cm used with both solid phantoms. The change in the response for both well-type chambers was only 0.1% for changes in the source position within +/-7.5 mm around the response peak. Good agreement was found between all dosimetry systems included in our study. Taking the HDR-1000 well-type chamber results as a reference, we observed percentage root mean square (RMS) values of 0.11% for the SDS well-type chamber, 0.44% for the cylindrical, and 0.60% for the plate phantom setup. A comparison of our results using the cylindrical phantom with those of the manufacturer showed a percentage RMS value of 3.3% with a percentage fractional error range of -13.0% to +6.0%. The comparison of our calibration results with those of PTB gave deviations less than 0.4% for all systems. CONCLUSIONS: Our results have shown that with careful use of all calibration system protocols an accurate determination of source strength can be obtained. However, the manufacturer's calibration is not accurate enough on its own, and it should be mandatory for clinics to always measure the source strength of newly delivered 192Ir brachytherapy sources. The influence of the applicator material, metal or plastic, should always be taken into account.