Lung tumor motion change during stereotactic body radiotherapy (SBRT): an evaluation using MRI.

The purpose of this study is to investigate changes in lung tumor internal target volume during stereotactic body radiotherapy treatment (SBRT) using magnetic resonance imaging (MRI). Ten lung cancer patients (13 tumors) undergoing SBRT (48 Gy over four consecutive days) were evaluated. Each patient underwent three lung MRI evaluations: before SBRT (MRI-1), after fraction 3 of SBRT (MRI-3), and three months after completion of SBRT (MRI-3m). Each MRI consisted of T1-weighted images in axial plane through the entire lung. A cone-beam CT (CBCT) was taken before each fraction. On MRI and CBCT taken before fractions 1 and 3, gross tumor volume (GTV) was contoured and differences between the two volumes were compared. Median tumor size on CBCT before fractions1 (CBCT-1) and 3 (CBCT-3) was 8.68 and 11.10 cm3, respectively. In 12 tumors, the GTV was larger on CBCT-3 compared to CBCT-1 (median enlargement, 1.56 cm3). Median tumor size on MRI-1, MRI-3, and MRI-3m was 7.91, 11.60, and 3.33 cm3, respectively. In all patients, the GTV was larger on MRI-3 compared to MRI-1 (median enlargement, 1.54 cm3). In all patients, GTV was smaller on MRI-3m compared to MRI-1 (median shrinkage, 5.44 cm3). On CBCT and MRI, all patients showed enlargement of the GTV during the treatment week of SBRT, except for one patient who showed minimal shrinkage (0.86 cm3). Changes in tumor volume are unpredictable; therefore, motion and breathing must be taken into account during treatment planning, and image-guided methods should be used, when treating with large fraction sizes.

Dosimetry tools and techniques for IMRT.

Intensity modulated radiation therapy (IMRT) poses a number of challenges for properly measuring commissioning data and quality assurance (QA) radiation dose distributions. This report provides a comprehensive overview of how dosimeters, phantoms, and dose distribution analysis techniques should be used to support the commissioning and quality assurance requirements of an IMRT program. The proper applications of each dosimeter are described along with the limitations of each system. Point detectors, arrays, film, and electronic portal imagers are discussed with respect to their proper use, along with potential applications of 3D dosimetry. Regardless of the IMRT technique utilized, some situations require the use of multiple detectors for the acquisition of accurate commissioning data. The overall goal of this task group report is to provide a document that aids the physicist in the proper selection and use of the dosimetry tools available for IMRT QA and to provide a resource for physicists that describes dosimetry measurement techniques for purposes of IMRT commissioning and measurement-based characterization or verification of IMRT treatment plans. This report is not intended to provide a comprehensive review of commissioning and QA procedures for IMRT. Instead, this report focuses on the aspects of metrology, particularly the practical aspects of measurements that are unique to IMRT. The metrology of IMRT concerns the application of measurement instruments and their suitability, calibration, and quality control of measurements. Each of the dosimetry measurement tools has limitations that need to be considered when incorporating them into a commissioning process or a comprehensive QA program. For example, routine quality assurance procedures require the use of robust field dosimetry systems. These often exhibit limitations with respect to spatial resolution or energy response and need to themselves be commissioned against more established dosimeters. A chain of dosimeters, from secondary standards to field instruments, is established to assure the quantitative nature of the tests. This report is intended to describe the characteristics of the components of these systems; dosimeters, phantoms, and dose evaluation algorithms. This work is the report of AAPM Task Group 120.

Assessment of the setup dependence of detector response functions for mega-voltage linear accelerators.

