The effect of integrating knowledge-based planning with multicriteria optimization in treatment planning for prostate SBRT.

Knowledge-based planning (KBP) and multicriteria optimization (MCO) are two powerful tools to assist treatment planners in achieving optimal target coverage and organ-at-risk (OAR) sparing. The purpose of this work is to investigate if integrating MCO with conventional KBP can further improve treatment plan quality for prostate cancer stereotactic body radiation therapy (SBRT). A two-phase study was designed to investigate the impact of MCO and KBP in prostate SBRT treatment planning. The first phase involved the creation of a KBP model based on thirty clinical SBRT plans, generated by manual optimization (KBP_M). A ten-patient validation cohort was used to compare manual, MCO, and KBP_M optimization techniques. The next phase involved replanning the original model cohort with additional tradeoff optimization via MCO to create a second model, KBP_MCO. Plans were then generated using linear integration (KBP_M+MCO), non-linear integration (KBP_MCO), and a combination of integration methods (KBP_MCO+MCO). All plans were analyzed for planning target volume (PTV) coverage, OAR constraints, and plan quality metrics. Comparisons were generated to evaluate plan and model quality. Phase 1 highlighted the necessity of KBP and MCO in treatment planning, as both optimization methods improved plan quality metrics (Conformity and Heterogeneity Indices) and reduced mean rectal dose by 2 Gy, as compared to manual planning. Integrating MCO with KBP did not further improve plan quality, as little significance was seen over KBP or MCO alone. Principal component score (PCS) fitting showed KBP_MCO improved bladder and rectum estimated and modeled dose correlation by 5% and 22%, respectively; however, model improvements did not significantly impact plan quality. KBP and MCO have shown to reduce OAR dose while maintaining desired PTV coverage in this study. Further integration of KBP and MCO did not show marked improvements in treatment plan quality while requiring increased time in model generation and optimization time.

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.

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.

A Method to Recognize Anatomical Site and Image Acquisition View in X-ray Images.

A method was developed to recognize anatomical site and image acquisition view automatically in 2D X-ray images that are used in image-guided radiation therapy. The purpose is to enable site and view dependent automation and optimization in the image processing tasks including 2D-2D image registration, 2D image contrast enhancement, and independent treatment site confirmation. The X-ray images for 180 patients of six disease sites (the brain, head-neck, breast, lung, abdomen, and pelvis) were included in this study with 30 patients each site and two images of orthogonal views each patient. A hierarchical multiclass recognition model was developed to recognize general site first and then specific site. Each node of the hierarchical model recognized the images using a feature extraction step based on principal component analysis followed by a binary classification step based on support vector machine. Given two images in known orthogonal views, the site recognition model achieved a 99% average F1 score across the six sites. If the views were unknown in the images, the average F1 score was 97%. If only one image was taken either with or without view information, the average F1 score was 94%. The accuracy of the site-specific view recognition models was 100%.

The role of oxygen in the photostimulation luminescence process of europium doped potassium chloride.

A recent suggestion that europium doped potassium chloride (KCl:Eu2+) has the potential to significantly advance the state-of-the-art in radiation therapy dosimetry has generated a renewed interest in a classic storage phosphor material. The purposes of this work are to investigate the role of oxygen in the photostimulation luminescence (PSL) process and to determine if both increased PSL yield and improved temporal stability could be realized in KCl:Eu2+ by incorporating oxygen in the material fabrication process. Regardless of synthesis atmosphere, air or pure nitrogen, PSL amplitude shows a maximum at 1.0 mol % Eu. Depending on europium concentration, dosimeters fabricated in air exhibit stronger PSL by a factor of 2 to 4 compared to those made in N2. There is no change in PSL stimulation spectrum while noticeable shifts in both photoluminescence and PSL emission spectra are observed for air versus nitrogen. Almost all charge-storage centers are spatially correlated, suggesting oxygen's stabilization role in the PSL process. However, oxygen alone does not improve material's temporal stability in the first few hours post irradiation at room temperature, probably because a significant portion of radiation-induced holes are stored in the Vk centers which are mobile.

