Engineering topics

Low, Daniel A.

Publications and source records attributed to Low, Daniel A..

At least 19 records

DOE OSTI2018

Using Big Data Analytics to Advance Precision Radiation Oncology

Big clinical data analytics as a primary component of precision medicine is discussed, identifying where these emerging tools fit in the spectrum of genomics and radiomics research. A learning health system (LHS) is conceptualized that uses clinically acquired data with machine learning to advance the initiatives of precision medicine. The LHS is comprehensive and can be used for clinical decision support, discovery, and hypothesis derivation. These developing uses can positively impact the ultimate management and therapeutic course for patients. The conceptual model for each use of clinical data, however, is different, and an overview of the implications is discussed. With advancements in technologies and culture to improve the efficiency, accuracy, and breadth of measurements of the patient condition, the concept of an LHS may be realized in precision radiation therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2016

Technical Note: Dosimetric effects of couch position variability on treatment plan quality with an MRI-guided Co-60 radiation therapy machine

Purpose: Magnetic resonance imaging (MRI) guidance in radiation therapy brings real-time imaging and adaptive planning into the treatment vault where it can account for interfraction and intrafraction movement of soft tissue. The only commercially available MRI-guided radiation therapy device is a three-head {sup 60}Co and MRI system with an integrated treatment planning system (TPS). Couch attenuation of the beam of up to 20% is well modeled in the TPS. Variations in the patient’s day-to-day position introduce discrepancies in the actual couch attenuation as modeled in the treatment plan. For this reason, the authors’ institution avoids plans with beams that pass through or near the couch edges. This study investigates the effects of differential beam attenuation by the couch due to couch shifts in order to determine whether couch edge avoidance restrictions can be lifted. Couch shifts were simulated using a Monte Carlo treatment planning system and ion chamber measurements performed for validation. Methods: A total of 27 plans from 23 patients were investigated. Couch shifts of 1 and 2 cm were introduced in combinations of lateral and vertical directions to simulate patient position variations giving 16 shifted plans per reference plan. The 1 and 2 cm shifts were based on shifts recorded in 320 treatment fractions. Results: Following TG176 recommendations for measurement methods, couch attenuation measurements agreed with TPS modeled attenuation to within 2.1%. Planning target volume D95 changed less than 1% for 1 and 2 cm couch shifts in only the x-direction and less than 3% for all directions. Conclusions: Dosimetry of all plans tested was robust to couch shifts up to ±2 cm. In general, couch shifts resulted in clinically insignificant dosimetric deviations. It is conceivable that in certain cases with large systematic couch shifts and plans that are particularly sensitive to shifts, dosimetric changes might rise to a clinically significant level.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

A comprehensive formulation for volumetric modulated arc therapy planning

Purpose: Volumetric modulated arc therapy (VMAT) is a widely employed radiation therapy technique, showing comparable dosimetry to static beam intensity modulated radiation therapy (IMRT) with reduced monitor units and treatment time. However, the current VMAT optimization has various greedy heuristics employed for an empirical solution, which jeopardizes plan consistency and quality. The authors introduce a novel direct aperture optimization method for VMAT to overcome these limitations. Methods: The comprehensive VMAT (comVMAT) planning was formulated as an optimization problem with an L2-norm fidelity term to penalize the difference between the optimized dose and the prescribed dose, as well as an anisotropic total variation term to promote piecewise continuity in the fluence maps, preparing it for direct aperture optimization. A level set function was used to describe the aperture shapes and the difference between aperture shapes at adjacent angles was penalized to control MLC motion range. A proximal-class optimization solver was adopted to solve the large scale optimization problem, and an alternating optimization strategy was implemented to solve the fluence intensity and aperture shapes simultaneously. Single arc comVMAT plans, utilizing 180 beams with 2° angular resolution, were generated for a glioblastoma multiforme case, a lung (LNG) case, and two head and neck cases—one with three PTVs (H&N{sub 3PTV}) and one with foue PTVs (H&N{sub 4PTV})—to test the efficacy. The plans were optimized using an alternating optimization strategy. The plans were compared against the clinical VMAT (clnVMAT) plans utilizing two overlapping coplanar arcs for treatment. Results: The optimization of the comVMAT plans had converged within 600 iterations of the block minimization algorithm. comVMAT plans were able to consistently reduce the dose to all organs-at-risk (OARs) as compared to the clnVMAT plans. On average, comVMAT plans reduced the max and mean OAR dose by 6.59% and 7.45%, respectively, of the prescription dose. Reductions in max dose and mean dose were as high as 14.5 Gy in the LNG case and 15.3 Gy in the H&N{sub 3PTV} case. PTV coverages measured by D95, D98, and D99 were within 0.25% of the prescription dose. By comprehensively optimizing all beams, the comVMAT optimizer gained the freedom to allow some selected beams to deliver higher intensities, yielding a dose distribution that resembles a static beam IMRT plan with beam orientation optimization. Conclusions: The novel nongreedy VMAT approach simultaneously optimizes all beams in an arc and then directly generates deliverable apertures. The single arc VMAT approach thus fully utilizes the digital Linac’s capability in dose rate and gantry rotation speed modulation. In practice, the new single VMAT algorithm generates plans superior to existing VMAT algorithms utilizing two arcs.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

