Direct aperture optimization: a turnkey solution for step-and-shoot IMRT

Simultaneous beam sampling and aperture shape optimization for SPORT

PURPOSE

Station parameter optimized radiation therapy (SPORT) was recently proposed to fully utilize the technical capability of emerging digital linear accelerators, in which the station parameters of a delivery system, such as aperture shape and weight, couch position/angle, gantry/collimator angle, can be optimized simultaneously. SPORT promises to deliver remarkable radiation dose distributions in an efficient manner, yet there exists no optimization algorithm for its implementation. The purpose of this work is to develop an algorithm to simultaneously optimize the beam sampling and aperture shapes.

METHODS

The authors build a mathematical model with the fundamental station point parameters as the decision variables. To solve the resulting large-scale optimization problem, the authors devise an effective algorithm by integrating three advanced optimization techniques: column generation, subgradient method, and pattern search. Column generation adds the most beneficial stations sequentially until the plan quality improvement saturates and provides a good starting point for the subsequent optimization. It also adds the new stations during the algorithm if beneficial. For each update resulted from column generation, the subgradient method improves the selected stations locally by reshaping the apertures and updating the beam angles toward a descent subgradient direction. The algorithm continues to improve the selected stations locally and globally by a pattern search algorithm to explore the part of search space not reachable by the subgradient method. By combining these three techniques together, all plausible combinations of station parameters are searched efficiently to yield the optimal solution.

RESULTS

A SPORT optimization framework with seamlessly integration of three complementary algorithms, column generation, subgradient method, and pattern search, was established. The proposed technique was applied to two previously treated clinical cases: a head and neck and a prostate case. It significantly improved the target conformality and at the same time critical structure sparing compared with conventional intensity modulated radiation therapy (IMRT). In the head and neck case, for example, the average PTV coverage D99% for two PTVs, cord and brainstem max doses, and right parotid gland mean dose were improved, respectively, by about 7%, 37%, 12%, and 16%.

CONCLUSIONS

The proposed method automatically determines the number of the stations required to generate a satisfactory plan and optimizes simultaneously the involved station parameters, leading to improved quality of the resultant treatment plans as compared with the conventional IMRT plans.

Total-variation regularization based inverse planning for intensity modulated arc therapy

Intensity modulated arc therapy (IMAT) delivers conformal dose distributions through continuous gantry rotation with constant or variable speed while modulating the field aperture shape and weight. The enlarged angular space and machine delivery constraints make inverse planning of IMAT more intractable as compared to its counterpart of fixed gantry IMRT. Currently, IMAT inverse planning is being done using two extreme methods: the first one computes in beamlet domain with a subsequent arc leaf sequencing, and the second proceeds in machine parameter domain with entire emphasis placed on a pre-determined delivery method without exploring potentially better alternative delivery schemes. Towards truly optimizing the IMAT treatment on a patient specific basis, in this work we propose a total-variation based inverse planning framework for IMAT, which takes advantage of the useful features of the above two existing approaches while avoiding their shortcomings. A quadratic optimization algorithm has been implemented to demonstrate the performance and advantage of the proposed approach. Applications of the technique to a prostate case and a head and neck case indicate that the algorithm is capable of generating IMAT plans with patient specific numbers of arcs efficiently. Superior dose distributions and delivery time are achieved with a maximum number of apertures of three for each field. As compared to conventional beamlet-based algorithms, our method regularizes the field modulation complexity during optimization, and permits us to obtain the best possible plan with a pre-set modulation complexity of fluences. As illustrated in both prostate and head-and-neck case studies, the proposed method produces more favorable dose distributions than the segment-based algorithms, by optimally accommodating the clinical need of intensity modulation levels for each individual field. On a more fundamental level, our formulation preserves the convexity of optimization and makes the search of the global optimal solution possible with a deterministic method.

Optimization approaches to volumetric modulated arc therapy planning

Volumetric modulated arc therapy (VMAT) has found widespread clinical application in recent years. A large number of treatment planning studies have evaluated the potential for VMAT for different disease sites based on the currently available commercial implementations of VMAT planning. In contrast, literature on the underlying mathematical optimization methods used in treatment planning is scarce. VMAT planning represents a challenging large scale optimization problem. In contrast to fluence map optimization in intensity-modulated radiotherapy planning for static beams, VMAT planning represents a nonconvex optimization problem. In this paper, the authors review the state-of-the-art in VMAT planning from an algorithmic perspective. Different approaches to VMAT optimization, including arc sequencing methods, extensions of direct aperture optimization, and direct optimization of leaf trajectories are reviewed. Their advantages and limitations are outlined and recommendations for improvements are discussed.

Shared data for intensity modulated radiation therapy (IMRT) optimization research: the CORT dataset

BACKGROUND

We provide common datasets (which we call the CORT dataset: common optimization for radiation therapy) that researchers can use when developing and contrasting radiation treatment planning optimization algorithms. The datasets allow researchers to make one-to-one comparisons of algorithms in order to solve various instances of the radiation therapy treatment planning problem in intensity modulated radiation therapy (IMRT), including beam angle optimization, volumetric modulated arc therapy and direct aperture optimization.

RESULTS

We provide datasets for a prostate case, a liver case, a head and neck case, and a standard IMRT phantom. We provide the dose-influence matrix from a variety of beam/couch angle pairs for each dataset. The dose-influence matrix is the main entity needed to perform optimizations: it contains the dose to each patient voxel from each pencil beam. In addition, the original Digital Imaging and Communications in Medicine (DICOM) computed tomography (CT) scan, as well as the DICOM structure file, are provided for each case.

CONCLUSIONS

Here we present an open dataset - the first of its kind - to the radiation oncology community, which will allow researchers to compare methods for optimizing radiation dose delivery.

Dosimetric Evaluation of Different Intensity-Modulated Radiotherapy Techniques for Breast Cancer After Conservative Surgery

Intensity-modulated radiotherapy (IMRT) potentially leads to a more favorite dose distribution compared to 3-dimensional or conventional tangential radiotherapy (RT) for breast cancer after conservative surgery or mastectomy. The aim of this study was to compare dosimetric parameters of the planning target volume (PTV) and organs at risk (OARs) among helical tomotherapy (HT), inverse-planned IMRT (IP-IMRT), and forward-planned field in field (FP-FIF) IMRT techniques after breast-conserving surgery. Computed tomography scans from 20 patients (12 left sided and 8 right sided) previously treated with T1N0 carcinoma were selected for this dosimetric planning study. We designed HT, IP-IMRT, and FP-FIF plans for each patient. Plans were compared according to dose-volume histogram analysis in terms of PTV homogeneity and conformity indices (HI and CI) as well as OARs dose and volume parameters. Both HI and CI of the PTV showed statistically significant difference among IP-IMRT, FP-FIF, and HT with those of HT were best (P < .05). Compared to FP-FIF, IP-IMRT showed smaller exposed volumes of ipsilateral lung, heart, contralateral lung, and breast, while HT indicated smaller exposed volumes of ipsilateral lung but larger exposed volumes of contralateral lung and breast as well as heart. In addition, HT demonstrated an increase in exposed volume of ipsilateral lung (except for fraction of lung volume receiving >30 Gy and 20 Gy), heart, contralateral lung, and breast compared with IP-IMRT. For breast cancer radiotherapy (RT) after conservative surgery, HT provides better dose homogeneity and conformity of PTV compared to IP-IMRT and FP-FIF techniques, especially for patients with supraclavicular lymph nodes involved. Meanwhile, HT decreases the OAR volumes receiving higher doses with an increase in the volumes receiving low doses, which is known to lead to an increased rate of radiation-induced secondary malignancies. Hence, composite factors including dosimetric advantage, clinical effect, and economic burden should be taken into comprehensive consideration when choosing an RT technique in clinical practice.

Direct leaf trajectory optimization for volumetric modulated arc therapy planning with sliding window delivery

PURPOSE

The authors propose a novel optimization model for volumetric modulated arc therapy (VMAT) planning that directly optimizes deliverable leaf trajectories in the treatment plan optimization problem, and eliminates the need for a separate arc-sequencing step.

METHODS

In this model, a 360° arc is divided into a given number of arc segments in which the leaves move unidirectionally. This facilitates an algorithm that determines the optimal piecewise linear leaf trajectories for each arc segment, which are deliverable in a given treatment time. Multileaf collimator constraints, including maximum leaf speed and interdigitation, are accounted for explicitly. The algorithm is customized to allow for VMAT delivery using constant gantry speed and dose rate, however, the algorithm generalizes to variable gantry speed if beneficial.

