Evaluation of diffusion-weighted MRI and geometric distortion on a 0.35T MR-LINAC at multiple gantry angles.

Diffusion-weighted imaging (DWI) provides a valuable diagnostic tool for tumor evaluation. Yet, it is difficult to acquire daily MRI data sets in the traditional radiotherapy clinical setting due to patient burden and limited resources. However, integrated MRI radiotherapy treatment systems facilitate daily functional MRI acquisitions like DWI during treatment exams. Before ADC values from MR-RT systems can be used clinically their reproducibility and accuracy must be quantified. This study used a NIST traceable DWI phantom to verify ADC values acquired on a 0.35 T MR-LINAC system at multiple gantry angles. A diffusion-weighted echo planar imaging sequence was used for all image acquisitions, with b-values of 0, 500, 900, 2000 s/mm2 for the 1.5 T and 3.0 T systems and 0, 200, 500, 800 s/mm2 for the 0.35 T system. Images were acquired at multiple gantry angles on the MR-LINAC system from 0° to 330° in 30° increments to assess the impact of gantry angle on geometric distortion and ADC values. CT images, and three fiducial markers were used as ground truth for geometric distortion measurements. The distance between fiducial markers increased by as much as 7.2 mm on the MR-LINAC at gantry angle 60°. ADC values of deionized water vials from the 1.5 T and 3.0 T systems were 8.30 × 10-6  mm2 /s and -0.85 × 10-6  mm2 /s off, respectively, from the expected value of 1127 × 10-6  mm2 /s. The MR-LINAC system provided an ADC value of the pure water vials that was -116.63 × 10-6  mm2 /s off from the expected value of 1127 × 10-6  mm2 /s. The MR-LINAC also showed a variation in ADC across all gantry angles of 33.72 × 10-6  mm2 /s and 20.41 × 10-6  mm2 /s for the vials with expected values of 1127 × 10-6  mm2 /s and 248 × 10-6  mm2 /s, respectively. This study showed that variation of the ADC values and geometric information on the 0.35 T MR-LINAC system was dependent on the gantry angle at acquisition.

Filmless quality assurance of a Leksell Gamma Knife® Icon™.

PURPOSE: The annual quality assurance (QA) of Leksell Gamma Knife® (LGK) systems are typically performed using films. Film is a good candidate for small field dosimetry due to its high spatial resolution and availability. However, there are multiple challenges with using film; film does not provide real-time measurement and requires batch-specific calibration. Our findings show that active detector-based QA can simplify the procedure and save time without loss of accuracy. METHODS: Annual QA tests for a LGK Icon™ system were performed using both film-based and filmless techniques. Output calibration, relative output factors (ROF), radiation profiles, sector uniformity/source counting, and verification of the unit center point (UCP) and radiation focal point (RFP) coincidence tests were performed. Radiochromic films, two ionization chambers, and a synthetic diamond detector were used for the measurements. Results were compared and verified with the treatment planning system (TPS). RESULTS: The measured dose rate of the LGK Icon was within 0.4% of the TPS value set at the time of commissioning using an ionization chamber. ROF for the 8 and 4-mm collimators were found to be 0.3% and 1.8% different from TPS values using the MicroDiamond detector and 2.6% and 1.9% different for film, respectively. Excellent agreement was found between TPS and measured dose profiles using the MicroDiamond detector which was within 1%/1 mm vs 2%/1 mm for film. Sector uniformity was found to be within 1% for all eight sectors measured using an ionization chamber. Verification of UCP and RFP coincidence using the MicroDiamond detector and pinprick film test was within 0.3 mm at isocenter for both. CONCLUSION: The annual QA of a LGK Icon was successfully performed by employing filmless techniques. Comparable results were obtained using radiochromic films. Utilizing active detectors instead of films simplifies the QA process and saves time without loss of accuracy.

Phase I/II Trial of Electrophysiology-Guided Noninvasive Cardiac Radioablation for Ventricular Tachycardia.