A treatment planning comparison between modulated tri-cobalt-60 teletherapy and linear accelerator–based stereotactic body radiotherapy for central early-stage non−small cell lung cancer

We evaluated the feasibility of planning stereotactic body radiotherapy (SBRT) for large central early-stage non−small cell lung cancer with a tri-cobalt-60 (tri-{sup 60}Co) system equipped with real-time magnetic resonance imaging (MRI) guidance, as compared to linear accelerator (LINAC)–based SBRT. In all, 20 patients with large central early-stage non−small cell lung cancer who were treated between 2010 and 2015 with LINAC-based SBRT were replanned using a tri-{sup 60}Co system for a prescription dose of 50 Gy in 4 fractions. Doses to organs at risk were evaluated based on established MD Anderson constraints for central lung SBRT. R{sub 100} values were calculated as the total tissue volume receiving 100% of the dose (V{sub 100}) divided by the planning target volume and compared to assess dose conformity. Dosimetric comparisons between LINAC-based and tri-{sup 60}Co SBRT plans were performed using Student's t-test and Wilcoxon Ranks test. Blinded reviews by radiation oncologists were performed to assess the suitability of both plans for clinical delivery. The mean planning target volume was 48.3 cc (range: 12.1 to 139.4 cc). Of the tri-{sup 60}Co SBRT plans, a mean 97.4% of dosimetric parameters per patient met MD Anderson dose constraints, whereas a mean 98.8% of dosimetric parameters per patient were met with LINAC-based SBRT planning (p = 0.056). R{sub 100} values were similar between both plans (1.20 vs 1.21, p = 0.79). Upon blinded review by 4 radiation oncologists, an average of 90% of the tri-{sup 60}Co SBRT plans were considered acceptable for clinical delivery compared with 100% of the corresponding LINAC-based SBRT plans (p = 0.17). SBRT planning using the tri-{sup 60}Co system with built-in MRI is feasible and achieves clinically acceptable plans for most central lung patients, with similar target dose conformity and organ at risk dosimetry. The added benefit of real-time MRI-guided therapy may further optimize tumor targeting while improving normal tissue sparing, which warrants further investigation in a prospective feasibility clinical trial.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2016

A novel software and conceptual design of the hardware platform for intensity modulated radiation therapy

Purpose: To deliver high quality intensity modulated radiotherapy (IMRT) using a novel generalized sparse orthogonal collimators (SOCs), the authors introduce a novel direct aperture optimization (DAO) approach based on discrete rectangular representation. Methods: A total of seven patients—two glioblastoma multiforme, three head & neck (including one with three prescription doses), and two lung—were included. 20 noncoplanar beams were selected using a column generation and pricing optimization method. The SOC is a generalized conventional orthogonal collimators with N leaves in each collimator bank, where N = 1, 2, or 4. SOC degenerates to conventional jaws when N = 1. For SOC-based IMRT, rectangular aperture optimization (RAO) was performed to optimize the fluence maps using rectangular representation, producing fluence maps that can be directly converted into a set of deliverable rectangular apertures. In order to optimize the dose distribution and minimize the number of apertures used, the overall objective was formulated to incorporate an L2 penalty reflecting the difference between the prescription and the projected doses, and an L1 sparsity regularization term to encourage a low number of nonzero rectangular basis coefficients. The optimization problem was solved using the Chambolle–Pock algorithm, a first-order primal–dual algorithm. Performance of RAO was compared to conventional two-step IMRT optimization including fluence map optimization and direct stratification for multileaf collimator (MLC) segmentation (DMS) using the same number of segments. For the RAO plans, segment travel time for SOC delivery was evaluated for the N = 1, N = 2, and N = 4 SOC designs to characterize the improvement in delivery efficiency as a function of N. Results: Comparable PTV dose homogeneity and coverage were observed between the RAO and the DMS plans. The RAO plans were slightly superior to the DMS plans in sparing critical structures. On average, the maximum and mean critical organ doses were reduced by 1.94% and 1.44% of the prescription dose. The average number of delivery segments was 12.68 segments per beam for both the RAO and DMS plans. The N = 2 and N = 4 SOC designs were, on average, 1.56 and 1.80 times more efficient than the N = 1 SOC design to deliver. The mean aperture size produced by the RAO plans was 3.9 times larger than that of the DMS plans. Conclusions: The DAO and dose domain optimization approach enabled high quality IMRT plans using a low-complexity collimator setup. The dosimetric quality is comparable or slightly superior to conventional MLC-based IMRT plans using the same number of delivery segments. The SOC IMRT delivery efficiency can be significantly improved by increasing the leaf numbers, but the number is still significantly lower than the number of leaves in a typical MLC.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2016