Temperature dependence of the photostimulated luminescence in KCl:Eu2+

The goal of this work is to understand the physical mechanism behind the signal stabilization process in KCl:Eu2+, a storage phosphor material that has generated renewed interest due to its potential in radiation therapy dosimetry application. The temperature dependency of the photostimulated luminescence (PSL) spectra and intensity vs. time post x ray irradiation was measured. Commercial BaFBr:Eu2+ materials were included in this study for comparison. Unlike BaFBr:Eu2+, broadening of the F(Cl-) stimulation band and red-shift of the peak were observed for KCl:Eu2+ with increasing temperature. For irradiations at temperatures lower than 200 K, PSL intensity of KCl:Eu2+ showed recuperation behavior in the first 2 hrs post-irradiation and stayed almost constant with time thereafter. Moreover, spatially-correlated storage centers increased from 24% for irradiation at 50 K to 31% at 195 K and almost 100% at room temperature. The data suggest that certain types of charge storage-centers were mobile and contribute to the fast fading in PSL.

X ray storage performance of KCl:Eu2+ with high cumulated dose.

The effects of high cumulative radiation dose on the luminescence properties of KCl:Eu2+ are investigated. Pellet samples of KCl:Eu2+ were given doses of up to 200 kGy at the Louisiana State University Synchrotron facility. After synchrotron irradiation, samples were optically bleached and given a clinical dose of 2 Gy from a 6 MV medical linear accelerator. Optical properties were evaluated using photostimulated luminescence (PSL), photoluminescence (PL), and temperature-dependent PSL measurements. For a cumulated dose of up to 5-10 kGy, the PSL emission intensity increased by 15% compared to the PSL signal with no radiation history. For doses higher than 10 kGy, the PSL emission intensity retained at least 70% of the original intensity. Spatial correlation of the charge storage centers increased for doses up to 5 kGy and then decreased for higher cumulative doses. Emission band at 975 nm was attributed to transitions of Eu1+. PL spectra showed an intense peak centered at 420 nm for all cumulative doses. The results of this work show that KCl:Eu2+ storage phosphors are excellent reusable materials for radiation therapy dosimetry.

Physical phantom studies of helical cone-beam CT with exact reconstruction.

PURPOSE: Onboard cone-beam computed tomography (CBCT) connected to radiotherapy linear accelerators suffers CT number consistency and uniformity problems in addition to limited longitudinal coverage. Such problems have prevented CBCT from being fully utilized in many quantitative applications including tumor response evaluation and daily radiation dose computation. This paper presents a feasibility study on the helical CBCT scan with exact reconstruction that could be a potential solution. METHODS: A Varian TrueBeam treatment machine was programmed in the research mode to accomplish helical scans that required synchronized gantry circular rotation and couch table linear motion. Two physical phantoms were scanned in both 360° and 720° helical trajectories. A Katsevich exact reconstruction algorithm was implemented and tested with digital phantom simulations. It was further optimized to account for mechanical instabilities of both gantry rotation and couch table motion from the physical phantom measurements. Preprocessing was employed to correct photon scattering, beam hardening, and bowtie filtration. The reconstructed images were compared to those reconstructed from the FDK approximate reconstruction algorithm using the same phantom projections. Comparisons have also been made with the clinical circular CBCT images and the diagnostic helical CT images of the same physical phantoms. RESULTS: Satisfactory reconstruction results were obtained for the Katsevich algorithm in digital phantom study. Physical phantom results demonstrated that a 360° helical scan could provide up to 19 cm longitudinal coverage, which could be increased to 54 cm with a 720° helical scan. Image spatial resolution and soft tissue contrast were sufficient. The Q-value, which combined the spatial frequency response (modulation transfer function) and the image noise, was calculated, and suggested that the Katsevich algorithm was superior to the FDK algorithm. CONCLUSIONS: A helical CBCT scan is useful to extend the longitudinal coverage. The Katsevich exact reconstruction algorithm could provide additional advantages in image qualities over the traditional FDK approximate algorithm. The combination of helical CBCT scan with exact reconstruction was proved feasible and would render CBCT more useful in image-guided radiation therapy.