RESULTS

The authors demonstrate the method for three different tumor sites: a head-and-neck case, a prostate case, and a paraspinal case. The authors first obtain a reference plan for intensity modulated radiotherapy (IMRT) using fluence map optimization and 20 intensity-modulated fields in equally spaced beam directions, which is beyond the standard of care. Modeling the typical clinical setup for the treatment sites considered, IMRT plans using seven or nine beams are also computed. Subsequently, VMAT plans are optimized by dividing the 360° arc into 20 corresponding arc segments. Assuming typical machine parameters (a dose rate of 600 MU/min, and a maximum leaf speed of 3 cm/s), it is demonstrated that the optimized VMAT plans with 2-3 min delivery time are of noticeably better quality than the 7-9 beam IMRT plans. The VMAT plan quality approaches the quality of the 20-beam IMRT benchmark plan for delivery times between 3 and 4 min.

CONCLUSIONS

The results indicate that high quality treatments can be delivered in a single arc with 20 arc segments if sufficient time is allowed for modulation in each segment.

MCTP system model based on linear programming optimization of apertures obtained from sequencing patient image data maps

PURPOSE

The authors present a hybrid direct multileaf collimator (MLC) aperture optimization model exclusively based on sequencing of patient imaging data to be implemented on a Monte Carlo treatment planning system (MC-TPS) to allow the explicit radiation transport simulation of advanced radiotherapy treatments with optimal results in efficient times for clinical practice.

METHODS

The planning system (called CARMEN) is a full MC-TPS, controlled through aMATLAB interface, which is based on the sequencing of a novel map, called "biophysical" map, which is generated from enhanced image data of patients to achieve a set of segments actually deliverable. In order to reduce the required computation time, the conventional fluence map has been replaced by the biophysical map which is sequenced to provide direct apertures that will later be weighted by means of an optimization algorithm based on linear programming. A ray-casting algorithm throughout the patient CT assembles information about the found structures, the mass thickness crossed, as well as PET values. Data are recorded to generate a biophysical map for each gantry angle. These maps are the input files for a home-made sequencer developed to take into account the interactions of photons and electrons with the MLC. For each linac (Axesse of Elekta and Primus of Siemens) and energy beam studied (6, 9, 12, 15 MeV and 6 MV), phase space files were simulated with the EGSnrc/BEAMnrc code. The dose calculation in patient was carried out with the BEAMDOSE code. This code is a modified version of EGSnrc/DOSXYZnrc able to calculate the beamlet dose in order to combine them with different weights during the optimization process.

RESULTS

Three complex radiotherapy treatments were selected to check the reliability of CARMEN in situations where the MC calculation can offer an added value: A head-and-neck case (Case I) with three targets delineated on PET/CT images and a demanding dose-escalation; a partial breast irradiation case (Case II) solved with photon and electron modulated beams (IMRT + MERT); and a prostatic bed case (Case III) with a pronounced concave-shaped PTV by using volumetric modulated arc therapy. In the three cases, the required target prescription doses and constraints on organs at risk were fulfilled in a short enough time to allow routine clinical implementation. The quality assurance protocol followed to check CARMEN system showed a high agreement with the experimental measurements.

CONCLUSIONS

A Monte Carlo treatment planning model exclusively based on maps performed from patient imaging data has been presented. The sequencing of these maps allows obtaining deliverable apertures which are weighted for modulation under a linear programming formulation. The model is able to solve complex radiotherapy treatments with high accuracy in an efficient computation time.

Investigation of pulsed IMRT and VMAT for re-irradiation treatments: dosimetric and delivery feasibilities

Many tumor cells demonstrate hyperradiosensitivity at doses below ~50 cGy. Together with the increased normal tissue repair under low dose rate, the pulsed low dose rate radiotherapy (PLDR), which separates a daily fractional dose of 200 cGy into 10 pulses with 3 min interval between pulses (~20 cGy/pulse and effective dose rate 6.7 cGy min−1), potentially reduces late normal tissue toxicity while still providing significant tumor control for re-irradiation treatments. This work investigates the dosimetric and technical feasibilities of intensity modulated radiotherapy (IMRT) and volumetric modulated arc therapy (VMAT)-based PLDR treatments using Varian Linacs. Twenty one cases (12 real re-irradiation cases) including treatment sites of pancreas, prostate, pelvis, lung, head-and-neck, and breast were recruited for this study. The lowest machine operation dose rate (100 MU min−1) was employed in the plan delivery. Ten-field step-and-shoot IMRT and dual-arc VMAT plans were generated using the Eclipse TPS with routine planning strategies. The dual-arc plans were delivered five times to achieve a 200 cGy daily dose (~20 cGy arc−1). The resulting plan quality was evaluated according to the heterogeneity and conformity indexes (HI and CI) of the planning target volume (PTV). The dosimetric feasibility of retaining the hyperradiosensitivity for PLDR was assessed based on the minimum and maximum dose in the target volume from each pulse. The delivery accuracy of VMAT and IMRT at the 100 MU min−1 machine operation dose rate was verified using a 2D diode array and ion chamber measurements. The delivery reproducibility was further investigated by analyzing the Dynalog files of repeated deliveries. A comparable plan quality was achieved by the IMRT (CI 1.10–1.38; HI 1.04–1.10) and the VMAT (CI 1.08–1.26; HI 1.05–1.10) techniques. The minimum/maximum PTV dose per pulse is 7.9 ± 5.1 cGy/33.7 ± 6.9 cGy for the IMRT and 12.3 ± 4.1 cGy/29.2 ± 4.7 cGy for the VMAT. Six out of the 186 IMRT pulses (fields) were found to exceed 50 cGy maximum PTV dose per pulse while the maximum PTV dose per pulse was within 40 cGy for all the VMAT pulses (arcs). However, for VMAT plans, the dosimetric quality of the entire treatment plan was less superior for the breast cases and large irregular targets. The gamma passing rates for both techniques at the 100 MU min−1 dose rate were at least 94.1% (3%/3 mm) and the point dose measurements agreed with the planned values to within 2.2%. The average root mean square error of the leaf position was 0.93 ± 0.83 mm for IMRT and 0.53 ± 0.48 mm for VMAT based on the Dynalog file analysis. The RMS error of the leaf position was nearly identical for the repeated deliveries of the same plans. In general, both techniques are feasible for PLDR treatments. VMAT was more advantageous for PLDR with more uniform target dose per pulse, especially for centrally located tumors. However, for large, irregular and/or peripheral tumors, IMRT could produce more favorable PLDR plans. By taking the biological benefit of PLDR delivery and the dosimetric benefit of IMRT and VMAT, the proposed methods have a great potential for those previously-irradiated recurrent patients.

VMAT monthly QA using two techniques: 2D ion chamber array with an isocentric gantry mount and an in vivo dosimetric device attached to gantry

Purpose

Varian RapidArc is a volumetric modulated arc therapy (VMAT) that obtains a conformal dose around the desired structure by employing variable gantry speed, dose rate and dynamic multileaf collimator (DMLC) speed as the gantry rotates about machine isocenter. This study is meant to build upon previous research by Ling et al. by completing the tests with an in vivo dosimetric device attached to the linac gantry and a 2D ionisation chamber array with an isocentric gantry mount.

Materials and methods

Two PTW detectors, seven29 array with gantry mount and DAVID, were attached to the linear accelerator gantry, allowing each device to remain perpendicular to the beam at all gantry angles. Three tests for RapidArc evaluation were performed on these devices including: dose rate and gantry speed variation, DMLC speed and dose rate variation and DMLC position accuracy. The reproducibility of the arc data was also reported.

Results

A picket fence plan varying dose rates (111 to 600 MU/minute) and gantry speeds (5·5 to 4·3°/second) was delivered consisting of seven sections of different combinations. These measurements were compared with static gantry, open field measurements and found to be within 2·39% for the DAVID device and 0·84% for the seven29. A four-section picket fence of varying DMLC speeds (0·46, 0·92, 1·84 and 2·76 cm/second) was similarly evaluated and found to be within 1·99% and 3·66% for the DAVID and seven29, respectively. For DMLC position accuracy, a picket fence arc plan was compared with a static picket fence and found to agree within 0.38% and 2.91%. Reproducibility for these three RapidArc plans was found to be within 0·30% and 2·70% for the DAVID and seven29.

