Adapting IMRT delivery fraction-by-fraction to cater for variable intrafraction motion

This paper presents a technique for coping with variable intrafraction organ motion when delivering intensity-modulated radiation therapy (IMRT). The strategy is an adaptive delivery in which the fluence delivered up to a particular fraction is subtracted from the required total-course planned fluence to create an adapted residual fluence for the next fraction. This requires that the fluence already delivered can be computed, knowing the intrafraction motion during each fraction. If the adaptation is unconstrained, as would be required for perfect delivery of the planned fluence, then the individual fractional fluences would become unphysical, with both negative components and spikes. Hence it is argued that constraints must be applied; first, positivity constraints and second, constraints to limit fluence spikes. Additionally, it is shown to be helpful to constrain other quantities which are explained. The power of the strategy is that it adapts to the (potentially variable) moving geometry during each fraction. It is not a perfect delivery but it is always better than making no adaptation. The fractionated nature of radiation therapy is thus exploited to advantage. The fluence adaptation method does not require re-planning at each fraction but this imposes limitations which are stated. The fuller theory of dose adaptation is also developed for intrafraction motion. The method is complementary to other adaptive strategies recently discussed with respect to interfraction motion.     

A new way of adapting IMRT delivery fraction-by-fraction to cater for variable intrafraction motion

In this paper a technique is presented for adaptive therapy to compensate for variable intrafraction tissue motion. So long as the motion can be measured or deduced for each fraction the technique modifies the fluence profile for the subsequent fractions in a repeatable cyclic way. The fluence modification is based on projecting the dose discrepancies between the cumulative delivered dose after each fraction and the expected planned dose at the same stage. It was shown that, in general, it is best to adapt the fluence profile to moving leaves that also have been modified to 'breathe' according to some regular default motion. However, it is important to point out that, if this regular default motion were to differ too much from the variable motion at each fraction, then the result can be worse than adapting to non-breathing leaves in a dynamic MLC technique. Furthermore, in general it should always be possible to improve results by starting the adaptation process with a constrained deconvolution of the regular default motion.     

Novel lung IMRT planning algorithms with nonuniform dose delivery strategy to account for respiratory motion

To effectively deliver radiation dose to lung tumors, respiratory motion has to be considered in treatment planning. In this paper we first present a new lung IMRT planning algorithm, referred as the dose shaping (DS) method, that shapes the dose distribution according to the probability distribution of the tumor over the breathing cycle to account for respiratory motion. In IMRT planning a dose-based convolution method was generally adopted to compensate for random organ motion by performing 4-D dose calculations using a tumor motion probability density function. We modified the CON-DOSE method to a dose volume histogram based convolution method (CON-DVH) that allows nonuniform dose distribution to account for respiratory motion. We implemented the two new planning algorithms on an in-house IMRT planning system that uses the Eclipse (Varian, Palo Alto, CA) planning workstation as the dose calculation engine. The new algorithms were compared with (1) the conventional margin extension approach in which margin is generated based on the extreme positions of the tumor, (2) the dose-based convolution method, and (3) gating with 3 mm residual motion. Dose volume histogram, tumor control probability, normal tissue complication probability, and mean lung dose were calculated and used to evaluate the relative performance of these approaches at the end-exhale phase of the respiratory cycle. We recruited six patients in our treatment planning study. The study demonstrated that the two new methods could significantly reduce the ipsilateral normal lung dose and outperformed the margin extension method and the dose-based convolution method. Compared with the gated approach that has the best performance in the low dose region, the two methods we proposed have similar potential to escalate tumor dose, but could be more efficient because dose is delivered continuously.     

The effects of tumor motion on planning and delivery of respiratory-gated IMRT

The purpose of this study is to investigate the effects of object motion on the planning and delivery of IMRT. Two phantoms containing objects were imaged using CT under a variety of motion conditions. The effects of object motion on axial CT acquisition with and without gating were assessed qualitatively and quantitatively. Measurements of effective slice width and position for the CT scans were made. Mutual information image fusion was adapted for use as a quantitative measure of object deformation in CT images. IMRT plans were generated on the CT scans of the moving and gated object images. These plans were delivered with motion, with and without gating, and the delivery error between the moving deliveries and a nonmoving delivery was assessed using a scalable vector-based index. Motion during CT acquisition produces motion artifact, object deformation, and object mispositioning, which can be substantially reduced with gating. Objects that vary in cross section in the direction of motion exhibit the most deformation in CT images. Mutual information provides a useful quantitative estimate of object deformation. The delivery of IMRT in the presence of target motion significantly alters the delivered dose distribution in relation to the planned distribution. The utilization of gating for IMRT treatment, including imaging, planning, and delivery, significantly reduces the errors introduced by object motion.     