Accuracy of UTE-MRI-based patient setup for brain cancer radiation therapy

Purpose: Radiation therapy simulations solely based on MRI have advantages compared to CT-based approaches. One feature readily available from computed tomography (CT) that would need to be reproduced with MR is the ability to compute digitally reconstructed radiographs (DRRs) for comparison against on-board radiographs commonly used for patient positioning. In this study, the authors generate MR-based bone images using a single ultrashort echo time (UTE) pulse sequence and quantify their 3D and 2D image registration accuracy to CT and radiographic images for treatments in the cranium. Methods: Seven brain cancer patients were scanned at 1.5 T using a radial UTE sequence. The sequence acquired two images at two different echo times. The two images were processed using an in-house software to generate the UTE bone images. The resultant bone images were rigidly registered to simulation CT data and the registration error was determined using manually annotated landmarks as references. DRRs were created based on UTE-MRI and registered to simulated on-board images (OBIs) and actual clinical 2D oblique images from ExacTrac™. Results: UTE-MRI resulted in well visualized cranial, facial, and vertebral bones that quantitatively matched the bones in the CT images with geometric measurement errors of less than 1 mm. The registration error between DRRs generated from 3D UTE-MRI and the simulated 2D OBIs or the clinical oblique x-ray images was also less than 1 mm for all patients. Conclusions: UTE-MRI-based DRRs appear to be promising for daily patient setup of brain cancer radiotherapy with kV on-board imaging.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2015

A Method for Assessing Ground-Truth Accuracy of the 5DCT Technique

Purpose: To develop a technique that assesses the accuracy of the breathing phase-specific volume image generation process by patient-specific breathing motion model using the original free-breathing computed tomographic (CT) scans as ground truths. Methods: Sixteen lung cancer patients underwent a previously published protocol in which 25 free-breathing fast helical CT scans were acquired with a simultaneous breathing surrogate. A patient-specific motion model was constructed based on the tissue displacements determined by a state-of-the-art deformable image registration. The first image was arbitrarily selected as the reference image. The motion model was used, along with the free-breathing phase information of the original 25 image datasets, to generate a set of deformation vector fields that mapped the reference image to the 24 nonreference images. The high-pitch helically acquired original scans served as ground truths because they captured the instantaneous tissue positions during free breathing. Image similarity between the simulated and the original scans was assessed using deformable registration that evaluated the pointwise discordance throughout the lungs. Results: Qualitative comparisons using image overlays showed excellent agreement between the simulated images and the original images. Even large 2-cm diaphragm displacements were very well modeled, as was sliding motion across the lung–chest wall boundary. The mean error across the patient cohort was 1.15 ± 0.37 mm, and the mean 95th percentile error was 2.47 ± 0.78 mm. Conclusion: The proposed ground truth–based technique provided voxel-by-voxel accuracy analysis that could identify organ-specific or tumor-specific motion modeling errors for treatment planning. Despite a large variety of breathing patterns and lung deformations during the free-breathing scanning session, the 5-dimensionl CT technique was able to accurately reproduce the original helical CT scans, suggesting its applicability to a wide range of patients.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2015

The development and verification of a highly accurate collision prediction model for automated noncoplanar plan delivery