Technical note: electronic chart checks in a paperless radiation therapy clinic.

PURPOSE: EcCk, which stands for Electronic Chart ChecK, is a computer software and database system. It was developed to improve quality and efficiency of patient chart checking in radiation oncology departments. The core concept is to automatically collect and analyze patient treatment data, and to report discrepancies and potential concerns. METHODS: EcCk consists of several different computer technologies, including relational database, DICOM, dynamic HTML, and image processing. Implemented in MATLAB and C#, EcCk processes patient data in DICOM, PDF, Microsoft Word, database, and Pinnacle native formats. Generated reports are stored on the storage server and indexed in the database. A standalone report-browser program is implemented to allow users to view reports on any computer in the department. Checks are performed according to predefined logical rules, and results are presented through color-coded reports in which discrepancies are summarized and highlighted. Users examine the reports and take appropriate actions. The core design is intended to automate human task and to improve the reliability of the performed tasks. The software is not intended to replace human audits but rather to aid as a decision support tool. RESULTS: The software was successfully implemented in the clinical environment and has demonstrated the feasibility of automation of this common task with modern clinical tools. The software integrates multiple disconnected systems and successfully supports analysis of data in diverse formats. CONCLUSIONS: While the human is the ultimate expert, EcCk has a significant potential to improve quality and efficiency of patient treatment record audits, and to allow verification of tasks that are not easily performed by humans. EcCk can potentially relieve human experts from simple and repetitive tasks, and allow them to work on other important tasks, and in the end to improve the quality and safety of radiation therapy treatments.

Radiation hardness of the storage phosphor europium doped potassium chloride for radiation therapy dosimetry.

PURPOSE: An important property of a reusable dosimeter is its radiation hardness, that is, its ability to retain its dosimetric merits after irradiation. The radiation hardness of europium doped potassium chloride (KC1:Eu2+), a storage phosphor material recently proposed for radiation therapy dosimetry, is examined in this study. METHODS: Pellet-style KCl:Eu2+ dosimeters, 6 mm in diameter, and 1 mm thick, were fabricated in-house for this study. The pellets were exposed by a 6 MV photon beam or in a high dose rate 137Cs irradiator. Macroscopic properties, such as radiation sensitivity, dose response linearity, and signal stability, were studied with a laboratory photostimulated luminescence (PSL) readout system. Since phosphor performance is related to the state of the storage centers and the activator, Eu2+, in the host lattice, spectroscopic and temporal measurements were carried out in order to explore radiation-induced changes at the microscopic level. RESULTS: KCl:Eu2+ dosimeters retained approximately 90% of their initial signal strength after a 5000 Gy dose history. Dose response was initially supralinear over the dose range of 100-700 cGy but became linear after 60 Gy. Linearity did not change significantly in the 0-5000 Gy dose history spanned in this study. Annealing high dose history chips resulted in a return of supralinearity and a recovery of sensitivity. There were no significant changes in the PSL stimulation spectra, PSL emission spectra, photoluminescence spectra, or luminescence lifetime, indicating that the PSL signal process remains intact after irradiation but at a reduced efficiency due to reparable radiation-induced perturbations in the crystal lattice. CONCLUSIONS: Systematic studies of KCl:Eu2+ material are important for understanding how the material can be optimized for radiation therapy dosimetry purposes. The data presented here indicate that KCl:Eu2+ exhibits strong radiation hardness and lends support for further investigations of this novel material.

Dual-storage phosphor proton therapy dosimetry: Simultaneous quantification of dose and linear energy transfer.