Quantifying the interplay effect in prostate IMRT delivery using a convolution-based method

The authors present a segment-based convolution method to account for the interplay effect between intrafraction organ motion and the multileaf collimator position for each particular segment in intensity modulated radiation therapy (IMRT) delivered in a step-and-shoot manner. In this method, the static dose distribution attributed to each segment is convolved with the probability density function (PDF) of motion during delivery of the segment, whereas in the conventional convolution method ("average-based convolution"), the static dose distribution is convolved with the PDF averaged over an entire fraction, an entire treatment course, or even an entire patient population. In the case of IMRT delivered in a step-and-shoot manner, the average-based convolution method assumes that in each segment the target volume experiences the same motion pattern (PDF) as that of population. In the segment-based convolution method, the dose during each segment is calculated by convolving the static dose with the motion PDF specific to that segment, allowing both intrafraction motion and the interplay effect to be accounted for in the dose calculation. Intrafraction prostate motion data from a population of 35 patients tracked using the Calypso system (Calypso Medical Technologies, Inc., Seattle, WA) was used to generate motion PDFs. These were then convolved with dose distributions from clinical prostate IMRT plans. For a single segment with a small number of monitor units, the interplay effect introduced errors of up to 25.9% in the mean CTV dose compared against the planned dose evaluated by using the PDF of the entire fraction. In contrast, the interplay effect reduced the minimum CTV dose by 4.4%, and the CTV generalized equivalent uniform dose by 1.3%, in single fraction plans. For entire treatment courses delivered in either a hypofractionated (five fractions) or conventional (> 30 fractions) regimen, the discrepancy in total dose due to interplay effect was negligible.     

Effects of intra-fraction motion on IMRT dose delivery: statistical analysis and simulation

There has been some concern that organ motion, especially intra-fraction organ motion due to breathing, can negate the potential merit of intensity-modulated radiotherapy (IMRT). We wanted to find out whether this concern is justified. Specifically, we wanted to investigate whether IMRT delivery techniques with moving parts, e.g., with a multileaf collimator (MLC), are particularly sensitive to organ motion due to the interplay between organ motion and leaf motion. We also wanted to know if, and by how much, fractionation of the treatment can reduce the effects. We performed a statistical analysis and calculated the expected dose values and dose variances for volume elements of organs that move during the delivery of the IMRT. We looked at the overall influence of organ motion during the course of a fractionated treatment. A linear-quadratic model was used to consider fractionation effects. Furthermore, we developed software to simulate motion effects for IMRT delivery with an MLC, with compensators, and with a scanning beam. For the simulation we assumed a sinusoidal motion in an isocentric plane. We found that the expected dose value is independent of the treatment technique. It is just a weighted average over the path of motion of the dose distribution without motion. If the treatment is delivered in several fractions, the distribution of the dose around the expected value is close to a Gaussian. For a typical treatment with 30 fractions, the standard deviation is generally within 1% of the expected value for MLC delivery if one assumes a typical motion amplitude of 5 mm (1 cm peak to peak). The standard deviation is generally even smaller for the compensator but bigger for scanning beam delivery. For the latter it can be reduced through multiple deliveries ('paintings') of the same field. In conclusion, the main effect of organ motion in IMRT is an averaging of the dose distribution without motion over the path of the motion. This is the same as for treatments with conventional beams. Additional effects that are specific to the IMRT delivery technique appear to be relatively small, except for the scanning beam.     

IMRT planning and delivery incorporating daily dose from mega-voltage cone-beam computed tomography imaging