PURPOSE: Accurate modeling of beam profiles is important for precise treatment planning dosimetry. Calculated beam profiles need to precisely replicate profiles measured during machine commissioning. Finite detector size introduces perturbations into the measured profiles, which, in turn, impact the resulting modeled profiles. The authors investigate a method for extracting the unperturbed beam profiles from those measured during linear accelerator commissioning. METHODS: In-plane and cross-plane data were collected for an Elekta Synergy linac at 6 MV using ionization chambers of volume 0.01, 0.04, 0.13, and 0.65 cm3 and a diode of surface area 0.64 mm2. The detectors were orientated with the stem perpendicular to the beam and pointing away from the gantry. Profiles were measured for a 10 x 10 cm2 field at depths ranging from 0.8 to 25.0 cm and SSDs from 90 to 110 cm. Shaping parameters of a Gaussian response function were obtained relative to the Edge detector. The Gaussian function was deconvolved from the measured ionization chamber data. The Edge detector profile was taken as an approximation to the true profile, to which deconvolved data were compared. Data were also collected with CC13 and Edge detectors for additional fields and energies on an Elekta Synergy, Varian Trilogy, and Siemens Oncor linear accelerator and response functions obtained. Response functions were compared as a function of depth, SSD, and detector scan direction. Variations in the shaping parameter were introduced and the effect on the resulting deconvolution profiles assessed. RESULTS: Up to 10% setup dependence in the Gaussian shaping parameter occurred, for each detector for a particular plane. This translated to less than a +/- 0.7 mm variation in the 80%-20% penumbral width. For large volume ionization chambers such as the FC65 Farmer type, where the cavity length to diameter ratio is far from 1, the scan direction produced up to a 40% difference in the shaping parameter between in-plane and cross-plane measurements. This is primarily due to the directional difference in penumbral width measured by the FC65 chamber, which can more than double in profiles obtained with the detector stem parallel compared to perpendicular to the scan direction. For the more symmetric CC13 chamber the variation was only 3% between in-plane and cross-plane measurements. CONCLUSIONS: The authors have shown that the detector response varies with detector type, depth, SSD, and detector scan direction. In-plane vs. cross-plane scanning can require calculation of a direction dependent response function. The effect of a 10% overall variation in the response function, for an ionization chamber, translates to a small deviation in the penumbra from that of the Edge detector measured profile when deconvolved. Due to the uncertainties introduced by deconvolution the Edge detector would be preferable in obtaining an approximation of the true profile, particularly for field sizes where the energy dependence of the diode can be neglected. However, an averaged response function could be utilized to provide a good approximation of the true profile for large ionization chambers and for larger fields for which diode detectors are not recommended.

Accurate heterogeneous dose calculation for lung cancer patients without high-resolution CT densities.

The aim of this study was to investigate the relative accuracy of megavoltage photon-beam dose calculations employing either 5 bulk densities or independent voxel densities determined by calibration of the CT Houndsfield number. Full-resolution CT and bulk density treatment plans were generated for 70 lung or esophageal cancer tumors (66 cases) using a commercial treatment planning system with an adaptive convolution dose calculation algorithm (Pinnacle3, Philips Medicals Systems). Bulk densities were applied to segmented regions. Individual and population average densities were compared to the full-resolution plan for each case. Monitor units were kept constant and no normalizations were employed. Dose volume histograms (DVH) and dose difference distributions were examined for all cases. The average densities of the segmented air, lung, fat, soft tissue, and bone for the entire set were found to be 0.14, 0.26, 0.89, 1.02, and 1.12 g/cc, respectively. In all cases, the normal tissue DVH agreed to better than 2% in dose. In 62 of 70 target DVHs, agreement to better than 3% in dose was observed. Six cases demonstrated emphysema, one with bullous formations and one with a hiatus hernia having a large volume of gas. They required the additional assignment of density to the emphysemic lung and inflammatory changes to the lung, the regions of collapsed lung, the bullous formations, and the hernia gas. Bulk tissue density dose calculation provides an accurate method of heterogeneous dose calculation. However, patients with advanced emphysema may require high-resolution CT studies for accurate treatment planning.

Evaluating target cold spots by the use of tail EUDs.

PURPOSE: To propose a new measure of target underdose that can be used in the evaluation and optimization of radiotherapy dose distributions. METHODS AND MATERIALS: We compare various formulations of the equivalent uniform dose (EUD) and introduce a modification of existing EUD definitions, which we call tail EUD. Tail EUD is a measure of "cold spots" below the prescription dose in the target dose distribution, using units of gray (Gy). We investigate the mathematical properties of various target EUD concepts, including tail EUD. We apply the tail EUD measure retrospectively to intensity modulated radiation therapy (IMRT) treatment plans from our plan database. We also use tail EUD as an optimization objective in the optimization of prostate, pancreas, and head-and-neck plans. RESULTS: Tail EUD has desirable mathematical properties. In particular, it is convex and it leads to convex level sets (i.e., no local minima) if the EUD from which it is derived is concave. The tail EUD value is correlated with the subjective degree of target coverage. Constraining tail EUDs to a certain level in plan optimization leads to comparable target coverage in different plans and treatment sites. CONCLUSIONS: The newly introduced concept of tail EUD appears to be useful for both plan evaluation and optimization. In addition it can potentially be applied in the design of new clinical protocols.