Purpose: Significant dosimetric benefits had been previously demonstrated in highly noncoplanar treatment plans. In this study, the authors developed and verified an individualized collision model for the purpose of delivering highly noncoplanar radiotherapy and tested the feasibility of total delivery automation with Varian TrueBeam developer mode. Methods: A hand-held 3D scanner was used to capture the surfaces of an anthropomorphic phantom and a human subject, which were positioned with a computer-aided design model of a TrueBeam machine to create a detailed virtual geometrical collision model. The collision model included gantry, collimator, and couch motion degrees of freedom. The accuracy of the 3D scanner was validated by scanning a rigid cubical phantom with known dimensions. The collision model was then validated by generating 300 linear accelerator orientations corresponding to 300 gantry-to-couch and gantry-to-phantom distances, and comparing the corresponding distance measurements to their corresponding models. The linear accelerator orientations reflected uniformly sampled noncoplanar beam angles to the head, lung, and prostate. The distance discrepancies between measurements on the physical and virtual systems were used to estimate treatment-site-specific safety buffer distances with 0.1%, 0.01%, and 0.001% probability of collision between the gantry and couch or phantom. Plans containing 20 noncoplanar beams to the brain, lung, and prostate optimized via an in-house noncoplanar radiotherapy platform were converted into XML script for automated delivery and the entire delivery was recorded and timed to demonstrate the feasibility of automated delivery. Results: The 3D scanner measured the dimension of the 14 cm cubic phantom within 0.5 mm. The maximal absolute discrepancy between machine and model measurements for gantry-to-couch and gantry-to-phantom was 0.95 and 2.97 cm, respectively. The reduced accuracy of gantry-to-phantom measurements was attributed to phantom setup errors due to the slightly deformable and flexible phantom extremities. The estimated site-specific safety buffer distance with 0.001% probability of collision for (gantry-to-couch, gantry-to-phantom) was (1.23 cm, 3.35 cm), (1.01 cm, 3.99 cm), and (2.19 cm, 5.73 cm) for treatment to the head, lung, and prostate, respectively. Automated delivery to all three treatment sites was completed in 15 min and collision free using a digital Linac. Conclusions: An individualized collision prediction model for the purpose of noncoplanar beam delivery was developed and verified. With the model, the study has demonstrated the feasibility of predicting deliverable beams for an individual patient and then guiding fully automated noncoplanar treatment delivery. This work motivates development of clinical workflows and quality assurance procedures to allow more extensive use and automation of noncoplanar beam geometries.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

Technology for Innovation in Radiation Oncology

Radiation therapy is an effective, personalized cancer treatment that has benefited from technological advances associated with the growing ability to identify and target tumors with accuracy and precision. Given that these advances have played a central role in the success of radiation therapy as a major component of comprehensive cancer care, the American Society for Radiation Oncology (ASTRO), the American Association of Physicists in Medicine (AAPM), and the National Cancer Institute (NCI) sponsored a workshop entitled “Technology for Innovation in Radiation Oncology,” which took place at the National Institutes of Health (NIH) in Bethesda, Maryland, on June 13 and 14, 2013. The purpose of this workshop was to discuss emerging technology for the field and to recognize areas for greater research investment. Expert clinicians and scientists discussed innovative technology in radiation oncology, in particular as to how these technologies are being developed and translated to clinical practice in the face of current and future challenges and opportunities. Technologies encompassed topics in functional imaging, treatment devices, nanotechnology, and information technology. The technical, quality, and safety performance of these technologies were also considered. A major theme of the workshop was the growing importance of innovation in the domain of process automation and oncology informatics. The technologically advanced nature of radiation therapy treatments predisposes radiation oncology research teams to take on informatics research initiatives. In addition, the discussion on technology development was balanced with a parallel conversation regarding the need for evidence of efficacy and effectiveness. The linkage between the need for evidence and the efforts in informatics research was clearly identified as synergistic.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2015

Near Real-Time Assessment of Anatomic and Dosimetric Variations for Head and Neck Radiation Therapy via Graphics Processing Unit–based Dose Deformation Framework

Purpose: The purpose of this study was to systematically monitor anatomic variations and their dosimetric consequences during intensity modulated radiation therapy (IMRT) for head and neck (H&N) cancer by using a graphics processing unit (GPU)-based deformable image registration (DIR) framework. Methods and Materials: Eleven IMRT H&N patients undergoing IMRT with daily megavoltage computed tomography (CT) and weekly kilovoltage CT (kVCT) scans were included in this analysis. Pretreatment kVCTs were automatically registered with their corresponding planning CTs through a GPU-based DIR framework. The deformation of each contoured structure in the H&N region was computed to account for nonrigid change in the patient setup. The Jacobian determinant of the planning target volumes and the surrounding critical structures were used to quantify anatomical volume changes. The actual delivered dose was calculated accounting for the organ deformation. The dose distribution uncertainties due to registration errors were estimated using a landmark-based gamma evaluation. Results: Dramatic interfractional anatomic changes were observed. During the treatment course of 6 to 7 weeks, the parotid gland volumes changed up to 34.7%, and the center-of-mass displacement of the 2 parotid glands varied in the range of 0.9 to 8.8 mm. For the primary treatment volume, the cumulative minimum and mean and equivalent uniform doses assessed by the weekly kVCTs were lower than the planned doses by up to 14.9% (P=.14), 2% (P=.39), and 7.3% (P=.05), respectively. The cumulative mean doses were significantly higher than the planned dose for the left parotid (P=.03) and right parotid glands (P=.006). The computation including DIR and dose accumulation was ultrafast (∼45 seconds) with registration accuracy at the subvoxel level. Conclusions: A systematic analysis of anatomic variations in the H&N region and their dosimetric consequences is critical in improving treatment efficacy. Nearly real-time assessment of anatomic and dosimetric variations is feasible using the GPU-based DIR framework. Clinical implementation of this technology may enable timely plan adaptation and improved outcome.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2015