The technology of online mega-voltage cone-beam (CB) computed tomography (MV-CBCT) imaging is currently used in many institutions to generate a 3D anatomical dataset of a patient in treatment position. It utilizes an accelerator therapy beam, delivered with 200 degrees gantry rotation, and captured by an electronic portal imager to account for organ motion and setup variations. Although the patient dose exposure from a single volumetric MV-CBCT imaging procedure is comparable to that from standard double-exposure orthogonal portal images, daily image localization procedures can result in a significant dose increase to healthy tissue. A technique to incorporate the daily dose, from a MV-CBCT imaging procedure, in the IMRT treatment planning optimization process was developed. A composite IMRT plan incorporating the total dose from the CB was optimized with the objective of ensuring uniform target coverage while sparing the surrounding normal tissue. One head and neck cancer patient and four prostate cancer patients were planned and treated using this technique. Dosimetric results from the prostate IMRT plans optimized with or without CB showed similar target coverage and comparable sparing of bladder and rectum volumes. Average mean doses were higher by 1.6 +/- 1.0 Gy for the bladder and comparable for the rectum (-0.3 +/- 1.4 Gy). In addition, an average mean dose increase of 1.9 +/- 0.8 Gy in the femoral heads and 1.7 +/- 0.6 Gy in irradiated tissue was observed. However, the V65 and V70 values for bladder and rectum were lower by 2.3 +/- 1.5% and 2.4 +/- 2.1% indicating better volume sparing at high doses with the optimized plans incorporating CB. For the head and neck case, identical target coverage was achieved, while a comparable sparing of the brain stem, optic chiasm, and optic nerves was observed. The technique of optimized planning incorporating doses from daily online MV-CBCT procedures provides an alternative method for imaging IMRT patients. It allows for daily treatment modifications where other volumetric tomographic imaging techniques may not be feasible and/or available and where accurate patient localization with a high degree of precision is required.     

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.     

Prostate intrafraction motion evaluation using kV fluoroscopy during treatment delivery: a feasibility and accuracy study

Margin reduction for prostate radiotherapy is limited by uncertainty in prostate localization during treatment. We investigated the feasibility and accuracy of measuring prostate intrafraction motion using kV fluoroscopy performed simultaneously with radiotherapy. Three gold coils used for target localization were implanted into the patient's prostate gland before undergoing hypofractionated online image-guided step-and-shoot intensity modulated radiation therapy (IMRT) on an Elekta Synergy linear accelerator. At each fraction, the patient was aligned using a cone-beam computed tomography (CBCT), after which the IMRT treatment delivery and fluoroscopy were performed simultaneously. In addition, a post-treatment CBCT was acquired with the patient still on the table. To measure the intrafraction motion, we developed an algorithm to register the fluoroscopy images to a reference image derived from the post-treatment CBCT, and we estimated coil motion in three-dimensional (3D) space by combining information from registrations at different gantry angles. We also detected the MV beam turning on and off using MV scatter incident in the same fluoroscopy images, and used this information to synchronize our intrafraction evaluation with the treatment delivery. In addition, we assessed the following: the method to synchronize with treatment delivery, the dose from kV imaging, the accuracy of the localization, and the error propagated into the 3D localization from motion between fluoroscopy acquisitions. With 0.16 mAs/frame and a bowtie filter implemented, the coils could be localized with the gantry at both 0 degrees and 270 degrees with the MV beam off, and at 270 degrees with the MV beam on when multiple fluoroscopy frames were averaged. The localization in two-dimensions for phantom and patient measurements was performed with submillimeter accuracy. After backprojection into 3D the patient localization error was (-0.04 +/- 0.30) mm, (0.09 +/- 0.36)mm, and (0.03 +/- 0.68)mm in the right-left (RL), anterior-posterior (AP), and superior-inferior (SI) axes, respectively. Simulations showed that while oscillating (stationary) motion cannot be effectively represented in 3D, linearly drifting (nonstationary) motion is detectable with good accuracy. These results show that measuring prostate intrafraction motion using a single kV imager during radiotherapy is feasible and can be performed with acceptable accuracy.     

The use of spatial dose gradients and probability density function to evaluate the effect of internal organ motion for prostate IMRT treatment planning