PURPOSE: To investigate the feasibility of using the high Zeff storage phosphor material BaFBrI:Eu2+ in conjunction with the low Zeff storage phosphor material KCl:Eu2+ for simultaneous proton dose and linear energy transfer (LET) measurements by (a) measuring the fundamental optical and dosimetric properties of BaFBrI:Eu2+ , (b) evaluating its compatibility in being readout simultaneously with KCl:Eu2+ dosimeters, and (c) modeling and validating its LET dependence under elevated proton LET irradiation. METHODS: A commercial BaFBrI:Eu2+ storage phosphor detector (Model ST-VI, Fujifilm) was characterized with energy dispersive x-ray spectroscopy (EDS) analysis to obtain its elemental composition. The dosimeters were irradiated using both a Mevion S250 proton therapy unit (at the center of a spread-out Bragg peak, SOBP) and a Varian Clinac iX linear accelerator with the latter being a low LET irradiation. The photostimulated luminescence (PSL) emission spectra, excitation spectra, and luminescent lifetimes of the detector were measured after proton and photon irradiations. Dosimetric properties including dose linearity, dose rate dependence, radiation hardness, temporal, and readout stabilities were studied using a laboratory optical reader after proton irradiations. In addition, its proton energy dependence was analytically modeled and experimentally validated by irradiating the detectors at various depths within the SOBP (Range: 15.0 g/cm2 , Modulation: 10.0 g/cm2 ). RESULTS: The active detector composition for the high Zeff storage phosphor detector was found to be BaFBr0.85 I0.15 :Eu2+ . The BaFBr0.85 I0.15 :Eu2+ material's excitation and emission spectra were in agreement under proton and photon irradiations, with peaks of 586 ± 1 nm and 400 ± 1 nm, respectively, with a full width at half maximum (FWHM) of 119 ± 3 nm and 30 ± 2 nm, respectively. As dosimeter response under photon irradiation is generally believed to be free from LET effect, these results suggest LET independence of charge storage center types resulted from ionizing radiations. There is sufficient spectral overlaps with KCl:Eu2+ dosimeters allowing both dosimeters to be readout under equivalent readout conditions, that is, 594 nm stimulation and 420 nm detection wavelengths. Its PSL characteristic lifetime was found to be less than 5 microseconds which would make it suitable for fast 2D readout post irradiation. Its 420 nm emission band intensity was found to be linear up to 10 Gy absolute proton dose under the same irradiation conditions, dose rate independent, stable in time and under multiple readouts, and with high radiation hardness under cumulative proton dose histories up to 200 Gy as tested in this study. BaFBr0.85 I0.15 :Eu2+ showed significant proton energy-dependent dose under-response in regions of high LET which could be modeled by stopping power ratio calculations with an accuracy of 3% in low LET regions and a distance-to-agreement (DTA) of 1 mm in high LET regions (>5 keV/µm). CONCLUSION: We have proven the feasibility of dual-storage phosphor proton dosimetry for simultaneous proton dose and LET measurements. BaFBr0.85 I0.15 :Eu2+ has shown equally excellent dosimetry performance as its low Zeff complement KCl:Eu2+ with distinctive LET dependence merely as a result of its higher Zeff . These promising results pave the way for future studies involving simultaneous proton dose and LET measurements using this novel approach.

Characterization of Dosimetric Differences in Strut-Adjusted Volume Implant Treatment Plans Calculated With TG-43 Formalism and a Model-Based Dose Calculation Algorithm.