High-Quality T2-Weighted 4-Dimensional Magnetic Resonance Imaging for Radiation Therapy Applications

Purpose: The purpose of this study was to improve triggering efficiency of the prospective respiratory amplitude-triggered 4-dimensional magnetic resonance imaging (4DMRI) method and to develop a 4DMRI imaging protocol that could offer T2 weighting for better tumor visualization, good spatial coverage and spatial resolution, and respiratory motion sampling within a reasonable amount of time for radiation therapy applications. Methods and Materials: The respiratory state splitting (RSS) and multi-shot acquisition (MSA) methods were analytically compared and validated in a simulation study by using the respiratory signals from 10 healthy human subjects. The RSS method was more effective in improving triggering efficiency. It was implemented in prospective respiratory amplitude-triggered 4DMRI. 4DMRI image datasets were acquired from 5 healthy human subjects. Liver motion was estimated using the acquired 4DMRI image datasets. Results: The simulation study showed the RSS method was more effective for improving triggering efficiency than the MSA method. The average reductions in 4DMRI acquisition times were 36% and 10% for the RSS and MSA methods, respectively. The human subject study showed that T2-weighted 4DMRI with 10 respiratory states, 60 slices at a spatial resolution of 1.5 × 1.5 × 3.0 mm{sup 3} could be acquired in 9 to 18 minutes, depending on the individual's breath pattern. Based on the acquired 4DMRI image datasets, the ranges of peak-to-peak liver displacements among 5 human subjects were 9.0 to 12.9 mm, 2.5 to 3.9 mm, and 0.5 to 2.3 mm in superior-inferior, anterior-posterior, and left-right directions, respectively. Conclusions: We demonstrated that with the RSS method, it was feasible to acquire high-quality T2-weighted 4DMRI within a reasonable amount of time for radiation therapy applications.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2015

Dose domain regularization of MLC leaf patterns for highly complex IMRT plans

Purpose: The advent of automated beam orientation and fluence optimization enables more complex intensity modulated radiation therapy (IMRT) planning using an increasing number of fields to exploit the expanded solution space. This has created a challenge in converting complex fluences to robust multileaf collimator (MLC) segments for delivery. A novel method to regularize the fluence map and simplify MLC segments is introduced to maximize delivery efficiency, accuracy, and plan quality. Methods: In this work, we implemented a novel approach to regularize optimized fluences in the dose domain. The treatment planning problem was formulated in an optimization framework to minimize the segmentation-induced dose distribution degradation subject to a total variation regularization to encourage piecewise smoothness in fluence maps. The optimization problem was solved using a first-order primal-dual algorithm known as the Chambolle-Pock algorithm. Plans for 2 GBM, 2 head and neck, and 2 lung patients were created using 20 automatically selected and optimized noncoplanar beams. The fluence was first regularized using Chambolle-Pock and then stratified into equal steps, and the MLC segments were calculated using a previously described level reducing method. Isolated apertures with sizes smaller than preset thresholds of 1–3 bixels, which are square units of an IMRT fluence map from MLC discretization, were removed from the MLC segments. Performance of the dose domain regularized (DDR) fluences was compared to direct stratification and direct MLC segmentation (DMS) of the fluences using level reduction without dose domain fluence regularization. Results: For all six cases, the DDR method increased the average planning target volume dose homogeneity (D95/D5) from 0.814 to 0.878 while maintaining equivalent dose to organs at risk (OARs). Regularized fluences were more robust to MLC sequencing, particularly to the stratification and small aperture removal. The maximum and mean aperture sizes using the DDR were consistently larger than those from DMS for all tested number of segments. Conclusions: The fluence map to MLC segmentation conversion problem was formulated as a secondary optimization problem in the dose domain to minimize the smoothness-regularized dose discrepancy. The large scale optimization problem was solved using a primal-dual algorithm that transformed complicated fluences into maps that were more robust to the MLC segmentation and sequencing, affording fewer and larger segments with minimal degradation to dose distribution.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

