O.L. Green

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 60Co 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.

Benchmark IMRT evaluation of a Co-60 MRI-guided radiation therapy system

A device for MRI-guided radiation therapy (MR-IGRT) that uses cobalt-60 sources to deliver intensity modulated radiation therapy is now commercially available. We investigated the performance of the treatment planning and delivery system against the benchmark recommended by the American Association of Physicists in Medicine (AAPM) Task Group 119 for IMRT commissioning and demonstrated that the device plans and delivers IMRT treatments within recommended confidence limits and with similar accuracy as linac IMRT.

Comparison of onboard low-field magnetic resonance imaging versus onboard computed tomography for anatomy visualization in radiotherapy.

ABSTRACT Background. Onboard magnetic resonance imaging (OB-MRI) for daily localization and adaptive radiotherapy has been under development by several groups. However, no clinical studies have evaluated whether OB-MRI improves visualization of the target and organs at risk (OARs) compared to standard onboard computed tomography (OB-CT). This study compared visualization of patient anatomy on images acquired on the MRI-60Co ViewRay system to those acquired with OB-CT. Material and methods. Fourteen patients enrolled on a protocol approved by the Institutional Review Board (IRB) and undergoing image-guided radiotherapy for cancer in the thorax (n = 2), pelvis (n = 6), abdomen (n = 3) or head and neck (n = 3) were imaged with OB-MRI and OB-CT. For each of the 14 patients, the OB-MRI and OB-CT datasets were displayed side-by-side and independently reviewed by three radiation oncologists. Each physician was asked to evaluate which dataset offered better visualization of the target and OARs. A quantitative contouring study was performed on two abdominal patients to assess if OB-MRI could offer improved inter-observer segmentation agreement for adaptive planning. Results. In total 221 OARs and 10 targets were compared for visualization on OB-MRI and OB-CT by each of the three physicians. The majority of physicians (two or more) evaluated visualization on MRI as better for 71% of structures, worse for 10% of structures, and equivalent for 14% of structures. 5% of structures were not visible on either. Physicians agreed unanimously for 74% and in majority for > 99% of structures. Targets were better visualized on MRI in 4/10 cases, and never on OB-CT. Conclusion. Low-field MR provides better anatomic visualization of many radiotherapy targets and most OARs as compared to OB-CT. Further studies with OB-MRI should be pursued.

Patient-Specific Quality Assurance for the Delivery of 60Co 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 ViewRays 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 ViewRays IMRT QA.

Characterization of the onboard imaging unit for the first clinical magnetic resonance image guided radiation therapy system.

PURPOSE To characterize the performance of the onboard imaging unit for the first clinical magnetic resonance image guided radiation therapy (MR-IGRT) system. METHODS The imaging performance characterization included four components: ACR (the American College of Radiology) phantom test, spatial integrity, coil signal to noise ratio (SNR) and uniformity, and magnetic field homogeneity. The ACR phantom test was performed in accordance with the ACR phantom test guidance. The spatial integrity test was evaluated using a 40.8 × 40.8 × 40.8 cm(3) spatial integrity phantom. MR and computed tomography (CT) images of the phantom were acquired and coregistered. Objects were identified around the surfaces of 20 and 35 cm diameters of spherical volume (DSVs) on both the MR and CT images. Geometric distortion was quantified using deviation in object location between the MR and CT images. The coil SNR test was performed according to the national electrical manufacturers association (NEMA) standards MS-1 and MS-9. The magnetic field homogeneity test was measured using field camera and spectral peak methods. RESULTS For the ACR tests, the slice position error was less than 0.10 cm, the slice thickness error was less than 0.05 cm, the resolved high-contrast spatial resolution was 0.09 cm, the resolved low-contrast spokes were more than 25, the image intensity uniformity was above 93%, and the percentage ghosting was less than 0.22%. All were within the ACR recommended specifications. The maximum geometric distortions within the 20 and 35 cm DSVs were 0.10 and 0.18 cm for high spatial resolution three-dimensional images and 0.08 and 0.20 cm for high temporal resolution two dimensional cine images based on the distance-to-phantom-center method. The average SNR was 12.0 for the body coil, 42.9 for the combined torso coil, and 44.0 for the combined head and neck coil. Magnetic field homogeneities at gantry angles of 0°, 30°, 60°, 90°, and 120° were 23.55, 20.43, 18.76, 19.11, and 22.22 ppm, respectively, using the field camera method over the 45 cm DSV. CONCLUSIONS The onboard imaging unit of the first commercial MR-IGRT system meets ACR, NEMA, and vendor specifications.

