Search PubMed⌕ Search

Biomedical subjects

A L Boyer

Publications and source records attributed to A L Boyer.

At least 19 recordsLinked to original sources

Breast-conserving radiation therapy using combined electron and intensity-modulated radiotherapy technique.

BACKGROUND AND PURPOSE: To explore the feasibility of a multi-modality breast-conserving radiation therapy treatment technique to reduce high dose to the ipsilateral lung and the heart when compared with the conventional treatment technique using two tangential fields. MATERIALS AND METHODS: An electron beam with appropriate energy was combined with four intensity modulated photon beams. The direction of the electron beam was chosen to be tilted 10-20 degrees laterally from the anteroposterior direction. Two of the intensity-modulated photon beams had the same gantry angles as the conventional tangential fields, whereas the other two beams were rotated 15-25 degrees toward the anteroposterior directions from the first two photon beams. An iterative algorithm was developed which optimizes the weight of the electron beam as well as the fluence profiles of the photon beams for a given patient. Two breast cancer patients with early-stage breast tumors were planned with the new technique and the results were compared with those from 3D planning using tangential fields as well as 9-field intensity-modulated radiotherapy (IMRT) techniques. RESULTS: The combined electron and IMRT plans showed better dose conformity to the target with significantly reduced dose to the ipsilateral lung and, in the case of the left-breast patient, reduced dose to the heart, than the tangential field plans. In both the right-sided and left-sided breast plans, the dose to other normal structures was similar to that from conventional plans and was much smaller than that from the 9-field IMRT plans. The optimized electron beam provided between 70 to 80% of the prescribed dose at the depth of maximum dose of the electron beam. CONCLUSIONS: The combined electron and IMRT technique showed improvement over the conventional treatment technique using tangential fields with reduced dose to the ipsilateral lung and the heart. The customized beam directions of the four IMRT fields also kept the dose to other critical structures to a minimum.

Breast Neoplasms↗

Dosimetric effects of patient displacement and collimator and gantry angle misalignment on intensity modulated radiation therapy.

PURPOSE AND OBJECTIVE: The primary goal of this study was to examine systematically the dosimetric effect of small patient movements and linear accelerator angular setting misalignments in the delivery of intensity modulated radiation therapy. We will also provide a method for estimating dosimetric errors for an arbitrary combination of these uncertainties. MATERIALS AND METHODS: Sites in two patients (lumbar-vertebra and nasopharynx) were studied. Optimized intensity modulated radiation therapy treatment plans were computed for each patient using a commercially available inverse planning system (CORVUS, NOMOS Corporation, Sewickley, PA). The plans used nine coplanar beams. For each patient the dose distributions and relevant dosimetric quantities were calculated, including the maximum, minimum, and average doses in targets and sensitive structures. The corresponding dose volumetric information was recalculated by purposely varying the collimator angle or gantry angle of an incident beam while keeping other beams unchanged. Similar calculations were carried out by varying the couch indices in either horizontal or vertical directions. The intensity maps of all the beams were kept the same as those in the optimized plan. The change of a dosimetric quantity, Q, for a combination of collimator and gantry angle misalignments and patient displacements was estimated using Delta=Sigma(DeltaQ/Deltax(i))Deltax(i). Here DeltaQ is the variation of Q due to Deltax(i), which is the change of the i-th variable (collimator angle, gantry angle, or couch indices), and DeltaQ/Deltax(i) is a quantity equivalent to the partial derivative of the dosimetric quantity Q with respect to x(i). RESULTS: While the change in dosimetric quantities was case dependent, it was found that the results were much more sensitive to small changes in the couch indices than to changes in the accelerator angular setting. For instance, in the first example in the paper, a 3-mm movement of the couch in the anterior-posterior direction can cause a 38% decrease in the minimum target dose or a 41% increase in the maximum cord dose, whereas a 5 degrees change in the θ(1)=20 degrees beam only gave rise to a 1.5% decrease in the target minimum or 5.1% in the cord maximum. The effect of systematic positioning uncertainties of the machine settings was more serious than random uncertainties, which tended to smear out the errors in dose distributions. CONCLUSIONS: The dose distribution of an intensity modulated radiation therapy (IMRT) plan changes with patient displacement and angular misalignment in a complex way. A method was proposed to estimate dosimetric errors for an arbitrary combination of uncertainties in these quantities. While it is important to eliminate the angular misalignment, it was found that the couch indices (or patient positioning) played a much more important role. Accurate patient set-up and patient immobilization is crucial in order to take advantage fully of the technological advances of IMRT. In practice, a sensitivity check should be useful to foresee potential IMRT treatment complications and a warning should be given if the sensitivity exceeds an empirical value. Quality assurance action levels for a given plan can be established out of the sensitivity calculation.