Accuracy of Routine Treatment Planning 4-Dimensional and Deep-Inspiration Breath-Hold Computed Tomography Delineation of the Left Anterior Descending Artery in Radiation Therapy

Purpose: To assess the feasibility of radiation therapy treatment planning 4-dimensional computed tomography (4DCT) and deep-inspiration breath-hold (DIBH) CT to accurately contour the left anterior descending artery (LAD), a primary indicator of radiation-induced cardiac toxicity for patients undergoing radiation therapy. Methods and Materials: Ten subjects were prospectively imaged with a cardiac-gated MRI protocol to determine cardiac motion effects, including the displacement of a region of interest comprising the LAD. A series of planar views were obtained and resampled to create a 3-dimensional (3D) volume. A 3D optical flow deformable image registration algorithm determined tissue displacement during the cardiac cycle. The measured motion was then used as a spatial boundary to characterize motion blurring of the radiologist-delineated LAD structure for a cohort of 10 consecutive patients enrolled prospectively on a breast study including 4DCT and DIBH scans. Coronary motion–induced blurring artifacts were quantified by applying an unsharp filter to accentuate the LAD structure despite the presence of motion blurring. The 4DCT maximum inhalation and exhalation respiratory phases were coregistered to determine the LAD displacement during tidal respiration, as visualized in 4DCT. Results: The average 90th percentile heart motion for the region of interest was 0.7 ± 0.1 mm (left–right [LR]), 1.3 ± 0.6 mm (superior–inferior [SI]), and 0.6 ± 0.2 mm (anterior–posterior [AP]) in the cardiac-gated MRI cohort. The average relative increase in the number of voxels comprising the LAD contour was 69.4% ± 4.5% for the DIBH. The LAD volume overestimation had the dosimetric impact of decreasing the reported mean LAD dose by 23% ± 9% on average in the DIBH. During tidal respiration the average relative LAD contour increase was 69.3% ± 5.9% and 67.9% ± 4.6% for inhalation and exhalation respiratory phases, respectively. The average 90th percentile LAD motion was 4.8 ± 1.1 mm (LR), 0.9 ± 0.4 mm (SI), and 1.9 ± 0.6 mm (AP) for the 4DCT cohort, in the absence of cardiac gating. Conclusions: An anisotropic margin of 2.7 mm (LR), 4.1 mm (SI), and 2.4 mm (AP) was quantitatively determined to account for motion blurring and patient setup error while placing minimum constraint on the plan optimization.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2015

Dependence of Achievable Plan Quality on Treatment Technique and Planning Goal Refinement: A Head-and-Neck Intensity Modulated Radiation Therapy Application

Purpose: To develop a practical workflow for retrospectively analyzing target and normal tissue dose–volume endpoints for various intensity modulated radiation therapy (IMRT) delivery techniques; to develop technique-specific planning goals to improve plan consistency and quality when feasible. Methods and Materials: A total of 165 consecutive head-and-neck patients from our patient registry were selected and retrospectively analyzed. All IMRT plans were generated using the same dose–volume guidelines for TomoTherapy (Tomo, Accuray), TrueBeam (TB, Varian) using fixed-field IMRT (TB-IMRT) or RAPIDARC (TB-RAPIDARC), or Siemens Oncor (Siemens-IMRT, Siemens). A MATLAB-based dose–volume extraction and analysis tool was developed to export dosimetric endpoints for each patient. With a fair stratification of patient cohort, the variation of achieved dosimetric endpoints was analyzed among different treatment techniques. Upon identification of statistically significant variations, technique-specific planning goals were derived from dynamically accumulated institutional data. Results: Retrospective analysis showed that although all techniques yielded comparable target coverage, the doses to the critical structures differed. The maximum cord doses were 34.1 ± 2.6, 42.7 ± 2.1, 43.3 ± 2.0, and 45.1 ± 1.6 Gy for Tomo, TB-IMRT, TB-RAPIDARC, and Siemens-IMRT plans, respectively. Analyses of variance showed significant differences for the maximum cord doses but no significant differences for other selected structures among the investigated IMRT delivery techniques. Subsequently, a refined technique-specific dose–volume guideline for maximum cord dose was derived at a confidence level of 95%. The dosimetric plans that failed the refined technique-specific planning goals were reoptimized according to the refined constraints. We observed better cord sparing with minimal variations for the target coverage and other organ at risk sparing for the Tomo cases, and higher parotid doses for C-arm linear accelerator–based IMRT and RAPIDARC plans. Conclusion: Patient registry–based processes allowed easy and systematic dosimetric assessment of treatment plan quality and consistency. Our analysis revealed the dependence of certain dosimetric endpoints on the treatment techniques. Technique-specific refinement of planning goals may lead to improvement in plan consistency and plan quality.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2015