The aim of this study is to investigate the effects of internal organ motion on IMRT treatment planning of prostate patients using a spatial dose gradient and probability density function. Spatial dose distributions were generated from a Pinnacle3 planning system using a co-planar, five-field intensity modulated radiation therapy (IMRT) technique. Five plans were created for each patient using equally spaced beams but shifting the angular displacement of the beam by 15 degree increments. Dose profiles taken through the isocentre in anterior-posterior (A-P), right-left (R-L) and superior-inferior (S-I) directions for IMRT plans were analysed by exporting RTOG file data from Pinnacle. The convolution of the 'static' dose distribution D0(x, y, z) and probability density function (PDF), denoted as P(x, y, z), was used to analyse the combined effect of repositioning error and internal organ motion. Organ motion leads to an enlarged beam penumbra. The amount of percentage mean dose deviation (PMDD) depends on the dose gradient and organ motion probability density function. Organ motion dose sensitivity was defined by the rate of change in PMDD with standard deviation of motion PDF and was found to increase with the maximum dose gradient in anterior, posterior, left and right directions. Due to common inferior and superior field borders of the field segments, the sharpest dose gradient will occur in the inferior or both superior and inferior penumbrae. Thus, prostate motion in the S-I direction produces the highest dose difference. The PMDD is within 2.5% when standard deviation is less than 5 mm, but the PMDD is over 2.5% in the inferior direction when standard deviation is higher than 5 mm in the inferior direction. Verification of prostate organ motion in the inferior directions is essential. The margin of the planning target volume (PTV) significantly impacts on the confidence of tumour control probability (TCP) and level of normal tissue complication probability (NTCP). Smaller margins help to reduce the dose to normal tissues, but may compromise the dose coverage of the PTV. Lower rectal NTCP can be achieved by either a smaller margin or a steeper dose gradient between PTV and rectum. With the same DVH control points, the rectum has lower complication in the seven-beam technique used in this study because of the steeper dose gradient between the target volume and rectum. The relationship between dose gradient and rectal complication can be used to evaluate IMRT treatment planning. The dose gradient analysis is a powerful tool to improve IMRT treatment plans and can be used for QA checking of treatment plans for prostate patients.     

Dosimetric and radiobiological impact of dose fractionation on respiratory motion induced IMRT delivery errors: a volumetric dose measurement study

Respiratory motion can introduce substantial dose errors during IMRT delivery. These errors are difficult to predict because of the nonsynchronous interplay between radiation beams and tissues. The present study investigates the impact of dose fractionation on respiratory motion induced dosimetric errors during IMRT delivery and their radiobiological implications by using measured 3D dose. We focused on IMRT delivery with dynamic multileaf collimation (DMLC-IMRT). IMRT plans using several beam arrangements were optimized for and delivered to a polystyrene phantom containing a simulated target and critical organs. The phantom was set in linear sinusoidal motion at a frequency of 15 cycles/min (0.25 Hz). The amplitude of the motion was +/- 0.75 cm in the longitudinal direction and +/- 0.25 cm in the lateral direction. Absolute doses were measured with a 0.125 cc ionization chamber while dose distributions were measured with transverse films spaced 6 mm apart. Measurements were performed for varying number of fractions with motion, with respiratory-gated motion, and without motion. A tumor control probability (TCP) model for an inhomogeneously irradiated tumor was used to calculate and compare TCPs for the measurements and the treatment plans. Equivalent uniform doses (EUD) were also computed. For individual fields, point measurements using an ionization chamber showed substantial dose deviations (-11.7% to 47.8%) for the moving phantom as compared to the stationary phantom. However, much smaller deviations (-1.7% to 3.5%) were observed for the composite dose of all fields. The dose distributions and DVHs of stationary and gated deliveries were in good agreement with those of treatment plans, while those of the nongated moving phantom showed substantial differences. Compared to the stationary phantom, the largest differences observed for the minimum and maximum target doses were -18.8% and +19.7%, respectively. Due to their random nature, these dose errors tended to average out over fractionated treatments. The results of five-fraction measurements showed significantly improved agreement between the moving and stationary phantom. The changes in TCP were less than 4.3% for a single fraction, and less than 2.3% for two or more fractions. Variation of average EUD per fraction was small (< 3.1 cGy for a fraction size of 200 cGy), even when the DVHs were noticeably different from that of the stationary tumor. In conclusion, IMRT treatment of sites affected by respiratory motion can introduce significant dose errors in individual field doses; however, these errors tend to cancel out between fields and average out over dose fractionation. 3D dose distributions, DVHs, TCPs, and EUDs for stationary and moving cases showed good agreement after two or more fractions, suggesting that tumors affected by respiration motion may be treated using IMRT without significant dosimetric and biological consequences.     

The susceptibility of IMRT dose distributions to intrafraction organ motion: an investigation into smoothing filters derived from four dimensional computed tomography data