Humans↗

Monitor unit calculation for an intensity modulated photon field by a simple scatter-summation algorithm.

An important issue in intensity modulated radiation therapy (IMRT) is the verification of the monitor unit (MU) calculation of the planning system using an independent procedure. Because of the intensity modulation and the dynamic nature of the delivery process, the problem becomes much more involved than that in conventional radiation therapy. In this work, a closed formula for MU calculation is derived. The approach is independent of the specific form of leaf sequence algorithms. It is straightforward to implement the procedure using a simple computer program. The approach is illustrated by a simplified example and is demonstrated by a few CORVUS (NOMOS Corporation, Sewickley, PA) treatment plans. The results indicate that it is robust and suitable for IMRT MU verification.

Algorithms↗

Beam orientation optimization in intensity-modulated radiation treatment planning.

Beam direction optimization is an important problem in radiation therapy. In intensity modulated radiation therapy (IMRT), the difficulty for computer optimization of the beam directions arises from the fact that they are coupled with the intensity profiles of the incident beams. In order to obtain the optimal incident beam directions using iterative or stochastic methods, the beam profiles ought to be optimized after every change of beam configuration. In this paper we report an effective algorithm to optimize gantry angles for IMRT. In our calculation the gantry angles and the beam profiles (beamlet weights) were treated as two separate groups of variables. The gantry angles were sampled according to a simulated annealing algorithm. For each sampled beam configuration, beam profile calculation was done using a fast filtered backprojection (FBP) method. Simulated annealing was also used for beam profile optimization to examine the performance of the FBP for beam orientation optimization. Relative importance factors were incorporated into the objective function to control the relative importance of the target and the sensitive structures. Minimization of the objective function resulted in the best possible beam orientations and beam profiles judged by the given objective function. The algorithm was applied to several model problems and the results showed that the approach has potential for IMRT applications.

Algorithms↗

A dosimetric comparison of fan-beam intensity modulated radiotherapy with Gamma Knife stereotactic radiosurgery for treating intermediate intracranial lesions.

PURPOSE: To compare and evaluate treatment plans for the fan-beam intensity modulated radiotherapy and the Gamma Knife radiosurgery for treating medium-size intracranial lesions (range 4-25 cm3). METHODS AND MATERIALS: Treatment plans were developed for the Leksell Gamma Knife and a fan-beam inverse treatment planning system for intensity modulated radiotherapy. Treatment plan comparisons were carried out using dose-volume histogram (DVH), tissue-volume ratio (TVR), and maximum dose to the prescription dose (MDPD) ratio. The study was carried out for both simulated targets and clinical targets with irregular shapes and at different locations. RESULTS: The MDPD ratio was significantly greater for the Gamma Knife plans than for the fan-beam IMRT plans. The Gamma Knife plans produced equivalent TVR values to the fan-beam IMRT plans. Based on the DVH comparison, the fan-beam IMRT delivered significantly more dose to the normal brain tissue than the Gamma Knife. The results of the comparison were found to be insensitive to the target locations. CONCLUSION: The Gamma Knife is better than the fan-beam IMRT in sparing normal brain tissue while producing equivalent tumor dose conformity for treating medium-size intracranial lesions. However, the target dose homogeneity is significantly better for the fan-beam IMRT than for the Gamma Knife.

Brain Neoplasms↗

Synchronizing dynamic multileaf collimators for producing two-dimensional intensity-modulated fields with minimum beam delivery time.

PURPOSE: Leaf motion synchronization of dynamic multileaf collimators (DMLC) for intensity-modulated radiotherapy (IMRT) is important in improving dose distribution and reducing "tongue-and-groove" effects for a prescribed intensity profile. Leaf synchronization could also be used in transforming a one-dimensional leaf-setting algorithm into a two-dimensional leaf-setting algorithm. In this work, we aim to develop a generalized leaf synchronization method for delivering IMRT with the minimized beam delivery time and the optimized subfield variations for a leaf-setting sequence. METHODS AND MATERIALS: With the leaf synchronization procedure, all active MLC leaf pairs start and finish off a leaf sequence simultaneously. In this work, the MLC leaf pairs were synchronized under the condition that the resulting leaf sequence produces the desired intensity profile with the minimum beam delivery time. The parameter of the leaf synchronization function was determined through the least-square minimization of the area variations of all subfields within a leaf sequence. The leaf synchronization and optimization procedure were applied and analyzed for clinical relevant intensity profiles for treating the head-and-neck cancer patients using IMRT. RESULTS: The total monitor units and the optimized beam delivery time of generating a two-dimensional intensity profile was proven through this work to be the global minimum of all leaf-setting sequences including the unsynchronized leaf-setting sequences. The optimized parameter for subfield variations of the synchronized leaf trajectories was found to be dependent on individual intensity profiles. For all our studied cases, the unsynchronized leaf trajectories always have significantly larger subfield variations than the synchronized leaf trajectories. CONCLUSION: It is important and also feasible to synchronize and optimize dynamic MLC leaf motions while still keeping the total beam delivery time minimum for delivering arbitrary two-dimensional intensity-modulated fields.