Conclusion

The DAVID and seven29 detectors were able to perform the RapidArc quality assurance tests efficiently and accurately and the results were reproducible. Periodic verification of DMLC movement, dose rate variation and gantry speed variation relating to RapidArc delivery can be completed in a timelier manner using this equipment.

Advantages and limitations of navigation-based multicriteria optimization (MCO) for localized prostate cancer IMRT planning

Efficacy of inverse planning is becoming increasingly important for advanced radiotherapy techniques. This study's aims were to validate multicriteria optimization (MCO) in RayStation (v2.4, RaySearch Laboratories, Sweden) against standard intensity-modulated radiation therapy (IMRT) optimization in Oncentra (v4.1, Nucletron BV, the Netherlands) and characterize dose differences due to conversion of navigated MCO plans into deliverable multileaf collimator apertures. Step-and-shoot IMRT plans were created for 10 patients with localized prostate cancer using both standard optimization and MCO. Acceptable standard IMRT plans with minimal average rectal dose were chosen for comparison with deliverable MCO plans. The trade-off was, for the MCO plans, managed through a user interface that permits continuous navigation between fluence-based plans. Navigated MCO plans were made deliverable at incremental steps along a trajectory between maximal target homogeneity and maximal rectal sparing. Dosimetric differences between navigated and deliverable MCO plans were also quantified. MCO plans, chosen as acceptable under navigated and deliverable conditions resulted in similar rectal sparing compared with standard optimization (33.7 ± 1.8 Gy vs 35.5 ± 4.2 Gy, p = 0.117). The dose differences between navigated and deliverable MCO plans increased as higher priority was placed on rectal avoidance. If the best possible deliverable MCO was chosen, a significant reduction in rectal dose was observed in comparison with standard optimization (30.6 ± 1.4 Gy vs 35.5 ± 4.2 Gy, p = 0.047). Improvements were, however, to some extent, at the expense of less conformal dose distributions, which resulted in significantly higher doses to the bladder for 2 of the 3 tolerance levels. In conclusion, similar IMRT plans can be created for patients with prostate cancer using MCO compared with standard optimization. Limitations exist within MCO regarding conversion of navigated plans to deliverable apertures, particularly for plans that emphasize avoidance of critical structures. Minimizing these differences would result in better quality treatments for patients with prostate cancer who were treated with radiotherapy using MCO plans.

Direct aperture optimization using an inverse form of back-projection

Direct aperture optimization (DAO) has been used to produce high dosimetric quality intensity-modulated radiotherapy (IMRT) treatment plans with fast treatment delivery by directly modeling the multileaf collimator segment shapes and weights. To improve plan quality and reduce treatment time for our in-house treatment planning system, we implemented a new DAO approach without using a global objective function (GFO). An index concept is introduced as an inverse form of back-projection used in the CT multiplicative algebraic reconstruction technique (MART). The index, introduced for IMRT optimization in this work, is analogous to the multiplicand in MART. The index is defined as the ratio of the optima over the current. It is assigned to each voxel and beamlet to optimize the fluence map. The indices for beamlets and segments are used to optimize multileaf collimator (MLC) segment shapes and segment weights, respectively. Preliminary data show that without sacrificing dosimetric quality, the implementation of the DAO reduced average IMRT treatment time from 13 min to 8 min for the prostate, and from 15 min to 9 min for the head and neck using our in-house treatment planning system PlanUNC. The DAO approach has also shown promise in optimizing rotational IMRT with burst mode in a head and neck test case.

Improving IMRT delivery efficiency with reweighted L1-minimization for inverse planning

PURPOSE

This study presents an improved technique to further simplify the fluence-map in intensity modulated radiation therapy (IMRT) inverse planning, thereby reducing plan complexity and improving delivery efficiency, while maintaining the plan quality.

METHODS

First-order total-variation (TV) minimization (min.) based on L1-norm has been proposed to reduce the complexity of fluence-map in IMRT by generating sparse fluence-map variations. However, with stronger dose sparing to the critical structures, the inevitable increase in the fluence-map complexity can lead to inefficient dose delivery. Theoretically, L0-min. is the ideal solution for the sparse signal recovery problem, yet practically intractable due to its nonconvexity of the objective function. As an alternative, the authors use the iteratively reweighted L1-min. technique to incorporate the benefits of the L0-norm into the tractability of L1-min. The weight multiplied to each element is inversely related to the magnitude of the corresponding element, which is iteratively updated by the reweighting process. The proposed penalizing process combined with TV min. further improves sparsity in the fluence-map variations, hence ultimately enhancing the delivery efficiency. To validate the proposed method, this work compares three treatment plans obtained from quadratic min. (generally used in clinic IMRT), conventional TV min., and our proposed reweighted TV min. techniques, implemented by a large-scale L1-solver (template for first-order conic solver), for five patient clinical data. Criteria such as conformation number (CN), modulation index (MI), and estimated treatment time are employed to assess the relationship between the plan quality and delivery efficiency.

RESULTS

The proposed method yields simpler fluence-maps than the quadratic and conventional TV based techniques. To attain a given CN and dose sparing to the critical organs for 5 clinical cases, the proposed method reduces the number of segments by 10-15 and 30-35, relative to TV min. and quadratic min. based plans, while MIs decreases by about 20%-30% and 40%-60% over the plans by two existing techniques, respectively. With such conditions, the total treatment time of the plans obtained from our proposed method can be reduced by 12-30 s and 30-80 s mainly due to greatly shorter multileaf collimator (MLC) traveling time in IMRT step-and-shoot delivery.

CONCLUSIONS

The reweighted L1-minimization technique provides a promising solution to simplify the fluence-map variations in IMRT inverse planning. It improves the delivery efficiency by reducing the entire segments and treatment time, while maintaining the plan quality in terms of target conformity and critical structure sparing.

Trajectory optimization for dynamic couch rotation during volumetric modulated arc radiotherapy

Non-coplanar radiation beams are often used in three-dimensional conformal and intensity modulated radiotherapy to reduce dose to organs at risk (OAR) by geometric avoidance. In volumetric modulated arc radiotherapy (VMAT) non-coplanar geometries are generally achieved by applying patient couch rotations to single or multiple full or partial arcs. This paper presents a trajectory optimization method for a non-coplanar technique, dynamic couch rotation during VMAT (DCR–VMAT), which combines ray tracing with a graph search algorithm. Four clinical test cases (partial breast, brain, prostate only, and prostate and pelvic nodes) were used to evaluate the potential OAR sparing for trajectory-optimized DCR–VMAT plans, compared with standard coplanar VMAT. In each case, ray tracing was performed and a cost map reflecting the number of OAR voxels intersected for each potential source position was generated. The least-cost path through the cost map, corresponding to an optimal DCR–VMAT trajectory, was determined using Dijkstra's algorithm. Results show that trajectory optimization can reduce dose to specified OARs for plans otherwise comparable to conventional coplanar VMAT techniques. For the partial breast case, the mean heart dose was reduced by 53%. In the brain case, the maximum lens doses were reduced by 61% (left) and 77% (right) and the globes by 37% (left) and 40% (right). Bowel mean dose was reduced by 15% in the prostate only case. For the prostate and pelvic nodes case, the bowel V50 Gy and V60 Gy were reduced by 9% and 45% respectively. Future work will involve further development of the algorithm and assessment of its performance over a larger number of cases in site-specific cohorts.

Deliverable navigation for multicriteria IMRT treatment planning by combining shared and individual apertures

We consider the problem of deliverable Pareto surface navigation for step-and-shoot intensity-modulated radiation therapy. This problem amounts to calculation of a collection of treatment plans with the property that convex combinations of plans are directly deliverable. Previous methods for deliverable navigation impose restrictions on the number of apertures of the individual plans, or require that all treatment plans have identical apertures. We introduce simultaneous direct step-and-shoot optimization of multiple plans subject to constraints that some of the apertures must be identical across all plans. This method generalizes previous methods for deliverable navigation to allow for treatment plans with some apertures from a collective pool and some apertures that are individual. The method can also be used as a post-processing step to previous methods for deliverable navigation in order to improve upon their plans. By applying the method to subsets of plans in the collection representing the Pareto set, we show how it can enable convergence toward the unrestricted (non-navigable) Pareto set where all apertures are individual.