This study investigated the sensitivity of static planning of intensity-modulated beams (IMBs) to intrafraction deformable organ motion and assessed whether smoothing of the IMBs at the treatment-planning stage can reduce this sensitivity. The study was performed with a 4D computed tomography (CT) data set for an IMRT treatment of a patient with liver cancer. Fluence profiles obtained from inverse-planning calculations on a standard reference CT scan were redelivered on a CT scan from the 4D data set at a different part of the breathing cycle. The use of a nonrigid registration model on the 4D data set additionally enabled detailed analysis of the overall intrafraction motion effects on the IMRT delivery during free breathing. Smoothing filters were then applied to the beam profiles within the optimization process to investigate whether this could reduce the sensitivity of IMBs to intrafraction organ motion. In addition, optimal fluence profiles from calculations on each individual phase of the breathing cycle were averaged to mimic the convolution of a static dose distribution with a motion probability kernel and assess its usefulness. Results from nonrigid registrations of the CT scan data showed a maximum liver motion of 7 mm in superior-inferior direction for this patient. Dose-volume histogram (DVH) comparison indicated a systematic shift when planning treatment on a motion-frozen, standard CT scan but delivering over a full breathing cycle. The ratio of the dose to 50% of the normal liver to 50% of the planning target volume (PTV) changed up to 28% between different phases. Smoothing beam profiles with a median-window filter did not overcome the substantial shift in dose due to a difference in breathing phase between planning and delivery of treatment. Averaging of optimal beam profiles at different phases of the breathing cycle mainly resulted in an increase in dose to the organs at risk (OAR) and did not seem beneficial to compensate for organ motion compared with using a large margin. Additionally, the results emphasized the need for 4D CT scans when aiming to reduce the internal margin (IM). Using only a single planning scan introduces a systematic shift in the dose distribution during delivery. Smoothing beam profiles either based on a single scan or over the different breathing phases was not beneficial for reducing this shift.     

A stochastic convolution/superposition method with isocenter sampling to evaluate intrafraction motion effects in IMRT

Current methods to calculate dose distributions with organ motion can be broadly classified as "dose convolution" and "fluence convolution" methods. In the former, a static dose distribution is convolved with the probability distribution function (PDF) that characterizes the motion. However, artifacts are produced near the surface and around inhomogeneities because the method assumes shift invariance. Fluence convolution avoids these artifacts by convolving the PDF with the incident fluence instead of the patient dose. In this paper we present an alternative method that improves the accuracy, generality as well as the speed of dose calculation with organ motion. The algorithm starts by sampling an isocenter point from a parametrically defined space curve corresponding to the patient-specific motion trajectory. Then a photon is sampled in the linac head and propagated through the three-dimensional (3-D) collimator structure corresponding to a particular MLC segment chosen randomly from the planned IMRT leaf sequence. The photon is then made to interact at a point in the CT-based simulation phantom. Randomly sampled monoenergetic kernel rays issued from this point are then made to deposit energy in the voxels. Our method explicitly accounts for MLC-specific effects (spectral hardening, tongue-and-groove, head scatter) as well as changes in SSD with isocentric displacement, assuming that the body moves rigidly with the isocenter. Since the positions are randomly sampled from a continuum, there is no motion discretization, and the computation takes no more time than a static calculation. To validate our method, we obtained ten separate film measurements of an IMRT plan delivered on a phantom moving sinusoidally, with each fraction starting with a random phase. For 2 cm motion amplitude, we found that a ten-fraction average of the film measurements gave an agreement with the calculated infinite fraction average to within 2 mm in the isodose curves. The results also corroborate the existing notion that the interfraction dose variability due to the interplay between the MLC motion and breathing motion averages out over typical multifraction treatments. Simulation with motion waveforms more representative of real breathing indicate that the motion can produce penumbral spreading asymmetric about the static dose distributions. Such calculations can help a clinician decide to use, for example, a larger margin in the superior direction than in the inferior direction. In the paper we demonstrate that a 15 min run on a single CPU can readily illustrate the effect of a patient-specific breathing waveform, and can guide the physician in making informed decisions about margin expansion and dose escalation.     

Reducing dose calculation time for accurate iterative IMRT planning