Algorithms↗

Clinical implementation of wedge filter optimization in three-dimensional radiotherapy treatment planning.

BACKGROUND AND PURPOSE: To describe a wedge filter optimization technique which automatically chooses the beam weights and wedge filters and to demonstrate the implementation of the algorithm in clinical three-dimensional (3D) radiotherapy treatment planning. MATERIAL AND METHODS: Given the incident directions and beam energies of J beams, the dose distribution is a function of the beam weights, wedge angles, and wedge orientations. Instead of decomposing an incident field into a superposition of an open and two nominal wedged fields and then optimizing their weights, the algorithm optimizes the objective function with respect to the beam weights, wedge angles and wedge orientations directly. A salient feature of the algorithm is that no planner intervention was required in the selection of wedge filters during the optimization process. A dose-based objective function which incorporated the relative importance of structures was adopted in this work. The objective function was minimized by the method of simulated annealing. The technique was demonstrated by using a phantom study and two clinical cases. RESULTS: For the phantom case, the classical wedge pair result was obtained, providing a useful test of the algorithm. Dose distributions and dose volume histograms for the target and surrounding organs were presented for the two clinical cases. It was also shown that dose homogeneity to the target could be compromised by increasing the relative importance factors to the surrounding organs. CONCLUSIONS: A 3D wedge filter optimization algorithm has been developed. The technique has the potential to fully automate the 3D radiotherapy treatment planning process. In addition, treatment planning time and efforts were significantly reduced.

Algorithms↗

Intensity-modulated radiation therapy with dynamic multileaf collimators.

Intensity-modulated radiotherapy (IMRT) has been considered as a means of providing dose distributions that conform to concave target volumes. For computer-controlled multileaf collimators (MLCs) to be used to modulate x-ray beams, some procedure must be used to determine the sequence of leaf positions used to produce the desired modulation. This article derives and compares four leaf-sequencing algorithms. MLC leaf sequencing can be accomplished by representing the areal intensity modulation of a beam with a series of beam profiles. A velocity-modulation equation for computing the modulation required for a one-dimensional profile, described originally using more extensive algebra, is derived using a graphic approach. The velocity-modulation approach is compared with an equal incremental step-and-shoot approach derived by Bortfeld and Boyer. An areal step-and-shoot technique derived by Xia and Verhey is introduced and compared with the profile-by-profile methods. Finally, an approach is considered using multiple repeated arcs developed by Yu. This wide variety of methods can yield an approach to IMRT that conforms to the engineering constraints imposed by the design of a particular linear accelerator.

Algorithms↗

Optimization of importance factors in inverse planning.

Inverse treatment planning starts with a treatment objective and obtains the solution by optimizing an objective function. The clinical objectives are usually multifaceted and potentially incompatible with one another. A set of importance factors is often incorporated in the objective function to parametrize trade-off strategies and to prioritize the dose conformality in different anatomical structures. Whereas the general formalism remains the same, different sets of importance factors characterize plans of obviously different flavour and thus critically determine the final plan. Up to now, the determination of these parameters has been a 'guessing' game based on empirical knowledge because the final dose distribution depends on the parameters in a complex and implicit way. The influence of these parameters is not known until the plan optimization is completed. In order to compromise properly the conflicting requirements of the target and sensitive structures, the parameters are usually adjusted through a trial-and-error process. In this paper, a method to estimate these parameters computationally is proposed and an iterative computer algorithm is described to determine these parameters numerically. The treatment plan selection is done in two steps. First, a set of importance factors are chosen and the corresponding beam parameters (e.g. beam profiles) are optimized under the guidance of a quadratic objective function using an iterative algorithm reported earlier. The 'optimal' plan is then evaluated by an additional scoring function. The importance factors in the objective function are accordingly adjusted to improve the ranking of the plan. For every change in the importance factors, the beam parameters need to be re-optimized. This process continues in an iterative fashion until the scoring function is saturated. The algorithm was applied to two clinical cases and the results demonstrated that it has the potential to improve significantly the existing method of inverse planning. It was noticed that near the final solution the plan became insensitive to small variations of the importance factors.