A prototype quantitative film scanner for radiochromic film dosimetry.

We have developed a high resolution, quantitative, two-dimensional optical film scanner for use with a commercial high sensitivity radiochromic film (RCF) for measuring single fraction external-beam radiotherapy dose distributions. The film scanner was designed to eliminate artifacts commonly observed in RCF dosimetry. The scanner employed a stationary light source and detector with a moving antireflective glass film platen attached to a high precision computerized X-Y translation stage. An ultrabright red light emitting diode (LED) with a peak output at 633 nm and full width at half maximum (FWHM) of 16 nm was selected as the scanner light source to match the RCF absorption peak. A dual detector system was created using two silicon photodiode detectors to simultaneously measure incident and transmitted light. The LED light output was focused to a submillimeter (FWHM 0.67 mm) spot size, which was determined from a scanning knife-edge technique for measuring Gaussian optical beams. Data acquisition was performed with a 16-bit A/D card in conjunction with commercial software. The linearity of the measured densities on the scanner was tested using a calibrated neutral-density step filter. Sensitometric curves and three IMRT field scans were acquired with a spatial resolution of 1 mm for both radiographic film and RCF. The results were compared with measurements taken with a commercial diode array under identical delivery conditions. The RCF was rotated by 90 deg and rescanned to study orientation effects. Comparison between the RCF and the diode array measurements using percent dose difference and distance-to-agreement criteria produced average passing rates of 99.0% using 3%/3 mm criteria and 96.7% using 2%/2 mm criteria. The same comparison between the radiographic film and diode array measurements resulted in average passing rates 96.6% and 91.6% for the above two criteria, respectively. No measurable light-scatter or interference scanner artifacts were observed. The RCF rotated by 90 deg showed no measurable orientation effect. A scan of a 15 x 15 cm2 area with 1 mm resolution required 22 min to acquire. The LED densitometer provides an accurate film dosimetry system with 1 mm or better resolution. The scanner eliminates the orientation dependence of RCF dosimetry that was previously reported with commercial flatbed scanners.

A computational implementation and comparison of several intensity modulated proton therapy treatment planning algorithms.

The authors present a comparative study of intensity modulated proton therapy (IMPT) treatment planning employing algorithms of three-dimensional (3D) modulation, and 2.5-dimensional (2.5D) modulation, and intensity modulated distal edge tracking (DET) [A. Lomax, Phys. Med. Biol. 44, 185-205 (1999)] applied to the treatment of head-and-neck cancer radiotherapy. These three approaches were also compared with 6 MV photon intensity modulated radiation therapy (IMRT). All algorithms were implemented in the University of Florida Optimized Radiation Therapy system using a finite sized pencil beam dose model and a convex fluence map optimization model. The 3D IMPT and the DET algorithms showed considerable advantages over the photon IMRT in terms of dose conformity and sparing of organs at risk when the beam number was not constrained. The 2.5D algorithm did not show an advantage over the photon IMRT except in the dose reduction to the distant healthy tissues, which is inherent in proton beam delivery. The influences of proton beam number and pencil beam size on the IMPT plan quality were also studied. Out of 24 cases studied, three cases could be adequately planned with one beam and 12 cases could be adequately planned with two beams, but the dose uniformity was often marginally acceptable. Adding one or two more beams in each case dramatically improved the dose uniformity. The finite pencil beam size had more influence on the plan quality of the 2.5D and DET algorithms than that of the 3D IMPT. To obtain a satisfactory plan quality, a 0.5 cm pencil beam size was required for the 3D IMPT and a 0.3 cm size was required for the 2.5D and the DET algorithms. Delivery of the IMPT plans produced in this study would require a proton beam spot scanning technique that has yet to be developed clinically.

Comparative analysis of 60Co intensity-modulated radiation therapy.