PURPOSE: To comprehensively characterize dosimetric differences between calculations with a commercial model-based dose calculation algorithm (MBDCA) and the TG-43 formalism in application to accelerated partial breast irradiation (APBI) with the strut-adjusted volume implant (SAVI) applicator. METHODS: Dose for 100 patients treated with the SAVI applicator was recalculated with an MBDCA for comparison to dose calculated via TG-43. For every pair of dose calculations, dose-volume histogram (DVH) metrics including V90%, V95%, V100%, V150%, and V200% for the PTV_EVAL were compared. Features were defined for each case including (1) applicator size, (2) ratio between PTV_EVAL contour and 1-cm rind surrounding SAVI applicator, (3) ratio between dwell time in central catheter and total dwell time, and (4) mean computed tomography (CT) number within the lumpectomy cavity. Wilcoxon rank sum tests were performed to test whether treatment plans could be stratified according to feature values into groups with statistically significant dosimetry differences between MBDCA and TG-43. RESULTS: For all DVH metrics, differences between TG-43 and MBDCA calculations were statistically significant (P < .05). Minimum (maximum) relative percent differences between the MBDCA and TG-43 for V90%, V95%, and V100% were -2.1% (0.1%), -3.1% (-0.1%), and -5.0% (-0.5%), respectively. The median relative percent difference in mean PTV_EVAL dose between the MBDCA and TG-43 was -3.9%, with minimum (maximum) difference of -6.5% (-1.8%). For V90%, V95%, and V100%, plan quality worsened beyond defined thresholds in 26, 23, and 31 cases with no instances of coverage improvement. Features 1, 2, and 4 were shown to be able to stratify treatment plans into groups with statistically significant differences in dosimetry metrics between MBDCA and TG-43. CONCLUSIONS: Investigated dose metrics for SAVI treatments were found to be systematically lower with MBDCA calculation in comparison to TG-43. Plans could be stratified according to several features by the magnitude of dosimetric differences between these calculations.

Target tracking using DMLC for volumetric modulated arc therapy: a simulation study.

PURPOSE: Target tracking using dynamic multileaf collimator (DMLC) is a promising approach for intrafraction motion management in radiation therapy. The purpose of this work is to develop a DMLC tracking algorithm capable of delivering volumetric-modulated arc therapy (VMAT) to the targets that experience two-dimensional (2D) rigid motion in the beam's eye view. METHODS: The problem of VMAT delivery to moving targets is formulated as a control problem with constraints. The relationships between gantry speed, gantry acceleration, MLC leaf-velocity, dose rate, and target motion are derived. An iterative search algorithm is developed to find numerical solutions for efficient delivery of a specific VMAT plan to the moving target using 2D DMLC tracking. The delivery of five VMAT lung plans is simulated. The planned and delivered fluence maps in the target-reference frame are calculated and compared. RESULTS: The simulation demonstrates that the 2D tracking algorithm is capable of delivering the VMAT plan to a moving target fast and accurately without violating the machine constraints and the integrity of the treatment plan. The average delivery time is only 29 s longer than that of no-tracking delivery, 101 versus 72 s, respectively. The fluence maps are normalized to 200 MU and the average root-mean-square error between the desired and the delivered fluence is 2.1 MU, compared to 14.8 MU for no-tracking and 3.6 MU for one-dimensional tracking. CONCLUSIONS: A locally optimal MLC tracking algorithm for VMAT delivery is proposed, aiming at shortest delivery time while maintaining treatment plan invariant. The inconsequential increase of treatment time due to DMLC tracking is clinically desirable, which makes VMAT with DMLC tracking attractive in treating moving tumors.

Theoretical and empirical investigations of KCl : Eu2+ for nearly water-equivalent radiotherapy dosimetry.