Relationship of segment area and monitor unit efficiency in aperture-based IMRT optimization

In step‐and‐shoot IMRT plans, aperture‐based optimization (or one‐step optimization) has been considered as a means of improving monitor unit (MU) efficiency compared to fluence‐based optimization (or two‐step optimization). However, the extent of improvement on MU efficiency varies, depending on the implementation and design of one‐step optimization. In this paper, we attempted to investigate MU efficiency issue in two methods of one‐step optimization implemented in two commercial treatment planning systems (TPSs). Five patients with nasopharyngeal cancer and five patients with advanced prostate cancer were selected for this study. For these patients, clinically used IMRT plans were generated using the Direct Machine Parameter Optimization (DMPO) in the Pinnacle TPS. New IMRT plans were created using the Direct Aperture Optimization (DAO) method in the Panther TPS. For the purpose of this study, we used the similar planning dose objectives and beam configurations with a similar total number of segments in each pair of DMPO and DAO plans. With similar plan quality, DMPO plans required more MUs than DAO plans. The average number of MUs (expressed in mean ±1 SD) for the DMPO and DAO plans was 1,169±186 and 671±135 for the nasopharynx cases, and 711±48 and 400±65 for the prostate cases, respectively. The average segment areas (expressed in mean ±1 SD) for the DMPO plans were smaller than those for the DAO plans: 46.0±7.6 cm2 vs. 100.9±32.3 cm2 for the nasopharynx cases, and 58.3±17.2 cm2 vs. 97.4±35.0 cm2 for the prostate cases, respectively. In conclusion, two one‐step optimization algorithms, DMPO and DAO, resulted in much different MU efficiency with the similar number of segments and optimization parameters. This MU difference is largely attributed to the fact that large area segments are used more often in DAO plans than in DMPO plans.

PACS number: 87.55.de

Segment edit and segment weight optimization: two techniques for intensity modulated radiation therapy and their application to the planning for nasopharyngeal carcinoma

The purpose of this study was to evaluate the two functions: segment weight optimization (SWO) and segment edit (SE) in the latest XiO 4.7 radiation treatment planning system and their effect on the planning of intensity modulated radiation therapy (IMRT) for Nasopharyngeal Carcinoma (NPC). SWO first appeared in XiO 4.5 and SE in XiO 4.7. Twelve patients with NPC were selected and there were three plans for each patient: the common step-and-shoot IMRT plan (C-IMRT); S-IMRT was based on the result of C-IMRT and the plan was further optimized with SWO; F-IMRT was based on S-IMRT and the segments were edited for lowering the dose received by normal tissues. The paired plans were analyzed by comparing the total number of segments, monitor units, the homogeneity index and conformity index of the target volumes and the dose delivered to organs at risk (OAR) including spinal cord, brain stem, optic nerves, chiasm, parotids and larynx. The study exhibited that the total number of segments and monitor units of S-IMRT and F-IMRT were around 25.3%, 3.4% less than those of C-IMRT respectively. The HI and CI indexes of target volumes among three kinds of plans did not show the significant difference. The doses received by spinal cord, brain stem, parotids, larynx were decreased at S-IMRT and F-IMRT as compared to C-IMRT; the highest doses delivered to chiasm and optic nerves were S-IMRT, the next C-IMRT, the lowest F-IMRT. This study showed that the SWO function could substantially reduce the total number of segments of step-and-shoot IMRT plans and the SE function had the incredible ability to lower the dose received by normal tissues.

Tomotherapy-like versus VMAT-like treatments: a multicriteria comparison for a prostate geometry

PURPOSE

To perform a methodological comparison of volumetric modulated arc therapy (VMAT)-like and tomotherapy-like techniques for a prostate geometry, exploring the dependence on machine, delivery, and optimization parameters of cost function values optimized for each technique.

METHODS

A gradient-descent algorithm is used to optimize tomotherapy-like treatments, while VMAT-like optimization is carried out using a direct-aperture simulated annealing algorithm with 180 control points equispaced at 2° angles. Dose distributions are linked to fluences via a three-dimensional double-gaussian pencil beam model. Plans are optimized for a prostate geometry, outlined according to the CHHiP protocol. The cost function used for optimization contains ten simple functions, each of which describes a single planning objective. These functions are split into three structure groups according to whether they are used to control PTV, rectal or bladder dose levels. Different optimizations have been performed by varying the relative weights of each of these structure groups, exploring in this way a three-dimensional Pareto front. Plan quality is studied according to the value of the optimized cost function and the relative Euclidean distance between the components of the cost function and those of the nearest plan lying on a reference Pareto front obtained for tomotherapy-like plans generated using a 1 cm fan-beam width and 1/3 pitch.

RESULTS

The quality of tomotherapy-like optimization depends on the fan-beam width, s, and rotation pitch, p, used to deliver the treatment. These values together define the effective longitudinal resolution with which fluence can be modulated, and lower cost function values are obtained for treatments optimized with tighter pitches and narrower fan-beam widths (higher modulation resolution). On the other hand, the cost function values of VMAT-like optimizations depends on the optimization running time, leaf displacement constraints, and number of arcs employed, as well as on the size of the beamlets used in the optimization (a change in leaf width from 5 to 10 mm clearly worsens the value of the objective function, but only a marginal improvement is observed when the leaf movement discretization step is reduced from 5 to 5/3 mm). However, for no combination of these parameter values did VMAT-like optimizations match the cost function values of optimized tomo-like plans obtained for s = 1 cm and p = 1∕3 (or 1/2). This is the case all across the Pareto front. On the other hand, cost function values of VMAT-like plans are generally lower than those of optimized tomotherapy-like plans obtained for s = 2.5 cm.

CONCLUSIONS

Tomotherapy-like plans created for the prostate geometry using a 1 cm fan-beam width and pitches of 1/3 or 1/2 have lower cost function values than VMAT-like plans, although the associated dosimetric improvements are quite small, both techniques generating very good dose distributions. When a 2.5 cm wide fan-beam is used for tomotherapy-like treatments the pattern is reversed, the tomotherapy-like plans having higher cost functions than the VMAT-like ones.

A column-generation-based method for multi-criteria direct aperture optimization

Navigation-based multi-criteria optimization has been introduced to radiotherapy planning in order to allow the interactive exploration of trade-offs between conflicting clinical goals. However, this has been mainly applied to fluence map optimization. The subsequent leaf sequencing step may cause dose discrepancy, leading to human iteration loops in the treatment planning process that multi-criteria methods were meant to avoid. To circumvent this issue, this paper investigates the application of direct aperture optimization methods in the context of multi-criteria optimization. We develop a solution method to directly obtain a collection of apertures that can adequately span the entire Pareto surface. To that end, we extend the column generation method for direct aperture optimization to a multi-criteria setting in which apertures that can improve the entire Pareto surface are sequentially identified and added to the treatment plan. Our proposed solution method can be embedded in a navigation-based multi-criteria optimization framework, in which the treatment planner explores the trade-off between treatment objectives directly in the space of deliverable apertures. Our solution method is demonstrated for a paraspinal case where the trade-off between target coverage and spinal-cord sparing is studied. The computational results validate that our proposed method obtains a balanced approximation of the Pareto surface over a wide range of clinically relevant plans.

Aperture shape optimization for IMRT treatment planning

We propose an algorithm for aperture shape optimization (ASO) for step and shoot delivery of intensity-modulated radiotherapy. The method is an approach to direct aperture optimization (DAO) that exploits gradient information to locally optimize the positions of the leafs of a multileaf collimator. Based on the dose-influence matrix, the dose distribution is locally approximated as a linear function of the leaf positions. Since this approximation is valid only in a small interval around the current leaf positions, we use a trust-region-like method to optimize the leaf positions: in one iteration, the leaf motion is confined to the beamlets where the leaf edges are currently positioned. This yields a well-behaved optimization problem for the leaf positions and the aperture weights, which can be solved efficiently. If, in one iteration, a leaf is moved to the edge of a beamlet, the leaf motion can be confined to the neighboring beamlet in the next iteration. This allows for large leaf position changes over the course of the algorithm. In this paper, the ASO algorithm is embedded into a column-generation approach to DAO. After a new aperture is added to the treatment plan, we use the ASO algorithm to simultaneously optimize aperture weights and leaf positions for the new set of apertures. We present results for a paraspinal tumor case, a prostate case and a head and neck case. The computational results indicate that, using this approach, treatment plans close to the ideal fluence map optimization solution can be obtained.

