Smoothing intensity-modulated treatment delivery under hardware constraints

Influence of segment width on plan quality for volumetric modulated arc based stereotactic body radiotherapy

AIM

To study the influence of segment width on plan quality for volumetric modulated arc based stereotactic body radiotherapy.

BACKGROUND

The redundancy of modulation for regularly shaped small volume tumors results in creation of many small segments and an increase of monitor units, with a consequent prolongation of treatment and uncertainty in treatment delivery.

MATERIALS AND METHODS

Six cases each in lung, abdomen and liver were taken for the study. For each case, three VMAT SBRT plans were generated with different penalties on minimum segment width of 0.5, 1.0 and 1.5 cm. A comparison was made on the metrics of dose volume histogram, dosimetric indices, monitor units (MUs) and delivery accuracy.

RESULTS

The mean reduction of total MUs when compared with 0.5 cm plan was observed as 12.7 ± 6.0% and 17.5 ± 7.2% for 1.0 cm and 1.5 cm of minimum segment width, respectively. The p value showed a significant degradation in dosimetric indices for 1.5 cm plans when compared with 0.5 cm and 1.0 cm plans. The average deviation of measured dose with TPS calculated was 3.0 ± 1.1%, 2.1 ± 0.84% and 1.8 ± 0.9% for 0.5, 1.0 and 1.5 cm, respectively. The calculated gamma index with pass criteria of 2% dose difference and 2 mm distance to agreement was 95.9 ± 2.8%, 96.5 ± 2.6% and 97.8 ± 1.6% as calculated for 0.5, 1.0 and 1.5 cm of penalties, respectively. In view of the trade off between delivery efficiency and plan quality, 1 cm minimum segment width plans showed an improvement.

CONCLUSIONS

VMAT SBRT plans with increased optimal value of minimum segment width showed better plan quality and delivery efficiency for stereotactic body radiotherapy.

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.

A hybrid algorithm for instant optimization of beam weights in anatomy-based intensity modulated radiotherapy: A performance evaluation study

The study aims to introduce a hybrid optimization algorithm for anatomy-based intensity modulated radiotherapy (AB-IMRT). Our proposal is that by integrating an exact optimization algorithm with a heuristic optimization algorithm, the advantages of both the algorithms can be combined, which will lead to an efficient global optimizer solving the problem at a very fast rate. Our hybrid approach combines Gaussian elimination algorithm (exact optimizer) with fast simulated annealing algorithm (a heuristic global optimizer) for the optimization of beam weights in AB-IMRT. The algorithm has been implemented using MATLAB software. The optimization efficiency of the hybrid algorithm is clarified by (i) analysis of the numerical characteristics of the algorithm and (ii) analysis of the clinical capabilities of the algorithm. The numerical and clinical characteristics of the hybrid algorithm are compared with Gaussian elimination method (GEM) and fast simulated annealing (FSA). The numerical characteristics include convergence, consistency, number of iterations and overall optimization speed, which were analyzed for the respective cases of 8 patients. The clinical capabilities of the hybrid algorithm are demonstrated in cases of (a) prostate and (b) brain. The analyses reveal that (i) the convergence speed of the hybrid algorithm is approximately three times higher than that of FSA algorithm; (ii) the convergence (percentage reduction in the cost function) in hybrid algorithm is about 20% improved as compared to that in GEM algorithm; (iii) the hybrid algorithm is capable of producing relatively better treatment plans in terms of Conformity Index (CI) [~ 2% - 5% improvement] and Homogeneity Index (HI) [~ 4% - 10% improvement] as compared to GEM and FSA algorithms; (iv) the sparing of organs at risk in hybrid algorithm-based plans is better than that in GEM-based plans and comparable to that in FSA-based plans; and (v) the beam weights resulting from the hybrid algorithm are about 20% smoother than those obtained in GEM and FSA algorithms. In summary, the study demonstrates that hybrid algorithms can be effectively used for fast optimization of beam weights in AB-IMRT.

Fast, accurate photon beam accelerator modeling using BEAMnrc: a systematic investigation of efficiency enhancing methods and cross-section data