Algorithms↗

Dosimetric verification of a commercial inverse treatment planning system.

A commercial three-dimensional (3D) inverse treatment planning system, Corvus (Nomos Corporation, Sewickley, PA), was recently made available. This paper reports our preliminary results and experience with commissioning this system for clinical implementation. This system uses a simulated annealing inverse planning algorithm to calculate intensity-modulated fields. The intensity-modulated fields are divided into beam profiles that can be delivered by means of a sequence of leaf settings by a multileaf collimator (MLC). The treatments are delivered using a computer-controlled MLC. To test the dose calculation algorithm used by the Corvus software, the dose distributions for single rectangularly shaped fields were compared with water phantom scan data. The dose distributions predicted to be delivered by multiple fields were measured using an ion chamber that could be positioned in a rotatable cylindrical water phantom. Integrated charge collected by the ion chamber was used to check the absolute dose of single- and multifield intensity modulated treatments at various spatial points. The measured and predicted doses were found to agree to within 4% at all measurement points. Another set of measurements used a cubic polystyrene phantom with radiographic film to record the radiation dose distribution. The films were calibrated and scanned to yield two-dimensional isodose distributions. Finally, a beam imaging system (BIS) was used to measure the intensity-modulated x-ray beam patterns in the beam's-eye view. The BIS-measured images were then compared with a theoretical calculation based on the MLC leaf sequence files to verify that the treatment would be executed accurately and without machine faults. Excellent correlation (correlation coefficients > or = 0.96) was found for all cases. Treatment plans generated using intensity-modulated beams appear to be suitable for treatment of irregularly shaped tumours adjacent to critical structures. The results indicated that the system has potential for clinical radiation treatment planning and delivery and may in the future reduce treatment complexity.

Algorithms↗

A comparison of Monte Carlo and analytic first scatter dose spread arrays.

We compare first scattered point dose spread arrays generated by Monte Carlo and an analytic method. The analytic method models energy deposition using Klein-Nishina cross sections for Compton scatter and approximations for electron transport. Assumptions in the analytic method are shown to be valid within a region of the point dose spread kernel in which meaningful comparisons can be made. Differences between the models are less than 10% for the forward scatter directions for radii greater than the electron range associated with the first scattered Compton photon. Differences in the backscatter region are discussed and indicate that the analytic model is useful for identifying large errors that might be present in numerically generated first scatter point dose spread arrays. The analytic method is simple and useful for validating first scatter kernels.

Monte Carlo Method↗

Clinical implementation of a Monte Carlo treatment planning system.

The purpose of this study was to implement the Monte Carlo method for clinical radiotherapy dose calculations. We used the EGS4/BEAM code to obtain the phase-space data for 6-20 MeV electron beams and 4, 6, and 15 MV photon beams for Varian Clinac 1800, 2100C, and 2300CD accelerators. A multiple-source model was used to reconstruct the phase-space data for both electron and photon beams, which retained the accuracy of the Monte Carlo beam data. The multiple-source model reduced the phase-space data storage requirement by a factor of 1000 and the accelerator simulation time by a factor of 10 or more. Agreement within 2% was achieved between the Monte Carlo calculations and measurements of the dose distributions in homogeneous and heterogeneous phantoms for various field sizes, source-surface distances, and beam modulations. The Monte Carlo calculated electron output factors were within 2% of the measured values for various treatment fields while the heterogeneity correction factors for various lung and bone phantoms were within 1% for photon beams and within 2% for electron beams. The EGS4/DOSXYZ Monte Carlo code was used for phantom and patient dose calculations. The results were compared to the dose distributions produced by a conventional treatment planning system and an intensity-modulated radiotherapy inverse-planning system. Significant differences (>5% in dose and >5 mm shift in isodose lines) were found between Monte Carlo calculations and the analytical calculations implemented in the commercial systems. Treatment sites showing the largest dose differences were for head and neck, lung, and breast cases.

Algorithms↗

Estimation theory and model parameter selection for therapeutic treatment plan optimization.