Deliverable navigation for multicriteria step and shoot IMRT treatment planning

We consider Pareto surface based multi-criteria optimization for step and shoot IMRT planning. By analyzing two navigation algorithms, we show both theoretically and in practice that the number of plans needed to form convex combinations of plans during navigation can be kept small (much less than the theoretical maximum number needed in general, which is equal to the number of objectives for on-surface Pareto navigation). Therefore a workable approach for directly deliverable navigation in this setting is to segment the underlying Pareto surface plans and then enforce the mild restriction that only a small number of these plans are active at any time during plan navigation, thus limiting the total number of segments used in the final plan.

Improving IMRT-plan quality with MLC leaf position refinement post plan optimization

PURPOSE

In intensity-modulated radiation therapy (IMRT) planning, reducing the pencil-beam size may lead to a significant improvement in dose conformity, but also increase the time needed for the dose calculation and plan optimization. The authors develop and evaluate a postoptimization refinement (POpR) method, which makes fine adjustments to the multileaf collimator (MLC) leaf positions after plan optimization, enhancing the spatial precision and improving the plan quality without a significant impact on the computational burden.

METHODS

The authors' POpR method is implemented using a commercial treatment planning system based on direct aperture optimization. After an IMRT plan is optimized using pencil beams with regular pencil-beam step size, a greedy search is conducted by looping through all of the involved MLC leaves to see if moving the MLC leaf in or out by half of a pencil-beam step size will improve the objective function value. The half-sized pencil beams, which are used for updating dose distribution in the greedy search, are derived from the existing full-sized pencil beams without need for further pencil-beam dose calculations. A benchmark phantom case and a head-and-neck (HN) case are studied for testing the authors' POpR method.

RESULTS

Using a benchmark phantom and a HN case, the authors have verified that their POpR method can be an efficient technique in the IMRT planning process. Effectiveness of POpR is confirmed by noting significant improvements in objective function values. Dosimetric benefits of POpR are comparable to those of using a finer pencil-beam size from the optimization start, but with far less computation and time.

CONCLUSIONS

The POpR is a feasible and practical method to significantly improve IMRT-plan quality without compromising the planning efficiency.

Dose optimization with first-order total-variation minimization for dense angularly sampled and sparse intensity modulated radiation therapy (DASSIM-RT)

PURPOSE

A new treatment scheme coined as dense angularly sampled and sparse intensity modulated radiation therapy (DASSIM-RT) has recently been proposed to bridge the gap between IMRT and VMAT. By increasing the angular sampling of radiation beams while eliminating dispensable segments of the incident fields, DASSIM-RT is capable of providing improved conformity in dose distributions while maintaining high delivery efficiency. The fact that DASSIM-RT utilizes a large number of incident beams represents a major computational challenge for the clinical applications of this powerful treatment scheme. The purpose of this work is to provide a practical solution to the DASSIM-RT inverse planning problem.

METHODS

The inverse planning problem is formulated as a fluence-map optimization problem with total-variation (TV) minimization. A newly released L1-solver, template for first-order conic solver (TFOCS), was adopted in this work. TFOCS achieves faster convergence with less memory usage as compared with conventional quadratic programming (QP) for the TV form through the effective use of conic forms, dual-variable updates, and optimal first-order approaches. As such, it is tailored to specifically address the computational challenges of large-scale optimization in DASSIM-RT inverse planning. Two clinical cases (a prostate and a head and neck case) are used to evaluate the effectiveness and efficiency of the proposed planning technique. DASSIM-RT plans with 15 and 30 beams are compared with conventional IMRT plans with 7 beams in terms of plan quality and delivery efficiency, which are quantified by conformation number (CN), the total number of segments and modulation index, respectively. For optimization efficiency, the QP-based approach was compared with the proposed algorithm for the DASSIM-RT plans with 15 beams for both cases.

RESULTS

Plan quality improves with an increasing number of incident beams, while the total number of segments is maintained to be about the same in both cases. For the prostate patient, the conformation number to the target was 0.7509, 0.7565, and 0.7611 with 80 segments for IMRT with 7 beams, and DASSIM-RT with 15 and 30 beams, respectively. For the head and neck (HN) patient with a complicated target shape, conformation numbers of the three treatment plans were 0.7554, 0.7758, and 0.7819 with 75 segments for all beam configurations. With respect to the dose sparing to the critical structures, the organs such as the femoral heads in the prostate case and the brainstem and spinal cord in the HN case were better protected with DASSIM-RT. For both cases, the delivery efficiency has been greatly improved as the beam angular sampling increases with the similar or better conformal dose distribution. Compared with conventional quadratic programming approaches, first-order TFOCS-based optimization achieves far faster convergence and smaller memory requirements in DASSIM-RT.

CONCLUSIONS

The new optimization algorithm TFOCS provides a practical and timely solution to the DASSIM-RT or other inverse planning problem requiring large memory space. The new treatment scheme is shown to outperform conventional IMRT in terms of dose conformity to both the targetand the critical structures, while maintaining high delivery efficiency.

Evaluation of the efficiency and effectiveness of independent dose calculation followed by machine log file analysis against conventional measurement based IMRT QA

Experimental methods are commonly used for patient-specific IMRT delivery verification. There are a variety of IMRT QA techniques which have been proposed and clinically used with a common understanding that not one single method can detect all possible errors. The aim of this work was to compare the efficiency and effectiveness of independent dose calculation followed by machine log file analysis to conventional measurement-based methods in detecting errors in IMRT delivery. Sixteen IMRT treatment plans (5 head-and-neck, 3 rectum, 3 breast, and 5 prostate plans) created with a commercial treatment planning system (TPS) were recalculated on a QA phantom. All treatment plans underwent ion chamber (IC) and 2D diode array measurements. The same set of plans was also recomputed with another commercial treatment planning system and the two sets of calculations were compared. The deviations between dosimetric measurements and independent dose calculation were evaluated. The comparisons included evaluations of DVHs and point doses calculated by the two TPS systems. Machine log files were captured during pretreatment composite point dose measurements and analyzed to verify data transfer and performance of the delivery machine. Average deviation between IC measurements and point dose calculations with the two TPSs for head-and-neck plans were 1.2 ± 1.3% and 1.4 ± 1.6%, respectively. For 2D diode array measurements, the mean gamma value with 3% dose difference and 3 mm distance-to-agreement was within 1.5% for 13 of 16 plans. The mean 3D dose differences calculated from two TPSs were within 3% for head-and-neck cases and within 2% for other plans. The machine log file analysis showed that the gantry angle, jaw position, collimator angle, and MUs were consistent as planned, and maximal MLC position error was less than 0.5 mm. The independent dose calculation followed by the machine log analysis takes an average 47 ± 6 minutes, while the experimental approach (using IC and 2D diode array measurements) takes an average about 2 hours in our clinic. Independent dose calculation followed by machine log file analysis can be a reliable tool to verify IMRT treatments. Additionally, independent dose calculations have the potential to identify several problems (heterogeneity calculations, data corruptions, system failures) with the primary TPS, which generally are not identifiable with a measurement-based approach. Additionally, machine log file analysis can identify many problems (gantry, collimator, jaw setting) which also may not be detected with a measurement-based approach. Machine log file analysis could also detect performance problems for individual MLC leaves which could be masked in the analysis of a measured fluence.

Sensitivity analysis for lexicographic ordering in radiation therapy treatment planning

PURPOSE

To introduce a method to efficiently identify and calculate meaningful tradeoffs between criteria in an interactive IMRT treatment planning procedure. The method provides a systematic approach to developing high-quality radiation therapy treatment plans.

METHODS

Treatment planners consider numerous dosimetric criteria of varying importance that, when optimized simultaneously through multicriteria optimization, yield a Pareto frontier which represents the set of Pareto-optimal treatment plans. However, generating and navigating this frontier is a time-consuming, nontrivial process. A lexicographic ordering (LO) approach to IMRT uses a physician's criteria preferences to partition the treatment planning decisions into a multistage treatment planning model. Because the relative importance of criteria optimized in the different stages may not necessarily constitute a strict prioritization, the authors introduce an interactive process, sensitivity analysis in lexicographic ordering (SALO), to allow the treatment planner control over the relative sequential-stage tradeoffs. By allowing this flexibility within a structured process, SALO implicitly restricts attention to and allows exploration of a subset of the Pareto efficient frontier that the physicians have deemed most important.