PURPOSE: The low effective atomic number, reusability, and other computed radiography-related advantages make europium doped potassium chloride (KCl : Eu2+) a promising dosimetry material. The purpose of this study is to model KCl : Eu2+ point dosimeters with a Monte Carlo (MC) method and, using this model, to investigate the dose responses of two-dimensional (2D) KCl : Eu2+ storage phosphor films (SPFs). METHODS: KCl : Eu2+ point dosimeters were irradiated using a 6 MV beam at four depths (5-20 cm) for each of five square field sizes (5 x 5-25 x 25 cm2). The dose measured by KCl : Eu2+ was compared to that measured by an ionization chamber to obtain the magnitude of energy dependent dose measurement artifact. The measurements were simulated using DOSXYZnrc with phase space files generated by BEAMnrcMP. Simulations were also performed for KCl : Eu2+ films with thicknesses ranging from 1 microm to 1 mm. The work function of the prototype KCl : Eu2+ material was determined by comparing the sensitivity of a 150 microm thick KCl : Eu2+ film to a commercial BaFBr0.85 I0.15 : Eu(2+)-based SPF with a known work function. The work function was then used to estimate the sensitivity of a 1 microm thick KCl : Eu2+ film. RESULTS: The simulated dose responses of prototype KCl : Eu2+ point dosimeters agree well with measurement data acquired by irradiating the dosimeters in the 6 MV beam with varying field size and depth. Furthermore, simulations with films demonstrate that an ultrathin KCl : Eu2+ film with thickness of the order of 1 microm would have nearly water-equivalent dose response. The simulation results can be understood using classic cavity theories. Finally, preliminary experiments and theoretical calculations show that ultrathin KCl : Eu2+ film could provide excellent signal in a 1 cGy dose-to-water irradiation. CONCLUSIONS: In conclusion, the authors demonstrate that KCl : Eu(2+)-based dosimeters can be accurately modeled by a MC method and that 2D KCl : Eu2+ films of the order of 1 microm thick would have minimal energy dependence. The data support the future research and development of a KCl : Eu2+ storage phosphor-based system for quantitative, high-resolution multidimensional radiation therapy dosimetry.

Dose uncertainty and resolution of polymer gel dosimetry using an MRI guided radiation therapy system's onboard 0.35 T scanner.

Magnetic Resonance Imaging (MRI) scanners are widely used for 3D gel dosimeters readout. However, limited access to MRI scanners is a challenge in MRI-based gel dosimetry. Recent clinical implementation of MRI-guided radiation therapy machines provides potential opportunities for onboard gel dosimetry using its MRI subsystem. The objective of this study was to investigate the feasibility of gel dosimetry using ViewRay's onboard 0.35 T MRI scanner. A BANG® polymer gel dosimeter was irradiated by three beams of 3 × 3 cm2 field size. The T2 relaxation rate (R2) of the irradiated gel was measured using a Philips 1.5 T Ingenia MRI and a ViewRay 0.35 T onboard MRI and spin-echo pulse sequences. The number of signal averages (NSA) was set to 16 for the ViewRay acquisitions and one for the Philips 1.5 T MRI to achieve similar signal-to-noise ratios. The in-plane spatial resolution was 1.5 × 1.5 mm2 and the slice thickness was 5 mm. The relative dose uncertainty was obtained using R2 versus dose curves to compare the performance of dosimetry using the two different MRIs and field strengths. The dose uncertainty decreased from 12% at 2 Gy to 3.5% at 7.5 Gy at 1.5 T. The dose uncertainty decreased from 13% at 2 Gy to 4% at 7.5 Gy with NSA = 16 and 3 × 3 mm2 pixel size, and from 10.5% at 2 Gy to 3.2% at 7.5 Gy with NSA = 16 and denoised R2 maps (1.5 × 1.5 mm2 pixel size) at 0.35 T. The mean of dose resolution was 0.4 Gy at 1.5 T while the mean of dose resolution was 0.8 Gy and 0.64 Gy at 0.35 T by downsampling and denoising the R2 map, respectively. Therefore, comparable dose uncertainty was achievable using the ViewRay's onboard 0.35 T and Philips 1.5 T MRI scanners. 3D gel dosimetry using onboard low-field MRI scanner provides ViewRay users a 3D high resolution dosimetry option besides film and ionization chamber.

Quantitative proton radiation therapy dosimetry using the storage phosphor europium-doped potassium chloride.