BACKGROUND: Case studies have suggested the efficacy of catheter-free, electrophysiology-guided noninvasive cardiac radioablation for ventricular tachycardia (VT) using stereotactic body radiation therapy, although prospective data are lacking. METHODS: We conducted a prospective phase I/II trial of noninvasive cardiac radioablation in adults with treatment-refractory episodes of VT or cardiomyopathy related to premature ventricular contractions (PVCs). Arrhythmogenic scar regions were targeted by combining noninvasive anatomic and electric cardiac imaging with a standard stereotactic body radiation therapy workflow followed by delivery of a single fraction of 25 Gy to the target. The primary safety end point was treatment-related serious adverse events in the first 90 days. The primary efficacy end point was any reduction in VT episodes (tracked by indwelling implantable cardioverter defibrillators) or any reduction in PVC burden (as measured by a 24-hour Holter monitor) comparing the 6 months before and after treatment (with a 6-week blanking window after treatment). Health-related quality of life was assessed using the Short Form-36 questionnaire. RESULTS: Nineteen patients were enrolled (17 for VT, 2 for PVC cardiomyopathy). Median noninvasive ablation time was 15.3 minutes (range, 5.4-32.3). In the first 90 days, 2/19 patients (10.5%) developed a treatment-related serious adverse event. The median number of VT episodes was reduced from 119 (range, 4-292) to 3 (range, 0-31; P<0.001). Reduction was observed for both implantable cardioverter defibrillator shocks and antitachycardia pacing. VT episodes or PVC burden were reduced in 17/18 evaluable patients (94%). The frequency of VT episodes or PVC burden was reduced by 75% in 89% of patients. Overall survival was 89% at 6 months and 72% at 12 months. Use of dual antiarrhythmic medications decreased from 59% to 12% ( P=0.008). Quality of life improved in 5 of 9 Short Form-36 domains at 6 months. CONCLUSIONS: Noninvasive electrophysiology-guided cardiac radioablation is associated with markedly reduced ventricular arrhythmia burden with modest short-term risks, reduction in antiarrhythmic drug use, and improvement in quality of life. CLINICAL TRIAL REGISTRATION: URL: https://www.clinicaltrials.gov/ . Unique identifier: NCT02919618.

Noninvasive Cardiac Radiation for Ablation of Ventricular Tachycardia.

BACKGROUND: Recent advances have enabled noninvasive mapping of cardiac arrhythmias with electrocardiographic imaging and noninvasive delivery of precise ablative radiation with stereotactic body radiation therapy (SBRT). We combined these techniques to perform catheter-free, electrophysiology-guided, noninvasive cardiac radioablation for ventricular tachycardia. METHODS: We targeted arrhythmogenic scar regions by combining anatomical imaging with noninvasive electrocardiographic imaging during ventricular tachycardia that was induced by means of an implantable cardioverter-defibrillator (ICD). SBRT simulation, planning, and treatments were performed with the use of standard techniques. Patients were treated with a single fraction of 25 Gy while awake. Efficacy was assessed by counting episodes of ventricular tachycardia, as recorded by ICDs. Safety was assessed by means of serial cardiac and thoracic imaging. RESULTS: From April through November 2015, five patients with high-risk, refractory ventricular tachycardia underwent treatment. The mean noninvasive ablation time was 14 minutes (range, 11 to 18). During the 3 months before treatment, the patients had a combined history of 6577 episodes of ventricular tachycardia. During a 6-week postablation "blanking period" (when arrhythmias may occur owing to postablation inflammation), there were 680 episodes of ventricular tachycardia. After the 6-week blanking period, there were 4 episodes of ventricular tachycardia over the next 46 patient-months, for a reduction from baseline of 99.9%. A reduction in episodes of ventricular tachycardia occurred in all five patients. The mean left ventricular ejection fraction did not decrease with treatment. At 3 months, adjacent lung showed opacities consistent with mild inflammatory changes, which had resolved by 1 year. CONCLUSIONS: In five patients with refractory ventricular tachycardia, noninvasive treatment with electrophysiology-guided cardiac radioablation markedly reduced the burden of ventricular tachycardia. (Funded by Barnes-Jewish Hospital Foundation and others.).

Direct tumor visual feedback during free breathing in 0.35T MRgRT.

To present a tumor motion control system during free breathing using direct tumor visual feedback to patients in 0.35 T magnetic resonance-guided radiotherapy (MRgRT). We present direct tumor visualization to patients by projecting real-time cine MR images on an MR-compatible display system inside a 0.35 T MRgRT bore. The direct tumor visualization included anatomical images with a target contour and an auto-segmented gating contour. In addition, a beam-status sign was added for patient guidance. The feasibility was investigated with a six-patient clinical evaluation of the system in terms of tumor motion range and beam-on time. Seven patients without visual guidance were used for comparison. Positions of the tumor and the auto-segmented gating contour from the cine MR images were used in probability analysis to evaluate tumor motion control. In addition, beam-on time was recorded to assess the efficacy of the visual feedback system. The direct tumor visualization system was developed and implemented in our clinic. The target contour extended 3 mm outside of the gating contour for 33.6 ± 24.9% of the time without visual guidance, and 37.2 ± 26.4% of the time with visual guidance. The average maximum motion outside of the gating contour was 14.4 ± 11.1 mm without and 13.0 ± 7.9 mm with visual guidance. Beam-on time as a percentage was 43.9 ± 15.3% without visual guidance, and 48.0 ± 21.2% with visual guidance, but was not significantly different (P = 0.34). We demonstrated the clinical feasibility and potential benefits of presenting direct tumor visual feedback to patients in MRgRT. The visual feedback allows patients to visualize and attempt to minimize tumor motion in free breathing. The proposed system and associated clinical workflow can be easily adapted for any type of MRgRT.