RESULTS

Improvements to treatment plans over a LO approach were found by implementing the SALO procedure on a brain case and a prostate case. In each stage, a physician assessed the tradeoff between previous stage and current stage criteria. The SALO method provided critical tradeoff information through curves approximating the relationship between criteria, which allowed the physician to determine the most desirable treatment plan.

CONCLUSIONS

The SALO procedure provides treatment planners with a directed, systematic process to treatment plan selection. By following a physician's prioritization, the treatment planner can avoid wasting effort considering clinically inferior treatment plans. The planner is guided by criteria importance, but given the information necessary to accurately adjust the relative importance at each stage. Through these attributes, the SALO procedure delivers an approach well balanced between efficiency and flexibility.

Inverse-optimized 3D conformal planning: minimizing complexity while achieving equivalence with beamlet IMRT in multiple clinical sites

PURPOSE

Inverse planned intensity modulated radiation therapy (IMRT) has helped many centers implement highly conformal treatment planning with beamlet-based techniques. The many comparisons between IMRT and 3D conformal (3DCRT) plans, however, have been limited because most 3DCRT plans are forward-planned while IMRT plans utilize inverse planning, meaning both optimization and delivery techniques are different. This work avoids that problem by comparing 3D plans generated with a unique inverse planning method for 3DCRT called inverse-optimized 3D (IO-3D) conformal planning. Since IO-3D and the beamlet IMRT to which it is compared use the same optimization techniques, cost functions, and plan evaluation tools, direct comparisons between IMRT and simple, optimized IO-3D plans are possible. Though IO-3D has some similarity to direct aperture optimization (DAO), since it directly optimizes the apertures used, IO-3D is specifically designed for 3DCRT fields (i.e., 1-2 apertures per beam) rather than starting with IMRT-like modulation and then optimizing aperture shapes. The two algorithms are very different in design, implementation, and use. The goals of this work include using IO-3D to evaluate how close simple but optimized IO-3D plans come to nonconstrained beamlet IMRT, showing that optimization, rather than modulation, may be the most important aspect of IMRT (for some sites).

METHODS

The IO-3D dose calculation and optimization functionality is integrated in the in-house 3D planning/optimization system. New features include random point dose calculation distributions, costlet and cost function capabilities, fast dose volume histogram (DVH) and plan evaluation tools, optimization search strategies designed for IO-3D, and an improved, reimplemented edge/octree calculation algorithm. The IO-3D optimization, in distinction to DAO, is designed to optimize 3D conformal plans (one to two segments per beam) and optimizes MLC segment shapes and weights with various user-controllable search strategies which optimize plans without beamlet or pencil beam approximations. IO-3D allows comparisons of beamlet, multisegment, and conformal plans optimized using the same cost functions, dose points, and plan evaluation metrics, so quantitative comparisons are straightforward. Here, comparisons of IO-3D and beamlet IMRT techniques are presented for breast, brain, liver, and lung plans.

RESULTS

IO-3D achieves high quality results comparable to beamlet IMRT, for many situations. Though the IO-3D plans have many fewer degrees of freedom for the optimization, this work finds that IO-3D plans with only one to two segments per beam are dosimetrically equivalent (or nearly so) to the beamlet IMRT plans, for several sites. IO-3D also reduces plan complexity significantly. Here, monitor units per fraction (MU/Fx) for IO-3D plans were 22%-68% less than that for the 1 cm × 1 cm beamlet IMRT plans and 72%-84% than the 0.5 cm × 0.5 cm beamlet IMRT plans.

CONCLUSIONS

The unique IO-3D algorithm illustrates that inverse planning can achieve high quality 3D conformal plans equivalent (or nearly so) to unconstrained beamlet IMRT plans, for many sites. IO-3D thus provides the potential to optimize flat or few-segment 3DCRT plans, creating less complex optimized plans which are efficient and simple to deliver. The less complex IO-3D plans have operational advantages for scenarios including adaptive replanning, cases with interfraction and intrafraction motion, and pediatric patients.

A new column-generation-based algorithm for VMAT treatment plan optimization

We study the treatment plan optimization problem for volumetric modulated arc therapy (VMAT). We propose a new column-generation-based algorithm that takes into account bounds on the gantry speed and dose rate, as well as an upper bound on the rate of change of the gantry speed, in addition to MLC constraints. The algorithm iteratively adds one aperture at each control point along the treatment arc. In each iteration, a restricted problem optimizing intensities at previously selected apertures is solved, and its solution is used to formulate a pricing problem, which selects an aperture at another control point that is compatible with previously selected apertures and leads to the largest rate of improvement in the objective function value of the restricted problem. Once a complete set of apertures is obtained, their intensities are optimized and the gantry speeds and dose rates are adjusted to minimize treatment time while satisfying all machine restrictions. Comparisons of treatment plans obtained by our algorithm to idealized IMRT plans of 177 beams on five clinical prostate cancer cases demonstrate high quality with respect to clinical dose-volume criteria. For all cases, our algorithm yields treatment plans that can be delivered in around 2 min. Implementation on a graphic processing unit enables us to finish the optimization of a VMAT plan in 25-55 s.

An Estimation of Radiobiological Parameters for Head-and-Neck Cancer Cells and the Clinical Implications

In vitro survival measurements using two human head-and-neck cancer (HNC) cell lines were performed. The specially designed split-dose surviving fraction was obtained and fitted to the linear-quadratic formalism. The repair halftime (Tr), the potential doubling time (Td), a/β and radiosensitivity a, were estimated. Other radiobiological models: EUD, BED, TCP, etc., were used to examine the potential treatment effectiveness of different IMRT techniques. Our data indicated the repair halftime of ~17 min based on two HNC cell lines. The combined a/β, a and Td are a/β = 8.1 ± 4.1 Gy, a = 0.22 ± 0.08 Gy-1, Td = 4.0 ± 1.8 day, respectively. The prolonged IMRT dose delivery for entire HNC treatment course could possibly result in the loss of biological effectiveness, i.e., the target EUDs decreased by 11% with fraction dose delivery time varying from 5 to 30 min. We determined the sublethal damage repair halftime and other radiobiological parameters for HNC cells, and to evaluate treatment effectiveness of the prolonged dose delivery times associated with different IMRT techniques. The estimated repair halftime for HNC is relatively short and may be comparable to the step-and-shoot IMRT fraction dose delivery time. The effectiveness of IMRT treatment may be improved by reducing the fraction delivery time for HNC treatment.

Generation of a novel phase-space-based cylindrical dose kernel for IMRT optimization

PURPOSE

Improving dose calculation accuracy is crucial in intensity-modulated radiation therapy (IMRT). We have developed a method for generating a phase-space-based dose kernel for IMRT planning of lung cancer patients.

METHODS

Particle transport in the linear accelerator treatment head of a 21EX, 6 MV photon beam (Varian Medical Systems, Palo Alto, CA) was simulated using the EGSnrc/BEAMnrc code system. The phase space information was recorded under the secondary jaws. Each particle in the phase space file was associated with a beamlet whose index was calculated and saved in the particle's LATCH variable. The DOSXYZnrc code was modified to accumulate the energy deposited by each particle based on its beamlet index. Furthermore, the central axis of each beamlet was calculated from the orientation of all the particles in this beamlet. A cylinder was then defined around the central axis so that only the energy deposited within the cylinder was counted. A look-up table was established for each cylinder during the tallying process. The efficiency and accuracy of the cylindrical beamlet energy deposition approach was evaluated using a treatment plan developed on a simulated lung phantom.

RESULTS

Profile and percentage depth doses computed in a water phantom for an open, square field size were within 1.5% of measurements. Dose optimized with the cylindrical dose kernel was found to be within 0.6% of that computed with the nontruncated 3D kernel. The cylindrical truncation reduced optimization time by approximately 80%.

CONCLUSIONS

A method for generating a phase-space-based dose kernel, using a truncated cylinder for scoring dose, in beamlet-based optimization of lung treatment planning was developed and found to be in good agreement with the standard, nontruncated scoring approach. Compared to previous techniques, our method significantly reduces computational time and memory requirements, which may be useful for Monte-Carlo-based 4D IMRT or IMAT treatment planning.

