Boyd McCurdy

Cone beam computerized tomography: the effect of calibration of the Hounsfield unit number to electron density on dose calculation accuracy for adaptive radiation therapy

The availability of cone beam computerized tomography (CBCT) images at the time of treatment has opened possibilities for dose calculations representing the delivered dose for adaptive radiation therapy. A significant component in the accuracy of dose calculation is the calibration of the Hounsfield unit (HU) number to electron density (ED). The aim of this work is to assess the impact of HU to ED calibration phantom insert composition and phantom volume on dose calculation accuracy for CBCT. CBCT HU to ED calibration curves for different commercial phantoms were measured and compared. The effect of the scattering volume of the phantom on the HU to ED calibration was examined as a function of phantom length and radial diameter. The resulting calibration curves were used at the treatment planning system to calculate doses for geometrically simple phantoms and a pelvic anatomical phantom to compare against measured doses. Three-dimensional dose distributions for the pelvis phantom were calculated using the HU to ED curves and compared using Chi comparisons. The HU to ED calibration curves for the commercial phantoms diverge at densities greater than that of water, depending on the elemental composition of the phantom insert. The effect of adding scatter material longitudinally, increasing the phantom length from 5 cm to 26 cm, was found to be up to 260 HU numbers for the high-density insert. The change in the HU value, by increasing the diameter of the phantom from 18 to 40 cm, was found to be up to 1200 HU for the high-density insert. The effect of phantom diameter on the HU to ED curve can lead to dose differences for 6 MV and 18 MV x-rays under bone inhomogeneities of up to 20% in extreme cases. These results show significant dosimetric differences when using a calibration phantom with materials which are not tissue equivalent. More importantly, the amount of scattering material used with the HU to ED calibration phantom has a significant effect on the dosimetric accuracy, particularly in the radial direction.

Measurement and modeling of the effect of support arm backscatter on dosimetry with a varian EPID.

PURPOSE Amorphous silicon EPIDs have been used for planar dose verification in IMRT treatments for many years. The support arm used to attach some types of EPIDs to linear accelerators can introduce inaccuracies to dosimetry measurements due to the presence of metallic parts in their structures. It is demonstrated that this uncertainty may be as large as approximately 6% of maximum image signal for large fields. In this study, a method has been described to quantify, model and correct for the effect of backscattered radiation from the EPID support arm (E-Arm type, Varian Medical Systems). METHODS Measurements of a support arm backscatter kernel were made using several 1 x 1 cm2 6 MV pencil beam irradiations at a sample of positions over the sensitive area of the EPID in standard clinical setup and repeated with the EPID removed from the support arm but at the same positions. A curve-fit to the subtraction of EPID response obtained on and off the arm was used to define the backscatter kernel. The measured kernel was compared with a backscatter kernel obtained by Monte Carlo simulations with EGS/BEAM code. A backscatter dose prediction using the measured backscatter kernel was added to an existing EPID dose prediction model. The improvement in the agreement of the modified model predictions with EPID measurements for a number of open fields and IMRT beams were investigated by comparison to the original model results. RESULTS Considering all functions tested to find the best functional fit to the data points, a broad Gaussian curve proved to be the optimum fit to the backscatter data. The best fit through the Monte Carlo simulated backscatter kernel was also found to be a Gaussian curve. The maximum decrease in normalized root mean squared deviation of the measured and modeled EPID image profiles for open fields was 13.7% for a 15 x 15 cm2 field with no decrease observed for a 3 x 3 cm2 (the smallest) field as it was not affected by the arm backscatter. Gamma evaluation (2%, 2 mm criteria) showed the improvement in agreement between the model and measurement results when the backscatter was incorporated. The average increase in Gamma pass rate was 2% for head and neck and 1.3% for prostate IMRT fields investigated in this study. CONCLUSIONS The application of the backscatter kernel determined in this study improved the accuracy of dosimetry using a Varian EPID with E-arm for open fields of different sizes: Eight head and neck and seven prostate IMRT fields. Further improvement in the agreement between the model predictions and EPID measurements requires more sophisticated modeling of the backscatter.

Detection and correction for EPID and gantry sag during arc delivery using cine EPID imaging

