A reconstruction approach for proton computed tomography by modeling the integral depth dose of the scanning proton pencil beam.

PURPOSE: To present a proton computed tomography (pCT) reconstruction approach that models the integral depth dose (IDD) of the clinical scanning proton beam into beamlets. Using a multilayer ionization chamber (MLIC) as the imager, the proposed pCT system and the reconstruction approach can minimize extra ambient neutron dose and simplify the beamline design by eliminating an additional collimator to confine the proton beam. METHODS: Monte Carlo simulation was applied to digitally simulate the IDDs of the exiting proton beams detected by the MLIC. A forward model was developed to model each IDD into a weighted sum of percentage depth doses of the constituent beamlets separated laterally by 1 mm. The water equivalent path lengths (WEPLs) of the beamlets were determined by iteratively minimizing the squared L2-norm between the forward projected and simulated IDDs. The final WEPL values were reconstructed to pCT images, that is, proton stopping power ratio (SPR) maps, through simultaneous algebraic reconstruction technique with total variation regularization. The reconstruction process was tested with a digital cylindrical water-based phantom and an ICRP adult reference computational phantom. The mean of SPR within regions of interest (ROIs) and the WEPL along a 4 mm-wide beam ( WEP L 4 mm ${\rm{WEP}}{{\rm{L}}_{4{\rm{mm}}}}$ ) were compared with the reference values. The spatial resolution was analyzed at the edge of a cortical insert of the cylindrical phantom. RESULTS: The percentage deviations from reference SPR were within ±1% in all selected ROIs. The mean absolute error of the reconstructed SPR was 0.33%, 0.19%, and 0.27% for the cylindrical phantom, the adult phantom at the head and lung region, respectively. The corresponding percentage deviations from reference WEP L 4 mm ${\rm{WEP}}{{\rm{L}}_{4{\rm{mm}}}}$ were 0.48 ± 0.64%, 0.28 ± 0.48%, and 0.22 ± 0.49%. The full width at half maximum of the line spread function (LSF) derived from the radial edge spread function (ESF) of a cortical insert was 0.13 cm. The frequency at 10% of the modulation transfer function (MTF) was 6.38 cm-1 . The mean signal-to-noise ratio (SNR) of all the inserts was 2.45. The mean imaging dose was 0.29 and 0.25 cGy at the head and lung region of the adult phantom, respectively. CONCLUSION: A new pCT reconstruction approach was developed by modeling the IDDs of the uncollimated scanning proton beams in the pencil beam geometry. SPR accuracy within ±1%, spatial resolution of better than 2 mm at 10% MTF, and imaging dose at the magnitude of mGy were achieved. Potential side effects caused by neutron dose were eliminated by removing the extra beam collimator.

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.

Development of a fast Monte Carlo dose calculation system for online adaptive radiation therapy quality assurance.

Online adaptive radiation therapy (ART) based on real-time magnetic resonance imaging represents a paradigm-changing treatment scheme. However, conventional quality assurance (QA) methods based on phantom measurements are not feasible with the patient on the treatment couch. The purpose of this work is to develop a fast Monte Carlo system for validating online re-optimized tri-60Co IMRT adaptive plans with both high accuracy and speed. The Monte Carlo system is based on dose planning method (DPM) code with further simplification of electron transport and consideration of external magnetic fields. A vendor-provided head model was incorporated into the code. Both GPU acceleration and variance reduction were implemented. Additionally, to facilitate real-time decision support, a C++ GUI was developed for visualizing 3D dose distributions and performing various analyses in an online adaptive setting. A thoroughly validated Monte Carlo code (gPENELOPE) was used to benchmark the new system, named GPU-accelerated DPM with variance reduction (gDPMvr). The comparison using 15 clinical IMRT plans demonstrated that gDPMvr typically runs 43 times faster with only 0.5% loss in accuracy. Moreover, gDPMvr reached 1% local dose uncertainty within 2.3 min on average, and thus is well-suited for ART QA.

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.

Two-dimensional high spatial-resolution dosimeter using europium doped potassium chloride: a feasibility study.