Comprehensive validation of halcyon 2.0 plans and the implementation of patient specific QA with multiple detector platforms.

PURPOSE: To perform a comprehensive validation of plans generated on a preconfigured Halcyon 2.0 with preloaded beam model, including evaluations of new features and implementing the patient specific quality assurance (PSQA) process with multiple detectors. METHODS: A total of 56 plans were generated in Eclipse V15.6 (Varian Medical System) with a preconfigured Halcyon treatment machine. Ten plans were developed via the AAPM TG-119 test suite with both IMRT and VMAT techniques. 34 clinically treated plans using C-arm LINAC from 24 patients were replanned on Halcyon using IMRT or VMAT techniques for a variety of sites including: brain, head and neck, lung, breast, abdomen, and pelvis. Six of those plans were breast VMAT plans utilizing the extended treatment field technique available with Halcyon 2.0. The dynamically flattened beam (DFB), another new feature on Halcyon 2.0, was also used for an AP/PA spine and four field box pelvis, as well as ten 3D breast plans. All 56 plans were measured with an ion chamber (IC), film, portal dosimetry (PD), ArcCHECK, and Delta4. Tolerance and action limits were calculated and compared to the recommendations of TG-218. RESULTS: TG-119 IC and film confidence limits met those set by the task group, except for IMRT target point dose. Forty-four of 46 clinical plans were within 3% for IC measurements. Average gamma passing rates with 3% dose difference and 2mm distance-to-agreement for IMRT/VMAT plans were: Film - 96.8%, PD - 99.9%, ArcCHECK - 99.1%, and Delta4 - 99.2%. Calculated action limits were: Film - 86.3%, PD - 98.4%, ArcCHECK - 96.1%, and Delta4 - 95.7%. Extended treatment field technique was fully validated and 3D plans with DFB had similar results to IMRT/VMAT plans. CONCLUSION: Halcyon plan deliveries were verified with multiple measurement devices. New features of Halcyon 2.0 were also validated. Traditional PSQA techniques and process specific tolerance and action limits were successfully implemented.

Automated and robust beam data validation of a preconfigured ring gantry linear accelerator using a 1D tank with synchronized beam delivery and couch motions.

PURPOSE: To develop an efficient and automated methodology for beam data validation for a preconfigured ring gantry linear accelerator using scripting and a one-dimensional (1D) tank with automated couch motions. MATERIALS AND METHODS: Using an application programming interface, a program was developed to allow the user to choose a set of beam data to validate with measurement. Once selected the program generates a set of instructions for radiation delivery with synchronized couch motions for the linear accelerator in the form of an extensible markup language (XML) file to be delivered on the ring gantry linear accelerator. The user then delivers these beams while measuring with the 1D tank and data logging electrometer. The program also automatically calculates this set of beams on the measurement geometry within the treatment planning system (TPS) and extracts the corresponding calculated dosimetric data for comparison to measurement. Once completed the program then returns a comparison of the measurement to the predicted result from the TPS to the user and prints a report. In this work lateral, longitudinal, and diagonal profiles were taken for fields sizes of 6 × 6, 8 × 8, 10 × 10, 20 × 20, and 28 × 28 cm2 at depths of 1.3, 5, 10, 20, and 30 cm. Depth dose profiles were taken for all field sizes. RESULTS: Using this methodology, the TPS was validated to agree with measurement. All compared points yielded a gamma value less than 1 for a 1.5%/1.5 mm criteria (100% passing rate). Off axis profiles had >98.5% of data points producing a gamma value <1 with a 1%/1 mm criteria. All depth profiles produced 100% of data points with a gamma value <1 with a 1%/1 mm criteria. All data points measured were within 1.5% or 2 mm distance to agreement. CONCLUSIONS: This methodology allows for an increase in automation in the beam data validation process. Leveraging the application program interface allows the user to use a single system to create the measurement files, predict the result, and then compare to actual measurement increasing efficiency and reducing the chance for user input errors.

Evaluation of a new secondary dose calculation software for Gamma Knife radiosurgery.