PURPOSE Electronic portal imaging devices (EPIDs) have been studied and used for pretreatment and in-vivo dosimetry applications for many years. The application of EPIDs for dosimetry in arc treatments requires accurate characterization of the mechanical sag of the EPID and gantry during rotation. Several studies have investigated the effects of gravity on the sag of these systems but each have limitations. In this study, an easy experiment setup and accurate algorithm have been introduced to characterize and correct for the effect of EPID and gantry sag during arc delivery. METHODS Three metallic ball bearings were used as markers in the beam: two of them fixed to the gantry head and the third positioned at the isocenter. EPID images were acquired during a 360° gantry rotation in cine imaging mode. The markers were tracked in EPID images and a robust in-house developed MATLAB code was used to analyse the images and find the EPID sag in three directions as well as the EPID + gantry sag by comparison to the reference gantry zero image. The algorithm results were then tested against independent methods. The method was applied to compare the effect in clockwise and counter clockwise gantry rotations and different source-to-detector distances (SDDs). The results were monitored for one linear accelerator over a course of 15 months and six other linear-accelerators from two treatment centers were also investigated using this method. The generalized shift patterns were derived from the data and used in an image registration algorithm to correct for the effect of the mechanical sag in the system. The Gamma evaluation (3%, 3 mm) technique was used to investigate the improvement in alignment of cine EPID images of a fixed field, by comparing both individual images and the sum of images in a series with the reference gantry zero image. RESULTS The mechanical sag during gantry rotation was dependent on the gantry angle and was larger in the in-plane direction, although the patterns were not identical for various linear-accelerators. The reproducibility of measurements was within 0.2 mm over a period of 15 months. The direction of gantry rotation and SDD did not affect the results by more than 0.3 mm. Results of independent tests agreed with the algorithm within the accuracy of the measurement tools. When comparing summed images, the percentage of points with Gamma index <1 increased from 85.4% to 94.1% after correcting for the EPID sag, and to 99.3% after correction for gantry + EPID sag. CONCLUSIONS The measurement method and algorithms introduced in this study use cine-images, are highly accurate, simple, fast, and reproducible. It tests all gantry angles and provides a suitable automatic analysis and correction tool to improve EPID dosimetry and perform comprehensive linac QA for arc treatments.

Gantry-angle resolved VMAT pretreatment verification using EPID image prediction

PURPOSE Pretreatment verification of volumetric modulated arc therapy (VMAT) dose delivery with electronic portal imaging device (EPID) uses images integrated over the entire delivery or over large subarcs. This work aims to develop a new method for gantry-angle-resolved verification of VMAT dose delivery using EPID. METHODS An EPID dose prediction model was used to calculate EPID images as a function of gantry angle for eight prostate patient deliveries. EPID image frames at 7.5 frames per second were acquired during delivery via a frame-grabber system. The gantry angle for each image was encoded in kV frames which were synchronized to the MV frames. Gamma analysis results as a function of gantry angle were assessed by integrating the frames over 2° subarcs with an angle-to-agreement tolerance of 0.5° about the measured image angle. RESULTS The model agreed with EPID images integrated over the entire delivery with average Gamma pass-rates at 2%, 2 mm of 99.7% (10% threshold). The accuracy of the kV derived gantry angle for each image was found to be 0.1° (1 SD) using a phantom test. For the gantry-resolved analysis all Gamma pass-rates were greater than 90% at 3%, 3 mm criteria (with only two exceptions), and more than 90% had a 95% pass-rate, with an average of 97.3%. The measured gantry angle lagged behind the predicted angle by a mean of 0.3°±0.3°, with a maximum lag of 1.3°. CONCLUSIONS The method provides a comprehensive and highly efficient pretreatment verification of VMAT delivery using EPID. Dose delivery accuracy is assessed as a function of gantry angle to ensure accurate treatment.

The influence of a novel transmission detector on 6 MV x-ray beam characteristics.

The purpose of this work was to investigate the influence of a new transmission detector on 6 MV x-ray beam properties. The device, COMPASS (IBA Dosimetry, Germany), contains 1600 plane parallel ionization chambers with a detector spacing of 6.5 mm and an active volume of 0.02 cm3. Surface dose measurements were carried out using a Markus chamber and radiochromic film for a range of field sizes and source-to-surface distances (SSDs). The surface dose and dose in the build-up region for COMPASS fields were compared to open fields. For moderately narrow beam geometric conditions, the increase in surface dose was small. For the largest field size investigated (20x20 cm2) at a 90 cm SSD, the surface dose with the detector was 34.9% versus 26.8% in the open field. However, the increase in surface dose in COMPASS fields was less than that observed with a standard block tray in the field (38.7% in the above example). It was found that beyond dmax, the difference in relative dose (profiles and PDDs) between open and COMPASS fields was insignificant. The mean transmission factor of the detector was 0.967 (standard deviation=0.002) measured over a range of field sizes from 3x3 to 20x20 cm2 at SSDs from 70 cm to 90 cm. In summary, the transmission detector was found to increase the relative dose in the buildup region but had a negligible effect on the beam parameters beyond dmax.

