Search PubMed⌕ Search

Biomedical subjects

A L Boyer

Publications and source records attributed to A L Boyer.

At least 19 recordsLinked to original sources

A quality assurance phantom for IMRT dose verification.

This paper investigates a quality assurance (QA) phantom specially designed to verify the accuracy of dose distributions and monitor units (MU) calculated by clinical treatment planning optimization systems and by the Monte Carlo method for intensity-modulated radiotherapy (IMRT). The QA phantom is a PMMA cylinder of 30 cm diameter and 40 cm length with various bone and lung inserts. A procedure (and formalism) has been developed to measure the absolute dose to water in the PMMA phantom. Another cylindrical phantom of the same dimensions, but made of water, was used to confirm the results obtained with the PMMA phantom. The PMMA phantom was irradiated by 4, 6 and 15 MV photon beams and the dose was measured using an ionization chamber and compared to the results calculated by a commercial inverse planning system (CORVUS, NOMOS, Sewickley, PA) and by the Monte Carlo method. The results show that the dose distributions calculated by both CORVUS and Monte Carlo agreed to within 2% of dose maximum with measured results in the uniform PMMA phantom for both open and intensity-modulated fields. Similar agreement was obtained between Monte Carlo calculations and measured results with the bone and lung heterogeneity inside the PMMA phantom while the CORVUS results were 4% different. The QA phantom has been integrated as a routine QA procedure for the patient's IMRT dose verification at Stanford since 1999.

Equipment Design↗

Role of beam orientation optimization in intensity-modulated radiation therapy.

PURPOSE: To investigate the role of beam orientation optimization in intensity-modulated radiation therapy (IMRT) and to examine the potential benefits of noncoplanar intensity-modulated beams. METHODS AND MATERIALS: A beam orientation optimization algorithm was implemented. For this purpose, system variables were divided into two groups: beam position (gantry and table angles) and beam profile (beamlet weights). Simulated annealing was used for beam orientation optimization and the simultaneous iterative inverse treatment planning algorithm (SIITP) for beam intensity profile optimization. Three clinical cases were studied: a localized prostate cancer, a nasopharyngeal cancer, and a paraspinal tumor. Nine fields were used for all treatments. For each case, 3 types of treatment plan optimization were performed: (1) beam intensity profiles were optimized for 9 equiangular spaced coplanar beams; (2) orientations and intensity profiles were optimized for 9 coplanar beams; (3) orientations and intensity profiles were optimized for 9 noncoplanar beams. RESULTS: For the localized prostate case, all 3 types of optimization described above resulted in dose distributions of a similar quality. For the nasopharynx case, optimized noncoplanar beams provided a significant gain in the gross tumor volume coverage. For the paraspinal case, orientation optimization using noncoplanar beams resulted in better kidney sparing and improved gross tumor volume coverage. CONCLUSION: The sensitivity of an IMRT treatment plan with respect to the selection of beam orientations varies from site to site. For some cases, the choice of beam orientations is important even when the number of beams is as large as 9. Noncoplanar beams provide an additional degree of freedom for IMRT treatment optimization and may allow for notable improvement in the quality of some complicated plans.

Algorithms↗

Modeling the extrafocal radiation and monitor chamber backscatter for photon beam dose calculation.

A simple analytical approach has been developed to model extrafocal radiation and monitor chamber backscatter for clinical photon beams. Model parameters for both the extrafocal source and monitor chamber backscatter are determined simultaneously using conventional measured data, i.e., in-air output factors for square and rectangular fields defined by the photon jaws. The model has been applied to 6 MV and 15 MV photon beams produced by a Varian Clinac 2300C/D accelerator. Contributions to the in-air output factor from the extrafocal radiation and monitor chamber backscatter, as predicted by the model, are in good agreement with the measurements. The model can be used to calculate the in-air output factors analytically, with an accuracy of 0.2% for symmetric or asymmetric rectangular fields defined by jaws when the calculation point is at the isocenter and 0.5% when the calculation point is at an extended SSD. For MLC-defined fields, with the jaws at the recommended positions, calculated in-air output factors agree with the measured data to within 0.3% at the isocenter and 0.7% at off-axis positions. The model has been incorporated into a Monte Carlo dose algorithm to calculate the absolute dose distributions in patients or phantoms. For three MLC-defined irregular fields (triangle shape, C-shape, and L-shape), the calculations agree with the measurements to about 1% even for points at off-axis positions. The model will be particularly useful for IMRT dose calculations because it accurately predicts beam output and penumbra dose.

Biophysical Phenomena↗

Verification of IMRT dose distributions using a water beam imaging system.

A water beam imaging system (WBIS) has been developed and used to verify dose distributions for intensity modulated radiotherapy using dynamic multileaf collimator. This system consisted of a water container, a scintillator screen, a charge-coupled device camera, and a portable personal computer. The scintillation image was captured by the camera. The pixel value in this image indicated the dose value in the scintillation screen. Images of radiation fields of known spatial distributions were used to calibrate the device. The verification was performed by comparing the image acquired from the measurement with a dose distribution from the IMRT plan. Because of light scattering in the scintillator screen, the image was blurred. A correction for this was developed by recognizing that the blur function could be fitted to a multiple Gaussian. The blur function was computed using the measured image of a 10 cm x 10 cm x-ray beam and the result of the dose distribution calculated using the Monte Carlo method. Based on the blur function derived using this method, an iterative reconstruction algorithm was applied to recover the dose distribution for an IMRT plan from the measured WBIS image. The reconstructed dose distribution was compared with Monte Carlo simulation result. Reasonable agreement was obtained from the comparison. The proposed approach makes it possible to carry out a real-time comparison of the dose distribution in a transverse plane between the measurement and the reference when we do an IMRT dose verification.