In this work, an investigation of efficiency enhancing methods and cross-section data in the BEAMnrc Monte Carlo (MC) code system is presented. Additionally, BEAMnrc was compared with VMC++, another special-purpose MC code system that has recently been enhanced for the simulation of the entire treatment head. BEAMnrc and VMC++ were used to simulate a 6 MV photon beam from a Siemens Primus linear accelerator (linac) and phase space (PHSP) files were generated at 100 cm source-to-surface distance for the 10 x 10 and 40 x 40 cm2 field sizes. The BEAMnrc parameters/techniques under investigation were grouped by (i) photon and bremsstrahlung cross sections, (ii) approximate efficiency improving techniques (AEITs), (iii) variance reduction techniques (VRTs), and (iv) a VRT (bremsstrahlung photon splitting) in combination with an AEIT (charged particle range rejection). The BEAMnrc PHSP file obtained without the efficiency enhancing techniques under study or, when not possible, with their default values (e.g., EXACT algorithm for the boundary crossing algorithm) and with the default cross-section data (PEGS4 and Bethe-Heitler) was used as the "base line" for accuracy verification of the PHSP files generated from the different groups described previously. Subsequently, a selection of the PHSP files was used as input for DOSXYZnrc-based water phantom dose calculations, which were verified against measurements. The performance of the different VRTs and AEITs available in BEAMnrc and of VMC++ was specified by the relative efficiency, i.e., by the efficiency of the MC simulation relative to that of the BEAMnrc base-line calculation. The highest relative efficiencies were approximately 935 (approximately 111 min on a single 2.6 GHz processor) and approximately 200 (approximately 45 min on a single processor) for the 10 x 10 field size with 50 million histories and 40 x 40 cm2 field size with 100 million histories, respectively, using the VRT directional bremsstrahlung splitting (DBS) with no electron splitting. When DBS was used with electron splitting and combined with augmented charged particle range rejection, a technique recently introduced in BEAMnrc, relative efficiencies were approximately 420 (approximately 253 min on a single processor) and approximately 175 (approximately 58 min on a single processor) for the 10 x 10 and 40 x 40 cm2 field sizes, respectively. Calculations of the Siemens Primus treatment head with VMC++ produced relative efficiencies of approximately 1400 (approximately 6 min on a single processor) and approximately 60 (approximately 4 min on a single processor) for the 10 x 10 and 40 x 40 cm2 field sizes, respectively. BEAMnrc PHSP calculations with DBS alone or DBS in combination with charged particle range rejection were more efficient than the other efficiency enhancing techniques used. Using VMC++, accurate simulations of the entire linac treatment head were performed within minutes on a single processor. Noteworthy differences (+/- 1%-3%) in the mean energy, planar fluence, and angular and spectral distributions were observed with the NIST bremsstrahlung cross sections compared with those of Bethe-Heitler (BEAMnrc default bremsstrahlung cross section). However, MC calculated dose distributions in water phantoms (using combinations of VRTs/AEITs and cross-section data) agreed within 2% of measurements. Furthermore, MC calculated dose distributions in a simulated water/air/water phantom, using NIST cross sections, were within 2% agreement with the BEAMnrc Bethe-Heitler default case.

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.

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

Currently, there are two types of treatment planning algorithms for intensity modulated radiation therapy (IMRT). The beamlet-based algorithm generates beamlet intensity maps with high complexity, resulting in large numbers of segments in the delivery after a leaf-sequencing algorithm is applied. The segment-based direct aperture optimization (DAO) algorithm includes the physical constraints of the deliverable apertures in the calculation, and achieves a conformal dose distribution using a small number of segments. However, the number of segments is pre-fixed in most of the DAO approaches, and the typical random search scheme in the optimization is computationally intensive. A regularization-based algorithm is proposed to overcome the drawbacks of the DAO method. Instead of smoothing the beamlet intensity maps as in many existing methods, we include a total-variation term in the optimization objective function to reduce the number of signal levels of the beam intensity maps. An aperture rectification algorithm is then applied to generate a significantly reduced number of deliverable apertures. As compared to the DAO algorithm, our method has an efficient form of quadratic optimization, with an additional advantage of optimizing field-specific numbers of segments based on the modulation complexity. The proposed approach is evaluated using two clinical cases. Under the condition that the clinical acceptance criteria of the treatment plan are satisfied, for the prostate patient, the total number of segments for five fields is reduced from 61 using the Eclipse planning system to 35 using the proposed algorithm; for the head and neck patient, the total number of segments for seven fields is reduced from 107 to 28. The head and neck result is also compared to that using an equal number of four segments for each field. The comparison shows that using field-specific numbers of segments achieves a much improved dose distribution.