Treatment optimization is usually formulated as an inverse problem, which starts with a prescribed dose distribution and obtains an optimized solution under the guidance of an objective function. The solution is a compromise between the conflicting requirements of the target and sensitive structures. In this paper, the treatment plan optimization is formulated as an estimation problem of a discrete and possibly nonconvex system. The concept of preference function is introduced. Instead of prescribing a dose to a structure (or a set of voxels), the approach prioritizes the doses with different preference levels and reduces the problem into selecting a solution with a suitable estimator. The preference function provides a foundation for statistical analysis of the system and allows us to apply various techniques developed in statistical analysis to plan optimization. It is shown that an optimization based on a quadratic objective function is a special case of the formalism. A general two-step method for using a computer to determine the values of the model parameters is proposed. The approach provides an efficient way to include prior knowledge into the optimization process. The method is illustrated using a simplified two-pixel system as well as two clinical cases. The generality of the approach, coupled with promising demonstrations, indicates that the method has broad implications for radiotherapy treatment plan optimization.

Algorithms↗

Matching photon and electron fields with dynamic intensity modulation.

A technique was developed to reduce the size and magnitude of the hot and cold spots in the abutting regions of photon and electron fields. The photon and electron fields were set up such that the photon field extended approximately 2 cm into the electron field in the abutting region. The region of the photon beam that overlapped the electron field was modulated using a multileaf collimator, effectively broadening the photon penumbra to make it complimentary to the electron penumbra. The computer calculations were verified using film measurements for abutting a 6 MV photon beam with a 9 MeV electron beam. A uniform dose was achieved at a prespecified depth of 2 cm, and dose uniformity was improved at the specified depth and beyond compared with unmodulated photon beams. A slight increase in dose inhomogeneity was seen at shallower depths. The overall areas of the hot and cold spots were significantly reduced. The technique also reduced the sensitivity of dose homogeneity to setup errors such that the magnitudes of the hot and cold spots were about half of those produced with unmodulated photon beam when an overlap or gap of 4 mm was introduced. The technique was applied to the treatment of a head and neck cancer and a lymphoma involving the right pleura with markedly reduced dose inhomogeneity in the abutting regions.

Computer Simulation↗

Using static MLC fields to replace partial transmission cerrobend blocks in treatment planning of rectal carcinoma cases.

A treatment planning technique has been developed for using static multileaf collimators to replace partial transmission blocks for treating rectal or cervix carcinoma. The static MLC fields were used to replace the partial transmission block of the anterior-posterior pelvis field for the so-called "thunderbird" technique. Treatment plans were developed and evaluated on a commercial three-dimensional treatment planning system (FOCUS, Computerized, Medical Systems, St. Louis, MO). The result of the treatment plan comparison indicates that the static MLC fields are capable of achieving the same target, inguinal and pelvic dose distribution as the partial transmission cerrobend blocked fields. The MLC fields are easy to modify particularly for the match line adjustments. In conclusion, it is efficient and effective to use static MLC fields to replace partial transmission blocks in the "thunderbird" technique for treating rectal or cervix carcinoma.

Humans↗

Monte Carlo calculations of electron beam output factors for a medical linear accelerator.

The purpose of this study was to investigate the application of the Monte Carlo technique to the calculation and analysis of output factors for electron beams used in radiotherapy. The code EGS4/BEAM was used to obtain phase-space files for 6, 12 and 20 MeV clinical electron beams from a scattering-foil linac (Varian Clinac 2100C) for a clinically representative range of applicator and square or rectangular insert combinations. The source-to-surface distance used was 100 cm. The field sizes ranged from 1 x 1 cm2 to 20 x 20 cm2. These phase-space files were analysed to study the intrinsic beam characteristics and used as source input for relative dose and output factor computations in homogeneous water phantoms using the code EGS4/DOSXYZ. The calculated relative central-axis depth-dose and transverse dose profiles at various depths of clinical interest agreed with the corresponding measured dose profiles to within 2% of the maximum dose. Calculated output factors for the fields studied agreed with measured output factors to about 2%. This demonstrated that for the Varian Clinac 2100C linear accelerator, electron beam dose calculations in homogeneous water phantoms can be performed accurately at the 2% level using Monte Carlo simulations.

Computer Simulation↗

Application of a video-optical beam imaging system for quality assurance of medical accelerators.

Method validation techniques were developed and experiments were carried out using a beam imaging system (BIS, Wellhöfer Dosimetrie, Schwarzenbruck, Germany) for routine quality assurance of medical accelerators. The routine quality assurance tasks include x-ray beam flatness and symmetry check, light/radiation field congruence test, beam energy constancy for electrons and mechanical checks for couch and collimator rotations. Comparisons were made between the BIS application and conventional quality assurance methods that use radiographic films or detector arrays. In this work, we have demonstrated efficiency and accuracy of the BIS to perform some of the routine quality assurance tasks for medical linear accelerators.

Electrons↗