In this study, we perform a scientific comparative analysis of using (60)Co beams in intensity-modulated radiation therapy (IMRT). In particular, we evaluate the treatment plan quality obtained with (i) 6 MV, 18 MV and (60)Co IMRT; (ii) different numbers of static multileaf collimator (MLC) delivered (60)Co beams and (iii) a helical tomotherapy (60)Co beam geometry. We employ a convex fluence map optimization (FMO) model, which allows for the comparison of plan quality between different beam energies and configurations for a given case. A total of 25 clinical patient cases that each contain volumetric CT studies, primary and secondary delineated targets, and contoured structures were studied: 5 head-and-neck (H&N), 5 prostate, 5 central nervous system (CNS), 5 breast and 5 lung cases. The DICOM plan data were anonymized and exported to the University of Florida optimized radiation therapy (UFORT) treatment planning system. The FMO problem was solved for each case for 5-71 equidistant beams as well as a helical geometry for H&N, prostate, CNS and lung cases, and for 3-7 equidistant beams in the upper hemisphere for breast cases, all with 6 MV, 18 MV and (60)Co dose models. In all cases, 95% of the target volumes received at least the prescribed dose with clinical sparing criteria for critical organs being met for all structures that were not wholly or partially contained within the target volume. Improvements in critical organ sparing were found with an increasing number of equidistant (60)Co beams, yet were marginal above 9 beams for H&N, prostate, CNS and lung. Breast cases produced similar plans for 3-7 beams. A helical (60)Co beam geometry achieved similar plan quality as static plans with 11 equidistant (60)Co beams. Furthermore, 18 MV plans were initially found not to provide the same target coverage as 6 MV and (60)Co plans; however, adjusting the trade-offs in the optimization model allowed equivalent target coverage for 18 MV. For plans with comparable target coverage, critical structure sparing was best achieved with 6 MV beams followed closely by (60)Co beams, with 18 MV beams requiring significantly increased dose to critical structures. In this paper, we report in detail on a representative set of results from these experiments. The results of the investigation demonstrate the potential for IMRT radiotherapy employing commercially available (60)Co sources and a double-focused MLC. Increasing the number of equidistant beams beyond 9 was not observed to significantly improve target coverage or critical organ sparing and static plans were found to produce comparable plans to those obtained using a helical tomotherapy treatment delivery when optimized using the same well-tuned convex FMO model. While previous studies have shown that 18 MV plans are equivalent to 6 MV for prostate IMRT, we found that the 18 MV beams actually required more fluence to provide similar quality target coverage.

TG-69: radiographic film for megavoltage beam dosimetry.

TG-69 is a task group report of the AAPM on the use of radiographic film for dosimetry. Radiographic films have been used for radiation dosimetry since the discovery of x-rays and have become an integral part of dose verification for both routine quality assurance and for complex treatments such as soft wedges (dynamic and virtual), intensity modulated radiation therapy (IMRT), image guided radiation therapy (IGRT), and small field dosimetry like stereotactic radiosurgery. Film is convenient to use, spatially accurate, and provides a permanent record of the integrated two dimensional dose distributions. However, there are several challenges to obtaining high quality dosimetric results with film, namely, the dependence of optical density on photon energy, field size, depth, film batch sensitivity differences, film orientation, processing conditions, and scanner performance. Prior to the clinical implementation of a film dosimetry program, the film, processor, and scanner need to be tested to characterize them with respect to these variables. Also, the physicist must understand the basic characteristics of all components of film dosimetry systems. The primary mission of this task group report is to provide guidelines for film selection, irradiation, processing, scanning, and interpretation to allow the physicist to accurately and precisely measure dose with film. Additionally, we present the basic principles and characteristics of film, processors, and scanners. Procedural recommendations are made for each of the steps required for film dosimetry and guidance is given regarding expected levels of accuracy. Finally, some clinical applications of film dosimetry are discussed.

Optimization of equivalent uniform dose using the L-curve criterion.

Optimization of equivalent uniform dose (EUD) in inverse planning for intensity-modulated radiation therapy (IMRT) prevents variation in radiobiological effect between different radiotherapy treatment plans, which is due to variation in the pattern of dose nonuniformity. For instance, the survival fraction of clonogens would be consistent with the prescription when the optimized EUD is equal to the prescribed EUD. One of the problems in the practical implementation of this approach is that the spatial dose distribution in EUD-based inverse planning would be underdetermined because an unlimited number of nonuniform dose distributions can be computed for a prescribed value of EUD. Together with ill-posedness of the underlying integral equation, this may significantly increase the dose nonuniformity. To optimize EUD and keep dose nonuniformity within reasonable limits, we implemented into an EUD-based objective function an additional criterion which ensures the smoothness of beam intensity functions. This approach is similar to the variational regularization technique which was previously studied for the dose-based least-squares optimization. We show that the variational regularization together with the L-curve criterion for the regularization parameter can significantly reduce dose nonuniformity in EUD-based inverse planning.

An exact approach to direct aperture optimization in IMRT treatment planning.

