Qiuwen Wu

Number and orientations of beams in intensity‐modulated radiation treatments

The fundamental question of how many equispaced coplanar intensity-modulated photon beams are required to obtain an optimum treatment plan is investigated in a dose escalation study for a typical prostate tumor. Furthermore, optimization of beam orientations to improve dose distributions is explored. A dose-based objective function and a fast gradient technique are employed for optimizing the intensity profiles (inverse planning). An exhaustive search and fast simulated annealing techniques (FSA) are used to optimize beam orientations. However, to keep computation times reasonable, the intensity profiles for each beam arrangement are still optimized using inverse planning. A pencil beam convolution algorithm is employed for dose calculation. All calculations are performed in three-dimensional (3D) geometry for 15 MV photons. DVHs, dose displays, TCP, NTCP, and biological score functions are used for evaluation of treatment plans. It is shown that for the prostate case presented here, the minimum required number of equiangular beams depends on the prescription dose level and ranges from three beams for 70 Gy plans to seven to nine beams for 81 Gy plans. For the highest dose level (81 Gy), beam orientations are optimized and compared to equiangular spaced arrangements. It is shown that (1) optimizing beam orientations is most valuable for a small numbers of beams (< or = 5) and the gain diminishes rapidly for higher numbers of beams; (2) if sensitive structures (for example rectum) are partially enclosed by the target volume, beams coming from their direction tend to be preferable, since they allow greater control over dose distributions; (3) while FSA and an exhaustive search lead to the same results, computation times using FSA are reduced by two orders of magnitude to clinically acceptable values. Moreover, characteristics of and demands on biology-based and dose-based objective functions for optimization of intensity-modulated treatments are discussed.

Algorithms and functionality of an intensity modulated radiotherapy optimization system

The main purpose of this paper is to describe formalisms, algorithms, and certain unique features of a system for optimization of intensity modulated radiotherapy (IMRT). The system is coupled to a commercial treatment planning system with an accurate dose calculation engine based on the kernel superposition algorithm. The system was designed for use for research as well as for routine clinical practice. It employs dose- and dose-volume-based objective functions. The system can optimize IMRT plans with multiple target volumes simultaneously. Each target volume may be assigned a different prescription dose with constraints on either underdosing, or overdosing, or both. For organs at risk more than one constraint may be applied. This feature allows simultaneous treatment of primary, regional disease and electively treated nodes. The system allows specification of constraints on logical combinations of anatomic structures, such as a region of overlap between the prostate planning target volume and rectum or the volume of lung excluding the tumor. The optimization may also be performed on plans which, in addition to intensity-modulated beams, include other modalities such as non-IMRT photon and electron beams and brachytherapy sources. The various features of the system are illustrated with one phantom example and two clinical examples: a brain stereotactic radiosurgery case and a nasopharynx case. In the cylindrical phantom example, the use of the system for overlap regions is demonstrated. The brain stereotactic radiosurgery example shows the improvement of IMRT plans over the conventional arcs based plan and the three-dimensional conformal plan with multiple fixed gantry angles and demonstrates the application of our system to cases where small grid sizes are important. The nasopharynx example shows the potential of IMRT to simultaneously treat large and boost fields. It also illustrates the power of IMRT to protect normal anatomic structures for highly complex situations and the efficiency in planning and delivery achievable with IMRT. The overall IMRT planning time is typically less than 2 h on a Sun Ultrasparc workstation, most of which is spent in repeated computation of dose distributions.

A comparison of three stereotactic radiotherapy techniques; arcs vs. noncoplanar fixed fields vs. intensity modulation