Adaptive diffusion smoothing: a diffusion-based method to reduce IMRT field complexity

Inverse-planned intensity modulated radiation therapy (IMRT) is often able to achieve complex treatment planning goals that are unattainable with forward three-dimensional (3D) conformal planning. However, the common use of IMRT has introduced several new challenges. The potentially high degree of modulation in IMRT beams risks the loss of some advantages of 3D planning, such as excellent target coverage and high delivery efficiency. Previous attempts to reduce beam complexity by smoothing often result in plan degradation because the smoothing algorithm cannot distinguish between areas of desirable and undesirable modulation. The purpose of this work is to introduce and evaluate adaptive diffusion smoothing (ADS), a novel procedure designed to preferentially reduce IMRT beam complexity. In this method, a discrete diffusion equation is used to smooth IMRT beams using diffusion coefficients, automatically defined for each beamlet, that dictate the degree of smoothing allowed for each beamlet. This yields a method that can distinguish between areas of desirable and undesirable modulation. The ADS method has been incorporated into our optimization system as a weighted cost function penalty, with two diffusion coefficient definitions designed to promote: (1) uniform smoothing everywhere or (2) smoothing based on cost function gradients with respect to the plan beamlet intensities. The ADS method (with both coefficient types) has been tested in a phantom and in two clinical examples (prostate and head/neck). Both types of diffusion coefficients produce plans with reduced modulation and minimal dosimetric impact, but the cost function gradient-based coefficients show more potential for reducing beam modulation without affecting dosimetric plan quality. In summary, adaptive diffusion smoothing is a promising tool for ensuring that only the necessary amount of beam modulation is used, promoting more efficient and accurate IMRT planning, QA, and delivery.

What is an acceptably smoothed fluence? Dosimetric and delivery considerations for dynamic sliding window IMRT

Background

The study summarised in this report aimed to investigate the interplay between fluence complexity, dose calculation algorithms, dose calculation spatial resolution and delivery characteristics (monitor units, effective field width and dose delivery against dose prediction agreement) was investigated. A sample set of complex planning cases was selected and tested using a commercial treatment planning system capable of inverse optimisation and equipped with tools to tune fluence smoothness.

Methods

A set of increasingly smoothed fluence patterns was correlated to a generalised expression of the Modulation Index (MI) concept, in nature independent from the specific planning system used that could therefore be recommended as a predictor to score fluence "quality" at a very early stage of the IMRT QA process. Fluence complexity was also correlated to delivery accuracy and characteristics in terms of number of MU, dynamic window width and agreement between calculation and measurement (expressed as percentage of field area with a γ > 1 (%FA)) when comparing calculated vs. delivered modulated dose maps. Different resolutions of the calculation grid and different photon dose algorithms (pencil beam and anisotropic analytical algorithm) were used for the investigations.

Results and Conclusion

i) MI can be used as a reliable parameter to test different approaches/algorithms to smooth fluences implemented in a TPS, and to identify the preferable default values for the smoothing parameters if appropriate tools are implemented; ii) a MI threshold set at MI < 19 could ensure that the planned beams are safely and accurately delivered within stringent quality criteria; iii) a reduction in fluence complexity is strictly correlated to a corresponding reduction in MUs, as well as to a decrease of the average sliding window width (for dynamic IMRT delivery); iv) a smoother fluence results in a reduction of dose in the healthy tissue with a potentially relevant clinical benefit; v) increasing the smoothing parameter s, MI decreases with %FA: fluence complexity has a significant impact on the accuracy of delivery and the agreement between calculation and measurements improves with the advanced algorithms.

Reduction of IMRT beam complexity through the use of beam modulation penalties in the objective function