We consider the problem of intensity-modulated radiation therapy (IMRT) treatment planning using direct aperture optimization. While this problem has been relatively well studied in recent years, most approaches employ a heuristic approach to the generation of apertures. In contrast, we use an exact approach that explicitly formulates the fluence map optimization (FMO) problem as a convex optimization problem in terms of all multileaf collimator (MLC) deliverable apertures and their associated intensities. However, the number of deliverable apertures, and therefore the number of decision variables and constraints in the new problem formulation, is typically enormous. To overcome this, we use an iterative approach that employs a subproblem whose optimal solution either provides a suitable aperture to add to a given pool of allowable apertures or concludes that the current solution is optimal. We are able to handle standard consecutiveness, interdigitation and connectedness constraints that may be imposed by the particular MLC system used, as well as jaws-only delivery. Our approach has the additional advantage that it can explicitly account for transmission of dose through the part of an aperture that is blocked by the MLC system, yielding a more precise assessment of the treatment plan than what is possible using a traditional beamlet-based FMO problem. Finally, we develop and test two stopping rules that can be used to identify treatment plans of high clinical quality that are deliverable very efficiently. Tests on clinical head-and-neck cancer cases showed the efficacy of our approach, yielding treatment plans comparable in quality to plans obtained by the traditional method with a reduction of more than 75% in the number of apertures and a reduction of more than 50% in beam-on time, with only a modest increase in computational effort. The results also show that delivery efficiency is very insensitive to the addition of traditional MLC constraints; however, jaws-only treatment requires about a doubling in beam-on time and number of apertures used. Finally, we showed the importance of accounting for transmission effects when assessing or, preferably, optimizing treatment plan quality.

The impact of heterogeneity correction on dosimetric parameters that predict for radiation pneumonitis.

PURPOSE: To determine if heterogeneity correction significantly affects commonly measured dosimetric parameters predicting pulmonary toxicity in patients receiving radiation for lung cancer. METHODS AND MATERIALS: Sixty-eight patients treated for lung cancer were evaluated. The conformal treatment technique mostly employed anteroposterior/posterior-anterior fields and off-cord obliques. The percent total lung volume receiving 20 Gy or higher (V20) and mean lung dose (MLD) were correlated with the incidence of radiation pneumonitis. Parameters from both heterogeneity-corrected and heterogeneity-uncorrected plans were used to assess this risk. RESULTS: Univariate analysis revealed a significant correlation between the development of radiation pneumonitis and both V20 and MLD. A best-fit line to a plot of V20 from the homogeneous plan against the corresponding V20 heterogeneous value produced a slope of 1.00 and zero offset, indicating no difference between the two parameters. For MLD, a similarly significant correlation is seen between the heterogeneous and homogeneous parameters, indicating a 4% difference when correcting for heterogeneity. A significant correlation was also observed between the MLD and V20 parameters (p < 0.0001). CONCLUSIONS: A high degree of correlation exists between heterogeneity-corrected and heterogeneity-uncorrected dosimetric parameters for lung and the risk of developing pneumonitis. Either V20 or MLD predicts the pneumonitis risk with similar effect.

A fourier analysis on the maximum acceptable grid size for discrete proton beam dose calculation.

We developed an analytical method for determining the maximum acceptable grid size for discrete dose calculation in proton therapy treatment plan optimization, so that the accuracy of the optimized dose distribution is guaranteed in the phase of dose sampling and the superfluous computational work is avoided. The accuracy of dose sampling was judged by the criterion that the continuous dose distribution could be reconstructed from the discrete dose within a 2% error limit. To keep the error caused by the discrete dose sampling under a 2% limit, the dose grid size cannot exceed a maximum acceptable value. The method was based on Fourier analysis and the Shannon-Nyquist sampling theorem as an extension of our previous analysis for photon beam intensity modulated radiation therapy [J. F. Dempsey, H. E. Romeijn, J. G. Li, D. A. Low, and J. R. Palta, Med. Phys. 32, 380-388 (2005)]. The proton beam model used for the analysis was a near monoenergetic (of width about 1% the incident energy) and monodirectional infinitesimal (nonintegrated) pencil beam in water medium. By monodirection, we mean that the proton particles are in the same direction before entering the water medium and the various scattering prior to entrance to water is not taken into account. In intensity modulated proton therapy, the elementary intensity modulation entity for proton therapy is either an infinitesimal or finite sized beamlet. Since a finite sized beamlet is the superposition of infinitesimal pencil beams, the result of the maximum acceptable grid size obtained with infinitesimal pencil beam also applies to finite sized beamlet. The analytic Bragg curve function proposed by Bortfeld [T. Bortfeld, Med. Phys. 24, 2024-2033 (1997)] was employed. The lateral profile was approximated by a depth dependent Gaussian distribution. The model included the spreads of the Bragg peak and the lateral profiles due to multiple Coulomb scattering. The dependence of the maximum acceptable dose grid size on the orientation of the beam with respect to the dose grid was also investigated. The maximum acceptable dose grid size depends on the gradient of dose profile and in turn the range of proton beam. In the case that only the phantom scattering was considered and that the beam was aligned with the dose grid, grid sizes from 0.4 to 6.8 mm were required for proton beams with ranges from 2 to 30 cm for 2% error limit at the Bragg peak point. A near linear relation between the maximum acceptable grid size and beam range was observed. For this analysis model, the resolution requirement was not significantly related to the orientation of the beam with respect to the grid.