Current available secondary dose calculation software for Gamma Knife radiosurgery falls short in situations where the target is shallow in depth or when the patient is positioned with a gamma angle other than 90°. In this work, we evaluate a new secondary calculation software which utilizes an innovative method to handle nonstandard gamma angles and image thresholding to render the skull for dose calculation. 800 treatment targets previously treated with our GammaKnife Icon system were imported from our treatment planning system (GammaPlan 11.0.3) and a secondary dose calculation was conducted. The agreement between the new calculations and the TPS were recorded and compared to the original secondary dose calculation agreement with the TPS using a Wilcoxon Signed Rank Test. Further comparisons using a Mann-Whitney test were made for targets treated at a 90° gamma angle against those treated with either a 70 or 110 gamma angle for both the new and commercial secondary dose calculation systems. Correlations between dose deviations from the treatment planning system against average target depth were evaluated using a Kendall's Tau correlation test for both programs. The Wilcoxon Signed Rank Test indicated a significant difference in the agreement between the two secondary calculations and the TPS, with a P-value < 0.0001. With respect to patients treated at nonstandard gamma angles, the new software was largely independent of patient setup, while the commercial software showed a significant dependence (P-value < 0.0001). The new secondary dose calculation software showed a moderate correlation with calculation depth, while the commercial software showed a weak correlation (Tau = -.322 and Tau = -.217 respectively). Overall, the new secondary software has better agreement with the TPS than the commercially available secondary calculation software over a range of diverse treatment geometries.

A Monte Carlo based analytic model of the in-room neutron ambient dose equivalent for a Mevion gantry-mounted passively scattered proton system.

The goal of this study was to develop a Monte Carlo (MC)-based analytical model that can predict the in-room ambient dose equivalent from a Mevion gantry-mounted passively scattered proton system. The Mevion S250 and treatment vault were simulated using the MCNPX MC code. The results of the in-room neutron dose measurements, using an FHT 762 WENDI-II detector, were employed to benchmark the MC-derived values. After tuning the MCNPX MC code, for the same beam delivery parameters, the code was used to calculate the neutron spectra and ambient dose equivalent in the vault and at varying angles from the isocenter. Then, based on the calculations, an analytical model was reconstructed and data were fitted to derive the model parameters at 95% confidence intervals (CI). The MCNPX codes were tuned to within about 19% of the measured values for most of the measurements in the vault. For the maze, up to 0.08 mSv Gy-1 discrepancies were found between the experimental measurements and MCNPX calculated results. The analytical model showed up to 18% discrepancy for distances between 100 and 600 cm from the isocenter compared to the MC calculations. The model may underestimate the neutron ambient dose equivalent up to 21% for distances less than 100 cm from the isocenter. The proposed analytical model can be used to estimate the contribution of the secondary neutron dose from the Mevion S250 for the design of local shielding inside the proton therapy treatment vault.

Direct comparison between surface imaging and orthogonal radiographic imaging for SRS localization in phantom.

PURPOSE: Surface imaging (SI) offers a nonionizing, near real time alternative to radiographic imaging for intrafraction radiosurgery localization. In this work, we systematically compared a commercial SI system vs a commercial room mounted x-ray localization system in phantom. METHODS: An anthropomorphic head phantom with fiducial markers was imaged with linear accelerator on-board x-ray imaging, SI, and room mounted x-ray imaging (RM) at ±45° and ±90° couch angles for three different head tilts and six different isocenters (72 total positions). The shifts generated by the three systems were compared as functions of couch angle, head tilt, and isocenter position with the on-board imaging shifts used as ground truth. Two sample Kolmogorov-Smirnov tests were used to evaluate equivalence of the groups. RESULTS: The magnitude of the displacement vectors for RM minus on-board imaging and SI minus on-board imaging over all 72 phantom positions were 0.7 ± 0.3 mm for both cases. The RM and SI showed no significant difference based on couch angle or isocenter position. Both systems showed decreasing accuracy with increasing couch angle, but both systems agreed with ground truth to <=1.1 mm at all couch angles. The exaggerated chin-up head orientation showed significantly different shifts for SI and RM based on increased variance in the SI measurements, although both had submillimeter accuracy on average. The standard deviation of the real time SI displacement vector was <0.06 mm over all measurements, during which the on-board imaging panels partially blocked the lateral camera pods for half the time. CONCLUSIONS: RM and SI showed similar accuracy over measurements at 72 different phantom positions. SI showed minimal performance loss with camera pods blocked. SI is a feasible option for intra-fraction radiosurgery localization based on these phantom measurements.

MRI quality control for low-field MR-IGRT systems: Lessons learned.