Calibration↗

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↗

Energy- and intensity-modulated electron beams for radiotherapy.

This work investigates the feasibility of optimizing energy- and intensity-modulated electron beams for radiation therapy. A multileaf collimator (MLC) specially designed for modulated electron radiotherapy (MERT) was investigated both experimentally and by Monte Carlo simulations. An inverse-planning system based on Monte Carlo dose calculations was developed to optimize electron beam energy and intensity to achieve dose conformity for target volumes near the surface. The results showed that an MLC with 5 mm leaf widths could produce complex field shapes for MERT. Electron intra- and inter-leaf leakage had negligible effects on the dose distributions delivered with the MLC, even at shallow depths. Focused leaf ends reduced the electron scattering contributions to the dose compared with straight leaf ends. As anticipated, moving the MLC position toward the patient surface reduced the penumbra significantly. There were significant differences in the beamlet distributions calculated by an analytic 3-D pencil beam algorithm and the Monte Carlo method. The Monte Carlo calculated beamlet distributions were essential to the accuracy of the MERT dose distribution in cases involving large air gaps, oblique incidence and heterogeneous treatment targets (at the tissue-bone and bone-lung interfaces). To demonstrate the potential of MERT for target dose coverage and normal tissue sparing for treatment of superficial targets, treatment plans for a hypothetical treatment were compared using photon beams and MERT.

Algorithms↗

Monte Carlo verification of IMRT dose distributions from a commercial treatment planning optimization system.

The purpose of this work was to use Monte Carlo simulations to verify the accuracy of the dose distributions from a commercial treatment planning optimization system (Corvus, Nomos Corp., Sewickley, PA) for intensity-modulated radiotherapy (IMRT). A Monte Carlo treatment planning system has been implemented clinically to improve and verify the accuracy of radiotherapy dose calculations. Further modifications to the system were made to compute the dose in a patient for multiple fixed-gantry IMRT fields. The dose distributions in the experimental phantoms and in the patients were calculated and used to verify the optimized treatment plans generated by the Corvus system. The Monte Carlo calculated IMRT dose distributions agreed with the measurements to within 2% of the maximum dose for all the beam energies and field sizes for both the homogeneous and heterogeneous phantoms. The dose distributions predicted by the Corvus system, which employs a finite-size pencil beam (FSPB) algorithm, agreed with the Monte Carlo simulations and measurements to within 4% in a cylindrical water phantom with various hypothetical target shapes. Discrepancies of more than 5% (relative to the prescribed target dose) in the target region and over 20% in the critical structures were found in some IMRT patient calculations. The FSPB algorithm as implemented in the Corvus system is adequate for homogeneous phantoms (such as prostate) but may result in significant under or over-estimation of the dose in some cases involving heterogeneities such as the air-tissue, lung-tissue and tissue-bone interfaces.

Computer Simulation↗

Calculation of x-ray transmission through a multileaf collimator.

A ray tracing based method has been developed to calculate the x-ray transmission through a multileaf collimator (MLC) for beam delivery verification and dose calculation in intensity modulated radiotherapy (IMRT). The path length of a ray line in the MLC is accurately calculated using the exact geometry of the MLC leaves. The fluence distribution of an IMRT field is calculated first using a point source. The fluence distribution for a realistic beam model is obtained, as an approximation, by convolving the point source fluence distribution with the distribution of source strength. Full ray tracing calculations are performed using analytic and Monte Carlo simulated beam models to verify the accuracy of the convolution method. The calculation is in better agreement with measurements using either film or a beam imaging system (BIS) than previous calculations for MLC transmission using a simplified model. This ray tracing calculation can be applied to the problem of verifying dynamic MLC leaf sequences as part of a patient-specific quality assurance process for IMRT.

Algorithms↗

A dose-volume histogram based optimization algorithm for ultrasound guided prostate implants.

The task of treatment planning for prostate implants is to find an optimal seed configuration, comprising the target coverage and dosimetric consideration of critical structures such as the rectum and urethra. An efficient method to accomplish this is to use an inverse planning technique that derives the optimized solution from a prescribed treatment goal. The goal can be specified in the voxel domain as the desired doses to the voxels of the target and critical structures, or in the dose volume representation as the desired dose volume histograms (DVHs) of the target and critical structures. The DVH based optimization has been successfully used in plan optimization for intensity-modulated radiation therapy (IMRT) but little attention has been paid to its application in prostate implants. Clinically, it has long been known that some normal structure tolerances are more accurately assessed by volumetric information. Dose-volume histograms are also widely used for plan evaluation. When working in the DVH domain for optimization one has more control over the final DVHs. We have constructed an objective function sensitive to the DVHs of the target and critical structures. The objective function is minimized using an iterative algorithm, starting from a randomly selected initial seed configuration. At each iteration step, a trial position is given to a randomly selected source and the trial position is accepted if the objective function is decreased. To avoid being trapped in a less optimal local minimum, the optimization process is repeated. The final plan is selected from a pool of optimized plans obtained from a series of randomized initial seed configurations.

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↗