Inverse planned intensity modulated radiation therapy (IMRT) has become commonplace in treatment centers across the world. Due to the implications of beam complexity on treatment planning, delivery, and quality assurance, several methods have been proposed to reduce the complexity. These methods include beamlet intensity restrictions, smoothing procedures, and direct aperture optimization. Many of these methods typically sacrifice target coverage and/or normal tissue sparing in return for increased beam smoothness and delivery efficiency. In the present work, we penalize beam modulation in the inverse planning cost function to reduce beam complexity and increase delivery efficiency, while maintaining dosimetric quality. Three modulation penalties were tested: two that penalized deviation from Savitzky-Golay filtered versions of the optimized beams, and one that penalized the plan intensity map variation (a measure of overall beam modulation). The modulation penalties were applied at varying weights in a weighted sum objective (or cost) function to investigate their ability to reduce beam complexity while preserving IMRT plan quality. The behavior of the penalties was characterized on a CT phantom, and then clinical optimization comparisons were performed in the brain, prostate, and head/neck. Comparisons were made between (i) plans with a baseline cost function (ii) plans with a baseline cost function employing maximum beamlet intensity limits, and (iii) plans with each of the modulation penalties added to the baseline cost function. Plan analysis was based upon dose-volume histograms, relevant dose metrics, beam modulation, and monitor units required for step and shoot delivery. Each of the techniques yielded improvements over a baseline cost function in terms of MU reduction. In most cases, this was achieved with minimal change to the plan DVHs and metrics. In all cases, an acceptable plan was reached with each of the methods while reducing MU substantially. Each individual method has merit as a tool for reducing IMRT beam complexity and could be easily applied in the clinic to improve overall inverse plan quality. However, the penalty based upon the plan intensity map variation consistently produced the most delivery-efficient plans with the fewest computations.

Elongated beamlets: a simple technique for segment and MU reduction for sMLC IMRT delivery on accelerators utilizing 5 mm leaf widths

The focus of this work is to demonstrate the effects of using an elongated beamlet to achieve similar dose conformity as achieved with a square beamlet while reducing the number of segments and subsequent MU required. A series of 10 patients were planned for IMRT delivery to the prostate using minimum beamlet sizes of 5x5 mm2 (default scheme), 10x5 mm2 with the short axis parallel to the prostate-rectum interface (scheme 1), and 10x5 mm2 with the short axis perpendicular to the prostate-rectum interface (scheme 2). All other parameters between plans were left unchanged. Plans were appropriately normalized and evaluated for R65, R40, conformity index, total number of segments and MU. All plans were generated using the Corvus inverse planning system. The average number of segments in this study decreased by approximately 49% for both schemes 1 and 2. The subsequent number of MU required decreased by approximately 34.6%. The resultant modified modulation scaling factor (MSFmod) decreased by approximately 34.3%. Additionally, we found that each isodose distribution using scheme 2 would still meet our clinical acceptance criteria with no visible degradation in the dose distribution as compared with the default scheme. In conclusion, we have demonstrated that it is possible to achieve similar results as those obtained using a 5x5 mm2 beamlet with respect to target coverage and critical structure sparing by using strategically oriented elongated beamlets. This technique directly translates to a decreased MSF(mod) allowing for decreased leakage dose to the patient, a decreased risk of exceeding secondary shielding limits in pre-existing vaults, and shorter treatment times.

Improving IMRT delivery efficiency using intensity limits during inverse planning

Inverse planned intensity modulated radiotherapy (IMRT) fields can be highly modulated due to the large number of degrees of freedom involved in the inverse planning process. Additional modulation typically results in a more optimal plan, although the clinical rewards may be small or offset by additional delivery complexity and/or increased dose from transmission and leakage. Increasing modulation decreases delivery efficiency, and may lead to plans that are more sensitive to geometrical uncertainties. The purpose of this work is to assess the use of maximum intensity limits in inverse IMRT planning as a simple way to increase delivery efficiency without significantly affecting plan quality. Nine clinical cases (three each for brain, prostate, and head/neck) were used to evaluate advantages and disadvantages of limiting maximum intensity to increase delivery efficiency. IMRT plans were generated using in-house protocol-based constraints and objectives for the brain and head/neck, and RTOG 9406 dose volume objectives in the prostate. Each case was optimized at a series of maximum intensity ratios (the product of the maximum intensity and the number of beams divided by the prescribed dose to the target volume), and evaluated in terms of clinical metrics, dose-volume histograms, monitor units (MU) required per fraction (SMLC and DMLC delivery), and intensity map variation (a measure of the beam modulation). In each site tested, it was possible to reduce total monitor units by constraining the maximum allowed intensity without compromising the clinical acceptability of the plan. Monitor unit reductions up to 38% were observed for SMLC delivery, while reductions up to 29% were achieved for DMLC delivery. In general, complicated geometries saw a smaller reduction in monitor units for both delivery types, although DMLC delivery required significantly more monitor units in all cases. Constraining the maximum intensity in an inverse IMRT plan is a simple way to improve delivery efficiency without compromising plan objectives.