PURPOSE: To present lessons learned from magnetic resonance imaging (MRI) quality control (QC) tests for low-field MRI-guided radiation therapy (MR-IGRT) systems. METHODS: MRI QC programs were established for low-field MRI-60 Co and MRI-Linac systems. A retrospective analysis of MRI subsystem performance covered system commissioning, operations, maintenance, and quality control. Performance issues were classified into three groups: (a) Image noise and artifact; (b) Magnetic field homogeneity and linearity; and (c) System reliability and stability. RESULTS: Image noise and artifacts were attributed to room noise sources, unsatisfactory system cabling, and broken RF receiver coils. Gantry angle-dependent magnetic field inhomogeneities were more prominent on the MRI-Linac due to the high volume of steel shielding in the gantry. B0 inhomogeneities measured in a 24-cm spherical phantom were <5 ppm for both MR-IGRT systems after using MRI gradient offset (MRI-GO) compensation on the MRI-Linac. However, significant signal dephasing occurred on the MRI-Linac while the gantry was rotating. Spatial integrity measurements were sensitive to gradient calibration and vulnerable to shimming. The most common causes of MR-IGRT system interruptions were software disconnects between the MRI and radiation therapy delivery subsystems caused by patient table, gantry, and multi-leaf collimator (MLC) faults. The standard deviation (SD) of the receiver coil signal-to-noise ratio was 1.83 for the MRI-60 Co and 1.53 for the MRI-Linac. The SD of the deviation from the mean for the Larmor frequency was 1.41 ppm for the MRI-60 Co and 1.54 ppm for the MRI-Linac. The SD of the deviation from the mean for the transmitter reference amplitude was 0.90% for the MRI-60 Co and 1.68% for the MRI-Linac. High SDs in image stability data corresponded to reports of spike noise. CONCLUSIONS: There are significant technological challenges associated with implementing and maintaining MR-IGRT systems. Most of the performance issues were identified and resolved during commissioning.

Field-in-field breast planning for a jawless, double-stack MLC LINAC using flattening-filter-free beams.

BACKGROUND: This study intends to develop an efficient field-in-field (FiF) planning technique with the Eclipse treatment planning system (TPS) to determine the feasibility of using the Halcyon treatment delivery system for 3D treatment of breast cancer. METHODS: Ten treatment plans were prepared on the Halcyon treatment planning system and compared to the same patients' clinically delivered TrueBeam plans which used flattened 6 MV and 10 MV beams. Patients selected for this study were treated via simple, tangential breast irradiation and did not receive radiotherapy of the supraclavicular or internal mammary lymph nodes. Planning target volumes (PTV) volumes ranged from 519 cc to 1211 cc with a mean target volume of 877 cc. Several planning techniques involving collimator, gantry rotation, and number of FiF segments were investigated as well as the use of the dynamically flattened beam (DFB) - a predefined MLC pattern that is designed to provide a flattened beam profile at 10 cm depth on a standard water phantom. For comparison, the clinically delivered TrueBeam plans remained unaltered except for normalization of the target coverage to more readily compare the two treatment delivery techniques. RESULTS: Using the physician defined PTV, normalized such that 98% of the volume was covered by 95% of the prescribed dose, the Halcyon plans were deemed clinically acceptable and comparable to the TrueBeam plans by the radiation oncologist. Resulting average global maximum doses in the test patients were identical between the TrueBeam and Halcyon plans (108% of Rx) and a mean PTV dose of 102.5% vs 101.6%, respectively. CONCLUSIONS: From this study a practical and efficient planning method for delivering 3D conformal breast radiotherapy using the Halcyon linear accelerator has been developed. When normalized to the clinically desired coverage, hot spots were maintained to acceptable levels and overall plan quality was comparable to plans delivered on conventional C-arm LINACs.

Characterization and validation of an intra-fraction motion management system for masked-based radiosurgery.

PURPOSE: Characterize the intra-fraction motion management (IFMM) system found on the Gamma Knife Icon (GKI), including spatial accuracy, latency, temporal performance, and overall effect on delivered dose. METHODS: A phantom was constructed, consisting of a three-axis translation mount, a remote motorized flipper, and a thermoplastic sphere surrounding a radiation detector. An infrared marker was placed on the translation mount secured to the flipper. The spatial accuracy of the IFMM was measured via the translation mount in all Cartesian planes. The detector was centered at the radiation focal point. A remote signal was used to move the marker out of the IFMM tolerance and pause the beam. A two-channel electrometer was used to record the signals from the detector and the flipper when motion was signaled. These signals determined the latency and temporal performance of the GKI. RESULTS: The spatial accuracy of the IFMM was found to be <0.1 mm. The measured latency was <200 ms. The dose difference with five interruptions was <0.5%. CONCLUSION: This work provides a quantitative characterization of the GKI IFMM system as required by the Nuclear Regulatory Commission. This provides a methodology for GKI users to satisfy these requirements using common laboratory equipment in lieu of a commercial solution.

