PubMed Health⌕ Search

Biomedical subjects

James M Galvin

Publications and source records attributed to James M Galvin.

12 recordsLinked to original sources

Inherent smoothness of intensity patterns for intensity modulated radiation therapy generated by simultaneous projection algorithms.

The efficient delivery of intensity modulated radiation therapy (IMRT) depends on finding optimized beam intensity patterns that produce dose distributions, which meet given constraints for the tumour as well as any critical organs to be spared. Many optimization algorithms that are used for beamlet-based inverse planning are susceptible to large variations of neighbouring intensities. Accurately delivering an intensity pattern with a large number of extrema can prove impossible given the mechanical limitations of standard multileaf collimator (MLC) delivery systems. In this study, we apply Cimmino's simultaneous projection algorithm to the beamlet-based inverse planning problem, modelled mathematically as a system of linear inequalities. We show that using this method allows us to arrive at a smoother intensity pattern. Including nonlinear terms in the simultaneous projection algorithm to deal with dose-volume histogram (DVH) constraints does not compromise this property from our experimental observation. The smoothness properties are compared with those from other optimization algorithms which include simulated annealing and the gradient descent method. The simultaneous property of these algorithms is ideally suited to parallel computing technologies.

Algorithms↗

The dose-volume constraint satisfaction problem for inverse treatment planning with field segments.

The prescribed goals of radiation treatment planning are often expressed in terms of dose-volume constraints. We present a novel formulation of a dose-volume constraint satisfaction search for the discretized radiation therapy model. This approach does not rely on any explicit cost function. Inverse treatment planning uses the aperture-based approach with predefined, according to geometric rules, segmental fields. The solver utilizes the simultaneous version of the cyclic subgradient projection algorithm. This is a deterministic iterative method designed for solving the convex feasibility problems. A prescription is expressed with the set of inequalities imposed on the dose at the voxel resolution. Additional constraint functions control the compliance with selected points of the expected cumulative dose-volume histograms. The performance of this method is tested on prostate and head-and-neck cases. The relationships with other models and algorithms of similar conceptual origin are discussed. The demonstrated advantages of the method are: the equivalence of the algorithmic and prescription parameters, the intuitive setup of free parameters, and the improved speed of the method as compared to similar iterative as well as other techniques. The technique reported here will deliver approximate solutions for inconsistent prescriptions.

Algorithms↗

A graph-searching method for MLC leaf sequencing under constraints.

A new leaf-sequencing algorithm for step-and-shoot IMRT that is based on a graph-searching technique is described. An iterative process guided by a quantitative measure for the complexity of the initial or residual intensity pattern is used to identify the field segments shaped by a multileaf collimator (MLC). Given a user selected number of intensity levels, the algorithm searches deliverable segment candidates considering all intensity levels and two collimator positions separated by 90 degrees. The candidates for each intensity level are obtained as the least number of segments to cover the areas with equal or higher intensity. The shape of a deliverable segment is adjusted by leaving out certain beam elements for later delivery if this results in a simpler residual intensity pattern and the segment is still deliverable. For a MLC design that does not allow leaf interdigitation, it is initially assumed that a single segment cannot cover two disjoined areas. Among all candidates the segment with the greatest reduction of the complexity of the residual intensity distribution is chosen for the current step of iteration. The iterative process generates a set of deliverable segments of simply connected areas. These segments are combined later under specific MLC constraints. Different orders of segment combination are considered for minimizing the beam-on time. The final segments are sequenced to minimize the leaf travel. This algorithm has been tested using randomly generated intensity distributions and clinical cases for the Varian, Siemens, and Elekta MLC systems. The results show that as the number of intensity levels is increased, the numbers of segments and MUs increase only modestly. Using two collimator angles results in decreases in the required number of segments and the number of monitor units that can be as much as 20%.

Algorithms↗

An optimization algorithm for intensity modulated radiotherapy--the simulated dynamics with dose-volume constraints.

We have developed a new method for optimization in intensity modulated radiation therapy (IMRT) that makes use of simulated dynamics in a classical system of interacting particles. An analogy is drawn between intensity profile optimization in IMRT and relaxation to the equilibrium configuration in a dynamic system. The intensities of beamlets are equivalent to the positions of the virtual particles. The potential energy of the system is defined by the objective function, which determines the equations of motion for the virtual particles. In this paper, we present the implementation of dose constraints and dose-volume constraints. Our strategy is to optimize the dose to the planned target volume (PTV) while keeping all constraints to the organs at risk (OARs) satisfied rigorously. A simple quadratic objective function is used that only includes terms for PTV voxels. By this approach, no additional parameters other than that for prescribing desired dose and constraints, such as importance factors, are needed. The hard constraints that require non-negative beamlet intensities and that the dose at any voxel in an OAR cannot exceed a maximum tolerance, are implemented as semi-transmittable potential barriers of infinite height. Dose-volume constraints are handled by placing hard constraints on partial volumes. Handling of the clinically applied constraints was tested using phantoms and clinical cases. Our results show that the dose-volume histogram (DVH) type of plan prescription can be fulfilled with satisfactory PTV coverage. In addition to the convenience of implementation, our method can achieve a high computational efficiency with the understanding of the dynamic behavior of the system.