Incorporating Cancer Stem Cells in Radiation Therapy Treatment Response Modeling and the Implication in Glioblastoma Multiforme Treatment Resistance

Purpose: To perform a preliminary exploration with a simplistic mathematical cancer stem cell (CSC) interaction model to determine whether the tumor-intrinsic heterogeneity and dynamic equilibrium between CSCs and differentiated cancer cells (DCCs) can better explain radiation therapy treatment response with a dual-compartment linear-quadratic (DLQ) model. Methods and Materials: The radiosensitivity parameters of CSCs and DCCs for cancer cell lines including glioblastoma multiforme (GBM), non–small cell lung cancer, melanoma, osteosarcoma, and prostate, cervical, and breast cancer were determined by performing robust least-square fitting using the DLQ model on published clonogenic survival data. Fitting performance was compared with the single-compartment LQ (SLQ) and universal survival curve models. The fitting results were then used in an ordinary differential equation describing the kinetics of DCCs and CSCs in response to 2- to 14.3-Gy fractionated treatments. The total dose to achieve tumor control and the fraction size that achieved the least normal biological equivalent dose were calculated. Results: Smaller cell survival fitting errors were observed using DLQ, with the exception of melanoma, which had a low α/β = 0.16 in SLQ. Ordinary differential equation simulation indicated lower normal tissue biological equivalent dose to achieve the same tumor control with a hypofractionated approach for 4 cell lines for the DLQ model, in contrast to SLQ, which favored 2 Gy per fraction for all cells except melanoma. The DLQ model indicated greater tumor radioresistance than SLQ, but the radioresistance was overcome by hypofractionation, other than the GBM cells, which responded poorly to all fractionations. Conclusion: The distinct radiosensitivity and dynamics between CSCs and DCCs in radiation therapy response could perhaps be one possible explanation for the heterogeneous intertumor response to hypofractionation and in some cases superior outcome from stereotactic ablative radiation therapy. The DLQ model also predicted the remarkable GBM radioresistance, a result that is highly consistent with clinical observations. The radioresistance putatively stemmed from accelerated DCC regrowth that rapidly restored compartmental equilibrium between CSCs and DCCs.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2015

Automated contouring error detection based on supervised geometric attribute distribution models for radiation therapy: A general strategy

Purpose: One of the most critical steps in radiation therapy treatment is accurate tumor and critical organ-at-risk (OAR) contouring. Both manual and automated contouring processes are prone to errors and to a large degree of inter- and intraobserver variability. These are often due to the limitations of imaging techniques in visualizing human anatomy as well as to inherent anatomical variability among individuals. Physicians/physicists have to reverify all the radiation therapy contours of every patient before using them for treatment planning, which is tedious, laborious, and still not an error-free process. In this study, the authors developed a general strategy based on novel geometric attribute distribution (GAD) models to automatically detect radiation therapy OAR contouring errors and facilitate the current clinical workflow. Methods: Considering the radiation therapy structures’ geometric attributes (centroid, volume, and shape), the spatial relationship of neighboring structures, as well as anatomical similarity of individual contours among patients, the authors established GAD models to characterize the interstructural centroid and volume variations, and the intrastructural shape variations of each individual structure. The GAD models are scalable and deformable, and constrained by their respective principal attribute variations calculated from training sets with verified OAR contours. A new iterative weighted GAD model-fitting algorithm was developed for contouring error detection. Receiver operating characteristic (ROC) analysis was employed in a unique way to optimize the model parameters to satisfy clinical requirements. A total of forty-four head-and-neck patient cases, each of which includes nine critical OAR contours, were utilized to demonstrate the proposed strategy. Twenty-nine out of these forty-four patient cases were utilized to train the inter- and intrastructural GAD models. These training data and the remaining fifteen testing data sets were separately employed to test the effectiveness of the proposed contouring error detection strategy. Results: An evaluation tool was implemented to illustrate how the proposed strategy automatically detects the radiation therapy contouring errors for a given patient and provides 3D graphical visualization of error detection results as well. The contouring error detection results were achieved with an average sensitivity of 0.954/0.906 and an average specificity of 0.901/0.909 on the centroid/volume related contouring errors of all the tested samples. As for the detection results on structural shape related contouring errors, an average sensitivity of 0.816 and an average specificity of 0.94 on all the tested samples were obtained. The promising results indicated the feasibility of the proposed strategy for the detection of contouring errors with low false detection rate. Conclusions: The proposed strategy can reliably identify contouring errors based upon inter- and intrastructural constraints derived from clinically approved contours. It holds great potential for improving the radiation therapy workflow. ROC and box plot analyses allow for analytically tuning of the system parameters to satisfy clinical requirements. Future work will focus on the improvement of strategy reliability by utilizing more training sets and additional geometric attribute constraints.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