Equivalency of beam scan data collection using a 1D tank and automated couch movements to traditional 3D tank measurements.

This work shows the feasibility of collecting linear accelerator beam data using just a 1-D water tank and automated couch movements with the goal to maximize the cost effectiveness in resource-limited clinical settings. Two commissioning datasets were acquired: (a) using a standard of practice 3D water tank scanning system (3DS) and (b) using a novel technique to translate a commercial TG-51 complaint 1D water tank via automated couch movements (1DS). The Extensible Markup Language (XML) was used to dynamically move the linear accelerator couch position (and thus the 1D tank) during radiation delivery for the acquisition of inline, crossline, and diagonal profiles. Both the 1DS and 3DS datasets were used to generate beam models (BM1 DS and BM3 DS ) in a commercial treatment planning system (TPS). 98.7% of 1DS measured points had a gamma value (2%/2 mm) < 1 when compared with the 3DS. Static jaw defined field and dynamic MLC field dose distribution comparisons for the TPS beam models BM1 DS and BM3 DS had 3D gamma values (2%/2 mm) < 1 for all 24,900,000 data points tested and >99.5% pass rate with gamma value (1%/1 mm) < 1. In conclusion, automated couch motions and a 1D scanning tank were used to collect commissioning beam data with accuracy comparable to traditionally acquired data using a 3D scanning system. TPS beam models generated directly from 1DS measured data were clinically equivalent to a model derived from 3DS data.

Normalize the response of EPID in pursuit of linear accelerator dosimetry standardization.

Normalize the response of electronic portal imaging device (EPID) is the first step toward an EPID-based standardization of Linear Accelerator (linac) dosimetry quality assurance. In this study, we described an approach to generate two-dimensional (2D) pixel sensitivity maps (PSM) for EPIDs response normalization utilizing an alternative beam and dark-field (ABDF) image acquisition technique and large overlapping field irradiations. The automated image acquisition was performed by XML-controlled machine operation and the PSM was generated based on a recursive calculation algorithm for Varian linacs equipped with aS1000 and aS1200 imager panels. Cross-comparisons of normalized beam profiles and 1.5%/1.5 mm 1D Gamma analysis was adopted to quantify the improvement of beam profile matching before and after PSM corrections. PSMs were derived for both photon (6, 10, 15 MV) and electron (6, 20 MeV) beams via proposed method. The PSM-corrected images reproduced a horn-shaped profile for photon beams and a relative uniform profiles for electrons. For dosimetrically matched linacs equipped with aS1000 panels, PSM-corrected images showed increased 1D-Gamma passing rates for all energies, with an average 10.5% improvement for crossline and 37% for inline beam profiles. Similar improvements in the phantom study were observed with a maximum improvement of 32% for 15 MV and 22% for 20 MeV. The PSM value showed no significant change for all energies over a 3-month period. In conclusion, the proposed approach correct EPID response for both aS1000 and aS1200 panels. This strategy enables the possibility to standardize linac dosimetry QA and to benchmark linac performance utilizing EPID as the common detector.

A method to reconstruct and apply 3D primary fluence for treatment delivery verification.

MOTIVATION: In this study, a method is reported to perform IMRT and VMAT treatment delivery verification using 3D volumetric primary beam fluences reconstructed directly from planned beam parameters and treatment delivery records. The goals of this paper are to demonstrate that 1) 3D beam fluences can be reconstructed efficiently, 2) quality assurance (QA) based on the reconstructed 3D fluences is capable of detecting additional treatment delivery errors, particularly for VMAT plans, beyond those identifiable by other existing treatment delivery verification methods, and 3) QA results based on 3D fluence calculation (3DFC) are correlated with QA results based on physical phantom measurements and radiation dose recalculations. METHODS: Using beam parameters extracted from DICOM plan files and treatment delivery log files, 3D volumetric primary fluences are reconstructed by forward-projecting the beam apertures, defined by the MLC leaf positions and modulated by beam MU values, at all gantry angles using first-order ray tracing. Treatment delivery verifications are performed by comparing 3D fluences reconstructed using beam parameters in delivery log files against those reconstructed from treatment plans. Passing rates are then determined using both voxel intensity differences and a 3D gamma analysis. QA sensitivity to various sources of errors is defined as the observed differences in passing rates. Correlations between passing rates obtained from QA derived from both 3D fluence calculations and physical measurements are investigated prospectively using 20 clinical treatment plans with artificially introduced machine delivery errors. RESULTS: Studies with artificially introduced errors show that common treatment delivery problems including gantry angle errors, MU errors, jaw position errors, collimator rotation errors, and MLC leaf position errors were detectable at less than normal machine tolerances. The reported 3DFC QA method has greater sensitivity than measurement-based QA methods. Statistical analysis-based Spearman's correlations shows that the 3DFC QA passing rates are significantly correlated with passing rates of physical phantom measurement-based QA methods. CONCLUSION: Among measurement-less treatment delivery verification methods, the reported 3DFC method is less demanding than those based on full dose re-calculations, and more comprehensive than those that solely checks beam parameters in treatment log files. With QA passing rates correlating to measurement-based passing rates, the 3DFC QA results could be useful for complementing the physical phantom measurements, or verifying treatment deliveries when physical measurements are not available. For the past 4+ years, the reported method has been implemented at authors' institution 1) as a complementary metric to physical phantom measurements for pretreatment, patient-specific QA of IMRT and VMAT plans, and 2) as an important part of the log file-based automated verification of daily patient treatment deliveries. It has been demonstrated to be useful in catching both treatment plan data transfer errors and treatment delivery problems.