PURPOSE: To (a) characterize the fundamental optical and dosimetric properties of the storage phosphor europium-doped potassium chloride for quantitative proton dosimetry, and (b) investigate if its dose radiation response can be described by an analytic radiation transport model. METHODS: Cylindrical KCl:Eu2+ dosimeters with dimensions of 6 mm diameter and 1 mm thickness were fabricated in-house. The dosimeters were irradiated using both a Mevion S250 passive scattering proton therapy system and a Varian Clinac iX linear accelerator. Photostimulated luminescence (PSL) emission spectra, excitation spectra, and luminescence lifetimes were measured for both proton and photon irradiations. Dosimetric properties including radiation hardness, dose linearity, signal stabilization, dose rate sensitivity, and energy dependence were studied using a laboratory optical reader after irradiations. The dosimeters were modeled using physical quantities including mass stopping powers in the storage phosphor and water for a given proton beam, and mass energy absorption coefficients and massing stopping powers in detector and water for a given photon beam. RESULTS: KCl:Eu2+ exhibited optical emission and stimulation peaks at 421 and 560 nm, respectively, for both proton and photon irradiations, enabling postirradiation readouts using a visible light source while detecting the PSL using a photomultiplier tube. KCl:Eu2+ showed a linear response from 0 to 8 Gy absorbed dose-to-water, a large dynamic range up to 60 Gy, dose-rate independence measured from 83 to 500 MU/min, and a PSL lifetime of <5 ms that is sufficiently short for supporting rapid scanning in a two-dimensional geometry. KCl:Eu2+ was highly reusable with only a slight signal decrease of ~3% at accumulated doses over 100 Gy, which could be managed by a periodic recalibration. The detected PSL signal strength of the dosimeter in the proton field had been calculated accurately to a maximum discrepancy of 2% using known physical quantities along with its prior signal strength as measured in a photon field at the same dose-to-water. This discrepancy might be attributed to an under-response due to linear energy transfer (LET) effect. However, comparisons of depth-dose measurements in a spread-out Bragg peak (SOBP) field with a parallel-plate ionization chamber showed no clear evidence of LET effects. Furthermore, range measurements agreed with ionization chamber measurements to within 1 mm. CONCLUSIONS: KCl:Eu2+ showed linear response over a large dynamic range for proton irradiations and reliably reproduced SOBP measurements as measured by ionization chambers. Its relatively low atomic number of 18 and near LET independence make it suited for quantitative proton dosimetry. In addition, its high radiation hardness means that it can be reused numerous times. Any potential measurement artifacts encountered in complex irradiation conditions should be able to be corrected for using known physical quantities.

Quantitative megavoltage radiation therapy dosimetry using the storage phosphor KCl: Eu2+.

This work, for the first time, reports the use of europium doped potassium chloride (KCl:Eu2+) storage phosphor for quantitative megavoltage radiation therapy dosimetry. In principle, KCl:Eu2+ functions using the same photostimulatated luminescence (PSL) mechanism as commercially available BaFBr0.85I0.15:Eu2+ material that is used for computed radiography (CR) but features a significantly smaller effective atomic number--18 versus 49--making it a potentially useful material for nearly tissue-equivalent radiation dosimetry. Cylindrical KCl:Eu2+ dosimeters, 7 mm in diameter and 1 mm thick, were fabricated in-house. Dosimetric properties, including radiation hardness, response linearity, signal fading, dose rate sensitivity, and energy dependence, were studied with a laboratory optical reader after irradiation by a linear accelerator. The overall experimental uncertainty was estimated to be within +/-2.5%. The findings were (1) KCl:Eu2+ showed satisfactory radiation hardness. There was no significant change in the stimulation spectra after irradiation up to 200 Gy when compared to a fresh dosimeter, indicating that this material could be reused at least 100 times if 2 Gy per use was assumed, e.g., for patient-specific IMRT QA. (2) KCl:Eu2+ exhibited supralinear response to dose after irradiation from 0 to 800 cGy. (3) After x ray irradiation, the PSL signal faded with time and eventually reached a fading rate of about 0.1 % /h after 12 h. (4) The sensitivity of the dosimeter was independent of the dose rate ranging from 15 to 1000 cGy/min. (5) The sensitivity showed no beam energy dependence for either open x ray or megavoltage electron fields. (6) Over-response to low-energy scattered photons was comparable to radiographic film, e.g., Kodak EDR2 film. By sandwiching dosimeters between low-energy photon filters (0.3 mm thick lead foils) during irradiation, the over-response was reduced. The authors have demonstrated that KCl:Eu2+ dosimeters have many desirable dosimetric characteristics that make the material conducive to radiation therapy dosimetry. In the future, a large-area KCl:Eu2+-based CR plate with a thickness of the order of a few microns, created using modern thin film techniques, could provide a reusable, quantitative, high-resolution two-dimensional dosimeter with minimal energy dependence.