An integrated model-driven method for in-treatment upper airway motion tracking using cine MRI in head and neck radiation therapy.

PURPOSE For the first time, MRI-guided radiation therapy systems can acquire cine images to dynamically monitor in-treatment internal organ motion. However, the complex head and neck (H&N) structures and low-contrast/resolution of on-board cine MRI images make automatic motion tracking a very challenging task. In this study, the authors proposed an integrated model-driven method to automatically track the in-treatment motion of the H&N upper airway, a complex and highly deformable region wherein internal motion often occurs in an either voluntary or involuntary manner, from cine MRI images for the analysis of H&N motion patterns. METHODS Considering the complex H&N structures and ensuring automatic and robust upper airway motion tracking, the authors firstly built a set of linked statistical shapes (including face, face-jaw, and face-jaw-palate) using principal component analysis from clinically approved contours delineated on a set of training data. The linked statistical shapes integrate explicit landmarks and implicit shape representation. Then, a hierarchical model-fitting algorithm was developed to align the linked shapes on the first image frame of a to-be-tracked cine sequence and to localize the upper airway region. Finally, a multifeature level set contour propagation scheme was performed to identify the upper airway shape change, frame-by-frame, on the entire image sequence. The multifeature fitting energy, including the information of intensity variations, edge saliency, curve geometry, and temporal shape continuity, was minimized to capture the details of moving airway boundaries. Sagittal cine MR image sequences acquired from three H&N cancer patients were utilized to demonstrate the performance of the proposed motion tracking method. RESULTS The tracking accuracy was validated by comparing the results to the average of two manual delineations in 50 randomly selected cine image frames from each patient. The resulting average dice similarity coefficient (93.28%  ±  1.46%) and margin error (0.49  ±  0.12 mm) showed good agreement between the automatic and manual results. The comparison with three other deformable model-based segmentation methods illustrated the superior shape tracking performance of the proposed method. Large interpatient variations of swallowing frequency, swallowing duration, and upper airway cross-sectional area were observed from the testing cine image sequences. CONCLUSIONS The proposed motion tracking method can provide accurate upper airway motion tracking results, and enable automatic and quantitative identification and analysis of in-treatment H&N upper airway motion. By integrating explicit and implicit linked-shape representations within a hierarchical model-fitting process, the proposed tracking method can process complex H&N structures and low-contrast/resolution cine MRI images. Future research will focus on the improvement of method reliability, patient motion pattern analysis for providing more information on patient-specific prediction of structure displacements, and motion effects on dosimetry for better H&N motion management in radiation therapy.

Technical Note: Magnetic field effects on Gafchromic‐film response in MR‐IGRT

PURPOSE Magnetokinetic changes may affect crystal orientation and polymerization within the active layer of radiochromic film (RCF). This effect is investigated in a magnetic resonance image-guided radiotherapy unit within the context of film dosimetry. METHODS Gafchromic EBT2 RCF was irradiated in a 30 × 30 × 30 cm3 solid water phantom using a Co-60 MRI guided radiotherapy system (B = 0.35 T) under normal operating conditions, and under the exact conditions and setup without a magnetic field. Fifteen 20.3 × 25.4 cm2 EBT2 film sheets were placed at three different depths (d = 0.5, 5, and 10 cm) using five different treatment plans. The plans were computed using the MRIdian (ViewRay, Inc.) treatment planning system to deliver doses between 0 and 17.6 Gy. Films were analyzed before and after irradiation to obtain the net optical density (netOD) for each color channel separately. Scanning electron microscope (SEM) images were obtained to compare the active layer of selected samples. RESULTS The results indicated that the red channel netOD decreased between 0.013 and 0.123 (average of 0.060 ± 0.033) for doses above 2.8 Gy, with a linear increase in this effect for higher doses. Green channel netOD showed similar results with a decrease between 0.012 and 0.105 (average of 0.041 ± 0.027) for doses above 3.5 Gy. The blue channel showed the weakest effect with a netOD decrease between 0.013 and 0.029 (average of 0.020 ± 0.006) for doses above 8.0 Gy. SEM images show changes in crystal orientation within active layer in RCF exposed in a magnetic field. CONCLUSIONS The presence of a magnetic field affects crystal orientation and polymerization during irradiation, where netOD decreased by an average of 8.7%, 8.0%, and 4.3% in the red, green, and blue channels, respectively. The under response was dependent on dose and differed by up to 15% at 17.6 Gy.

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 MRIdians 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.

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 commercial 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. PACS number(s): 87.55.kmThe 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 commercial 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. PACS number(s): 87.55.km.