Practical considerations for noise power spectra estimation for clinical CT scanners.

Local noise power spectra (NPS) have been commonly calculated to represent the noise properties of CT imaging systems, but their properties are significantly affected by the utilized calculation schemes. In this study, the effects of varied calculation parameters on the local NPS were analyzed, and practical suggestions were provided regarding the estimation of local NPS for clinical CT scanners. The uniformity module of a Catphan phantom was scanned with a Philips Brilliance 64 slice CT simulator with varied scanning protocols. Images were reconstructed using FBP and iDose4 iterative reconstruction with noise reduction levels 1, 3, and 6. Local NPS were calculated and compared for varied region of interest (ROI) locations and sizes, image background removal methods, and window functions. Additionally, with a predetermined NPS as a ground truth, local NPS calculation accuracy was compared for computer simulated ROIs, varying the aforementioned parameters in addition to ROI number. An analysis of the effects of these varied calculation parameters on the magnitude and shape of the NPS was conducted. The local NPS varied depending on calculation parameters, particularly at low spatial frequencies below ~ 0.15 mm-1. For the simulation study, NPS calculation error decreased exponentially as ROI number increased. For the Catphan study the NPS magnitude varied as a function of ROI location, which was better observed when using smaller ROI sizes. The image subtraction method for background removal was the most effective at reducing low-frequency background noise, and produced similar results no matter which ROI size or window function was used. The PCA background removal method with a Hann window function produced the closest match to image subtraction, with an average percent difference of 17.5%. Image noise should be analyzed locally by calculating the NPS for small ROI sizes. A minimum ROI size is recommended based on the chosen radial bin size and image pixel dimensions. As the ROI size decreases, the NPS becomes more dependent on the choice of background removal method and window function. The image subtraction method is most accurate, but other methods can achieve similar accuracy if certain window functions are applied. All dependencies should be analyzed and taken into account when considering the interpretation of the NPS for task-based image quality assessment.

A comparative study based on image quality and clinical task performance for CT reconstruction algorithms in radiotherapy.

CT image reconstruction is typically evaluated based on the ability to reduce the radiation dose to as-low-as-reasonably-achievable (ALARA) while maintaining acceptable image quality. However, the determination of common image quality metrics, such as noise, contrast, and contrast-to-noise ratio, is often insufficient for describing clinical radiotherapy task performance. In this study we designed and implemented a new comparative analysis method associating image quality, radiation dose, and patient size with radiotherapy task performance, with the purpose of guiding the clinical radiotherapy usage of CT reconstruction algorithms. The iDose4 iterative reconstruction algorithm was selected as the target for comparison, wherein filtered back-projection (FBP) reconstruction was regarded as the baseline. Both phantom and patient images were analyzed. A layer-adjustable anthropomorphic pelvis phantom capable of mimicking 38-58 cm lateral diameter-sized patients was imaged and reconstructed by the FBP and iDose4 algorithms with varying noise-reduction-levels, respectively. The resulting image sets were quantitatively assessed by two image quality indices, noise and contrast-to-noise ratio, and two clinical task-based indices, target CT Hounsfield number (for electron density determination) and structure contouring accuracy (for dose-volume calculations). Additionally, CT images of 34 patients reconstructed with iDose4 with six noise reduction levels were qualitatively evaluated by two radiation oncologists using a five-point scoring mechanism. For the phantom experiments, iDose4 achieved noise reduction up to 66.1% and CNR improvement up to 53.2%, compared to FBP without considering the changes of spatial resolution among images and the clinical acceptance of reconstructed images. Such improvements consistently appeared across different iDose4 noise reduction levels, exhibiting limited interlevel noise (< 5 HU) and target CT number variations (< 1 HU). The radiation dose required to achieve similar contouring accuracy decreased when using iDose4 in place of FBP, up to 32%. Contouring accuracy improvement for iDose4 images, when compared to FBP, was greater in larger patients than smaller-sized patients. Overall, the iDose4 algorithm provided superior radiation dose control while maintaining or improving task performance, when compared to FBP. The reader study on image quality improvement of patient cases shows that physicians preferred iDose4-reconstructed images on all cases compared to those from FBP algorithm with overall quality score: 1.21 vs. 3.15, p = 0.0022. However, qualitative evaluation strongly indicated that the radiation oncologists chose iDose4 noise reduction levels of 3-4 with additional consideration of task performance, instead of image quality metrics alone. Although higher iDose4 noise reduction levels improved the CNR through the further reduction of noise, there was pixelization of anatomical/tumor structures. Very-low-dose scans yielded severe photon starvation artifacts, which decreased target visualization on both FBP and iDose4 reconstructions, especially for the 58 cm phantom size. The iDose4 algorithm with a moderate noise reduction level is hence suggested for CT simulation and treatment planning. Quantitative task-based image quality metrics should be further investigated to accommodate additional clinical applications.