Algorithms↗

Guidance document on delivery, treatment planning, and clinical implementation of IMRT: report of the IMRT Subcommittee of the AAPM Radiation Therapy Committee.

Intensity-modulated radiation therapy (IMRT) represents one of the most significant technical advances in radiation therapy since the advent of the medical linear accelerator. It allows the clinical implementation of highly conformal nonconvex dose distributions. This complex but promising treatment modality is rapidly proliferating in both academic and community practice settings. However, these advances do not come without a risk. IMRT is not just an add-on to the current radiation therapy process; it represents a new paradigm that requires the knowledge of multimodality imaging, setup uncertainties and internal organ motion, tumor control probabilities, normal tissue complication probabilities, three-dimensional (3-D) dose calculation and optimization, and dynamic beam delivery of nonuniform beam intensities. Therefore, the purpose of this report is to guide and assist the clinical medical physicist in developing and implementing a viable and safe IMRT program. The scope of the IMRT program is quite broad, encompassing multileaf-collimator-based IMRT delivery systems, goal-based inverse treatment planning, and clinical implementation of IMRT with patient-specific quality assurance. This report, while not prescribing specific procedures, provides the framework and guidance to allow clinical radiation oncology physicists to make judicious decisions in implementing a safe and efficient IMRT program in their clinics.

Algorithms↗

Beam orientation optimization for IMRT by a hybrid method of the genetic algorithm and the simulated dynamics.

We have developed a new method for beam orientation optimization in intensity-modulated radiation therapy (IMRT). The problem of beam orientation optimization in IMRT is solved by a decoupled two-step iterative process: (1) optimization of the intensity profiles for given beam configurations; (2) selection of optimal beam configurations based on the ranking by an objective function score for the results of the intensity profile optimization. The simulated dynamics algorithm is used for the intensity profile optimization. This algorithm enforces both the hard constraints and dose-volume constraints. A genetic algorithm is used to select beam orientation configurations. The method has been tested for both a simulated and clinical case, and the results show that beam orientation optimization significantly improved IMRT plans within a time period that is clinically acceptable. The results also show the dependence of the optimal orientation configurations on the prescribed constraints. In addition, beam orientation optimization by the method described here can provide multiple plans with similar dose distributions. This degeneracy characteristic can be exploited to our advantage in introducing additional planning objectives, e.g., the smoothness of intensity profiles, for the selection of the optimal plan among the degenerate configurations for treatment delivery.

Algorithms↗

The use of mixed-integer programming for inverse treatment planning with pre-defined field segments.

Complex intensity patterns generated by traditional beamlet-based inverse treatment plans are often very difficult to deliver. In the approach presented in this work the intensity maps are controlled by pre-defining field segments to be used for dose optimization. A set of simple rules was used to define a pool of allowable delivery segments and the mixed-integer programming (MIP) method was used to optimize segment weights. The optimization problem was formulated by combining real variables describing segment, weights with a set of binary variables, used to enumerate voxels in targets and critical structures. The MIP method was compared to the previously used Cimmino projection algorithm. The field segmentation approach was compared to an inverse planning system with a traditional beamlet-based beam intensity optimization. In four complex cases of oropharyngeal cancer the segmental inverse planning produced treatment plans, which competed with traditional beamlet-based IMRT plans. The mixed-integer programming provided mechanism for imposition of dose-volume constraints and allowed for identification of the optimal solution for feasible problems. Additional advantages of the segmental technique presented here are: simplified dosimetry, quality assurance and treatment delivery.

Algorithms↗

A dosimetric comparison of various multileaf collimators.