Recent research has shown that KCl:Eu²âº has great potential for use in megavoltage radiation therapy dosimetry because this material exhibits excellent storage performance and is reusable due to strong radiation hardness. This work reports the authors' attempts to fabricate 2D KCl:Eu²âº storage phosphor films (SPFs) using both a physical vapor deposition (PVD) method and a tape casting method. X-ray diffraction analysis showed that a 10 µm thick PVD sample was composed of highly crystalline KCl. No additional phases were observed, suggesting that the europium activator had been completely incorporated into the KCl matrix. Photostimulated luminescence and photoluminescence spectra suggested that F (Cl(-)) centers were the electron storage centers post x-ray irradiation and that Eu²âº cations acted as luminescence centers in the photostimulation process. The 150 µm thick casted KCl:Eu²âº SPF showed sub-millimeter spatial-resolution. Monte Carlo simulations further demonstrated that the admixture of 20% KCl:Eu²âº and 80% low Z polymer binder exhibited almost no energy-dependence in a 6 MV beam. KCl:Eu²âº pellet samples showed a large dynamic range from 0.01 cGy to 60 Gy dose-to-water, and saturated at approximately 500 Gy as a result of KCl's intrinsic high radiation hardness. Taken together, this work provides strong evidence that KCl:Eu²âº-based SPF with associated readout apparatus could result in a novel electronic film system that has all the desirable features associated with classic radiographic film and, importantly, water equivalence and the capability of permanent identification of each detector.

Image-based modeling and precision medicine: patient-specific carotid and coronary plaque assessment and predictions.

Atherosclerotic plaques may rupture without warning and cause acute cardiovascular events such as heart attack and stroke. Current clinical screening tools are insufficient to identify those patients with risks early and prevent the adverse events from happening. Medical imaging and image-based modeling have made considerable progress in recent years in identifying plaque morphological and mechanical risk factors which may be used in developing improved patient screening strategies. The key steps and factors in image-based models for human carotid and coronary plaques were illustrated, as well as grand challenges facing the researchers in the field to develop more accurate screening tools.

Automated radiotherapy treatment plan integrity verification.

PURPOSE: In our clinic, physicists spend from 15 to 60 min to verify the physical and dosimetric integrity of radiotherapy plans before presentation to radiation oncology physicians for approval. The purpose of this study was to design and implement a framework to automate as many elements of this quality control (QC) step as possible. METHODS: A comprehensive computer application was developed to carry out a majority of these verification tasks in the Philips PINNACLE treatment planning system (TPS). This QC tool functions based on both PINNACLE scripting elements and PERL sub-routines. The core of this technique is the method of dynamic scripting, which involves a PERL programming module that is flexible and powerful for treatment plan data handling. Run-time plan data are collected, saved into temporary files, and analyzed against standard values and predefined logical rules. The results were summarized in a hypertext markup language (HTML) report that is displayed to the user. RESULTS: This tool has been in clinical use for over a year. The occurrence frequency of technical problems, which would cause delays and suboptimal plans, has been reduced since clinical implementation. CONCLUSIONS: In addition to drastically reducing the set of human-driven logical comparisons, this QC tool also accomplished some tasks that are otherwise either quite laborious or impractical for humans to verify, e.g., identifying conflicts amongst IMRT optimization objectives.

A fast inverse consistent deformable image registration method based on symmetric optical flow computation.

Deformable image registration is widely used in various radiation therapy applications including daily treatment planning adaptation to map planned tissue or dose to changing anatomy. In this work, a simple and efficient inverse consistency deformable registration method is proposed with aims of higher registration accuracy and faster convergence speed. Instead of registering image I to a second image J, the two images are symmetrically deformed toward one another in multiple passes, until both deformed images are matched and correct registration is therefore achieved. In each pass, a delta motion field is computed by minimizing a symmetric optical flow system cost function using modified optical flow algorithms. The images are then further deformed with the delta motion field in the positive and negative directions respectively, and then used for the next pass. The magnitude of the delta motion field is forced to be less than 0.4 voxel for every pass in order to guarantee smoothness and invertibility for the two overall motion fields that are accumulating the delta motion fields in both positive and negative directions, respectively. The final motion fields to register the original images I and J, in either direction, are calculated by inverting one overall motion field and combining the inversion result with the other overall motion field. The final motion fields are inversely consistent and this is ensured by the symmetric way that registration is carried out. The proposed method is demonstrated with phantom images, artificially deformed patient images and 4D-CT images. Our results suggest that the proposed method is able to improve the overall accuracy (reducing registration error by 30% or more, compared to the original and inversely inconsistent optical flow algorithms), reduce the inverse consistency error (by 95% or more) and increase the convergence rate (by 100% or more). The overall computation speed may slightly decrease, or increase in most cases because the new method converges faster. Compared to previously reported inverse consistency algorithms, the proposed method is simpler, easier to implement and more efficient.