A time-consuming component of IMRT optimization is the dose computation required in each iteration for the evaluation of the objective function. Accurate superposition/convolution (SC) and Monte Carlo (MC) dose calculations are currently considered too time-consuming for iterative IMRT dose calculation. Thus, fast, but less accurate algorithms such as pencil beam (PB) algorithms are typically used in most current IMRT systems. This paper describes two hybrid methods that utilize the speed of fast PB algorithms yet achieve the accuracy of optimizing based upon SC algorithms via the application of dose correction matrices. In one method, the ratio method, an infrequently computed voxel-by-voxel dose ratio matrix (R = D(SC)/D(PB)) is applied for each beam to the dose distributions calculated with the PB method during the optimization. That is, D(PB) x R is used for the dose calculation during the optimization. The optimization proceeds until both the IMRT beam intensities and the dose correction ratio matrix converge. In the second method, the correction method, a periodically computed voxel-by-voxel correction matrix for each beam, defined to be the difference between the SC and PB dose computations, is used to correct PB dose distributions. To validate the methods, IMRT treatment plans developed with the hybrid methods are compared with those obtained when the SC algorithm is used for all optimization iterations and with those obtained when PB-based optimization is followed by SC-based optimization. In the 12 patient cases studied, no clinically significant differences exist in the final treatment plans developed with each of the dose computation methodologies. However, the number of time-consuming SC iterations is reduced from 6-32 for pure SC optimization to four or less for the ratio matrix method and five or less for the correction method. Because the PB algorithm is faster at computing dose, this reduces the inverse planning optimization time for our implementation by a factor of 2 to 8 compared with pure SC optimization, without compromising the quality or accuracy of the final treatment plan.     

基于Halcyon加速器的乳腺癌IMRT计划设计方法

目的:探索一种基于Halcyon平台的快速简洁的乳腺癌调强放射治疗（IMRT）计划设计方法,旨在降低乳腺癌IMRT计划设计难度。方法:回顾性选择10例包含锁骨上靶区的乳腺癌病例,在Halcyon上进行固定模式12野IMRT计划设计,勾画特定的辅助结构,使用基本固定的优化条件得到乳腺癌IMRT计划（HAL_IMRT）。并将该计划同TrueBeam临床使用的容积调强弧形治疗（TrueBeam_VMAT）计划进行比较。结果:HAL_IMRT计划患侧肺[V5]、[V10]、[V20]、[Dmean],健侧肺[V5],健侧乳腺[V5]及[Dmean]显著低于TrueBeam_VMAT计划（P<0.05）;靶区方面CI、HI、[D2]、[D98]均不存在统计学差异（P>0.05）,但HAL_IMRT计划靶区最大点剂量[Dmax]高于TrueBeam_VMAT计划（P<0.05）。结论:使用Halcyon固定模式的IMRT计划方法在满足临床剂量要求下,能够降低危及器官剂量,并且能够明显降低乳腺癌IMRT计划设计的难度。     

Dosimetric effect of respiration-gated beam on IMRT delivery

Intensity modulated radiation therapy (IMRT) with a dynamic multileaf collimator (DMLC) requires synchronization of DMLC leaf motion with dose delivery. A delay in DMLC communication is known to cause leaf lag and lead to dosimetric errors. The errors may be exacerbated by gated operation. The purpose of this study was to investigate the effect of leaf lag on the accuracy of doses delivered in gated IMRT. We first determined the effective leaf delay time by measuring the dose in a stationary phantom delivered by wedge-shaped fields. The wedge fields were generated by a DMLC at various dose rates. The so determined delay varied from 88.3 to 90.5 ms. The dosimetric effect of this delay on gated IMRT was studied by delivering wedge-shaped and clinical IMRT fields to moving and stationary phantoms at dose rates ranging from 100 to 600 MU/min, with and without gating. Respiratory motion was simulated by a linear sinusoidal motion of the phantom. An ionization chamber and films were employed for absolute dose and 2-D dose distribution measurements. Discrepancies between gated and nongated delivery to the stationary phantom were observed in both absolute dose and 2-D dose distribution measurements. These discrepancies increased monotonically with dose rate and frequency of beam interruptions, and could reach 3.7% of the total dose delivered to a 0.6 cm3 ion chamber. Isodose lines could be shifted by as much as 3 mm. The results are consistent with the explanation that beam hold-offs in gated delivery allowed the lagging leaves to catch up with the delivered monitor units each time that the beam was interrupted. Low dose rates, slow leaf speeds and low frequencies of beam interruptions reduce the effect of this delay-and-catch-up cycle. For gated IMRT it is therefore important to find a good balance between the conflicting requirements of rapid dose delivery and delivery accuracy.     