A system for EPID-based real-time treatment delivery verification during dynamic IMRT treatment

PURPOSE To design and develop a real-time electronic portal imaging device (EPID)-based delivery verification system for dynamic intensity modulated radiation therapy (IMRT) which enables detection of gross treatment delivery errors before delivery of substantial radiation to the patient. METHODS The system utilizes a comprehensive physics-based model to generate a series of predicted transit EPID image frames as a reference dataset and compares these to measured EPID frames acquired during treatment. The two datasets are using MLC aperture comparison and cumulative signal checking techniques. The system operation in real-time was simulated offline using previously acquired images for 19 IMRT patient deliveries with both frame-by-frame comparison and cumulative frame comparison. Simulated error case studies were used to demonstrate the system sensitivity and performance. RESULTS The accuracy of the synchronization method was shown to agree within two control points which corresponds to approximately ∼1% of the total MU to be delivered for dynamic IMRT. The system achieved mean real-time gamma results for frame-by-frame analysis of 86.6% and 89.0% for 3%, 3 mm and 4%, 4 mm criteria, respectively, and 97.9% and 98.6% for cumulative gamma analysis. The system can detect a 10% MU error using 3%, 3 mm criteria within approximately 10 s. The EPID-based real-time delivery verification system successfully detected simulated gross errors introduced into patient plan deliveries in near real-time (within 0.1 s). CONCLUSIONS A real-time radiation delivery verification system for dynamic IMRT has been demonstrated that is designed to prevent major mistreatments in modern radiation therapy.

Model-based prediction of portal dose images during patient treatment

PURPOSE Dosimetric verification of radiation therapy is crucial when delivering complex treatments like intensity modulated radiation therapy (IMRT) or volumetric modulated arc therapy (VMAT. Pretreatment verification, characterized by methods applied without the patient present and before the treatment start date, is typically carried out at most centers. In vivo dosimetric verification, characterized by methods applied with the patient present, is not commonly carried out in the clinic. This work presents a novel, model-based EPID dosimetry method that could be used for routine clinical in vivo patient treatment verification. METHODS The authors integrated a detailed fluence model with a patient scatter prediction model that uses a superposition of scatter energy fluence kernels, generated via Monte Carlo techniques, to determine patient scatter fluence delivered to the EPID. The total dose to the EPID was calculated using the sum of convolutions of the calculated energy fluence distribution entering the EPID with monoenergetic dose kernels, specific to the a-Si EPID. Measured images with simple, square fields delivered to slab phantoms were validated against predicted images. Measured and predicted images acquired during the delivery of IMRT fields to slabs and an anthropomorphic phantom were compared using the χ-comparison for 3% dose difference and 3 mm distance-to-agreement criteria. RESULTS Predicted and measured images of the square fields with slabs in the field agreed within 2.5%. Predicted portal dose images of clinical IMRT fields delivered to slabs and an anthropomorphic phantom agreed with measured images within 3% and 3 mm for an average of at least 97% of the infield pixels (defined as >10% maximum field dose) for each case, over all fields. CONCLUSIONS This work presents the first validation of the integration of a comprehensive fluence model with a patient and EPID radiation transport model that accounts for patient transmission, including complex factors such as patient scatter and the energy response of the a-Si detector. The portal dose image prediction model satisfies the 3% and 3 mm criteria for IMRT fields delivered to slab phantoms and could be used for patient treatment verification.

A feasibility study to investigate the use of thin-plate splines to account for prostate deformation

Image registration is an important step in the radiotherapy treatment planning process. It provides a method of fusing different types of diagnostic imaging information. One such application is to combine magnetic resonance spectroscopic images (MRSI) of the prostate with anatomical MRI and/or computed tomography images that are routinely used in the radiation treatment planning of prostate cancer. MRSI provides in vivo information related to the underlying metabolic activity of tissues, and can be related to the presence of cancer. However, the inflated endorectal coil required during MRS imaging poses a potential problem by deforming the prostate when it is filled with approximately 100 cm3 of air during image acquisition. This pushes the prostate superiorly/anteriorly, deforming the prostate and consequently the spectroscopic imaging data in a nonlinear manner. In this application, the coil-deformed MRS images are warped back to a non-deformed state, using a single data set. A nonlinear warping algorithm is presented to achieve this. Results indicate that the algorithm attains an accuracy of 97% (4 cm3 difference) when reproducing the total prostate volume compared to a Radiation Oncologist defined prostate volume. This difference is slightly smaller than the measured intra-operator variance of +/-1.5 cm3 (deflated coil) and the measured algorithm variance of +/-1.0 cm3. Additionally, intraprostatic nodules were used to assess the accuracy of the warping algorithm in regions inside the prostate. While choosing anatomical tie points along the external prostate surface, analysis of the nodules revealed the algorithm accuracy reduced to 63-93%.

First Experience With Real-Time EPID-Based Delivery Verification During IMRT and VMAT Sessions

PURPOSE Gantry-mounted megavoltage electronic portal imaging devices (EPIDs) have become ubiquitous on linear accelerators. WatchDog is a novel application of EPIDs, in which the image frames acquired during treatment are used to monitor treatment delivery in real time. We report on the preliminary use of WatchDog in a prospective study of cancer patients undergoing intensity modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT) and identify the challenges of clinical adoption. METHODS AND MATERIALS At the time of submission, 28 cancer patients (head and neck, pelvis, and prostate) undergoing fractionated external beam radiation therapy (24 IMRT, 4 VMAT) had ≥1 treatment fraction verified in real time (131 fractions or 881 fields). EPID images acquired continuously during treatment were synchronized and compared with model-generated transit EPID images within a frame time (∼0.1 s). A χ comparison was performed to cumulative frames to gauge the overall delivery quality, and the resulting pass rates were reported graphically during treatment delivery. Every frame acquired (500-1500 per fraction) was saved for postprocessing and analysis. RESULTS The system reported the mean ± standard deviation in real time χ 91.1% ± 11.5% (83.6% ± 13.2%) for cumulative frame χ analysis with 4%, 4 mm (3%, 3 mm) criteria, global over the integrated image. CONCLUSIONS A real-time EPID-based radiation delivery verification system for IMRT and VMAT has been demonstrated that aims to prevent major mistreatments in radiation therapy.