PURPOSE Linac arc based stereotactic radiotherapy is being used with increasing frequency to treat brain tumors. This approach can be used for single or fractionated treatments, and is typically carried out with circular collimators which are optimal for small, spherical targets. Treatment planning using fixed noncoplanar beams or intensity-modulated beams may enhance the ability to conform to irregularly shaped and/or large tumors, especially when combined with stereotactic localization. We compare the dose conformity and normal brain dose characteristics of three stereotactic techniques for various nonspherical target shapes. METHODS AND MATERIALS Three intracranial test targets were constructed using a 3D treatment planning system after a patient underwent CT simulation. The targets included an ellipsoid with major axis dimensions of 4.0, 2.0, and 2.0 cm, a hemisphere with a diameter of 4.0 cm, and an irregularly shaped patient tumor with a maximum dimension of 5.3 cm. The following stereotactic techniques were compared for each target: a) 5 arcs as used in traditional linac radiosurgery/radiotherapy (noncoplanar arcs [ARCS]), b) 6 fixed noncoplanar custom blocked fields (3D), c) intensity modulation using 6 noncoplanar beams and a mini-multileaf collimator (intensity-modulated radiation therapy [IMRTI). Dose volume histograms were performed for each target/technique combination. RESULTS For the ellipsoid, dose conformity is similar for all three techniques and normal brain isodose distributions are more favorable with the ARCS plan. For the hemisphere and irregular tumor targets, dose conformity and high/low isodose normal brain volumes are more favorable with the IMRT technique. CONCLUSIONS For the targets described above, the intensity-modulated technique results in improved dose conformity and decreased dose to nontarget brain in high and low isodose regions as compared to the standard noncoplanar arc technique or noncoplanar fixed fields for the hemisphere and tumor targets. Intensity-modulated treatment delivery may allow for an increase in the therapeutic ratio for treating stereotactically defined large and/or irregularly shaped intracranial targets.

Adaptive Replanning Strategies Accounting for Shrinkage in Head and Neck IMRT

PURPOSE Significant anatomic and volumetric changes occur in head and neck cancer patients during fractionated radiotherapy, and the actual dose can be considerably different from the original plan. The purposes of this study were (1) to evaluate the differences between planned and delivered dose, (2) to investigate margins required for anatomic changes, and (3) to find optimal replanning strategies. METHODS AND MATERIALS Eleven patients, each with one planning and six weekly helical CTs, were included. Intensity-modulated radiotherapy plans were generated using the simultaneous integrated boost technique. Weekly CTs were rigidly registered to planning CT before deformable registration was performed. The following replanning strategies were investigated with different margins (0, 3, 5 mm): midcourse (one replan), every other week (two replans), and every week (six replans). Doses were accumulated on the planning CT for comparison of various dose indices for target and critical structures. RESULTS The cumulative doses to targets were preserved even at the 0-mm margin. Doses to cord, brainstem, and mandible were unchanged. Significant increases in parotid doses were observed. Margin reduction from 5 to 0 mm led to a 22% improvement in parotid mean dose. Parotid sparing could be preserved with replanning. More frequent replanning led to better preservation; replanning more than once a week is unnecessary. CONCLUSION Shrinkage does not result in significant dosimetric difference in targets and critical structures, except for the parotid gland, for which the mean dose increases by approximately 10%. The benefit of replanning is improved sparing of the parotid. The combination of replanning and reduced margins can provide up to a 30% difference in parotid dose.

Quantifying the effect of intrafraction motion during breast IMRT planning and dose delivery

Respiratory motion during intensity modulated radiation therapy (IMRT) causes two types of problems. First, the clinical target volume (CTV) to planning target volume (PTV) margin needed to account for respiratory motion means that the lung and heart dose is higher than would occur in the absence of such motion. Second, because respiratory motion is not synchronized with multileaf collimator (MLC) motion, the delivered dose is not the same as the planned dose. The aims of this work were to evaluate these problems to determine (a) the effects of respiratory motion and setup error during breast IMRT treatment planning, (b) the effects of the interplay between respiratory motion and multileaf collimator (MLC) motion during breast IMRT delivery, and (c) the potential benefits of breast IMRT using breath-hold, respiratory gated, and 4D techniques. Seven early stage breast cancer patient data sets were planned for IMRT delivered with a dynamic MLC (DMLC). For each patient case, eight IMRT plans with varying respiratory motion magnitudes and setup errors (and hence CTV to PTV margins) were created. The effects of respiratory motion and setup error on the treatment plan were determined by comparing the eight dose distributions. For each fraction of these plans, the effect of the interplay between respiratory motion and MLC motion during IMRT delivery was simulated by superimposing the respiratory trace on the planned DMLC leaf motion, facilitating comparisons between the planned and expected dose distributions. When considering respiratory motion in the CTV-PTV expansion during breast IMRT planning, our results show that PTV dose heterogeneity increases with respiratory motion. Lung and heart doses also increase with respiratory motion. Due to the interplay between respiratory motion and MLC motion during IMRT delivery, the planned and expected dose distributions differ. This difference increases with respiratory motion. The expected dose varies from fraction to fraction. However, for the seven patients studied and respiratory trace used, for no breathing, shallow breathing, and normal breathing, there were no statistically significant differences between the planned and expected dose distributions. Thus, for breast IMRT, intrafraction motion degrades treatment plans predominantly by the necessary addition of a larger CTV to PTV margin than would be required in the absence of such motion. This motion can be limited by breath-hold, respiratory gated, or 4D techniques.

Intensity-modulated stereotactic radiosurgery using dynamic micro-multileaf collimation.

PURPOSE The implementation of dynamic leaf motion on a micro-multileaf collimator system provides the capability for intensity-modulated stereotactic radiosurgery (IMSRS), and the consequent potential for improved dose distributions for irregularly shaped tumor volumes adjacent to critical organs. This study explores the use of IMSRS to provide improved tumor coverage and normal tissue sparing for small cranial tumors relative to plans based on multiple fixed uniform-intensity beams or traditional circular collimator arc-based stereotactic techniques. METHODS AND MATERIALS Four patient cases involving small brain lesions are presented and analyzed. The cases were chosen to include a representative selection of target shapes, number of targets, and adjacent critical areas. Patient plans generated for these comparisons include standard arcs with multiple circular collimators, and fixed noncoplanar static fields with uniform-intensity beams and IMSRS. Parameters used for evaluation of the plans include the percentage of irradiated volume to tumor volume (PITV), normal tissue dose-volume histograms, and dose-homogeneity ratios. All IMSRS plans were computed using previously established IMRT techniques adapted for use with the BrainLAB M3 micro-multileaf collimator. The algorithms comprising the IMRT system for optimization of intensity distributions and conversion into leaf trajectories of the BrainLab M3 were developed at our institution. The ADAC Pinnacle(3) radiation treatment-planning system was used for dose calculations and for input of contours for target volumes and normal critical structures. RESULTS For all cases, the IMSRS plans showed a high degree of conformity of the dose distribution with the target shape. The IMSRS plans provided either (1) a smaller volume of normal tissue irradiated to significant dose levels, generally taken as doses greater than 50% of the prescription, or (2) a lower dose to an important adjacent critical organ. The reduction in volume of normal tissue irradiated in the IMSRS plans ranged from 10% to 50% relative to the other arc and uniform fixed-field plans. CONCLUSION The case studies presented for IMSRS demonstrate significant dosimetric improvements for small, irregularly shaped lesions of the brain when compared to treatments using multiple static fields or standard SRS arc techniques with circular collimators. For all cases, the IMSRS plan yielded a smaller volume of normal tissue irradiated, and/or a reduction in the volume of an adjacent critical organ (i.e., brainstem) irradiated to significant dose levels.

Algorithm and performance of a clinical IMRT beam-angle optimization system

This paper describes the algorithm and examines the performance of an intensity-modulated radiation therapy (IMRT) beam-angle optimization (BAO) system. In this algorithm successive sets of beam angles are selected from a set of predefined directions using a fast simulated annealing (FSA) algorithm. An IMRT beam-profile optimization is performed on each generated set of beams. The IMRT optimization is accelerated by using a fast dose calculation method that utilizes a precomputed dose kernel. A compact kernel is constructed for each of the predefined beams prior to starting the FSA algorithm. The IMRT optimizations during the BAO are then performed using these kernels in a fast dose calculation engine. This technique allows the IMRT optimization to be performed more than two orders of magnitude faster than a similar optimization that uses a convolution dose calculation engine. Any type of optimization criterion present in the IMRT system can be used in this BAO system. An objective function based on clinically-relevant dose-volume (DV) criteria is used in this study. This facilitates the comparison between a BAO plan and the corresponding plan produced by a planner since the latter is usually optimized using a DV-based objective function. A simple prostate case and a complex head-and-neck (HN) case were used to evaluate the usefulness and performance of this BAO method. For the prostate case we compared the BAO results for three, five and seven coplanar beams with those of the same number of equispaced coplanar beams. For the HN case we compare the BAO results for seven and nine non-coplanar beams with that for nine equispaced coplanar beams. In each case the BAO algorithm was allowed to search up to 1000 different sets of beams. The BAO for the prostate cases were finished in about 1-2 h on a moderate 400 MHz workstation while that for the head-and-neck cases were completed in 13-17 h on a 750 MHz machine. No a priori beam-selection criteria have been used in achieving this performance. In both the prostate and the head-and-neck cases, BAO is shown to provide improvements in plan quality over that of the equispaced beams. The use of DV-based objective function also allows us to study the dependence of the improvement of plan quality offered by BAO on the DV criteria used in the optimization. We found that BAO is especially useful for cases that require strong DV criteria. The main advantages of this BAO system are its speed and its direct link to a clinical IMRT system.

A method for determining multileaf collimator transmission and scatter for dynamic intensity modulated radiotherapy

The main purpose of this work is to demonstrate a practical means of determining the leaf transmission and scatter characteristics of a multileaf collimator (MLC) pertinent to the commissioning of dynamic intensity modulated radiotherapy, especially for the sweeping window technique. The data are necessary for the conversion of intensity distributions produced by intensity-modulated radiotherapy optimization systems into trajectories of MLC leaves for dynamic delivery. Measurements are described for two, tungsten alloy MLCs: a Mark II 80-leaf MLC on a Varian 2100C accelerator and a Millenium 120-leaf MLC on a Varian 2100EX accelerator. MLC leakage was measured by film for a series of field sizes. Measured MLC leakage was 1.68% for a 10 x 10 cm2 field for both 6 and 18 MV for the 80-leaf MLC. For the 6 MV field, the 1.68% leakage consisted of 1.48% direct transmission and 0.20% leaf scatter. Direct transmission through the 80-leaf MLC, including the rounded leaf tip, was calculated analytically taking into account the detailed leaf geometry and a Monte Carlo-generated energy spectrum of the accelerator. The integrated fluence under the leaf tip was equivalent to an inward shift of 0.06 cm of a hypothetical leaf with a flat, focused tip. Monte Carlo calculations of the dose to phantom beyond a closed 80-leaf MLC showed excellent agreement with the analytic results. The transmission depends on the density of the MLC alloy, which may differ among individual MLCs. Thus, it is important to measure the transmission of any particular MLC. Calculated doses for a series of uniform fields produced by dynamic sweeping windows of various widths agree with measurements within 2%.

Multiple local minima in IMRT optimization based on dose-volume criteria

Multiple local minima traps are known to exist in dose-volume and dose-response objective functions. Nevertheless, their presence and consequences are not considered impediments in finding satisfactory solutions in routine optimization of IMRT plans using gradient methods. However, there is often a concern that a significantly superior solution may exist unbeknownst to the planner and that the optimization process may not be able to reach it. We have investigated the soundness of the assumption that the presence of multiple minima traps can be ignored. To find local minima, we start the optimization process a large number of times with random initial intensities. We investigated whether the occurrence of local minima depends upon the choice of the objective function parameters and the number of variables and whether their existence is an impediment in finding a satisfactory solution. To learn about the behavior of multiple minima, we first used a symmetric cubic phantom containing a cubic target and an organ-at-risk surrounding it to optimize the beam weights of two pairs of parallel-opposed beams using a gradient technique. The phantom studies also served to test our software. Objective function parameters were chosen to ensure that multiple minima would exist. Data for 500 plans, optimized with random initial beam weights, were analyzed. The search process did succeed in finding the local minima and showed that the number of minima depends on the parameters of the objective functions. It was also found that the consequences of local minima depended on the number of beams. We further searched for the multiple minima in intensity-modulated treatment plans for a head-and-neck case and a lung case. In addition to the treatment plan scores and the dose-volume histograms, we examined the dose distributions and intensity patterns. We did not find any evidence that multiple local minima affect the outcome of optimization using gradient techniques in any clinically significant way. Our study supports the notion that multiple minima should not be an impediment to finding a good solution when gradient-based optimization techniques are employed. Changing the parameters for the objective function had no observable effect on our findings.

Application of dose compensation in image-guided radiotherapy of prostate cancer

In image-guided radiation therapy (IGRT), volumetric information on patient anatomy at treatment conditions is made available with in-room imaging devices capable of cone-beam CT. Setup error and inter-fraction rigid motion can be corrected online. The planning margin can therefore be reduced significantly. However, to compensate for uncertainties including organ deformation and intra-fraction motion, offline evaluation and replanning are necessary. The purpose of this study is to investigate the use of an offline dose compensation technique to further reduce the margin safely. In IGRT, online CT scan, rigid image registration and setup correction are performed at each fraction. Later the regions of interest are registered offline between treatment and planning CTs using a finite element method to account for non-rigid organ motion. Cumulative dose distribution is calculated and compared with the prescription dose. The discrepancy, if found significant, is repaired using the dose compensation technique, in which the cumulative dose distribution is incorporated in adaptive IMRT planning for future fractions. Two compensation schedules were tested in this study: single compensation at the end of the treatment course and compensation performed weekly. One patient with one planning CT and 16 treatment CTs were used in this simulation study. Due to the aggressive smaller planning margin used, severe underdose was observed in the clinical target volume. The size and magnitude of the underdose were reduced substantially with online guidance but were still significant. Both dose compensation strategies were able to reduce the dose deficit to an acceptable level without additional planning margin. Weekly compensation is more biologically beneficial and can spread the execution error into multiple fractions. The offline dose compensation technique allows further margin reduction and can complement the online guidance by compensating for uncertainties that cannot be reduced online, thereby increasing the confidence in IGRT delivery.