Fast voxel and polygon ray-tracing algorithms in intensity modulated radiation therapy treatment planning.

We present work on combining three algorithms to improve ray-tracing efficiency in radiation therapy dose computation. The three algorithms include: An improved point-in-polygon algorithm, incremental voxel ray tracing algorithm, and stereographic projection of beamlets for voxel truncation. The point-in-polygon and incremental voxel ray-tracing algorithms have been used in computer graphics and nuclear medicine applications while the stereographic projection algorithm was developed by our group. These algorithms demonstrate significant improvements over the current standard algorithms in peer reviewed literature, i.e., the polygon and voxel ray-tracing algorithms of Siddon for voxel classification (point-in-polygon testing) and dose computation, respectively, and radius testing for voxel truncation. The presented polygon ray-tracing technique was tested on 10 intensity modulated radiation therapy (IMRT) treatment planning cases that required the classification of between 0.58 and 2.0 million voxels on a 2.5 mm isotropic dose grid into 1-4 targets and 5-14 structures represented as extruded polygons (a.k.a. Siddon prisms). Incremental voxel ray tracing and voxel truncation employing virtual stereographic projection was tested on the same IMRT treatment planning cases where voxel dose was required for 230-2400 beamlets using a finite-size pencil-beam algorithm. Between a 100 and 360 fold cpu time improvement over Siddon's method was observed for the polygon ray-tracing algorithm to perform classification of voxels for target and structure membership. Between a 2.6 and 3.1 fold reduction in cpu time over current algorithms was found for the implementation of incremental ray tracing. Additionally, voxel truncation via stereographic projection was observed to be 11-25 times faster than the radial-testing beamlet extent approach and was further improved 1.7-2.0 fold through point-classification using the method of translation over the cross product technique.

A high-speed scintillation-based electronic portal imaging device to quantitatively characterize IMRT delivery.

We have developed an electronic portal imaging device (EPID) employing a fast scintillator and a high-speed camera. The device is designed to accurately and independently characterize the fluence delivered by a linear accelerator during intensity modulated radiation therapy (IMRT) with either step-and-shoot or dynamic multileaf collimator (MLC) delivery. Our aim is to accurately obtain the beam shape and fluence of all segments delivered during IMRT, in order to study the nature of discrepancies between the plan and the delivered doses. A commercial high-speed camera was combined with a terbium-doped gadolinium-oxy-sulfide (Gd2O2S:Tb) scintillator to form an EPID for the unaliased capture of two-dimensional fluence distributions of each beam in an IMRT delivery. The high speed EPID was synchronized to the accelerator pulse-forming network and gated to capture every possible pulse emitted from the accelerator, with an approximate frame rate of 360 frames-per-second (fps). A 62-segment beam from a head-and-neck IMRT treatment plan requiring 68 s to deliver was recorded with our high speed EPID producing approximately 6 Gbytes of imaging data. The EPID data were compared with the MLC instruction files and the MLC controller log files. The frames were binned to provide a frame rate of 72 fps with a signal-to-noise ratio that was sufficient to resolve leaf positions and segment fluence. The fractional fluence from the log files and EPID data agreed well. An ambiguity in the motion of the MLC during beam on was resolved. The log files reported leaf motions at the end of 33 of the 42 segments, while the EPID observed leaf motions in only 7 of the 42 segments. The static IMRT segment shapes observed by the high speed EPID were in good agreement with the shapes reported in the log files. The leaf motions observed during beam-on for step-and-shoot delivery were not temporally resolved by the log files.

Prospects for quantitative two-dimensional radiochromic film dosimetry for low dose-rate brachytherapy sources.

