Assessment of a novel commercial large field of view phantom for comprehensive MR imaging quality assurance of a 0.35T MRgRT system.

Consistent quality assurance (QA) programs are vital to MR-guided radiotherapy (MRgRT), for ensuring treatment is delivered accurately and the onboard MRI system is providing the expected image quality. However, daily imaging QA with a dedicated phantom is not common at many MRgRT centers, especially with large phantoms that cover a field of view (FOV), similar to the human torso. This work presents the first clinical experience with a purpose-built phantom for large FOV daily and periodic comprehensive quality assurance (QUASAR™ MRgRT Insight Phantom (beta)) from Modus Medical Devices Inc. (Modus QA) on an MRgRT system. A monthly American College of Radiology (ACR) QA phantom was also imaged for reference. Both phantoms were imaged on a 0.35T MR-Linac, a 1.5T Philips wide bore MRI, and a 3.0T Siemens MRI, with T1-weighted and T2-weighted acquisitions. The Insight phantom was imaged in axial and sagittal orientations. Image quality tests including geometric accuracy, spatial resolution accuracy, slice thickness accuracy, slice position accuracy, and image intensity uniformity were performed on each phantom, following their respective instruction manuals. The geometric distortion test showed similar distortions of -1.7 mm and -1.9 mm across a 190 mm and a 283 mm lengths for the ACR and MRgRT Insight phantoms, respectively. The MRgRT Insight phantom utilized a modulation transform function (MTF) for spatial resolution evaluation, which showed decreased performance on the lower B0 strength MRIs, as expected, and could provide a good daily indicator of machine performance. Both the Insight and ACR phantoms showed a match with scan parameters for slice thickness analysis. During the imaging and analysis of this novel MRgRT Insight phantom the authors found setup to be straightforward allowing for easy acquisition each day, and useful image analysis parameters for tracking MRI performance.

Optical imaging method to quantify spatial dose variation due to the electron return effect in an MR-linac.

PURPOSE: Treatment planning systems (TPSs) for MR-linacs must employ Monte Carlo-based simulations of dose deposition to model the effects of the primary magnetic field on dose. However, the accuracy of these simulations, especially for areas of tissue-air interfaces where the electron return effect (ERE) is expected, is difficult to validate due to physical constraints and magnetic field compatibility of available detectors. This study employs a novel dosimetric method based on remotely captured, real-time optical Cherenkov and scintillation imaging to visualize and quantify the ERE. METHODS: An intensified CMOS camera was used to image two phantoms with designed ERE cavities. Phantom A was a 40 cm × 10 cm × 10 cm clear acrylic block drilled with five holes of increasing diameters (0.5, 1, 2, 3, 4 cm). Phantom B was a clear acrylic block (25 cm × 20 cm × 5 cm) with three cavities of increasing diameter (3, 2, 1 cm) split into two halves in the transverse plane to accommodate radiochromic film. Both phantoms were imaged while being irradiated by 6 MV flattening filter free (FFF) beams within a MRIdian Viewray (Viewray, Cleveland, OH) MR-linac (0.34 T primary field). Phantom A was imaged while being irradiated by 6 MV FFF beams on a conventional linac (TrueBeam, Varian Medical Systems, San Jose, CA) to serve as a control. Images were post processed in Matlab (Mathworks Inc., Natick, MA) and compared to TPS dose volumes. RESULTS: Control imaging of Phantom A without the presence of a magnetic field supports the validity of the optical image data to a depth of 6 cm. In the presence of the magnetic field, the optical data shows deviations from the commissioned TPS dose in both intensity and localization. The largest air cavity examined (3 cm) indicated the largest dose differences, which were above 20% at some locations. Experiments with Phantom B illustrated similar agreement between optical and film dosimetry comparisons with TPS data in areas not affected by ERE. CONCLUSION: There are some appreciable differences in dose intensity and spatial dose distribution observed between the novel experimental data set and the dose models produced by the current clinically implemented MR-IGRT TPS.

Remote Cherenkov imaging-based quality assurance of a magnetic resonance image-guided radiotherapy system.