PARETO: A novel evolutionary optimization approach to multiobjective IMRT planning

PURPOSE In radiation therapy treatment planning, the clinical objectives of uniform high dose to the planning target volume (PTV) and low dose to the organs-at-risk (OARs) are invariably in conflict, often requiring compromises to be made between them when selecting the best treatment plan for a particular patient. In this work, the authors introduce Pareto-Aware Radiotherapy Evolutionary Treatment Optimization (pareto), a multiobjective optimization tool to solve for beam angles and fluence patterns in intensity-modulated radiation therapy (IMRT) treatment planning. METHODS pareto is built around a powerful multiobjective genetic algorithm (GA), which allows us to treat the problem of IMRT treatment plan optimization as a combined monolithic problem, where all beam fluence and angle parameters are treated equally during the optimization. We have employed a simple parameterized beam fluence representation with a realistic dose calculation approach, incorporating patient scatter effects, to demonstrate feasibility of the proposed approach on two phantoms. The first phantom is a simple cylindrical phantom containing a target surrounded by three OARs, while the second phantom is more complex and represents a paraspinal patient. RESULTS pareto results in a large database of Pareto nondominated solutions that represent the necessary trade-offs between objectives. The solution quality was examined for several PTV and OAR fitness functions. The combination of a conformity-based PTV fitness function and a dose-volume histogram (DVH) or equivalent uniform dose (EUD) -based fitness function for the OAR produced relatively uniform and conformal PTV doses, with well-spaced beams. A penalty function added to the fitness functions eliminates hotspots. Comparison of resulting DVHs to those from treatment plans developed with a single-objective fluence optimizer (from a commercial treatment planning system) showed good correlation. Results also indicated that pareto shows promise in optimizing the number of beams. CONCLUSIONS This initial evaluation of the evolutionary optimization software tool pareto for IMRT treatment planning demonstrates feasibility and provides motivation for continued development. Advantages of this approach over current commercial methods for treatment planning are many, including: (1) fully automated optimization that avoids human controlled iterative optimization and potentially improves overall process efficiency, (2) formulation of the problem as a true multiobjective one, which provides an optimized set of Pareto nondominated solutions refined over hundreds of generations and compiled from thousands of parameter sets explored during the run, and (3) rapid exploration of the final nondominated set accomplished by a graphical interface used to select the best treatment option for the patient.