Radiochromic film (RCF) has been shown to be a precise and accurate two-dimensional dosimeter for acute exposure radiation fields. However, "temporal history" mismatch between calibration and brachytherapy films due to RCF dose-rate effects could introduce potentially large uncertainties in low dose-rate (LDR) brachytherapy absolute dose measurement. This article presents a quantitative evaluation of the precision and accuracy of a laser scanner-based RCF-dosimetry system and the effect of the temporal history mismatch in LDR absolute dose measurement. MD-55-2 RCF was used to measure absolute dose for a low dose-rate 137Cs brachytherapy source using both single- and double-exposure techniques. Dose-measurement accuracy was evaluated by comparing RCF to Monte Carlo photon-transport simulation. The temporal history mismatch effect was investigated by examining dependence of RCF accuracy on irradiation-to-densitometry time interval. The predictions of the empirical cumulative dose superposition model (CDSM) were compared with measurements. For the double-exposure technique, the agreement between measurement and Monte Carlo simulation was better than 4% in the 3-60 Gy dose range with measurement precisions (coverage factor k = 1) of <2% and <6% for the doses greater or less than 3 Gy, respectively. The overall uncertainty (k = 1) of dose rate/air-kerma strength measurements achievable by this dosimetry system for a spatial resolution of 0.1 mm is less than 4% for doses greater than 5 Gy. The measured temporal history mismatch systematic error is about 1.8% for a 48 h postexposure time when using the double exposure technique and agrees with CDSM's prediction qualitatively. This work demonstrates that the model MD-55-2 RCF detector has the potential to support quantitative dose measurements about LDR brachytherapy sources with precision and accuracy better than that of previously described dosimeters. The impacts of this work on the future use of new type of RCF were also discussed.

Important considerations for radiochromic film dosimetry with flatbed CCD scanners and EBT GAFCHROMIC film.

In this study, we present three significant artifacts that have the potential to negatively impact the accuracy and precision of film dosimetry measurements made using GAFCHROMIC EBT radiochromic film when read out with CCD flatbed scanners. Films were scanned using three commonly employed instruments: a Macbeth TD932 spot densitometer, an Epson Expression 1680 CCD array scanner, and a Microtek ScanMaker i900 CCD array scanner. For the two scanners we assessed the variation in optical density (OD) of GAFCHROMIC EBT film with scanning bed position, angular rotation of the film with respect to the scan line direction, and temperature inside the scanner due to repeated scanning. Scanning uniform radiochromic films demonstrated a distinct bowing effect in profiles in the direction of the CCD array with a nonuniformity of up to 17%. Profiles along a direction orthogonal to the CCD array demonstrated a 7% variation. A strong angular dependence was found in measurements made with the flatbed scanners; the effect could not be reproduced with the spot densitometer. An IMRT quality assurance film was scanned twice rotating the film 90' between the scans. For films scanned on the Epson scanner, up to 12% variation was observed in unirradiated EBT films rotated between 0 degrees and 90 degrees, which decreased to approximately 8% for EBT films irradiated to 300 cGy. Variations of up to 80% were observed for films scanned with the Microtek scanner. The scanners were found to significantly increase the film temperature with repeated scanning. Film temperature between 18 and 33 degrees C caused OD changes of approximately 7%. Considering these effects, we recommend adherence to a strict scanning protocol that includes: maintaining the orientation of films scanned on flatbed scanners, limiting scanning to the central portion of the scanner bed, and limiting the number of consecutive scans to minimize changes in OD caused by film heating.

Novel breathing motion model for radiotherapy.