PURPOSE: Tools to perform regular quality assurance of magnetic resonance image-guided radiotherapy (MRIgRT) systems should ideally be independent of interference from the magnetic fields. Remotely acquired optical Cherenkov imaging-based dosimetry measurements in water were investigated for this purpose, comparing measures of dose accuracy, temporal dynamics, and overall integrated IMRT delivery. METHODS: A 40 × 30.5 × 37.5 cm3 water tank doped with 1 g/L of quinine sulfate was imaged using an intensified charge-coupled device (ICCD) to capture the Cherenkov emission while being irradiated by a commercial MRIgRT system (ViewRay™). The ICCD was placed down-bore at the end of the couch, 4 m from treatment isocenter and behind the 5-Gauss line of the 0.35-T MRI. After establishing optimal camera acquisition settings, square beams of increasing size (4.2 × 4.2 cm2 , 10.5 × 10.5 cm2 , and 14.7 × 14.7 cm2 ) were imaged at 0.93 frames per second, from an individual cobalt-60 treatment head, to develop projection measures related to percent depth dose (PDD) curves and cross beam profiles (CPB). These Cherenkov-derived measurements were compared to ionization chamber (IC) and radiographic film dosimetry data, as well as simulation data from the treatment planning system (TPS). An intensity-modulated radiotherapy (IMRT) commissioning plan from AAPM TG-119 (C4:C-Shape) was also imaged at 2.1 frames per second, and the single linear sum image from 509 s of plan delivery was compared to the dose volume prediction generated by the TPS using gamma index analysis. RESULTS: Analysis of standardized test target images (1024 × 1024 pixels) yielded a pixel resolution of 0.37 mm/pixel. The beam width measured from the Cherenkov image-generated projection CBPs was within 1 mm accuracy when compared to film measurements for all beams. The 502 point measurements (i.e., pixels) of the Cherenkov image-based projection percent depth dose curves (pPDDs) were compared to pPDDs simulated by the treatment planning system (TPS), with an overall average error of 0.60%, 0.56%, and 0.65% for the 4.2, 10.5, and 14.7 cm square beams, respectively. The relationships between pPDDs and central axis PDDs derived from the TPS were used to apply a weighting factor to the Cherenkov pPDD, so that the Cherenkov data could be directly compared to IC PDDs (average error of -0.07%, 0.10%, and -0.01% for the same sized beams, respectively). Finally, the composite image of the TG-119 C4 treatment plan achieved a 95.1% passing rate using 4%/4 mm gamma index agreement criteria between Cherenkov intensity and TPS dose volume data. CONCLUSIONS: This is the first examination of Cherenkov-generated pPDDs and pCBPs in an MR-IGRT system. Cherenkov imaging measurements were fast to acquire, and minimal error was observed overall. Cherenkov imaging also provided novel real-time data for IMRT QA. The strengths of this imaging are the rapid data capture ability providing real-time, high spatial resolution data, combined with the remote, noncontact nature of imaging. The biggest limitation of this method is the two-dimensional (2D) projection-based imaging of three-dimensional (3D) dose distributions through the transparent water tank.

Performance of a multi leaf collimator system for MR-guided radiation therapy.

PURPOSE: The purpose of this study was to investigate and characterize the performance of a Multi Leaf Collimator (MLC) designed for Cobalt-60 based MR-guided radiation therapy system in a 0.35 T magnetic field. METHODS: The MLC design and unique assembly features in the ViewRay MRIdian system were first reviewed. The RF cage shielding of MLC motor and cables were evaluated using ACR phantoms with real-time imaging and quantified by signal-to-noise ratio. The dosimetric characterizations, including the leaf transmission, leaf penumbra, tongue-and-groove effect, were investigated using radiosensitive films. The output factor of MLC-defined fields was measured with ionization chambers for both symmetric fields from 2.1 × 2.1 cm2 to 27.3 × 27.3 cm2 and asymmetric fields from 10.5 × 10.5 cm2 to 10.5 × 2.0 cm2 . Multi leaf collimator (MLC) positional accuracy was assessed by delivering either a picket fence (PF) style pattern on radiochromic films with wire-jig phantom or double and triple-rectangular patterns on ArcCheck-MR (Sun Nuclear, Melbourne, FL, USA) with gamma analysis as the pass/fail indicator. Leaf speed tests were performed to assess the capability of full range leaf travel within manufacture's specifications. Multi leaf collimator plan delivery reproducibility was tested by repeatedly delivering both open fields and fields with irregular shaped segments over 1-month period. RESULTS: Comparable SNRs within 4% were observed for MLC moving and stationary plans on vendor-reconstructed images, and the direct k-space reconstructed images showed that the three SNRs are within 1%. The maximum leaf transmission for all three MLCs was less than 0.35% and the average leakage was 0.153 ± 0.006%, 0.151 ± 0.008%, and 0.159 ± 0.015% for head 1, 2, and 3, respectively. Both the leaf edge and leaf end penumbra showed comparable values within 0.05 cm, and the measured values are within 0.1 cm with TPS values. The leaf edge TG effect indicated 10% underdose and the leaf end TG showed a shifted dose distribution with 0.3 cm offset. The leaf positioning test showed a 0.2 cm accuracy in the PF style test, and a gamma passing rate above 96% was observed with a 3%/2 mm criteria when comparing the measured double/triple-rectangular pattern fluence with TPS calculated fluence. The average leaf speed when executing the test plan fell in a range from 1.86 to 1.95 cm/s. The measured and TPS calculated output factors were within 2% for squared fields and within 3% for rectangular fields. The reproducibility test showed the deviation of output factors were well within 2% for square fields and the gamma passing rate within 1.5% for fields with irregular segments. The Monte Carlo predicted output factors were within 2% compared to TPS values. 15 out of the 16 IMRT plans have gamma passing rate more than 98% compared to the TPS fluence with an average passing rate of 99.1 ± 0.6%. CONCLUSION: The MRIdian MLC has a good RF noise shielding design, low radiation leakage, good positioning accuracy, comparable TG effect, and can be modeled by an independent Monte Carlo calculation platform.

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.