Direct aperture optimization of breast IMRT and the dosimetric impact of respiration motion

We have studied the application of direct aperture optimization (DAO) as an inverse planning tool for breast IMRT. Additionally, we have analysed the impact of respiratory motion on the quality of the delivered dose distribution. From this analysis, we have developed guidelines for balancing the desire for a high-quality optimized plan with the need to create a plan that will not degrade significantly in the presence of respiratory motion. For a DAO optimized breast IMRT plan, the tangential fields incorporate a flash field to cover the range of respiratory motion. The inverse planning algorithm then optimizes the shapes and weights of additional segments that are delivered in combination with the open fields. IMRT plans were generated using DAO with the relative weights of the open segments varied from 0% to 95%. To assess the impact of breathing motion, the dose distribution for the optimized IMRT plan was recalculated with the isocentre sampled from a predefined distribution in a Monte Carlo convolution/superposition dose engine with the breast simulated as a rigid object. The motion amplitudes applied in this study ranged from 0.5 to 2.0 cm. For a range of weighting levels assigned to the open field, comparisons were made between the static plans and the plans recalculated with motion. For the static plans, we found that uniform dose distributions could be generated with relative weights for the open segments equal to and below 80% and unacceptable levels of underdosage were observed with the weights larger than 80%. When simulated breathing motion was incorporated into the dose calculation, we observed a loss in dose uniformity as the weight of the open field was decreased to below 65%. More quantitatively, for each 1% decrease in the weight, the per cent volume of the target covered by at least 95% of the prescribed dose decreased by approximately 0.10% and 0.16% for motion amplitudes equal to 1.5 cm and 2.0 cm, respectively. When taking into account the motion effects, the most uniform and conformal dose distributions were achieved when the open segment weights were in the range of 65-80%. Within this range, high-quality IMRT plans were produced for each case. The study demonstrates that DAO with tangential fields provides a robust and efficient technique for breast IMRT planning and delivery when the open segment weight is selected between 65% and 80%.     

A generalized a priori dose uncertainty model of IMRT delivery.

Multileaf collimator-based intensity modulated radiation therapy (IMRT) is complex because each intensity modulated field consists of hundreds of subfields, each of which is associated with an intricate interplay of uncertainties. In this study, the authors have revised the previously introduced uncertainty model to provide an a priori accurate prediction of dose uncertainty during treatment planning in IMRT. In the previous model, the dose uncertainties were categorized into space-oriented dose uncertainty (SOU) and nonspace-oriented dose uncertainty (NOU). The revised model further divided the uncertainty sources into planning and delivery. SOU and NOU associated with a planning system were defined as inherent dose uncertainty. A convolution method with seven degrees of freedom was also newly applied to generalize the model for practical clinical cases. The model parameters were quantified through a set of measurements, accumulated routine quality assurance (QA) data, and peer-reviewed publications. The predicted uncertainty maps were compared with dose difference distributions between computations and 108 simple open-field measurements using a two-dimensional diode array detector to verify the validity of the model parameters and robustness of the generalized model. To examine the applicability of the model to overall dose uncertainty prediction in IMRT, a retrospective analysis of QA measurements using the diode array detector for 32 clinical IM fields was also performed. A scatter diagram and a correlation coefficient were employed to investigate a correlation of the predicted dose uncertainty distribution with the dose discrepancy distribution between calculation and delivery. In addition, a gamma test was performed to correlate failed regions in dose verification with the dose uncertainty map. The quantified model parameters well correlated the predicted dose uncertainty with the probable dose difference between calculations and measurements. It was visually validated with the scatter diagrams. The average correlation coefficient between uncertainty and dose difference of 108 verification measurements was 0.80 +/- 0.04, indicating a strong linear correlation. In the clinical IM field studies, the dose uncertainty map mimicked the probable dose difference distribution. The average correlation coefficient between the overall dose uncertainty and the dose difference of 32 QA measurements (total 13 184 comparison points) was 0.75 +/- 0.07, which also indicated a strong linear correlation between them. The failed regions of the gamma test remarkably corresponded to relatively high dose uncertainty. In conclusion, the dose uncertainty map was able to highlight high dose uncertainty regions, where more care should be taken during the treatment plan. The a priori accurate prediction of dose uncertainty in IMRT will significantly improve the treatment plan evaluation process, thus improving the quality of radiation treatments.