The dosimetric characteristics of three multileaf collimator (MLC) systems (Elekta, Siemens and Varian) having 10 mm leaf width are compared. A 6 MV photon beam was used from each unit for measurements. Film dosimetry was performed for the measurements and the analysis techniques were exactly duplicated in each system. Two of the collimators have rounded leaf ends (Elekta and Varian) and the third (Siemens) has a flat end that follows beam divergence. A scanning densitometer (Wellhöfer with 0.45 mm spot and 0.5 mm step size) was used for film analysis. The dosimetric characteristics studied include: penumbra width (80-20%) as a function of position of the leaf end in the field, inter- and intra-leaf radiation leakage, dose distribution of the tongue and groove, and isodose curves for stepped leaves forming 45 degrees angle beam edge. Results show that MLC designs with divergent and non-divergent leaves produce penumbra (80-20%) widths that are within 2.0 mm of each other. However, the distance of the collimator from the x-ray target plays an important role, and the smallest penumbra width was noted for the Varian MLC despite its rounded leaf-end design. Compared to the other systems, this collimator is positioned about 15 cm closer to the patient which affects the skin dose. The MLC with flat leaf end, although closer to the target, showed slightly poorer penumbra width. Inter-leaf leakage through the leaves is 1.3% for two of the collimators (Elekta and Varian) with the backup jaws and is nearly 1% for the third system (Siemens). The Siemens MLC produces reduced tongue-and-groove effect compared to the other two collimators (Elekta and Varian). The isodose undulation for a stepped edge is found to be significant for the collimator closest to the patient (Varian) and does not depend on the leaf-end shape. There is no perfect MLC system that can be recommended, rather each one has unique advantages and disadvantages that should be weighed with comfort, ease and cost effectiveness for clinical use.

Radiometry↗

A deterministic iterative least-squares algorithm for beam weight optimization in conformal radiotherapy.

Currently, inverse treatment planning in conformal radiotherapy is, in part, a trial-and-error process due to the interplay of many competing criteria for obtaining a clinically acceptable dose distribution. A new method is developed for beam weight optimization that incorporates clinically relevant nonlinear and linear constraints. The process is driven by a nonlinear, quasi-quadratic objective function and the solution space is defined by a set of linear constraints. At each step of iteration, the optimization problem is linearized by a self-consistent approximation that is local to the existing dose distribution. The dose distribution is then improved by solving a series of constrained least-squares problems using an established method until all prescribed constraints are satisfied. This differs from the current approaches in that it does not rely on the search for the global minimum of a specific objective function. Essentially, our proposed objective function can be construed as a functional that comprises a class of dose-based quadratic objective functions. Empirical adjustment for appropriate model parameters in the construction of objective function is minimized, since these parameters are in effect adaptively adjusted during optimization. The method is robust in solving difficult clinical cases using either aperture or pencil beam based planning techniques for intensity-modulated radiation therapy.

Algorithms↗

Comparison of characteristics of photon and electron beams generated by Philips/Elekta and Varian linear accelerators.

Commissioning data of the common energy of 6 MV photon beams are analyzed for seven Philips/ Elekta linacs of five different models: SL 75/5, SL 15, SL 25, SLi Precise, and SLi. The observed beam quality %dd(10)x of 67.7% +/- 0.3 agrees with the manufacturer's specification of 67.5% +/- 2% but with a much smaller variation, suggesting the possibility of a more accurate beam quality guideline for acceptance commissioning of these linac models. Field size factors are less uniform, especially for the biggest field size of 40 x 40 cm2, but all agree within 1.1% between the five models. These 6 MV photon beams are also shown to be 0.4 MV harder than the Varian 2100C model. For commissioning electron beams generated by the Philips/Elekta models SL 15, SLi Precise and SLi, guidelines can be found from this work for the following parameters: actual field size at 100 SSD (up to +/- 0.6%) (different from manufacturer specified values at 95 SSD), and virtual source distance (up to +/- 0.7%). Significant differences are seen for R50 (up to 3.7 mm), sigma-theta-x, effective source distance, and cone size factor. However, except for R50 where Varian 2100C datasets are not available for comparison, the differences between the different Philips/Elekta models are comparable to those between the same Varian 2100C model linacs measured by us and Watts.

Electrons↗

Output variation from an intensity modulating dynamic collimator.

Intensity modulated radiation therapy (IMRT) offers a method of delivering radiation dose conforming to the shape of targets while minimizing the dose to the surrounding tissue and nearby critical organs. One popular device is the NOMOS MIMiC Collimator coupled to the CORVUS treatment planning system. The MIMiC collimator, mounted on a linac head, opens and closes one or more of its 40 small leaves as determined by the planner while the linac delivers radiation and the gantry rotates. This dynamic IMRT allows the intensity to be modulated yielding a highly conformal dose distribution. However, the dose output becomes a function of the detailed manner in which the leaves open and close, since the opening and closing are not instantaneous. We investigate the effect of switch rates and delay in the open/close events on the output profiles. The output is enhanced as the switch rate increases. The enhancement factor at any point of measurement is dependent on its distance from the central plane. We interpret these variations in terms of a simple model, which includes the effect of leaf travel time during the process of opening and closing. We also include the time delay in establishing the specified pressure in the pneumonic device, which controls the opening and closing of the leaves. The information presented here offers a means for incorporating these output changes into the planning system. This may avoid the current situation where many patient plans need to be renormalized based on the actual measurement taken during the delivery of the specified intensity pattern to a phantom.

Equipment Failure Analysis↗