A comparative study of automatic image segmentation algorithms for target tracking in MR-IGRT.

On-board magnetic resonance (MR) image guidance during radiation therapy offers the potential for more accurate treatment delivery. To utilize the real-time image information, a crucial prerequisite is the ability to successfully segment and track regions of interest (ROI). The purpose of this work is to evaluate the performance of different segmentation algorithms using motion images (4 frames per second) acquired using a MR image-guided radiotherapy (MR-IGRT) system. Manual con-tours of the kidney, bladder, duodenum, and a liver tumor by an experienced radiation oncologist were used as the ground truth for performance evaluation. Besides the manual segmentation, images were automatically segmented using thresholding, fuzzy k-means (FKM), k-harmonic means (KHM), and reaction-diffusion level set evolution (RD-LSE) algorithms, as well as the tissue tracking algorithm provided by the ViewRay treatment planning and delivery system (VR-TPDS). The performance of the five algorithms was evaluated quantitatively by comparing with the manual segmentation using the Dice coefficient and target registration error (TRE) measured as the distance between the centroid of the manual ROI and the centroid of the automatically segmented ROI. All methods were able to successfully segment the bladder and the kidney, but only FKM, KHM, and VR-TPDS were able to segment the liver tumor and the duodenum. The performance of the thresholding, FKM, KHM, and RD-LSE algorithms degraded as the local image contrast decreased, whereas the performance of the VP-TPDS method was nearly independent of local image contrast due to the reference registration algorithm. For segmenting high-contrast images (i.e., kidney), the thresholding method provided the best speed (< 1 ms) with a satisfying accuracy (Dice = 0.95). When the image contrast was low, the VR-TPDS method had the best automatic contour. Results suggest an image quality determination procedure before segmentation and a combination of different methods for optimal segmentation with the on-board MR-IGRT system.

A software tool to automatically assure and report daily treatment deliveries by a cobalt-60 radiation therapy device.

The aims of this study were to develop a method for automatic and immediate verification of treatment delivery after each treatment fraction in order to detect and correct errors, and to develop a comprehensive daily report which includes delivery verification results, daily image-guided radiation therapy (IGRT) review, and information for weekly physics reviews. After systematically analyzing the requirements for treatment delivery verification and understanding the available information from a commercial MRI-guided radiotherapy treatment machine, we designed a procedure to use 1) treatment plan files, 2) delivery log files, and 3) beam output information to verify the accuracy and completeness of each daily treatment delivery. The procedure verifies the correctness of delivered treatment plan parameters including beams, beam segments and, for each segment, the beam-on time and MLC leaf positions. For each beam, composite primary fluence maps are calculated from the MLC leaf positions and segment beam-on time. Error statistics are calculated on the fluence difference maps between the plan and the delivery. A daily treatment delivery report is designed to include all required information for IGRT and weekly physics reviews including the plan and treatment fraction information, daily beam output information, and the treatment delivery verification results. A computer program was developed to implement the proposed procedure of the automatic delivery verification and daily report generation for an MRI guided radiation therapy system. The program was clinically commissioned. Sensitivity was measured with simulated errors. The final version has been integrated into the com-mercial version of the treatment delivery system. The method automatically verifies the EBRT treatment deliveries and generates the daily treatment reports. Already in clinical use for over one year, it is useful to facilitate delivery error detection, and to expedite physician daily IGRT review and physicist weekly chart review.