PURPOSE: An accurate model of breathing motion under quiet respiration is desirable to obtain the most accurate and conformal dose distributions for mobile lung cancer lesions. On the basis of recent lung motion measurements and the physiologic functioning of the lungs, we have determined that the motion of lung and lung tumor tissues can be modeled as a function of five degrees of freedom, the position of the tissues at a user-specified reference breathing phase, tidal volume and its temporal derivative airflow (tidal volume phase space). Time is an implicit variable in this model. METHODS AND MATERIALS: To test this hypothesis, a mathematical model of motion was developed that described the motion of objects p in the lungs as linear functions of tidal volume and airflow. The position of an object was described relative to its position -->P0 at the reference tidal volume and zero airflow, and the motion of the object was referenced to this position. Hysteresis behavior was hypothesized to be caused by pressure imbalances in the lung during breathing and was, in this model, a function of airflow. The motion was modeled as independent tidal volume and airflow displacement vectors, with the position of the object at time t equal to the vector sum -->rP(t) = -->rv(t) + -->rf(t) where -->rv(t) and -->rf(t) were displacement vectors with magnitudes approximated by linear functions of the tidal volume and airflow. To test this model, we analyzed five-dimensional CT scans (CT scans acquired with simultaneous real-time monitoring of the tidal volume) of 4 patients. The scans were acquired throughout the lungs, but the trajectories were analyzed in the couch positions near the diaphragm. A template-matching algorithm was implemented to identify the positions of the points throughout the 15 scans. In total, 76 points throughout the 4 patients were tracked. The lateral motion of these points was minimal; thus, the model was described in two spatial dimensions, with a total of six parameters necessary to describe the 30 degrees of freedom inherent in the 15 positions. RESULTS: For the 76 evaluated points, the average discrepancy (the distance between the measured and prediction positions) of the 15 locations for each tracked point was 0.75 +/- 0.25 mm, with an average maximal discrepancy of 1.55 +/- 0.54 mm. The average discrepancy was also tabulated as a fraction of the breathing motion. Discrepancies of <10% and 15% of the overall motion occurred in 73% and 95% of the tracked points, respectively. CONCLUSION: The motion tracking algorithms are being improved and automated to provide more motion data to test the models. This may allow a measurement of the motion-fitting parameters throughout the lungs. If the parameters vary smoothly, interpolation may be possible, yielding a continuous mathematical model of the breathing motion throughout the lungs. The utility of the model will depend on its stability as a function of time. If the model is only robust during the measurement session, it may be useful for determining lung function. If it is robust for weeks, it may be useful for treatment planning and gating of lung treatments. The use of tidal volume phase space for characterizing breathing motion appears to have provided, for the first time, the potential for a patient-specific mathematical model of breathing motion.

A fourier analysis of the dose grid resolution required for accurate IMRT fluence map optimization.

We present a theoretical and empirical analysis of the errors associated with the spatial discretization of the dose grid employed in optimized intensity modulated radiation therapy (IMRT) treatment plans. An information theory based Fourier analysis of the accuracy of discrete representations of three-dimensional dose distributions is presented. When applied to beamlet-based IMRT dose distributions, the theory produces analytic integrals that can bound worst case aliasing errors that can occur regardless of the location and orientation of the dose grid. The predictions of this theory are compared to empirical results obtained by solving a linear-programming based fluence-map optimization model to global optimality. A reasonable agreement between worst case estimates and the empirical results is attributed to the fact that the optimization takes advantage of aliasing to produce an optimal plan. We predicted and empirically demonstrated that an isotropic dose grid with <2.5 mm spacing is sufficient to prevent dose errors larger than a percent. However, we noted that in practice this resolution is mostly needed in high-dose target regions. Finally, a multiresolution 2-4-6 mm spacing model was developed and empirically tested where these spacings were applied to targets, structures, and tissue, respectively.

Evaluation of surface and build-up region dose for intensity-modulated radiation therapy in head and neck cancer.

Despite much development, there remains dosimetric uncertainty in the surface and build-up regions in intensity-modulated radiation therapy treatment plans for head and neck cancers. Experiments were performed to determine the dosimetric discrepancies in the surface and build-up region between the treatment planning system (TPS) prediction and experimental measurement using radiochromic film. A head and neck compression film phantom was constructed from two semicylindrical solid water slabs. Treatment plans were generated using two commercial TPSs (PINNACLE3 and CORVUS) for two cases, one with a shallow (approximately 0.5 cm depth) target and another with a deep (approximately 6 cm depth) target. The plans were evaluated for a 54 Gy prescribed dose. For each case, two pieces of radiochromic film were used for dose measurement. A small piece of film strip was placed on the surface and another was inserted within the phantom. Overall, both TPSs showed good agreement with the measurement. For the shallow target case, the dose differences were within +/- 300 cGy (5.6% with respect to the prescribed dose) for PINNACLE3 and +/- 240 cGy (4.4%) for CORVUS in 90% of the region of interest. For the deep target case, the dose differences were +/- 350 (6.5%) for PINNACLE3 and +/- 260 cGy (4.8%) for CORVUS in 90% of the region of interest. However, it was found that there were significant discrepancies from the surface to about 0.2 cm in depth for both the shallow and deep target cases. It was concluded that both TPSs overestimated the surface dose for both shallow and deep target cases. The amount of overestimation ranges from 400 to 1000 cGy (approximately 7.4% to 18.5% with respect to the prescribed dose, 5400 cGy).