4π Noncoplanar Stereotactic Body Radiation Therapy for Head-and-Neck Cancer: Potential to Improve Tumor Control and Late Toxicity

Purpose: To evaluate the potential benefit of 4π radiation therapy in recurrent, locally advanced, or metastatic head-and-neck cancer treated with stereotactic body radiation therapy (SBRT). Methods and Materials: Twenty-seven patients with 29 tumors who were treated using SBRT were included. In recurrent disease (n=26), SBRT was delivered with a median 44 Gy (range, 35-44 Gy) in 5 fractions. Three patients with sinonasal mucosal melanoma, metastatic breast cancer, and primary undifferentiated carcinoma received 35 Gy, 22.5 Gy, and 40 Gy in 5 fractions, respectively. Novel 4π treatment plans were created for each patient to meet the objective that 95% of the planning target volume was covered by 100% of the prescription dose. Doses to organs at risk (OARs) and 50% dose spillage volumes were compared against the delivered clinical SBRT plans. Local control (LC), late toxicity, tumor control probability (TCP), and normal tissue complication probability were determined. Results: Using 4π plans, mean/maximum doses to all OARs were reduced by 22% to 89%/10% to 86%. With 4π plans, the 50% dose spillage volume was decreased by 33%. Planning target volume prescription dose escalation by 10 Gy and 20 Gy were achieved while keeping doses to OARs significantly improved or unchanged from clinical plans, except for the carotid artery maximum dose at 20-Gy escalation. At a median follow-up of 10 months (range, 1-41 months), crude LC was 52%. The 2-year LC of 39.2% approximated the predicted mean TCP of 42.2%, which increased to 45.9% with 4π plans. For 10-Gy and 20-Gy dose escalation, 4π plans increased TCP from 80.1% and 88.1% to 85.5% and 91.4%, respectively. The 7.4% rate of grade ≥3 late toxicity was comparable to the predicted 5.6% mean normal tissue complication probability for OARs, which was significantly reduced by 4π planning at the prescribed and escalated doses. Conclusions: 4π plans may allow dose escalation with significant and consistent improvements in critical organ sparing, tumor control, and coverage.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

Feasibility of using intermediate x-ray energies for highly conformal extracranial radiotherapy

Purpose: To investigate the feasibility of using intermediate energy 2 MV x-rays for extracranial robotic intensity modulated radiation therapy. Methods: Two megavolts flattening filter free x-rays were simulated using the Monte Carlo code MCNP (v4c). A convolution/superposition dose calculation program was tuned to match the Monte Carlo calculation. The modeled 2 MV x-rays and actual 6 MV flattened x-rays from existing Varian Linacs were used in integrated beam orientation and fluence optimization for a head and neck, a liver, a lung, and a partial breast treatment. A column generation algorithm was used for the intensity modulation and beam orientation optimization. Identical optimization parameters were applied in three different planning modes for each site: 2, 6 MV, and dual energy 2/6 MV. Results: Excellent agreement was observed between the convolution/superposition and the Monte Carlo calculated percent depth dose profiles. For the patient plans, overall, the 2/6 MV x-ray plans had the best dosimetry followed by 2 MV only and 6 MV only plans. Between the two single energy plans, the PTV coverage was equivalent but 2 MV x-rays improved organs-at-risk sparing. For the head and neck case, the 2MV plan reduced lips, mandible, tongue, oral cavity, brain, larynx, left and right parotid gland mean doses by 14%, 8%, 4%, 14%, 24%, 6%, 30% and 16%, respectively. For the liver case, the 2 MV plan reduced the liver and body mean doses by 17% and 18%, respectively. For the lung case, lung V20, V10, and V5 were reduced by 13%, 25%, and 30%, respectively. V10 of heart with 2 MV plan was reduced by 59%. For the partial breast treatment, the 2 MV plan reduced the mean dose to the ipsilateral and contralateral lungs by 27% and 47%, respectively. The mean body dose was reduced by 16%. Conclusions: The authors showed the feasibility of using flattening filter free 2 MV x-rays for extracranial treatments as evidenced by equivalent or superior dosimetry compared to 6 MV plans using the same inverse noncoplanar intensity modulated planning method.

62 RADIOLOGY AND NUCLEAR MEDICINE↗