A graph-searching method for MLC leaf sequencing under constraints

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

A method to improve spatial resolution and smoothness of intensity profiles in IMRT treatment planning

In IMRT optimization, the size of beamlets used for optimizing the beam intensity distributions is a planner-selected parameter. The appropriate setting for the beamlet size is critical to the outcome of IMRT planning. With too small beamlets the dose calculation can be inaccurate, and the resulting intensity profiles can be unnecessarily complex and difficult to generate. Relatively simple intensity profiles can be obtained using large beamlets. However, this may compromise the conformity of the dose distribution. In this paper we present a method, in which multiple beamlet matrices displaced from each other by a shift in MLC leaf travel direction are used instead of the single beamlet matrix per beam in a conventional method, to achieve finer spatial resolution for the intensity distribution than the given beamlet size. Two test cases were used to assess the method by the resultant DVHs and dose distributions and characteristic indices of the intensity profiles. The results show that this method can produce optimized dose distributions that are similar to those produced by the conventional inverse planning method with the benefit of smoother intensity profiles that are easier to deliver with a computer controlled MLC.

Simultaneous minimizing monitor units and number of segments without leaf end abutment for segmental intensity modulated radiation therapy delivery

Leaf end abutment is seldom studied when delivering segmental intensity modulated radiation therapy (IMRT) fields. We developed an efficient leaf sequencing method to eliminate leaf end abutment for segmental IMRT delivery. Our method uses simple matrix and sorting operations to obtain a solution that simultaneously minimizes total monitor units and number of segments without leaf end abutment between segments. We implemented and demonstrated our method for multiple clinical cases. We compared the results of our method with the results from exhaustive search method. We found that our solution without leaf end abutment produced equivalent results to the unconstrained solutions in terms of minimum total monitor units and minimum number of leaf segments. We conclude that the leaf end abutment fields can be avoided without affecting the efficiency of segmental IMRT delivery. The major strength of our method is its simplicity and high computing speed. This potentially provides a useful means for generating segmental IMRT fields that require high spatial resolution or complex intensity distributions.

The physical basis of IMRT and inverse planning

Intensity-modulated radiation therapy (IMRT) can sculpt the high-dose volume around the site of disease with hitherto unachievable precision. Conformal avoidance of normal tissues goes hand in hand with this. Inhomogeneous dose painting is possible. The technique has become a clinical reality and is likely to be the dominant approach this decade for improving the clinical practice of photon therapy. This Series will explore all aspects of the "IMRT chain". Only 15 years ago just a handful of physicists were working on this subject. IMRT has developed so rapidly that its recent past is also its ancient history. This article will review the history of IMRT with just a glance at precursors. The physical basis of IMRT is then described including an attempt to introduce the concepts of convex and concave dose distributions, ill-conditioning, inverse-problem degeneracy, cost functions and complex solutions all with a minimum of technical jargon or mathematics. The many techniques for inverse planning are described and the review concludes with a look forward to the future of image-guided IMRT (IG-IMRT).

A quality assurance method for analyzing and verifying intensity modulated fields

A quality assurance method is developed for measuring, verifying and analyzing intensity modulated radiation fields. It is applicable for rotational and fixed-beam intensity modulated radiation therapy (IMRT) treatments. A gantry-mount device was constructed to measure the transmission dose of an IMRT field using radiographic films. A double-exposure technique with optimal kernel estimate method was developed to minimize the errors from measurements. A chi2 confidence level test method was developed to detect the discrepancies between measured and prescribed IMRT fluence distributions. Our method was tested for rotational and fixed-beam IMRT treatment verifications. The method was found insensitive to the hardware-related parameters for rotational and fixed-beam IMRT deliveries. The chi2 confidence level test was found to be more sensitive than linear correlation method in detecting relative small errors for cases with a few segments or narrow regions of interest. In conclusion, we demonstrated a quantitative method for verifying and analyzing IMRT treatment deliveries.