A GPU-accelerated Monte Carlo dose calculation platform and its application toward validating an MRI-guided radiation therapy beam model.

PURPOSE: The clinical commissioning of IMRT subject to a magnetic field is challenging. The purpose of this work is to develop a GPU-accelerated Monte Carlo dose calculation platform based on penelope and then use the platform to validate a vendor-provided MRIdian head model toward quality assurance of clinical IMRT treatment plans subject to a 0.35 T magnetic field. METHODS: penelope was first translated from fortran to c++ and the result was confirmed to produce equivalent results to the original code. The c++ code was then adapted to cuda in a workflow optimized for GPU architecture. The original code was expanded to include voxelized transport with Woodcock tracking, faster electron/positron propagation in a magnetic field, and several features that make gpenelope highly user-friendly. Moreover, the vendor-provided MRIdian head model was incorporated into the code in an effort to apply gpenelope as both an accurate and rapid dose validation system. A set of experimental measurements were performed on the MRIdian system to examine the accuracy of both the head model and gpenelope. Ultimately, gpenelope was applied toward independent validation of patient doses calculated by MRIdian's kmc. RESULTS: An acceleration factor of 152 was achieved in comparison to the original single-thread fortran implementation with the original accuracy being preserved. For 16 treatment plans including stomach (4), lung (2), liver (3), adrenal gland (2), pancreas (2), spleen(1), mediastinum (1), and breast (1), the MRIdian dose calculation engine agrees with gpenelope with a mean gamma passing rate of 99.1% ± 0.6% (2%/2 mm). CONCLUSIONS: A Monte Carlo simulation platform was developed based on a GPU- accelerated version of penelope. This platform was used to validate that both the vendor-provided head model and fast Monte Carlo engine used by the MRIdian system are accurate in modeling radiation transport in a patient using 2%/2 mm gamma criteria. Future applications of this platform will include dose validation and accumulation, IMRT optimization, and dosimetry system modeling for next generation MR-IGRT systems.

Online Magnetic Resonance Image Guided Adaptive Radiation Therapy: First Clinical Applications.

PURPOSE: To demonstrate the feasibility of online adaptive magnetic resonance (MR) image guided radiation therapy (MR-IGRT) through reporting of our initial clinical experience and workflow considerations. METHODS AND MATERIALS: The first clinically deployed online adaptive MR-IGRT system consisted of a split 0.35T MR scanner straddling a ring gantry with 3 multileaf collimator-equipped (60)Co heads. The unit is supported by a Monte Carlo-based treatment planning system that allows real-time adaptive planning with the patient on the table. All patients undergo computed tomography and MR imaging (MRI) simulation for initial treatment planning. A volumetric MRI scan is acquired for each patient at the daily treatment setup. Deformable registration is performed using the planning computed tomography data set, which allows for the transfer of the initial contours and the electron density map to the daily MRI scan. The deformed electron density map is then used to recalculate the original plan on the daily MRI scan for physician evaluation. Recontouring and plan reoptimization are performed when required, and patient-specific quality assurance (QA) is performed using an independent in-house software system. RESULTS: The first online adaptive MR-IGRT treatments consisted of 5 patients with abdominopelvic malignancies. The clinical setting included neoadjuvant colorectal (n=3), unresectable gastric (n=1), and unresectable pheochromocytoma (n=1). Recontouring and reoptimization were deemed necessary for 3 of 5 patients, and the initial plan was deemed sufficient for 2 of the 5 patients. The reasons for plan adaptation included tumor progression or regression and a change in small bowel anatomy. In a subsequently expanded cohort of 170 fractions (20 patients), 52 fractions (30.6%) were reoptimized online, and 92 fractions (54.1%) were treated with an online-adapted or previously adapted plan. The median time for recontouring, reoptimization, and QA was 26 minutes. CONCLUSION: Online adaptive MR-IGRT has been successfully implemented with planning and QA workflow suitable for routine clinical application. Clinical trials are in development to formally evaluate adaptive treatments for a variety of disease sites.