Radiation therapy for nasopharyngeal carcinoma using simultaneously integrated boost (SIB) protocol: a comparison planning study between intensity modulated arc radiotherapy vs. intensity modulated radiotherapy

The aim of this paper is to compare the dosimetric difference between intensity-modulated arc therapy (IMAT) and conventional intensity-modulated radiation therapy (IMRT) for radiotherapy of nasopharyngeal carcinoma (NPC) using simultaneously integrated boost (SIB) protocol. Ten patients with nasopharyngeal carcinoma underwent SIB protocol were retrospectively studied. The plan target volume (PTV) of NPC contained nasopharynx gross target volume and the positive neck lymph nodes, PTV1 contained the high-risk sites of microscopic extension and the whole nasopharynx and PTV2 contained the low-risk sites. The prescription dose of PTV was 66 Gy/30 fractions, and for PTV1 60 Gy/30 fractions and for PTV2 54 Gy/30 fractions. IMAT (two 358° arcs) and IMRT (7 fields) plans were designed for each patients using SIB strategies. The monitor unit (MU), treatment time (T) and dosimetric difference between IMRT and IMAT were compared. IMAT can achieve better conformal index (CI) than IMRT (P < 0.05) for all PTVs, while no significant difference were found in homogeneity index (HI) (P > 0.05). There's no significant difference found in radiation dose of brain stem, lenses and parotids, while the maximum dose of spinal cord of IMAT was higher than IMRT (P < 0.05). The monitor unit of IMRT (1308 ± 213) was more than IMAT (606 ± 96) (P < 0.05), while the treatment time of IMRT (540 ± 160S) was more than IMAT (160 ± 10S). This study shows that IMAT using SIB strategies for NPC radiotherapy can achieve similar target coverage with better conformity with less MU and delivery time comparing to IMRT.

Dosimetric evaluation of conventional radiotherapy, 3-D conformal radiotherapy and direct machine parameter optimisation intensity-modulated radiotherapy for breast cancer after conservative surgery

INTRODUCTION

The use of conservative surgery combined with whole-breast irradiation (WBI) has been established as a valid alternative to mastectomy for the management of early-stage breast cancer. The aim of this study was to compare dosimetric parameters of the planning target volume(PTV) and organs at risk (OARs) between conventional radiation therapy (CR), 3-D conformal radiation therapy (3DCRT), and direct machine parameter optimisation intensity-modulated radiation therapy (DMPO-IMRT) after breast-conserving surgery.

METHODS AND MATERIALS

Computed tomography (CT) scans from 20 patients (13 left-sided and 7 right-sided) previously treated with T1N0 or ductal carcinoma were selected for this dosimetric planning study. We designed CR, 3DCRT and DMPO-IMRT plans for each patient. The prescribed dose was 50 Gy/2 Gy/25 f, 95% of PTV received the prescription dose. Doses were computed with a commercially available treatment planning system using convolution/superimposition (CS) algorithm. Plans were compared according to dose-volume histogram (DVH) analysis in terms of PTV homogeneity and conformity indices (HI and CI) as well as OARs dose and volume parameters.

RESULTS

Both the HI and CI of the PTV showed statistically significant difference between CR, 3DCRT and DMPO-IMRT with those of DMPO-IMRT were best (P < 0.05). Compared with CR, 3DCRT showed smaller exposed volumes of ipsilateral lung, contralateral breast and heart while DMPO-IMRT indicated larger exposed volumes of ipsilateral lung (except for V20 and V30), contralateral breast and heart. In addition, DMPO-IMRT demonstrated an increase of exposed volume of ipsilateral lung (except for V30), contralateral breast and heart compared with 3DCRT.

CONCLUSIONS

In WBI of breast cancer after conservative surgery, 3DCRT and DMPO-IMRT improved the homogeneity and conformity of the PTV compared with CR. Meanwhile, 3DCRT reduced the irradiated volumes of OARs at all dose levels listed in our study while DMPO-IMRT reduced the irradiated volumes of OARs in high-dose areas but increased the irradiated volumes of OARs in low-dose areas.

From analytic inversion to contemporary IMRT optimization: radiation therapy planning revisited from a mathematical perspective

In this paper we look at the development of radiation therapy treatment planning from a mathematical point of view. Historically, planning for Intensity-Modulated Radiation Therapy (IMRT) has been considered as an inverse problem. We discuss first the two fundamental approaches that have been investigated to solve this inverse problem: Continuous analytic inversion techniques on one hand, and fully-discretized algebraic methods on the other hand. In the second part of the paper, we review another fundamental question which has been subject to debate from the beginning of IMRT until the present day: The rotation therapy approach versus fixed angle IMRT. This builds a bridge from historic work on IMRT planning to contemporary research in the context of Intensity-Modulated Arc Therapy (IMAT).

The impact of direct aperture optimization on plan quality and efficiency in complex head and neck IMRT

Background

Conventional step&shoot intensity modulated radio therapy (IMRT) approaches potentially lead to treatment plans with high numbers of segments and monitor units (MU) and, therefore, could be time consuming at the linear accelerator. Direct optimization methods are able to reduce the complexity without degrading the quality of the plan. The aim of this study is the evaluation of different IMRT approaches at standardized conditions for head and neck tumors.

Method

For 27 patients with carcinomas in the head and neck region a planning study with a 2-step-IMRT system (KonRad), a direct optimization system (Panther DAO) and a mixture of both approaches (MasterPlan DSS) was created. In order to avoid different prescription doses for boost volumes a simple standardization was realized. The dose was downscaled to 50 Gy to the planning target volume (PTV) which included the primary tumor as well as the bilateral lymphatic drainage (cervical and supraclavicular). Dose restrictions for the organs at risk (OAR) were downscaled to this prescription from high dose concepts up to 72 Gy. Those limits were defined as planning objectives while reaching definable PTV coverage with a standardized field setup. The parameters were evaluated from the corresponding dose volume histogram (DVH). Special attention was paid to the efficiency of the method, measured by means of calculated MU and required segments. Statistical tests of significance were applied to quantify the differences between the evaluated systems.

Results

PTV coverage for all systems in terms of V_90% and V_95% fell short of the requested 100% and 95%, respectively, but were still acceptable (range: 98.7% to 99.1% and 94.2% to 94.7%). Overall for OAR sparing and the burden of healthy tissue with low doses no technique was superior for all evaluated parameters. Differences were found for the number of segments where the direct optimization systems generated less segments. Lowest average numbers of MU were 308 by Panther DAO calculated for 2 Gy fractions. Based on these findings the treatment time at the linear accelerator is the lowest for Panther DAO.

Conclusions

All IMRT approaches implemented in the different treatment planning systems (TPS) generated clinically acceptable and comparable plans. No superior system in terms of PTV coverage and OAR sparing was found. Major differences in efficiency of the method in terms of calculated MU and treatment times were found.

3D-conformal RT, fixed-field IMRT and RapidArc, which one is better for esophageal carcinoma treated with elective nodal irradiation

The purpose of this study is to compare the characteristics of 3D-conformal radiotherapy (3D-CRT), fixed-field intensity-modulated radiotherapy (IMRT) and RapidArc for esophageal squamous cell carcinoma (ESCC) treated with elective nodal irradiation (ENI). CT datasets of 20 patients with ESCC were included and plans for single and double arcs of RapidArc (RA1 and RA2), 7-field IMRT and 3D-CRT were created and optimized for each patient. The goal was to deliver 59.6 Gy to ≥95% of the planning target volume (40 Gy to electively irradiated lymph nodal regions) while meeting the same normal-tissue dose constraints. The plans were compared based on dosimetric characteristics of target and organs at risk (OARs), monitor units (MUs), and appraised beam-on time. Both RA2 and IMRT resulted in similar target coverage (V95%, 97.84±1.50% for RA2 versus 96.96±1.15% for IMRT), homogeneity index (HI, 0.11±0.02 for RA2 versus 0.10±0.01 for IMRT) and conformity index (CI, 0.81±0.03 for RA2 versus 0.79±0.04 for IMRT), which displayed slightly better than single arc (V95%=94.55±1.50%, HI=0.12±0.02, CI=0.80±0.02) and much better than 3D-CRT (V95%=91.17±2.89%, HI=0.15±0.03, CI=0.60±0.07). The total lung V20, V30 was reduced approximately from 31%, 16% (3D-CRT) to 22%, 13% (IMRT) and 20%, 12% (RA2); the heart V30, V40 from 29%, 21% (3D-CRT) to 28%, 20% (IMRT) and 27%, 18% (RA2). The maximum dose to the spinal cord was 44.26±2.60 Gy for 3D-CRT, 42.47±2.40 Gy for IMRT, and 42.79±1.81 Gy for RA2. The number of MUs per fraction reduced from 990±165 (IMRT) to 503±70 (3D-CRT) and 502±79 (RA2). Appraised beam-on time of RapidArc was 1.2-2.4 min, which was lower than IMRT with 5.4 min by average. RapidArc, especially for double arcs plan could provide slight improvements in OARs sparing and lower MUs without compromised target qualities compared with IMRT, which was much better than 3D-CRT for ESCC treated with ENI.

