Using total-variation regularization for intensity modulated radiation therapy inverse planning with field-specific numbers of segments

Built-in wavelet-induced smoothness to reduce plan complexity in intensity modulated radiation therapy (IMRT)

Objective.Reducing plan complexity in intensity modulated radiation therapy (IMRT) to ensure dosimetric accuracy and delivery efficiency of the radiation treatment plans. We propose a novel approach by representing the beamlet intensities using an incomplete wavelet basis that explicitly excludes fluctuating intensity maps from the decision space (explicit hard constraint). This technique provides a built-in wavelet-induced smoothness that improves both dosimetric plan quality and delivery efficiency.Approach.The beamlet intensity maps need to be especially smooth in the leaf travel direction (referred to as theX-direction). We treat the intensity map of each beam as a 2D image and represent it using the wavelets corresponding to low-frequency changes in theX-direction (i.e. approximation and horizontal). The absence of wavelets corresponding to high-frequency changes (i.e. vertical and diagonal) induces built-in smoothness. We still utilize a regularization term in the objective function to promote smoothness in theY-direction (perpendicular to theX-direction) and further possible smoothness in theX-direction. This technique has been tested on three patient cases of different disease sites (paraspinal, lung, prostate) and all final evaluations and comparisons have been performed on an FDA-approved commercial treatment planning system (Varian Eclipse^TM).Main results.Wavelet-induced smoothness reduced monitor units by about 10%, 45%, and 14% for paraspinal, lung, and prostate cases, respectively. It also improved organ at risk sparing, especially on the complex paraspinal case where it resulted in about 7%, 13%, and 14% less mean dose to esophagus, lung, and cord, respectively. Moreover, built-in wavelet-induced smoothness desensitizes the results to changing the weight associated to the regularization term, and thereby mitigates the weight fine-tuning difficulty.Significance.Fluence smoothness is often achieved by smoothing the beamlet intensity maps using a proper regularization term in the objective function aiming at disincentivizing fluctuation in the beamlet intensities (implicit soft constraint). This work reports a novel application of wavelets in imposing an explicit smoothness hard constraint in the search space using an incomplete wavelet basis. This idea has been successfully applied to exclude complex and clinically irrelevant radiation plans from the search space. The code and pertained models along with a sample dataset are released on our LowDimRT GitHub (https://github.com/PortPy-Project/LowDimRT).

Description and evaluation of a new volumetric-modulated arc therapy plan complexity metric

This study describes a new plan complexity metric for volumetric-modulated arc therapy (VMAT) and evaluates the relationship of this metric with the VMAT dosimetric accuracy. The new modulation complexity score for VMAT (NMCSv) that is based on the aperture shape and multi-leaf collimator (MLC) leaf travel is described. Its performance is evaluated through correlation and receiver operating characteristic (ROC) analyses with patient-specific gamma passing rates using 2 3-dimensional diode arrays. For comparison, the following metrics are evaluated using the same correlation analyses: average field width, average leaf travel, modulation complexity score, and leaf travel modulation complexity score. Spearman's rank correlation analysis is performed to examine any relationships between the complexity metrics and the patient-specific gamma passing rates. ROC curves are used to assess the performance of the plan metrics using a gamma passing rate of 3%/3 mm criterion with a 95% tolerance level. In both the diode arrays, the gamma passing rates (3%/3 mm and 2%/2 mm) for patient-specific dosimetric verification of VMAT plans are moderately or weakly correlated to all the complexity metrics. NMCSv demonstrates the highest correlation with the passing rates (r = 0.652, p < 0.001 for Delta4 and r = 0.499, p < 0.001 for ArcCheck) and the highest area under the curve value (0.809, p < 0.01 for Delta4 and 0.734, p < 0.01 for ArcCheck). While using the Delta4 system, NMCSv exhibits an excellent classification performance with area under the curves of 0.926 (sensitivity: 0.913; specificity: 0.860; p < 0.01) and 0.918 (sensitivity: 0.943; specificity: 0.720; p < 0.01) for rectal and cervical cancer plans, respectively. NMCSv as a novel potential clinical plan complexity metric is moderately correlated with the gamma passing rate. It demonstrates the best performance with respect to distinguishing the dosimetric accuracy of VMAT plans among the evaluated metrics. The classification performance of complexity metrics can be affected by various dosimetry verification devices and treatment sites.

Many-isocenter optimization for robotic radiotherapy

Despite significant dosimetric gains, clinical implementation of the 4π non-coplanar radiotherapy on the widely available C-arm gantry system is hindered by limited clearance, and the need to perform complex coordinated gantry and couch motion. A robotic radiotherapy platform would be conducive to such treatment but a new conflict between field size and MLC modulation resolution needs to be managed for versatile applications. This study investigates the dosimetry and delivery efficiency of purposefully creating many isocenters to achieve simultaneously high MLC modulation resolution and large tumor coverage. An integrated optimization framework was proposed for simultaneous beam orientation optimization (BOO), isocenter selection, and fluence map optimization (FMO). The framework includes a least-square dose fidelity objective, a total variation term for regularizing the fluence smoothness, and a group sparsity term for beam selection. A minimal number of isocenters were identified for efficient target coverage. Colliding beams excluded, high-resolution small-field 4π intensity-modulated radiotherapy (IMRT) treatment plans with 50 cm source-to-isocenter distance (SID-50) on 10 Head and Neck (H&N) cancer patients were compared with low-resolution large-field plans with 100 cm SID (SID-100). With the same or better target coverage, the average reduction of [Dmean, Dmax] of 20-beam SID-50 plans from 20-beam SID-100 plans were [2.09 Gy, 1.19 Gy] for organs at risk (OARs) overall, [3.05 Gy, 0.04 Gy] for parotid gland, [3.62 Gy, 5.19 Gy] for larynx, and [3.27 Gy, 1.10 Gy] for mandible. R50 and integral dose were reduced by 5.3% and 9.6%, respectively. Wilcoxon signed-rank test showed significant difference (p   <  0.05) in planning target volume (PTV) homogeneity, PTV Dmax, R50, Integral dose, and OAR Dmean and Dmax. The estimated delivery time of 20-beam [SID-50, SID-100] plans were [19, 18] min and [14, 9] min, assuming 5 fractions and 30 fractions, respectively. With clinically acceptable delivery efficiency, many-isocenter optimization is dosimetrically desirable for treating large targets with high modulation resolution on the robotic platform.

Monte Carlo tree search -based non-coplanar trajectory design for station parameter optimized radiation therapy (SPORT)

An important yet challenging problem in LINAC-based rotational arc radiation therapy is the design of beam trajectory, which requires simultaneous consideration of delivery efficiency and final dose distribution. In this work, we propose a novel trajectory selection strategy by developing a Monte Carlo tree search (MCTS) algorithm during the beam trajectory selection process. To search through the vast number of possible trajectories, the MCTS algorithm was implemented. In this approach, a candidate trajectory is explored by starting from a leaf node and sequentially examining the next level of linked nodes with consideration of geometric and physical constraints. The maximum Upper Confidence Bounds for Trees, which is a function of average objective function value and the number of times the node under testing has been visited, was employed to intelligently select the trajectory. For each candidate trajectory, we run an inverse fluence map optimization with an infinity norm regularization. The ranking of the plan as measured by the corresponding objective function value was then fed back to update the statistics of the nodes on the trajectory. The method was evaluated with a chest wall and a brain case, and the results were compared with the coplanar and noncoplanar 4pi beam configurations. For both clinical cases, the MCTS method found effective and easy-to-deliver trajectories within an hour. As compared with the coplanar plans, it offers much better sparing of the OARs while maintaining the PTV coverage. The quality of the MCTS-generated plan is found to be comparable to the 4pi plans. Artificial intelligence based on MCTS is valuable to facilitate the design of beam trajectory and paves the way for future clinical use of non-coplanar treatment delivery.

A novel optimization framework for VMAT with dynamic gantry couch rotation

Existing volumetric modulated arc therapy (VMAT) optimization using coplanar arcs is highly efficient but usually dosimetrically inferior to intensity modulated radiation therapy (IMRT) with optimized non-coplanar beams. To achieve both dosimetric quality and delivery efficiency, we proposed in this study, a novel integrated optimization method for non-coplanar VMAT (4πVMAT). 4πVMAT with direct aperture optimization (DAO) was achieved by utilizing a least square dose fidelity objective, along with an anisotropic total variation term for regularizing the fluence smoothness, a single segment term for imposing simple apertures, and a group sparsity term for selecting beam angles. Continuous gantry/couch angle trajectories were selected using the Dijkstra's algorithm, where the edge and node costs were determined based on the maximal gantry rotation speed and the estimated fluence map at the current iteration, respectively. The couch-gantry-patient collision space was calculated based on actual machine geometry and a human subject 3D surface. Beams leading to collision are excluded from the DAO and beam trajectory selection (BTS). An alternating optimization strategy was implemented to solve the integrated DAO and BTS problem. The feasibility of 4πVMAT using one full-arc or two full-arcs was tested on nine patients with brain, lung, or prostate cancer. The plan was compared against a coplanar VMAT (2πVMAT) plan using one additional arc and collimator rotation. Compared to 2πVMAT, 4πVMAT reduced the average maximum and mean organs-at-risk dose by 9.63% and 3.08% of the prescription dose with the same target coverage. R50 was reduced by 23.0%. Maximum doses to the dose limiting organs, such as the brainstem, the major vessels, and the proximal bronchus, were reduced by 8.1 Gy (64.8%), 16.3 Gy (41.5%), and 19.83 Gy (55.5%), respectively. The novel 4πVMAT approach affords efficient delivery of non-coplanar arc trajectories that lead to dosimetric improvements compared with coplanar VMAT using more arcs.

Isodose feature-preserving voxelization (IFPV) for radiation therapy treatment planning

PURPOSE

Inverse planning involves iterative optimization of a large number of parameters and is known to be a labor-intensive procedure. To reduce the scale of computation and improve characterization of isodose plan, this paper presents an isodose feature-preserving voxelization (IFPV) framework for radiation therapy applications and demonstrates an implementation of inverse planning in the IFPV domain.

METHODS

A dose distribution in IFPV scheme is characterized by partitioning the voxels into subgroups according to their geometric and dosimetric values. Computationally, the isodose feature-preserving (IFP) clustering combines the conventional voxels that are spatially and dosimetrically close into physically meaningful clusters. A K-means algorithm and support vector machine (SVM) runs sequentially to group the voxels into IFP clusters. The former generates initial clusters according to the geometric and dosimetric information of the voxels and SVM is invoked to improve the connectivity of the IFP clusters. To illustrate the utility of the formalism, an inverse planning framework in the IFPV domain is implemented, and the resultant plans of three prostate IMRT and one head-and-neck cases are compared quantitatively with that obtained using conventional inverse planning technique.

RESULTS

The IFPV generates models with significant dimensionality reduction without compromising the spatial resolution seen in traditional downsampling schemes. The implementation of inverse planning in IFPV domain is demonstrated. In addition to the improved computational efficiency, it is found that, for the cases studied here, the IFPV-domain inverse planning yields better treatment plans than that of DVH-based planning, primarily because of more effective use of both geometric and dose information of the system during plan optimization.

CONCLUSIONS

The proposed IFPV provides a low parametric representation of isodose plan without compromising the essential characteristics of the plan, thus providing a practically valuable framework for various applications in radiation therapy.

VMAT optimization with dynamic collimator rotation

PURPOSE

Although collimator rotation is an optimization variable that can be exploited for dosimetric advantages, existing Volumetric Modulated Arc Therapy (VMAT) optimization uses a fixed collimator angle in each arc and only rotates the collimator between arcs. In this study, we develop a novel integrated optimization method for VMAT, accounting for dynamic collimator angles during the arc motion.

METHODS

Direct Aperture Optimization (DAO) for Dynamic Collimator in VMAT (DC-VMAT) was achieved by adding to the existing dose fidelity objective an anisotropic total variation term for regulating the fluence smoothness, a binary variable for forming simple apertures, and a group sparsity term for controlling collimator rotation. The optimal collimator angle for each beam angle was selected using the Dijkstra's algorithm, where the node costs depend on the estimated fluence map at the current iteration and the edge costs account for the mechanical constraints of multi-leaf collimator (MLC). An alternating optimization strategy was implemented to solve the DAO and collimator angle selection (CAS). Feasibility of DC-VMAT using one full-arc with dynamic collimator rotation was tested on a phantom with two small spherical targets, a brain, a lung and a prostate cancer patient. The plan was compared against a static collimator VMAT (SC-VMAT) plan using three full arcs with 60 degrees of collimator angle separation in patient studies.

RESULTS

With the same target coverage, DC-VMAT achieved 20.3% reduction of R50 in the phantom study, and reduced the average max and mean OAR dose by 4.49% and 2.53% of the prescription dose in patient studies, as compared with SC-VMAT. The collimator rotation co-ordinated with the gantry rotation in DC-VMAT plans for deliverability. There were 13 beam angles in the single-arc DC-VMAT plan in patient studies that requires slower gantry rotation to accommodate multiple collimator angles.

CONCLUSIONS

The novel DC-VMAT approach utilizes the dynamic collimator rotation during arc delivery. In doing so, DC-VMAT affords more sophisticated intensity modulation, alleviating the limitation previously imposed by the square beamlet from the MLC leaf thickness and achieves higher effective modulation resolution. Consequently, DC-VMAT with a single arc manages to achieve superior dosimetry than SC-VMAT with three full arcs.

Fraction-variant beam orientation optimization for non-coplanar IMRT

Conventional beam orientation optimization (BOO) algorithms for IMRT assume that the same set of beam angles is used for all treatment fractions. In this paper we present a BOO formulation based on group sparsity that simultaneously optimizes non-coplanar beam angles for all fractions, yielding a fraction-variant (FV) treatment plan. Beam angles are selected by solving a multi-fraction fluence map optimization problem involving 500-700 candidate beams per fraction, with an additional group sparsity term that encourages most candidate beams to be inactive. The optimization problem is solved using the fast iterative shrinkage-thresholding algorithm. Our FV BOO algorithm is used to create five-fraction treatment plans for digital phantom, prostate, and lung cases as well as a 30-fraction plan for a head and neck case. A homogeneous PTV dose coverage is maintained in all fractions. The treatment plans are compared with fraction-invariant plans that use a fixed set of beam angles for all fractions. The FV plans reduced OAR mean dose and D 2 values on average by 3.3% and 3.8% of the prescription dose, respectively. Notably, mean OAR dose was reduced by 14.3% of prescription dose (rectum), 11.6% (penile bulb), 10.7% (seminal vesicle), 5.5% (right femur), 3.5% (bladder), 4.0% (normal left lung), 15.5% (cochleas), and 5.2% (chiasm). D 2 was reduced by 14.9% of prescription dose (right femur), 8.2% (penile bulb), 12.7% (proximal bronchus), 4.1% (normal left lung), 15.2% (cochleas), 10.1% (orbits), 9.1% (chiasm), 8.7% (brainstem), and 7.1% (parotids). Meanwhile, PTV homogeneity defined as D 95/D 5 improved from .92 to .95 (digital phantom), from .95 to .98 (prostate case), and from .94 to .97 (lung case), and remained constant for the head and neck case. Moreover, the FV plans are dosimetrically similar to conventional plans that use twice as many beams per fraction. Thus, FV BOO offers the potential to reduce delivery time for non-coplanar IMRT.

Deterministic direct aperture optimization using multiphase piecewise constant segmentation

PURPOSE

Direct aperture optimization (DAO) attempts to incorporate machine constraints in the inverse optimization to eliminate the post-processing steps in fluence map optimization (FMO) that degrade plan quality. Current commercial DAO methods utilize a stochastic or greedy approach to search a small aperture solution space. In this study, we propose a novel deterministic direct aperture optimization that integrates the segmentation of fluence map in the optimization problem using the multiphase piecewise constant Mumford-Shah formulation.

METHODS

The Mumford-Shah based direct aperture optimization problem was formulated to include an L2-norm dose fidelity term to penalize differences between the projected dose and the prescribed dose, an anisotropic total variation term to promote piecewise continuity in the fluence maps, and the multiphase piecewise constant Mumford-Shah function to partition the fluence into pairwise discrete segments. A proximal-class, first-order primal-dual solver was implemented to solve the large scale optimization problem, and an alternating module strategy was implemented to update fluence and delivery segments. Three patients of varying complexity-one glioblastoma multiforme (GBM) patient, one lung (LNG) patient, and one bilateral head and neck (H&N) patient with 3 PTVs-were selected to test the new DAO method. For each patient, 20 non-coplanar beams were first selected using column generation, followed by the Mumford-Shah based DAO (DAO_MS ). For comparison, a popular and successful approach to DAO known as simulated annealing-a stochastic approach-was replicated. The simulated annealing DAO (DAO_SA ) plans were then created using the same beam angles and maximum number of segments per beam. PTV coverage, PTV homogeneity D95D5, and OAR sparing were assessed for each plan. In addition, high dose spillage, defined as the 50% isodose volume divided by the tumor volume, as well as conformity, defined as the van't Riet conformation number, were evaluated.

RESULTS

DAO_MS achieved essentially the same OAR doses compared with the DAO_SA plans for the GBM case. The average difference of OAR D_max and D_mean between the two plans were within 0.05% of the plan prescription dose. The lung case showed slightly improved critical structure sparing using the DAO_MS approach, where the average OAR D_max and D_mean were reduced by 3.67% and 1.08%, respectively, of the prescription dose. The DAO_MS plan substantially improved OAR dose sparing for the H&N patient, where the average OAR D_max and D_mean were reduced by over 10% of the prescription dose. The DAO_MS and DAO_SA plans were comparable for the GBM and LNG PTV coverage, while the DAO_MS plan substantially improved the H&N PTV coverage, increasing D99 by 6.98% of the prescription dose. For the GBM and LNG patients, the DAO_MS and DAO_SA plans had comparable high dose spillage but slightly worse conformity with the DAO_MS approach. For the H&N plan, DAO_MS was considerably superior in high dose spillage and conformity to the DAO_SA . The deterministic approach is able to solve the DAO problem substantially faster than the simulated annealing approach, with a 9.5- to 40-fold decrease in total solve time, depending on the patient case.

CONCLUSIONS

A novel deterministic direct aperture optimization formulation was developed and evaluated. It combines fluence map optimization and the multiphase piecewise constant Mumford-Shah segmentation into a unified framework, and the resulting optimization problem can be solved efficiently. Compared to the widely and commercially used simulated annealing DAO approach, it showed comparable dosimetry behavior for simple plans, and substantially improved OAR sparing, PTV coverage, PTV homogeneity, high dose spillage, and conformity for the more complex head and neck plan.

A new sparse optimization scheme for simultaneous beam angle and fluence map optimization in radiotherapy planning

[Formula: see text]-minimization-based sparse optimization was employed to solve the beam angle optimization (BAO) in intensity-modulated radiation therapy (IMRT) planning. The technique approximates the exact BAO formulation with efficiently computable convex surrogates, leading to plans that are inferior to those attainable with recently proposed gradient-based greedy schemes. In this paper, we alleviate/reduce the nontrivial inconsistencies between the [Formula: see text]-based formulations and the exact BAO model by proposing a new sparse optimization framework based on the most recent developments in group variable selection. We propose the incorporation of the group-folded concave penalty (gFCP) as a substitution to the [Formula: see text]-minimization framework. The new formulation is then solved by a variation of an existing gradient method. The performance of the proposed scheme is evaluated by both plan quality and the computational efficiency using three IMRT cases: a coplanar prostate case, a coplanar head-and-neck case, and a noncoplanar liver case. Involved in the evaluation are two alternative schemes: the [Formula: see text]-minimization approach and the gradient norm method (GNM). The gFCP-based scheme outperforms both counterpart approaches. In particular, gFCP generates better plans than those obtained using the [Formula: see text]-minimization for all three cases with a comparable computation time. As compared to the GNM, the gFCP improves both the plan quality and computational efficiency. The proposed gFCP-based scheme provides a promising framework for BAO and promises to improve both planning time and plan quality.

A comprehensive formulation for volumetric modulated arc therapy planning

PURPOSE

Volumetric modulated arc therapy (VMAT) is a widely employed radiation therapy technique, showing comparable dosimetry to static beam intensity modulated radiation therapy (IMRT) with reduced monitor units and treatment time. However, the current VMAT optimization has various greedy heuristics employed for an empirical solution, which jeopardizes plan consistency and quality. The authors introduce a novel direct aperture optimization method for VMAT to overcome these limitations.

METHODS

The comprehensive VMAT (comVMAT) planning was formulated as an optimization problem with an L2-norm fidelity term to penalize the difference between the optimized dose and the prescribed dose, as well as an anisotropic total variation term to promote piecewise continuity in the fluence maps, preparing it for direct aperture optimization. A level set function was used to describe the aperture shapes and the difference between aperture shapes at adjacent angles was penalized to control MLC motion range. A proximal-class optimization solver was adopted to solve the large scale optimization problem, and an alternating optimization strategy was implemented to solve the fluence intensity and aperture shapes simultaneously. Single arc comVMAT plans, utilizing 180 beams with 2° angular resolution, were generated for a glioblastoma multiforme case, a lung (LNG) case, and two head and neck cases-one with three PTVs (H&N3PTV) and one with foue PTVs (H&N4PTV)-to test the efficacy. The plans were optimized using an alternating optimization strategy. The plans were compared against the clinical VMAT (clnVMAT) plans utilizing two overlapping coplanar arcs for treatment.

RESULTS

The optimization of the comVMAT plans had converged within 600 iterations of the block minimization algorithm. comVMAT plans were able to consistently reduce the dose to all organs-at-risk (OARs) as compared to the clnVMAT plans. On average, comVMAT plans reduced the max and mean OAR dose by 6.59% and 7.45%, respectively, of the prescription dose. Reductions in max dose and mean dose were as high as 14.5 Gy in the LNG case and 15.3 Gy in the H&N3PTV case. PTV coverages measured by D95, D98, and D99 were within 0.25% of the prescription dose. By comprehensively optimizing all beams, the comVMAT optimizer gained the freedom to allow some selected beams to deliver higher intensities, yielding a dose distribution that resembles a static beam IMRT plan with beam orientation optimization.

CONCLUSIONS

The novel nongreedy VMAT approach simultaneously optimizes all beams in an arc and then directly generates deliverable apertures. The single arc VMAT approach thus fully utilizes the digital Linac's capability in dose rate and gantry rotation speed modulation. In practice, the new single VMAT algorithm generates plans superior to existing VMAT algorithms utilizing two arcs.

Beam’s-eye-view dosimetrics (BEVD) guided rotational station parameter optimized radiation therapy (SPORT) planning based on reweighted total-variation minimization

Conventional VMAT optimizes aperture shapes and weights at uniformly sampled stations, which is a generalization of the concept of a control point. Recently, rotational station parameter optimized radiation therapy (SPORT) has been proposed to improve the plan quality by inserting beams to the regions that demand additional intensity modulations, thus formulating nonuniform beam sampling. This work presents a new rotational SPORT planning strategy based on reweighted total-variation (TV) minimization (min.), using beam’s-eye-view dosimetrics (BEVD) guided beam selection. The convex programming based reweighted TV min. assures the simplified fluence-map, which facilitates single-aperture selection at each station for single-arc delivery. For the rotational arc treatment planning and non-uniform beam angle setting, the mathematical model needs to be modified by additional penalty term describing the fluence-map similarity and by determination of appropriate angular weighting factors. The proposed algorithm with additional penalty term is capable of achieving more efficient and deliverable plans adaptive to the conventional VMAT and SPORT planning schemes by reducing the dose delivery time about 5 to 10 s in three clinical cases (one prostate and two head-and-neck (HN) cases with a single and multiple targets). The BEVD guided beam selection provides effective and yet easy calculating methodology to select angles for denser, non-uniform angular sampling in SPORT planning. Our BEVD guided SPORT treatment schemes improve the dose sparing to femoral heads in the prostate and brainstem, parotid glands and oral cavity in the two HN cases, where the mean dose reduction of those organs ranges from 0.5 to 2.5 Gy. Also, it increases the conformation number assessing the dose conformity to the target from 0.84, 0.75 and 0.74 to 0.86, 0.79 and 0.80 in the prostate and two HN cases, while preserving the delivery efficiency, relative to conventional single-arc VMAT 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.

High-density 3D single molecular analysis based on compressed sensing

Single molecule fitting-based superresolution microscopy achieves sub-diffraction-limit image resolution but suffers from a need for long acquisition times to gather enough molecules. Several methods have recently been developed that analyze high molecule density images but most are only applicable to two dimensions. In this study, we implemented a high-density superresolution localization algorithm based on compressed sensing and a biplane approach that provides three-dimensional information about molecules, achieving super-resolution imaging at higher molecule densities than those achieved using the conventional single molecule fitting method.

Nonisocentric treatment strategy for breast radiation therapy: a proof of concept study

PURPOSE

To propose a nonisocentric treatment strategy as a special form of station parameter optimized radiation therapy, to improve sparing of critical structures while preserving target coverage in breast radiation therapy.

METHODS AND MATERIALS

To minimize the volume of exposed lung and heart in breast irradiation, we propose a novel nonisocentric treatment scheme by strategically placing nonconverging beams with multiple isocenters. As its name suggests, the central axes of these beams do not intersect at a single isocenter as in conventional breast treatment planning. Rather, the isocenter locations and beam directions are carefully selected, in that each beam is only responsible for a certain subvolume of the target, so as to minimize the volume of irradiated normal tissue. When put together, the beams will provide an adequate coverage of the target and expose only a minimal amount of normal tissue to radiation. We apply the nonisocentric planning technique to 2 previously treated clinical cases (breast and chest wall).

RESULTS

The proposed nonisocentric technique substantially improved sparing of the ipsilateral lung. Compared with conventional isocentric plans using 2 tangential beams, the mean lung dose was reduced by 38% and 50% using the proposed technique, and the volume of the ipsilateral lung receiving ≥ 20 Gy was reduced by a factor of approximately 2 and 3 for the breast and chest wall cases, respectively. The improvement in lung sparing is even greater compared with volumetric modulated arc therapy.

CONCLUSIONS

A nonisocentric implementation of station parameter optimized radiation therapy has been proposed for breast radiation therapy. The new treatment scheme overcomes the limitations of existing approaches and affords a useful tool for conformal breast radiation therapy, especially in cases with extreme chest wall curvature.

Toward a planning scheme for emission guided radiation therapy (EGRT): FDG based tumor tracking in a metastatic breast cancer patient

PURPOSE

Emission guided radiation therapy (EGRT) is a new modality that uses PET emissions in real-time for direct tumor tracking during radiation delivery. Radiation beamlets are delivered along positron emission tomography (PET) lines of response (LORs) by a fast rotating ring therapy unit consisting of a linear accelerator (Linac) and PET detectors. The feasibility of tumor tracking and a primitive modulation method to compensate for attenuation have been demonstrated using a 4D digital phantom in our prior work. However, the essential capability of achieving dose modulation as in conventional intensity modulated radiation therapy (IMRT) treatments remains absent. In this work, the authors develop a planning scheme for EGRT to accomplish sophisticated intensity modulation based on an IMRT plan while preserving tumor tracking.

METHODS

The planning scheme utilizes a precomputed LOR response probability distribution to achieve desired IMRT planning modulation with effects of inhomogeneous attenuation and nonuniform background activity distribution accounted for. Evaluation studies are performed on a 4D digital patient with a simulated lung tumor and a clinical patient who has a moving breast cancer metastasis in the lung. The Linac dose delivery is simulated using a voxel-based Monte Carlo algorithm. The IMRT plan is optimized for a planning target volume (PTV) that encompasses the tumor motion using the MOSEK package and a Pinnacle3™ workstation (Philips Healthcare, Fitchburg, WI) for digital and clinical patients, respectively. To obtain the emission data for both patients, the Geant4 application for tomographic emission (GATE) package and a commercial PET scanner are used. As a comparison, 3D and helical IMRT treatments covering the same PTV based on the same IMRT plan are simulated.

RESULTS

3D and helical IMRT treatments show similar dose distribution. In the digital patient case, compared with the 3D IMRT treatment, EGRT achieves a 15.1% relative increase in dose to 95% of the gross tumor volume (GTV) and a 31.8% increase to 50% of the GTV. In the patient case, EGRT yields a 15.2% relative increase in dose to 95% of the GTV and a 20.7% increase to 50% of the GTV. The organs at risk (OARs) doses are kept similar or lower for EGRT in both cases. Tumor tracking is observed in the presence of planning modulation in all EGRT treatments.

CONCLUSIONS

As compared to conventional IMRT treatments, the proposed EGRT planning scheme allows an escalated target dose while keeping dose to the OARs within the same planning limits. With the capabilities of incorporating planning modulation and accurate tumor tracking, EGRT has the potential to greatly improve targeting in radiation therapy and enable a practical and effective implementation of 4D radiation therapy for planning and delivery.

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.

Inverse planning optimization method for intensity modulated radiation therapy

In order to facilitate the leaf sequencing process in intensity modulated radiation therapy (IMRT), and design of a practical leaf sequencing algorithm, it is an important issue to smooth the planned fluence maps. The objective is to achieve both high-efficiency and high-precision dose delivering by considering characteristics of leaf sequencing process. The key factor which affects total number of monitor units for the leaf sequencing optimization process is the max flow value of the digraph which formulated from the fluence maps. Therefore, we believe that one strategy for compromising dose conformity and total number of monitor units in dose delivery is to balance the dose distribution function and the max flow value mentioned above. However, there are too many paths in the digraph, and we don't know the flow value of which path is the maximum. The maximum flow value among the horizontal paths was selected and used in the objective function of the fluence map optimization to formulate the model. The model is a traditional linear constrained quadratic optimization model which can be solved by interior point method easily. We believe that the smoothed maps from this model are more suitable for leaf sequencing optimization process than other smoothing models. A clinical head-neck case and a prostate case were tested and compared using our proposed model and the smoothing model which is based on the minimization of total variance. The optimization results with the same level of total number of monitor units (TNMU) show that the fluence maps obtained from our model have much better dose performance for the target/non-target region than the maps from total variance based on the smoothing model. This indicates that our model achieves better dose distribution when the algorithm suppresses the TNMU at the same level. Although we have just used the max flow value of the horizontal paths in the diagraph in the objective function, a good balance has been achieved between the dose conformity and the total number of monitor units. This idea can be extended to other fluence map optimization model, and we believe it can also achieve good performance.

Multicriteria optimization for volumetric-modulated arc therapy by decomposition into a fluence-based relaxation and a segment weight-based restriction

PURPOSE

To develop a method for inverse volumetric-modulated arc therapy (VMAT) planning that combines multicriteria optimization (MCO) with direct machine parameter optimization. The ultimate goal is to provide an efficient and intuitive method for generating high quality VMAT plans.

METHODS

Multicriteria radiation therapy treatment planning amounts to approximating the relevant treatment options by a discrete set of plans, and selecting the combination thereof that strikes the best possible balance between conflicting objectives. This approach is applied to two decompositions of the inverse VMAT planning problem: a fluence-based relaxation considered at a coarsened gantry angle spacing and under a regularizing penalty on fluence modulation, and a segment weight-based restriction in a neighborhood of the solution to the relaxed problem. The two considered variable domains are interconnected by direct machine parameter optimization toward reproducing the dose-volume histogram of the fluence-based solution.

RESULTS

The dose distribution quality of plans generated by the proposed MCO method was assessed by direct comparison with benchmark plans generated by a conventional VMAT planning method. The results for four patient cases (prostate, pancreas, lung, and head and neck) are highly comparable between the MCO plans and the benchmark plans: Discrepancies between studied dose-volume statistics for organs at risk were-with the exception of the kidneys of the pancreas case-within 1 Gy or 1 percentage point. Target coverage of the MCO plans was comparable with that of the benchmark plans, but with a small tendency toward a shift from conformity to homogeneity.

CONCLUSIONS

MCO allows tradeoffs between conflicting objectives encountered in VMAT planning to be explored in an interactive manner through search over a continuous representation of the relevant treatment options. Treatment plans selected from such a representation are of comparable dose distribution quality to conventionally optimized VMAT plans.

Emission guided radiation therapy for lung and prostate cancers: a feasibility study on a digital patient

PURPOSE

Accurate tumor tracking remains a challenge in current radiation therapy. Many strategies including image guided radiation therapy alleviate the problem to certain extents. The authors propose a new modality called emission guided radiation therapy (EGRT) to accurately and directly track the tumor based on its biological signature. This work is to demonstrate the feasibility of EGRT under two clinical scenarios using a 4D digital patient model.

METHODS

EGRT uses lines of response (LOR's) from positron emission events to direct beamlets of therapeutic radiation through the emission sites inside a tumor. This is accomplished by a radiation delivery system consisting of a Linac and positron emission tomography (PET) detectors on a fast rotating closed-ring gantry. During the treatment of radiotracer-administrated cancer patients, PET detectors collect LOR's from tumor uptake sites and the Linac responds in nearly real-time with beamlets of radiation along the same LOR paths. Moving tumors are therefore treated with a high targeting accuracy. Based on the EGRT concept, the authors design a treatment method with additional modulation algorithms including attenuation correction and an integrated boost scheme. Performance is evaluated using simulations of a lung tumor case with 3D motion and a prostate tumor case with setup errors. The emission process is simulated by Geant4 Application for Tomographic Emission package (GATE) and Linac dose delivery is simulated using a voxel-based Monte Carlo algorithm (VMC++).

RESULTS

In the lung case with attenuation correction, compared to a conventional helical treatment, EGRT achieves a 41% relative increase in dose to 95% of the gross tumor volume (GTV) and a 55% increase to 50% of the GTV. All dose distributions are normalized for the same dose to the lung. In the prostate case with the integrated boost and no setup error, EGRT yields a 19% and 55% relative dose increase to 95% and 50% of the GTV, respectively, when all methods are normalized for the same dose to the rectum. In the prostate case with integrated boost where setup error is present, EGRT contributes a 21% and 52% relative dose increase to 95% and 50% of the GTV, respectively.

CONCLUSIONS

As a new radiation therapy modality with inherent tumor tracking, EGRT has the potential to substantially improve targeting in radiation therapy in the presence of intrafractional and interfractional motion.

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.

Accelerated barrier optimization compressed sensing (ABOCS) reconstruction for cone-beam CT: phantom studies

PURPOSE

Recent advances in compressed sensing (CS) enable accurate CT image reconstruction from highly undersampled and noisy projection measurements, due to the sparsifiable feature of most CT images using total variation (TV). These novel reconstruction methods have demonstrated advantages in clinical applications where radiation dose reduction is critical, such as onboard cone-beam CT (CBCT) imaging in radiation therapy. The image reconstruction using CS is formulated as either a constrained problem to minimize the TV objective within a small and fixed data fidelity error, or an unconstrained problem to minimize the data fidelity error with TV regularization. However, the conventional solutions to the above two formulations are either computationally inefficient or involved with inconsistent regularization parameter tuning, which significantly limit the clinical use of CS-based iterative reconstruction. In this paper, we propose an optimization algorithm for CS reconstruction which overcomes the above two drawbacks.

METHODS

The data fidelity tolerance of CS reconstruction can be well estimated based on the measured data, as most of the projection errors are from Poisson noise after effective data correction for scatter and beam-hardening effects. We therefore adopt the TV optimization framework with a data fidelity constraint. To accelerate the convergence, we first convert such a constrained optimization using a logarithmic barrier method into a form similar to that of the conventional TV regularization based reconstruction but with an automatically adjusted penalty weight. The problem is then solved efficiently by gradient projection with an adaptive Barzilai-Borwein step-size selection scheme. The proposed algorithm is referred to as accelerated barrier optimization for CS (ABOCS), and evaluated using both digital and physical phantom studies.

RESULTS

ABOCS directly estimates the data fidelity tolerance from the raw projection data. Therefore, as demonstrated in both digital Shepp-Logan and physical head phantom studies, consistent reconstruction performances are achieved using the same algorithm parameters on scans with different noise levels and∕or on different objects. On the contrary, the penalty weight in a TV regularization based method needs to be fine-tuned in a large range (up to seven times) to maintain the reconstructed image quality. The improvement of ABOCS on computational efficiency is demonstrated in the comparisons with adaptive-steepest-descent-projection-onto-convex-sets (ASD-POCS), an existing CS reconstruction algorithm also using constrained optimization. ASD-POCS alternatively minimizes the TV objective using adaptive steepest descent (ASD) and the data fidelity error using projection onto convex sets (POCS). For similar image qualities of the Shepp-Logan phantom, ABOCS requires less computation time than ASD-POCS in MATLAB by more than one order of magnitude.

CONCLUSIONS

We propose ABOCS for CBCT reconstruction. As compared to other published CS-based algorithms, our method has attractive features of fast convergence and consistent parameter settings for different datasets. These advantages have been demonstrated on phantom studies.

Faster STORM using compressed sensing

In super-resolution microscopy methods based on single-molecule switching, the rate of accumulating single-molecule activation events often limits the time resolution. Here we developed a sparse-signal recovery technique using compressed sensing to analyze images with highly overlapping fluorescent spots. This method allows an activated fluorophore density an order of magnitude higher than what conventional single-molecule fitting methods can handle. Using this method, we demonstrated imaging microtubule dynamics in living cells with a time resolution of 3 s.

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 4D IMRT planning method using deformable image registration to improve normal tissue sparing with contemporary delivery techniques

We propose a planning method to design true 4-dimensional (4D) intensity-modulated radiotherapy (IMRT) plans, called the t4Dplan method, in which the planning target volume (PTV) of the individual phases of the 4D computed tomography (CT) and the conventional PTV receive non-uniform doses but the cumulative dose to the PTV of each phase, computed using deformable image registration (DIR), are uniform. The non-uniform dose prescription for the conventional PTV was obtained by solving linear equations that required motion-convolved 4D dose to be uniform to the PTV for the end-exhalation phase (PTV50) and by constraining maximum inhomogeneity to 20%. A plug-in code to the treatment planning system was developed to perform the IMRT optimization based on this non-uniform PTV dose prescription. The 4D dose was obtained by summing the mapped doses from individual phases of the 4D CT using DIR. This 4D dose distribution was compared with that of the internal target volume (ITV) method. The robustness of the 4D plans over the course of radiotherapy was evaluated by computing the 4D dose distributions on repeat 4D CT datasets. Three patients with lung tumors were selected to demonstrate the advantages of the t4Dplan method compared with the commonly used ITV method. The 4D dose distribution using the t4Dplan method resulted in greater normal tissue sparing (such as lung, stomach, liver and heart) than did plans designed using the ITV method. The dose volume histograms of cumulative 4D doses to the PTV50, clinical target volume, lung, spinal cord, liver, and heart on the 4D repeat CTs for the two patients were similar to those for the 4D dose at the time of original planning.

Inverse planning for IMRT with nonuniform beam profiles using total-variation regularization (TVR)

PURPOSE

Radiation therapy with high dose rate and flattening filter-free (FFF) beams has the potential advantage of greatly reduced treatment time and out-of-field dose. Current inverse planning algorithms are, however, not customized for beams with nonuniform incident profiles and the resultant IMRT plans are often inefficient in delivery. The authors propose a total-variation regularization (TVR)-based formalism by taking the inherent shapes of incident beam profiles into account.

METHODS

A novel TVR-based inverse planning formalism is established for IMRT with nonuniform beam profiles. The authors introduce a TVR term into the objective function, which encourages piecewise constant fluence in the nonuniform FFF fluence domain. The proposed algorithm is applied to lung and prostate and head and neck cases and its performance is evaluated by comparing the resulting plans to those obtained using a conventional beamlet-based optimization (BBO).

RESULTS

For the prostate case, the authors' algorithm produces acceptable dose distributions with only 21 segments, while the conventional BBO requires 114 segments. For the lung case and the head and neck case, the proposed method generates similar coverage of target volume and sparing of the organs-at-risk as compared to BBO, but with a markedly reduced segment number.

CONCLUSIONS

TVR-based optimization in nonflat beam domain provides an effective way to maximally leverage the technical capacity of radiation therapy with FFF fields. The technique can generate effective IMRT plans with improved dose delivery efficiency without significant deterioration of the dose distribution.

Toward truly optimal IMRT dose distribution: inverse planning with voxel-specific penalty

PURPOSE

To establish an inverse planning framework with adjustable voxel penalty for more conformal IMRT dose distribution as well as improved interactive controllability over the regional dose distribution of the resultant plan.

MATERIALS AND METHOD

In the proposed coarse-to-fine planning scheme, a conventional inverse planning with organ specific parameters is first performed. The voxel penalty scheme is then "switched on" by allowing the prescription dose to change on an individual voxel scale according to the deviation of the actual voxel dose from the ideally desired dose. The rationale here is intuitive: when the dose at a voxel does not meet its ideal dose, it simply implies that this voxel is not competitive enough when compared with the ones that have met their planning goal. In this case, increasing the penalty of the voxel by varying the prescription can boost its competitiveness and thus improve its dose. After the prescription adjustment, the plan is re-optimized. The dose adjustment/re-optimization procedure is repeated until the resultant dose distribution cannot be improved anymore. The prescription adjustment on a finer scale can be accomplished either automatically or manually. In the latter case, the regions/voxels where a dose improvement is needed are selected visually, unlike in the automatic case where the selection is done purely based on the difference of the actual dose at a given voxel and its ideal prescription. The performance of the proposed method is evaluated using a head and neck and a prostate case.

RESULTS

An inverse planning framework with the voxel-specific penalty is established. By adjusting voxel prescriptions iteratively to boost the region where large mismatch between the actual calculated and desired doses occurs, substantial improvements can be achieved in the final dose distribution. The proposed method is applied to a head and neck case and a prostate case. For the former case, a significant reduction in the maximum dose to the brainstem is achieved while the PTV dose coverage is greatly improved. The doses to other organs at risk are also reduced, ranging from 10% to 30%. For the prostate case, the use of the voxel penalty scheme also results in vast improvements to the final dose distribution. The PTV experiences improved dose uniformity and the mean dose to the rectum and bladder is reduced by as much as 15%.

CONCLUSION

Introduction of the spatially non-uniform and adjustable prescription provides room for further improvements of currently achievable dose distributions and equips the planner with an effective tool to modify IMRT dose distributions interactively. The technique is easily implementable in any existing inverse planning platform, which should facilitate clinical IMRT planning process and, in future, off-line/on-line adaptive IMRT.

Clinical development of a failure detection-based online repositioning strategy for prostate IMRT--experiments, simulation, and dosimetry study

PURPOSE

To implement and evaluate clinic-ready adaptive imaging protocols for online patient repositioning (motion tracking) during prostate IMRT using treatment beam imaging supplemented by minimal, as-needed use of on-board kV.

METHODS

The authors examine the two-step decision-making strategy: (1) Use cine-MV imaging and online-updated characterization of prostate motion to detect target motion that is potentially beyond a predefined threshold and (2) use paired MV-kV 3D localization to determine overthreshold displacement and, if needed, reposition the patient. Two levels of clinical implementation were evaluated: (1) Field-by-field based motion correction for present-day linacs and (2) instantaneous repositioning for new-generation linacs with capabilities of simultaneous MV-kV imaging and remote automatic couch control during treatment delivery. Experiments were performed on a Varian Trilogy linac in clinical mode using a 4D motion phantom programed with prostate motion trajectories taken from patient data. Dosimetric impact was examined using a 2D ion chamber array. Simulations were done for 536 trajectories from 17 patients.

RESULTS

Despite the loss of marker detection efficiency caused by the MLC leaves sometimes obscuring the field at the marker's projected position on the MV imager, the field-by-field correction halved (from 23% to 10%) the mean percentage of time that target displacement exceeded a 3 mm threshold, as compared to no intervention. This was achieved at minimal cost in additional imaging (average of one MV-kV pair per two to three treatment fractions) and with a very small number of repositionings (once every four to five fractions). Also with low kV usage (approximation 2/fraction), the instantaneous repositioning approach reduced overthreshold time by more than 75% (23% to 5%) even with severe MLC blockage as often encountered in current IMRT and could reduce the overthreshold time tenfold (to < 2%) if the MLC blockage problem were relieved. The information acquired for repositioning using combined MV-kV images was found to have submillimeter accuracy.

CONCLUSIONS

This work demonstrated with a current clinical setup that substantial reduction of adverse targeting effects of intrafraction prostate motion can be realized. The proposed adaptive imaging strategy incurs minimal imaging dose to the patient as compared to other stereoscopic imaging techniques.

A unified framework for 3D radiation therapy and IMRT planning: plan optimization in the beamlet domain by constraining or regularizing the fluence map variations

The purpose of this work is to demonstrate that physical constraints on fluence gradients in 3D radiation therapy (RT) planning can be incorporated into beamlet optimization explicitly by direct constraint on the spatial variation of the fluence maps or implicitly by using total-variation regularization (TVR). The former method forces the fluence to vary in accordance with the known form of a wedged field and latter encourages the fluence to take the known form of the wedged field by requiring the derivatives of the fluence maps to be piece-wise constant. The performances of the proposed methods are evaluated by using a brain cancer case and a head and neck case. It is found that both approaches are capable of providing clinically sensible 3D RT solutions with monotonically varying fluence maps. For currently available 3D RT delivery schemes based on the use of customized physical or dynamic wedges, constrained optimization seems to be more useful because the optimized fields are directly deliverable. Working in the beamlet domain provides a natural way to model the spatial variation of the beam fluence. The proposed methods take advantage of the fact that 3D RT is a special form of intensity-modulated radiation therapy (IMRT) and finds the optimal plan by searching for fields with a certain type of spatial variation. The approach provides a unified framework for 3D CRT and IMRT plan optimization.

Search for IMRT inverse plans with piecewise constant fluence maps using compressed sensing techniques

An intensity-modulated radiation therapy (IMRT) field is composed of a series of segmented beams. It is practically important to reduce the number of segments while maintaining the conformality of the final dose distribution. In this article, the authors quantify the complexity of an IMRT fluence map by introducing the concept of sparsity of fluence maps and formulate the inverse planning problem into a framework of compressing sensing. In this approach, the treatment planning is modeled as a multiobjective optimization problem, with one objective on the dose performance and the other on the sparsity of the resultant fluence maps. A Pareto frontier is calculated, and the achieved dose distributions associated with the Pareto efficient points are evaluated using clinical acceptance criteria. The clinically acceptable dose distribution with the smallest number of segments is chosen as the final solution. The method is demonstrated in the application of fixed-gantry IMRT on a prostate patient. The result shows that the total number of segments is greatly reduced while a satisfactory dose distribution is still achieved. With the focus on the sparsity of the optimal solution, the proposed method is distinct from the existing beamlet- or segment-based optimization algorithms.