Development and Validation of a Bayesian Network Method to Detect External Beam Radiation Therapy Physician Order Errors.

PURPOSE: To investigate a Bayesian network (BN)-based method to detect errors in external beam radiation therapy physician orders. METHODS AND MATERIALS: A total of 4431 external beam radiation therapy orders from 2008 to 2017 at the authors' institution were obtained from clinical treatment management systems and divided into 3 groups: single prescription, concurrent boost, and sequential boost. Multiple BNs were developed for each group to detect errors in new orders using joint posterior probabilities of the order parameters, given disease information. Each BN was trained with a group of orders using a Bayesian learning algorithm. A procedure was developed to select the optimal BN for each treatment site in each group and to determine site-specific parameters and error detection thresholds. Potential clinical errors, created both manually and automatically, were applied to test error detection performance. RESULTS: The average true-positive rate (TPR) and false-positive rate (FPR) of error detection were 95.72% and 1.99%, respectively, for the single-prescription cohort with 9 treatment sites. For the concurrent-boost cohort, the TPR and FPR were 92.94% and 14.53%, respectively. For the sequential-boost cohort, the TPR and FPR were 100% and 9.48%, respectively, for the prescribed dose values and 100% and 4.34%, respectively, for the remaining order parameters. For the patient simulation and imaging parameters for 9 treatment sites, the TPR and FPR were 100% and 4.96%, respectively. CONCLUSIONS: The probabilistic BN method was able to perform physician order error detection at a higher accuracy than previously reported in a variety of complex prescription instances, thus warranting further development in incorporating BNs into clinical error detection tools to assist manual physician order checks.

An efficient model to guide prospective T2-weighted 4D magnetic resonance imaging acquisition.

PURPOSE: To establish a mathematical model to guide prospective T2-weighted four-dimensional magnetic resonance imaging (4DMRI) acquisition and to propose an efficient solution to expedite prospective T2-weighted 4DMRI acquisition. METHODS: Prospective T2-weighted 4DMRI acquisition was characterized by a mathematical model with 4DMRI acquisition time as the objective function and completeness of the image set, acquisition timing, image contrast, and image artifacts as constraints. Given the irregular nature of human respiration, an efficient solution based on the greedy strategy (ESGS) was proposed. The efficiency of the ESGS method was validated using healthy human subjects. Comparisons were made with the previous 4DMRI method incorporating the prefixed-order respiratory state splitting (PO-RSS) technique. RESULTS: 4DMRI image sets acquired using the ESGS and PO-RSS methods had similar image quality. The average time to acquire a 4DMRI image set covering 60 slices at 10 respiratory states was reduced by 30%, from 13.1 min using the PO-RSS method to 9.0 min using the ESGS method. It was demonstrated that high-quality T2-weighted 4DMRI could be obtained within a reasonable amount of time and all slices within each of the three-dimensional volumes were indeed acquired at the same respiratory state. CONCLUSIONS: The ESGS method substantially reduces the acquisition time for T2-weighted 4DMRI, making it ready for clinical evaluation to obtain abdominal tumor motion for radiotherapy treatment planning.