Effect of photon-beam energy on VMAT and IMRT treatment plan quality and dosimetric accuracy for advanced prostate cancer

PURPOSE

The goal of the research was to evaluate treatment plan quality and dosimetric accuracy of volumetric modulated arc therapy (VMAT) and intensity-modulated radiotherapy (IMRT) plans using 6, 10, and 15 MV photon beams for prostate cancer including lymph nodes.

METHODS

In this retrospective study, VMAT and IMRT plans were generated with the Pinnacle© treatment planning system (TPS) (V9.0) for 10 prostate cancer cases. Each plan consisted of two target volumes: PTV(B) included the prostate bed, PTV(PC+LN) contained PTV(B) and lymph nodes. For plan evaluation statistics, the homogeneity index, conformity index, mean doses, and near-max doses to organs at risk (OAR) were analyzed. Treatment time and number of monitor units were assessed to compare delivery efficiency. Dosimetric plan verification was performed with a 2D ionization chamber array placed in a full scatter phantom.

RESULTS

No differences were found for target and OAR parameters in low and high energy photon beam plans for both VMAT and IMRT. A slightly higher low dose volume was detected for 6 MV VMAT plans (normal tissue: D(mean) = 16.47 Gy) compared to 10 and 15 MV VMAT plans (D(mean) = 15.90 Gy and 15.74 Gy, respectively), similar to the findings in IMRT. In VMAT, > 96% of detector points passed the 3%/ 3 mm γ criterion; marginally better accuracy was found in IMRT (> 97%).

CONCLUSION

For static and rotational IMRT, 15 MV photons did not show advantages over 6 and 10 MV high energy photon beams in large volume pelvic plans. For the investigated TPS and linac combination, 10 MV photon beams can be used as the general purpose energy for intensity modulation.

The effect of gantry spacing resolution on plan quality in a single modulated arc optimization

Volumetric‐modulated arc technique (VMAT) is an efficient form of IMRT delivery. It is advantageous over conventional IMRT in terms of treatment delivery time. This study investigates the relation between the number of segments and plan quality in VMAT optimization for a single modulated arc. Five prostate, five lung, and five head‐and‐neck (HN) patient plans were studied retrospectively. For each case, four VMAT plans were generated. The plans differed only in the number of control points used in the optimization process. The control points were spaced 2°, 3°, 4°, and 6° apart, respectively. All of the optimization parameters were the same among the four schemes. The 2° spacing plan was used as a reference to which the other three plans were compared. The plan quality was assessed by comparison of dose indices (DIs) and generalized equivalent uniform doses (gEUDs) for targets and critical structures. All optimization schemes generated clinically acceptable plans. The differences between the majority of reference and compared DIs and gEUDs were within 3%. DIs and gEUDs which differed in excess of 3% corresponded to dose levels well below the organ tolerances. The DI and the gEUD differences increased with an increase in plan complexity from prostates to HNs. Optimization with gantry spacing resolution of 4° seems to be a very balanced alternative between plan quality and plan complexity.

PACS number: 87.55.de

Bridging the gap between IMRT and VMAT: dense angularly sampled and sparse intensity modulated radiation therapy

PURPOSE

To propose an alternative radiation therapy (RT) planning and delivery scheme with optimal angular beam sampling and intrabeam modulation for improved dose distribution while maintaining high delivery efficiency.

METHODS

In the proposed approach, coined as dense angularly sampled and sparse intensity modulated RT (DASSIM-RT), a large number of beam angles are used to increase the angular sampling, leading to potentially more conformal dose distributions as compared to conventional IMRT. At the same time, intensity modulation of the incident beams is simplified to eliminate the dispensable segments, compensating the increase in delivery time caused by the increased number of beams and facilitating the plan delivery. In a sense, the proposed approach shifts and transforms, in an optimal fashion, some of the beam segments in conventional IMRT to the added beams. For newly available digital accelerators, the DASSIM-RT delivery can be made very efficient by concatenating the beams so that they can be delivered sequentially without operator's intervention. Different from VMAT, the level of intensity modulation in DASSIS-RT is field specific and optimized to meet the need of each beam direction. Three clinical cases (a head and neck (HN) case, a pancreas case, and a lung case) are used to evaluate the proposed RT scheme. DASSIM-RT, VMAT, and conventional IMRT plans are compared quantitatively in terms of the conformality index (CI) and delivery efficiency.

RESULTS

Plan quality improves generally with the number and intensity modulation of the incident beams. For a fixed number of beams or fixed level of intensity modulation, the improvement saturates after the intensity modulation or number of beams reaches to a certain level. An interplay between the two variables is observed and the saturation point depends on the values of both variables. For all the cases studied here, the CI of DASSIM-RT with 15 beams and 5 intensity levels (0.90, 0.79, and 0.84 for the HN, pancreas, and lung cases, respectively) is similar with that of conventional IMRT with seven beams and ten intensity levels (0.88, 0.79, and 0.83) and is higher than that of single-arc VMAT (0.75, 0.75, and 0.82). It is also found that the DASSIM-RT plans generally have better sparing of organs-at-risk than IMRT plans. It is estimated that the dose delivery time of DASSIM-RT with 15 beams and 5 intensity levels is about 4.5, 4.4, and 4.2 min for the HN, pancreas, and lung case, respectively, similar to that of IMRT plans with 7 beams and 10 intensity levels.

CONCLUSION

DASSIS-RT bridges the gap between IMRT and VMAT and allows optimal sampling of angular space and intrabeam modulation, thus it provides improved conformity in dose distributions while maintaining high delivery efficiency.

A planning and delivery study of a rotational IMRT technique with burst delivery

PURPOSE

A novel rotational IMRT (rIMRT) technique using burst delivery (continuous gantry rotation with beam off during MLC repositioning) is investigated. The authors evaluate the plan quality and delivery efficiency and accuracy of this dynamic technique with a conventional flat 6 MV photon beam.

METHODS

Burst-delivery rIMRT was implemented in a planning system and delivered with a 160-MLC linac. Ten rIMRT plans were generated for five anonymized patient cases encompassing head and neck, brain, prostate, and prone breast. All plans were analyzed retrospectively and not used for treatment. Among the varied plan parameters were the number of optimization points, number of arcs, gantry speed, and gantry angle range (alpha) over which the beam is turned on at each optimization point. Combined rotational/step-and-shoot rIMRT plans were also created by superimposing multiple-segment static fields at several optimization points. The rIMRT trial plans were compared with each other and with plans generated using helical tomotherapy and VMAT. Burst-mode rotational IMRT plans were delivered and verified using a diode array, ionization chambers, thermoluminescent dosimeters, and film.

RESULTS

Burst-mode rIMRT can achieve plan quality comparable to helical tomotherapy, while the former may lead to slightly better OAR sparing for certain cases and the latter generally achieves slightly lower hot spots. Few instances were found in which increasing the number of optimization points above 36, or superimposing step-and-shoot IMRT segments, led to statistically significant improvements in OAR sparing. Using an additional rIMRT partial arc yielded substantial OAR dose improvements for the brain case. Measured doses from the rIMRT plan delivery were within 4% of the plan calculation in low dose gradient regions. Delivery time range was 228-375 s for single-arc rIMRT 200-cGy prescription with a 300 MU/min dose rate, comparable to tomotherapy and VMAT.

CONCLUSIONS

Rotational IMRT with burst delivery, whether combined with static fields or not, yields clinically acceptable and deliverable treatment plans.