Quality of Intensity Modulated Radiation Therapy Treatment Plans Using a 6°Co Magnetic Resonance Image Guidance Radiation Therapy System.

PURPOSE: This work describes a commercial treatment planning system, its technical features, and its capabilities for creating (60)Co intensity modulated radiation therapy (IMRT) treatment plans for a magnetic resonance image guidance radiation therapy (MR-IGRT) system. METHODS AND MATERIALS: The ViewRay treatment planning system (Oakwood Village, OH) was used to create (60)Co IMRT treatment plans for 33 cancer patients with disease in the abdominal, pelvic, thorax, and head and neck regions using physician-specified patient-specific target coverage and organ at risk (OAR) objectives. Backup plans using a third-party linear accelerator (linac)-based planning system were also created. Plans were evaluated by attending physicians and approved for treatment. The (60)Co and linac plans were compared by evaluating conformity numbers (CN) with 100% and 95% of prescription reference doses and heterogeneity indices (HI) for planning target volumes (PTVs) and maximum, mean, and dose-volume histogram (DVH) values for OARs. RESULTS: All (60)Co IMRT plans achieved PTV coverage and OAR sparing that were similar to linac plans. PTV conformity for (60)Co was within <1% and 3% of linac plans for 100% and 95% prescription reference isodoses, respectively, and heterogeneity was on average 4% greater. Comparisons of OAR mean dose showed generally better sparing with linac plans in the low-dose range <20 Gy, but comparable sparing for organs with mean doses >20 Gy. The mean doses for all (60)Co plan OARs were within clinical tolerances. CONCLUSIONS: A commercial (60)Co MR-IGRT device can produce highly conformal IMRT treatment plans similar in quality to linac IMRT for a variety of disease sites. Additional work is in progress to evaluate the clinical benefit of other novel features of this MR-IGRT system.

Patient-specific quality assurance for the delivery of (60)Co intensity modulated radiation therapy subject to a 0.35-T lateral magnetic field.

PURPOSE: This work describes a patient-specific dosimetry quality assurance (QA) program for intensity modulated radiation therapy (IMRT) using ViewRay, the first commercial magnetic resonance imaging-guided RT device. METHODS AND MATERIALS: The program consisted of: (1) a 1-dimensional multipoint ionization chamber measurement using a customized 15-cm(3) cube-shaped phantom; (2) 2-dimensional (2D) radiographic film measurement using a 30- × 30- × 20-cm(3) phantom with multiple inserted ionization chambers; (3) quasi-3D diode array (ArcCHECK) measurement with a centrally inserted ionization chamber; (4) 2D fluence verification using machine delivery log files; and (5) 3D Monte Carlo (MC) dose reconstruction with machine delivery files and phantom CT. RESULTS: Ionization chamber measurements agreed well with treatment planning system (TPS)-computed doses in all phantom geometries where the mean ± SD difference was 0.0% ± 1.3% (n=102; range, -3.0%-2.9%). Film measurements also showed excellent agreement with the TPS-computed 2D dose distributions where the mean passing rate using 3% relative/3 mm gamma criteria was 94.6% ± 3.4% (n=30; range, 87.4%-100%). For ArcCHECK measurements, the mean ± SD passing rate using 3% relative/3 mm gamma criteria was 98.9% ± 1.1% (n=34; range, 95.8%-100%). 2D fluence maps with a resolution of 1 × 1 mm(2) showed 100% passing rates for all plan deliveries (n=34). The MC reconstructed doses to the phantom agreed well with planned 3D doses where the mean passing rate using 3% absolute/3 mm gamma criteria was 99.0% ± 1.0% (n=18; range, 97.0%-100%), demonstrating the feasibility of evaluating the QA results in the patient geometry. CONCLUSIONS: We developed a dosimetry program for ViewRay's patient-specific IMRT QA. The methodology will be useful for other ViewRay users. The QA results presented here can assist the RT community to establish appropriate tolerance and action limits for ViewRay's IMRT QA.