Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Histogram”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

MO-DE-207A-06: ECG-Gated CT Reconstruction for a C-Arm Inverse Geometry X-Ray System

Purpose: To obtain ECG-gated CT images from truncated projection data acquired with a C-arm based inverse geometry fluoroscopy system, for the purpose of cardiac chamber mapping in interventional procedures. Methods: Scanning-beam digital x-ray (SBDX) is an inverse geometry fluoroscopy system with a scanned multisource x-ray tube and a photon-counting detector mounted to a C-arm. In the proposed method, SBDX short-scan rotational acquisition is performed followed by inverse geometry CT (IGCT) reconstruction and segmentation of contrast-enhanced objects. The prior image constrained compressed sensing (PICCS) framework was adapted for IGCT reconstruction to mitigate artifacts arising from data truncation and angular undersampling due to cardiac gating. The performance of the reconstruction algorithm was evaluated in numerical simulations of truncated and non-truncated thorax phantoms containing a dynamic ellipsoid to represent a moving cardiac chamber. The eccentricity of the ellipsoid was varied at frequencies from 1–1.5 Hz. Projection data were retrospectively sorted into 13 cardiac phases. Each phase was reconstructed using IGCT-PICCS, with a nongated gridded FBP (gFBP) prior image. Surface accuracy was determined using Dice similarity coefficient and a histogram of the point distances between the segmented surface and ground truth surface. Results: The gated IGCT-PICCS algorithm improved surface accuracy and reduced streaking and truncation artifacts when compared to nongated gFBP. For the non-truncated thorax with 1.25 Hz motion, 99% of segmented surface points were within 0.3 mm of the 15 mm diameter ground truth ellipse, versus 1.0 mm for gFBP. For the truncated thorax phantom with a 40 mm diameter ellipse, IGCT-PICCS surface accuracy measured 0.3 mm versus 7.8 mm for gFBP. Dice similarity coefficient was 0.99–1.00 (IGCT-PICCS) versus 0.63–0.75 (gFBP) for intensity-based segmentation thresholds ranging from 25–75% maximum contrast. Conclusions: The PICCS algorithm was successfully applied to reconstruct truncated IGCT projection data with angular undersampling resulting from simulated cardiac gating. Research supported by the National Heart, Lung, and Blood Institute of the NIH under award number R01HL084022. The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

MO-DE-207A-09: Low-Dose CT Image Reconstruction Via Learning From Different Patient Normal-Dose Images

Purpose: To investigate a novel low-dose CT (LdCT) image reconstruction strategy for lung CT imaging in radiation therapy. Methods: The proposed approach consists of four steps: (1) use the traditional filtered back-projection (FBP) method to reconstruct the LdCT image; (2) calculate structure similarity (SSIM) index between the FBP-reconstructed LdCT image and a set of normal-dose CT (NdCT) images, and select the NdCT image with the highest SSIM as the learning source; (3) segment the NdCT source image into lung and outside tissue regions via simple thresholding, and adopt multiple linear regression to learn high-order Markov random field (MRF) pattern for each tissue region in the NdCT source image; (4) segment the FBP-reconstructed LdCT image into lung and outside regions as well, and apply the learnt MRF prior in each tissue region for statistical iterative reconstruction of the LdCT image following the penalized weighted least squares (PWLS) framework. Quantitative evaluation of the reconstructed images was based on the signal-to-noise ratio (SNR), local binary pattern (LBP) and histogram of oriented gradients (HOG) metrics. Results: It was observed that lung and outside tissue regions have different MRF patterns predicted from the NdCT. Visual inspection showed that our method obviously outperformed the traditional FBP method. Comparing with the region-smoothing PWLS method, our method has, in average, 13% increase in SNR, 15% decrease in LBP difference, and 12% decrease in HOG difference from reference standard for all regions of interest, which indicated the superior performance of the proposed method in terms of image resolution and texture preservation. Conclusion: We proposed a novel LdCT image reconstruction method by learning similar image characteristics from a set of NdCT images, and the to-be-learnt NdCT image does not need to be scans from the same subject. This approach is particularly important for enhancing image quality in radiation therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

MO-DE-207B-10: Impact of Morphologic Characteristics On Radiomics Features From Contast-Enhanced CT for Primary Lung Tumors

Purpose: Determine the impact of morphologic characteristics (e.g. necrosis, vascular enhancement, and cavitation) on radiomic features from contrast enhanced CT (CE-CT) in primary lung tumors. Methods: We developed an auto-segmentation algorithm to separate lung tumors on contrast-enhanced CT into cavitation (air), necrosis, tissue, and enhancing vessels using a combination of thresholding and region-growing. An auto-segmentation algorithm was also designed to identify necrosis on FDG-PET scans. Wilcoxon rank-sum tests were used to determine if significant differences existed in radiomics features (histogram-uniformity and Laplacian-of-Gaussian average) from 249 patients, found to prognostic in previous work, based on the presence/absence of morphologic features. Feature values were also compared between the original tumor contours and contours excluding a specific morphologic feature. Comparison of necrosis segmentation on CE-CT versus FDG-PET was performed in 78 patients to assess for agreement using the concordance correlation coefficient (CCC). Results: Tumors with cavitation and enhancing vasculature had lower uniformity values (p = 0.001 and p = 0.03, respectively). Tumors with enhancing vasculature and necrosis had higher Laplacian-of-Gaussian average values (measure of “edges” within the tumor) (p < 0.001). Removing these tissue types from regions-of-interest did not drastically alter either radiomic feature value (all scenarios had R{sup 2} > 0.8). This suggests there may be interactions between morphologic characteristics and the radiomic feature value of tumor tissue. Comparison of necrosis volume and percent necrosis volume of tumor were found to have CCC values of 0.85 and 0.76, respectively between CE-CT and FDG-PET segmentation methods. Conclusions: Tumors with enhancing vasculature, necrosis, and cavitation have higher radiomic feature values that are associated with poor prognosis than tumors without these features. Removing these tissue types from quantitative assessment did not drastically impact radiomic feature values. High reproducibility of CE-CT segmented necrosis compared to FDG-PET segmented necrosis provides a reasonable validation of segmentation accuracy on CE-CT.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

MO-FG-BRA-04: A Novel Time Weighted Density Correction for Stereotactic Lung Radiotherapy: A Phantom Study

Purpose: Conventional treatment plans for lung radiotherapy are created using either the free breathing (FB) scheme which represents the tumor at an arbitrary breathing phase of the patient’s respiratory cycle, or the average computed tomography (ACT) intensity projection over 10-binned phases. Neither method is entirely accurate because of the absence of time dependence of tumor movement. In the present “Hybrid” method, the HU of tumor in 3D space is determined by relative weighting of the HU of the tumor and lung in proportion to the time they spend at that location during the entire breathing cycle. Methods: A Quasar respiratory motion phantom was employed to simulate lung tumor movement. Utilizing 4DCT image scans, volumetric modulated arc therapy (VMAT) plans were generated for three treatment planning scenarios which included conventional FB and ACT schemes, along with a third alternative Hybrid approach. Our internal target volume (ITV) hybrid structure was created using Boolean operation in Eclipse (ver. 11) treatment planning system, where independent sub-regions created by the gross tumor volume (GTV) overlap from the 10 motion phases were each assigned a time weighted CT value. The dose-volume-histograms (DVH) for each scheme were compared and analyzed. Results: Using our hybrid technique, we have demonstrated a reduction of 1.9% – 3.4% in total monitor units with respect to conventional treatment planning strategies, along with a 6 fold improvement in high dose spillage over the FB plan. The higher density ACT and Hybrid schemes also produced a slight enhancement in target conformity and reduction in low dose spillage. Conclusion: All treatment plans created in this study exceeded RTOG protocol criteria. Our results determine the free breathing approach yields an inaccurate account of the target treatment density. A significant decrease in unnecessary lung irradiation can be achieved by implementing Hybrid HU method with ACT method second best.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

MO-G-201-01: A Multi-Institutional Study Investigating the Performance of a Knowledge-Based Planning System Against Pinnacle Auto-Planning Engine in SIB-IMRT for the Head-And-Neck Cancer

Purpose: Knowledge-based Planning (KBP) founded on prior planning experience and Auto-Planning Engine (APE; commercialized in Pinnacle v9.10 TPS) based on progressive optimization algorithm both aim to eliminate the trial-and-error process in radiotherapy inverse planning. This study investigates the performance of the approaches in a multi-institutional setting to evaluate their functionalities in oropharyngeal cancer and offers suggestions how they can be implemented in the clinic. Methods: Radboud University Medical Center (RUMC) provided 35 oropharyngeal cancer patients (SIB-IMRT with two-dose-level prescription: 68 Gy to PTV68 and 50.3 Gy to PTV50.3) with corresponding comparative APE plans. Johns Hopkins University (JHU) contributed to a three-dose-level (70 Gy 63 Gy and 58.1 Gy) plan library for RUMC’s patient KBP generation. MedStar Georgetown University Hospital (MGUH) contributed to a KBP approach employing overlap-volume histogram (OVH-KBP) for generating RUMC’s patient KBP plans using JHU’s plan library. Since both approaches need their own user-defined parameters as initial inputs the first 10 patients were set aside as training set to finalize them. Meanwhile cross-institutional comparisons and adjustments were implemented for investigating institutions’ protocol discrepancies and the approaches’ user-defined parameters were updated accordingly. The finalized parameters were then applied to the remaining 25 patients for OVH-KBP and APE generation. A Wilcoxon rank-sum test was used for statistical comparison with significance level of p<0.05. Results: On average PTV68’s V95 was 96.5% in APE plans vs. 97% in OVH-KBP plans (p=0.36); PTV50.3’s V95 in APE plans was 97.8% vs.97.6% in OVH-KBP plans (p=0.6); cord’s D0.1 cc was 38.6 Gy in OVH-KBP plans vs. 43.7 Gy in APE plans (p=0.0001); mean doses to larynxes oral cavities parotids and submandibular glands were similar with p>0.2. Conclusions: The study demonstrates that KBP and APE can generate plans of comparable quality in a multi-institutional setting. Variations in clinical protocols can be effectively addressed for cross-institutional adaptations. Binbin Wu and Todd McNutt are the co-inventors of a patent associated with the proposed knowledge-based planning system which was licensed to Varian Medical Systems in 2015; This research was in part supported by Philips Radiation Oncology Systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-AB-BRA-03: Atlas-Based Algorithms with Local Registration-Goodness Weighting for MRI-Driven Electron Density Mapping

Purpose: To develop image-analysis algorithms to synthesize CT with accurate electron densities for MR-only radiotherapy of head & neck (H&N) and pelvis anatomies. Methods: CT and 3T-MRI (Philips, mDixon sequence) scans were randomly selected from a pool of H&N (n=11) and pelvis (n=12) anatomies to form an atlas. All MRIs were pre-processed to eliminate scanner and patient-induced intensity inhomogeneities and standardize their intensity histograms. CT and MRI for each patient were then co-registered to construct CT-MRI atlases. For more accurate CT-MR fusion, bone intensities in CT were suppressed to improve the similarity between CT and MRI. For a new patient, all CT-MRI atlases are deformed onto the new patients’ MRI initially. A newly-developed generalized registration error (GRE) metric was then calculated as a measure of local registration accuracy. The synthetic CT value at each point is a 1/GRE-weighted average of CTs from all CT-MR atlases. For evaluation, the mean absolute error (MAE) between the original and synthetic CT (generated in a leave-one-out scheme) was computed. The planning dose from the original and synthetic CT was also compared. Results: For H&N patients, MAE was 67±9, 114±22, and 116±9 HU over the entire-CT, air and bone regions, respectively. For pelvis anatomy, MAE was 47±5 and 146±14 for the entire and bone regions. In comparison with MIRADA medical, an FDA-approved registration tool, we found that our proposed registration strategy reduces MAE by ∼30% and ∼50% over the entire and bone regions, respectively. GRE-weighted strategy further lowers MAE by ∼15% to ∼40%. Our primary dose calculation also showed highly consistent results between the original and synthetic CT. Conclusion: We’ve developed a novel image-analysis technique to synthesize CT for H&N and pelvis anatomies. Our proposed image fusion strategy and GRE metric help generate more accurate synthetic CT using locally more similar atlases (Support: Philips Healthcare). The research is supported by Philips HealthCare.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-D-201-06: HDR Plan Prechecks Using Eclipse Scripting API

Purpose: Automate brachytherapy treatment plan quality check using Eclipse v13.6 scripting API based on pre-configured rules to minimize human error and maximize efficiency. Methods: The HDR Precheck system is developed based on a rules-driven approach using Eclipse scripting API. This system checks for critical plan parameters like channel length, first source position, source step size and channel mapping. The planned treatment time is verified independently based on analytical methods. For interstitial or SAVI APBI treatment plans, a Patterson-Parker system calculation is performed to verify the planned treatment time. For endobronchial treatments, an analytical formula from TG-59 is used. Acceptable tolerances were defined based on clinical experiences in our department. The system was designed to show PASS/FAIL status levels. Additional information, if necessary, is indicated appropriately in a separate comments field in the user interface. Results: The HDR Precheck system has been developed and tested to verify the treatment plan parameters that are routinely checked by the clinical physicist. The report also serves as a reminder or checklist for the planner to perform any additional critical checks such as applicator digitization or scenarios where the channel mapping was intentionally changed. It is expected to reduce the current manual plan check time from 15 minutes to <1 minute. Conclusion: Automating brachytherapy plan prechecks significantly reduces treatment plan precheck time and reduces human errors. When fully developed, this system will be able to perform TG-43 based second check of the treatment planning system’s dose calculation using random points in the target and critical structures. A histogram will be generated along with tabulated mean and standard deviation values for each structure. A knowledge database will also be developed for Brachyvision plans which will then be used for knowledge-based plan quality checks to further reduce treatment planning errors and increase confidence in the planned treatment.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

TU-D-207B-02: Delta-Radiomics: The Prognostic Value of Therapy-Induced Changes in Radiomics Features for Stage III Non-Small Cell Lung Cancer Patients

Purpose: To determine how radiomics features change during radiation therapy and whether those changes (delta-radiomics features) can improve prognostic models built with clinical factors. Methods: 62 radiomics features, including histogram, co-occurrence, run-length, gray-tone difference, and shape features, were calculated from pretreatment and weekly intra-treatment CTs for 107 stage III NSCLC patients (5–9 images per patient). Image preprocessing for each feature was determined using the set of pretreatment images: bit-depth resample and/or a smoothing filter were tested for their impact on volume-correlation and significance of each feature in univariate cox regression models to maximize their information content. Next, the optimized features were calculated from the intratreatment images and tested in linear mixed-effects models to determine which features changed significantly with dose-fraction. The slopes in these significant features were defined as delta-radiomics features. To test their prognostic potential multivariate cox regression models were fitted, first using only clinical features and then clinical+delta-radiomics features for overall-survival, local-recurrence, and distant-metastases. Leave-one-out cross validation was used for model-fitting and patient predictions. Concordance indices(c-index) and p-values for the log-rank test with patients stratified at the median were calculated. Results: Approximately one-half of the 62 optimized features required no preprocessing, one-fourth required smoothing, and one-fourth required smoothing and resampling. From these, 54 changed significantly during treatment. For overall-survival, the c-index improved from 0.52 for clinical factors alone to 0.62 for clinical+delta-radiomics features. For distant-metastases, the c-index improved from 0.53 to 0.58, while for local-recurrence it did not improve. Patient stratification significantly improved (p-value<0.05) for overallsurvival and distant-metastases when delta-radiomics features were included. The delta-radiomics versions of autocorrelation, kurtosis, and compactness were selected most frequently in leave-one-out iterations. Conclusion: Weekly changes in radiomics features can potentially be used to evaluate treatment response and predict patient outcomes. High-risk patients could be recommended for dose escalation or consolidation chemotherapy. This project was funded in part by grants from the National Cancer Institute (NCI) and the Cancer Prevention Research Institute of Texas (CPRIT).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-H-CAMPUS-JeP2-05: Can Automatic Delineation of Cardiac Substructures On Noncontrast CT Be Used for Cardiac Toxicity Analysis?

Purpose: To evaluate the feasibility of using an automatic segmentation tool to delineate cardiac substructures from computed tomography (CT) images for cardiac toxicity analysis for non-small cell lung cancer (NSCLC) patients after radiotherapy. Methods: A multi-atlas segmentation tool developed in-house was used to delineate eleven cardiac substructures including the whole heart, four heart chambers, and six greater vessels automatically from the averaged 4DCT planning images for 49 NSCLC patients. The automatic segmented contours were edited appropriately by two experienced radiation oncologists. The modified contours were compared with the auto-segmented contours using Dice similarity coefficient (DSC) and mean surface distance (MSD) to evaluate how much modification was needed. In addition, the dose volume histogram (DVH) of the modified contours were compared with that of the auto-segmented contours to evaluate the dosimetric difference between modified and auto-segmented contours. Results: Of the eleven structures, the averaged DSC values ranged from 0.73 ± 0.08 to 0.95 ± 0.04 and the averaged MSD values ranged from 1.3 ± 0.6 mm to 2.9 ± 5.1mm for the 49 patients. Overall, the modification is small. The pulmonary vein (PV) and the inferior vena cava required the most modifications. The V30 (volume receiving 30 Gy or above) for the whole heart and the mean dose to the whole heart and four heart chambers did not show statistically significant difference between modified and auto-segmented contours. The maximum dose to the greater vessels did not show statistically significant difference except for the PV. Conclusion: The automatic segmentation of the cardiac substructures did not require substantial modification. The dosimetric evaluation showed no statistically significant difference between auto-segmented and modified contours except for the PV, which suggests that auto-segmented contours for the cardiac dose response study are feasible in the clinical practice with a minor modification to the PV vessel.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-AB-BRA-05: PET-Guided Delivery Quality Evaluation of Yttrium-90 Microsphere Radioembolizaton for Hepatocellular Carcinoma Patients: The Optimal Sequence of Radioembolizaton and Chemoembolization Treatments

Purpose: Selective-internal-radiation-therapy (SIRT) and transarterial-chemoembolization (TACE) are commonly used for treatment of liver tumors. The use of TACE, which is macroembolic, prior to SIRT may cause hemodynamic changes in tumor vasculature that impair yttrium-90 (90Y) microsphere delivery to the targeted lesions. This work aims to quantify dosimetric tumor coverage using 90Y positron emission tomography (PET) dosimetry after SIRT alone compared to TACE followed by SIRT. Methods: A total of 40 consecutive hepatocellular carcinoma (HCC) SIRT patients who had a post-SIRT 90Y PET/CT scan were evaluated. The patient-specific-3D-dose was reconstructed from the PET images. Patients were categorized into two groups: patients received TACE prior SIRT procedure (n=18) and patient received SIRT alone (n=22). The lesions and liver were delineated by a senior radiologist. We evaluated both the lesion-specific dose-volume-histogram (DVH) and the selectivity index (SI) defined as the ratio of the average dose inside the total lesion(s) and the average dose of the normal liver. The SI values of patients were compared based on whether TACE was previously used. Results: A wide spectrum was observed in the lesion-specific DVH-evaluation and SI appeared to be suitable of evaluating the quality of each SIRT infusion. The average SI of the entire patient group was 3.0, i.e. targeted lesion receiving three times higher dose than normal liver. The average SI was 1.8 for patients who had prior TACE and 3.9 for patients who did not have prior TACE (p=0.008). 85% of the patients with prior TACE demonstrated poor 90Y-microsphere delivery (SI <2) while none demonstrated excellent delivery (SI >4). On the other hand, the incidence SI >4 among patients with no prior TACE was 37%. Conclusion: 3D dose evaluation using post-SIRT PET suggests that 90Y microsphere delivery to liver tumors is impaired among patients who received prior TACE compared to those who receive SIRT alone.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-AB-202-03: Quantifying Ventilation Change Due to Radiation Therapy Using 4DCT Jacobian Calculations

Purpose: Four-dimensional computed tomography (4DCT) and image registration can be used to determine regional lung ventilation changes after radiation therapy (RT). This study aimed to determine if lung ventilation change following radiation therapy was affected by the pre-RT ventilation of the lung. Methods: 13 subjects had three 4DCT scans: two repeat scans acquired before RT and one three months after RT. Regional ventilation was computed using Jacobian determinant calculations on the registered 4DCT images. The post-RT ventilation map was divided by the pre-RT ventilation map to get a voxel-by-voxel Jacobian ratio map depicting ventilation change over the course of RT. Jacobian ratio change was compared over the range of delivered doses. The first pre-RT ventilation image was divided by the second to establish a control for Jacobian ratio change without radiation delivered. The functional change between scans was assessed using histograms of the Jacobian ratios. Results: There were significantly (p < 0.05) more voxels that had a large decrease in Jacobian ratio in the post-RT divided by pre-RT map (15.6%) than the control (13.2%). There were also significantly (p < .01) more voxels that had a large increase in Jacobian ratio (16.2%) when compared to control (13.3%). Lung regions with low function (<10% expansion by Jacobian) showed a slight linear reduction in expansion (0.2%/10 Gy delivered), while high function regions (>10% expansion) showed a greater response (1.2% reduction/10 Gy). Contiguous high function regions > 1 liter occurred in 11 of 13 subjects. Conclusion: There is a significant change in regional ventilation following a course of radiation therapy. The change in Jacobian following RT is dependent both on the delivered dose and the initial ventilation of the lung tissue: high functioning lung has greater ventilation loss for equivalent radiation doses. Substantial regions of high function lung tissue are prevalent. Research support from NIH grants CA166119 and CA166703, a gift from Roger Koch, and a Pilot Grant from University of Iowa Carver College of Medicine.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-AB-207B-01: Dose Tolerance for SBRT/SABR

Purpose: Stereotactic body radiation therapy (SBRT) / stereotactic ablative body radiotherapy (SABR) is gaining popularity, but quantitative dose tolerance has still been lacking. To improve this, the April 2016 issue of Seminars in Radiation Oncology will have normal tissue complication probability (NTCP) models for 10 critical structures: optic pathway, cochlea, oral mucosa, esophagus, chestwall, aorta, bronchi, duodenum, small bowel, and spinal cord. Methods: The project included more than 1500 treatments in 1–5 fractions using CyberKnife, Gamma Knife, or LINAC, with 60 authors from 15 institutions. NTCP models were constructed from the 97 grade 2–3 complications, predominantly scored using the common terminology criteria for adverse events (CTCAEv4). Dose volume histogram (DVH) data from each institutional dataset was loaded into the DVH Evaluator software (DiversiLabs, LLC, Huntingdon Valley, Pa) for modeling. The current state of the literature for the critical structures was depicted using DVH Risk Maps: comparative graphs of dose tolerance limits that can include estimated risk levels, reported complications, DVH data for study patients, as well as high- and low-risk dose tolerance limits. Results: For relatively acceptable toxicity like grade 1–3 rib fractures and chestwall pain, the high-risk limits have 50% risk and the low-risk limits have 5% risk. Emami et al (IJROBP 1991 May 15;21(1):109–22) used 50% and 5% risk levels for all structures, whereas this effort used clinically acceptable ranges for each: in structures like aorta or spinal cord where complications must be avoided, the high- and low-risk limits have about 3% and 1% risk, respectively, in this issue of Seminars. These statistically based guidelines can help ensure plan quality for each patient. Conclusion: NTCP for SBRT is now becoming available. Hypofractionated dose tolerance can be dramatically different than extrapolations of conventional fractionation so NTCP analysis of the SBRT/SBRT data is important to ensure safe clinical practice. Dr. Grimm, designed and holds intellectual property rights to the DVH Evaluator software tool which is an FDA-cleared product in commercial use, and was used to analyze the data.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-AB-207B-06: Dose and Biological Uncertainties in Sarcoma

Purpose: To understand the clinical impact of key uncertainties in proton therapy potentially affecting the analysis of clinical trials, namely the assumption of using a constant relative biological effectiveness (RBE) of 1.1 compared to variable RBE for proton therapy and the use of analytical dose calculation (ADC) methods. Methods: Proton dose distributions were compared for analytical and Monte Carlo (TOPAS) dose calculations. In addition, differences between using a constant RBE of 1.1 (RBE-constant) were compared with four different RBE models (to assess model variations). 10 patients were selected from an ongoing clinical trial on IMRT versus scanned protons for sarcoma. Comparisons were performed using dosimetric indices based on dose-volume histogram analyses and γ-index analyses. Results: For three of the RBE-models the mean dose, D95, D50 and D02 (dose values covering 95%, 50% and 2% of the target volume, respectively) were up to 5% lower than for RBE-constant. The dosimetric indices for one of the RBE-models were around 9% lower than for the RBE-constant model. The differences for V90 (the percentage of the target volume covered by 90% of the prescription dose) were up to 40% for three RBE-models, whereas for one the difference was around 95%. All ADC dosimetric indices were up to 5% larger than for RBE-constant. The γ-index passing rate for the target volume with a 3%/3mm criterion was above 97% for all models except for one, which was below 24%. Conclusion: Interpretation of clinical trials on sarcoma may depend on dose calculation uncertainties (as assessed by Monte Carlo). In addition, the biological dose distribution depends notably on which RBE model is utilized. The current practice of using a constant RBE of 1.1 may overestimate the target dose by as much as 5% for biological dose calculations. Performing an RBE uncertainty analysis is recommended for trial analysis. U19 projects - U19 CA 021239. PI: Delaney.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-DE-201-01: BEST IN PHYSICS (THERAPY): A Fast Multi-Target Inverse Treatment Planning Strategy Optimizing Dosimetric Measures for High-Dose-Rate (HDR) Brachytherapy

Purpose: Inverse treatment planning (ITP) for interstitial HDR brachytherapy of gynecologic cancers seeks to maximize coverage of the clinical target volumes (tumor and vagina) while respecting dose-volume-histogram related dosimetric measures (DMs) for organs at risk (OARs). Commercially available ITP tools do not support DM-based planning because it is computationally too expensive to solve. In this study we present a novel approach that allows fast ITP for gynecologic cancers based on DMs for the first time. Methods: This novel strategy is an optimization model based on a smooth DM-based objective function. The smooth approximation is achieved by utilizing a logistic function for the evaluation of DMs. The resulting nonconvex and constrained optimization problem is then optimized with a BFGS algorithm. The model was evaluated using the implant geometry extracted from 20 patient treatment plans under an IRB-approved retrospective study. For each plan, the final DMs were evaluated and compared to the original clinical plans. The CTVs were the contoured tumor volume and the contoured surface of the vagina. Statistical significance was evaluated with a one-sided paired Wilcoxon signed-rank test. Results: As did the clinical plans, all generated plans fulfilled the defined DMs for OARs. The proposed strategy showed a statistically significant improvement (p<0.001) in coverage of the tumor and vagina, with absolute improvements of related DMs of (6.9 +/− 7.9)% and (28.2 +/− 12.0)%, respectively. This was achieved with a statistically significant (p<0.01) decrease of the high-dose-related DM for the tumor. The runtime of the optimization was (2.3 +/− 2.0) seconds. Conclusion: We demonstrated using clinical data that our novel approach allows rapid DM-based optimization with improved coverage of CTVs with fewer hot spots. Being up to three orders of magnitude faster than the current clinical practice, the method dramatically shortens planning time.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-FG-207B-03: Multi-Energy CT Reconstruction with Spatial Spectral Nonlocal Means Regularization

Purpose: Multi-energy computed tomography (MECT) is an emerging application in medical imaging due to its ability of material differentiation and potential for molecular imaging. In MECT, image correlations at different spatial and channels exist. It is desirable to incorporate these correlations in reconstruction to improve image quality. For this purpose, this study proposes a MECT reconstruction technique that employes spatial spectral non-local means (ssNLM) regularization. Methods: We consider a kVp-switching scanning method in which source energy is rapidly switched during data acquisition. For each energy channel, this yields projection data acquired at a number of angles, whereas projection angles among channels are different. We formulate the reconstruction task as an optimziation problem. A least square term enfores data fidelity. A ssNLM term is used as regularization to encourage similarities among image patches at different spatial locations and channels. When comparing image patches at different channels, intensity difference were corrected by a transformation estimated via histogram equalization during the reconstruction process. Results: We tested our method in a simulation study with a NCAT phantom and an experimental study with a Gammex phantom. For comparison purpose, we also performed reconstructions using conjugate-gradient least square (CGLS) method and conventional NLM method that only considers spatial correlation in an image. ssNLM is able to better suppress streak artifacts. The streaks are along different projection directions in images at different channels. ssNLM discourages this dissimilarity and hence removes them. True image structures are preserved in this process. Measurements in regions of interests yield 1.1 to 3.2 and 1.5 to 1.8 times higher contrast to noise ratio than the NLM approach. Improvements over CGLS is even more profound due to lack of regularization in the CGLS method and hence amplified noise. Conclusion: The proposed ssNLM method for kVp-switching MECT reconstruction can achieve high quality MECT images.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TH-AB-201-11: Small-Field Beam Quality Corrections and Their Clinical Impact for Flattened and Unflattened Fields of the TrueBeam

Purpose: To determine small-field beam quality correction factors for several detectors used with flattened and unflattened beams of the Varian TrueBeam accelerator, and to quantify the clinical impact of these corrections through a treatment planning study. Methods: Small-field beam quality correction factors were measured and simulated in EGSnrc for several microchambers and diodes. The Exradin W1 plastic scintillation detector was used as a reference for the measured results, and its response was validated using alanine measurements traceable to a primary standard. Symmetric, jaw-defined field sizes ranging from 6 mm to 10 cm were used. Twelve hypofractionated clinical cases were then recalculated and/or reoptimized using beam models that were modified based on the resulting output factors. Three-dimensional dose distributions and cumulative dose-volume-histogram plots were compared between the models. The impact of these measurements on the flattening-filter-scatter-source (FFSS) model, and its subsequent effect on treatment plans, was also investigated. Results: The diodes exhibited over-response in small fields, whereas the microchambers tended to under-respond as field size decreased. Field-size-dependent effects were less significant for the flattening-filter-free beams (i.e., corrections were closer to unity) than for flattened fields of the same size. Recalculated and reoptimized plans displayed varying levels of deviation from the original dose distributions, with larger discrepancies at effective jaw settings of 3×3 cm{sup 2} or smaller. Dose distributions calculated using uncorrected microchamber output factors and/or modified FFSS showed large discrepancies in the low-dose regions, which manifested as DVH changes for the organs at risk (OAR). Conclusion: Field-size-dependent beam quality corrections were measured and simulated for ionization chambers and diodes typically used in small fields. Plans calculated with models based on these measured data showed dosimetric deviations mainly out of field, which affected OAR dose metrics. Accordingly, these beam quality corrections should be accounted for during commissioning to ensure greater treatment planning accuracy.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

TH-AB-BRA-03: Backscatter Dose Factors Re-Evaluated for Inhomogeneities in the Presence of a 1.5 T Magnetic Field Using the GPUMCD Monte Carlo Algorithm

Purpose: To quantify the backscatter dose factors near the interfaces for clinically relevant high atomic number materials using GPUMCD for the Elekta MRI Linac. Methods: Backscatter dose factors (BSDF) were calculated as the ratio of the dose with and without the presence of the heterogeneity. The BSDF’s were calculated either in the absence or presence of an orthogonal 1.5 T magnetic field. Doses were scored in small voxels of side 1 mm in a water phantom with dimensions of 20×20×20 cm using GPUMCD (Elekta). The minimum uncertainty in dose calculations was kept to 0.5%. A slab of thickness 2 cm, representing the inhomogeneity, was placed inside the phantom with variable position from the surface of the phantom. The slab was filled with either bone, aluminum, titanium, stainless steel, or dental amalgam. The phantom was irradiated using particles sampled from a histogram which represented the true MRI Linac spectrum. Results: With the application of the 1.5 T magnetic field (B-On), all of the BSDF’s were reduced by at least 8% compared to the no magnetic field (B-Off) cases. For the B-Off case, the BSDF decreases exponentially with the upstream distance away from the interface. With B-On, the BSDF decreases exponentially for titanium, SS, and amalgam. However, it remains constant for Aluminum. In the case of bone, the BSDF increases up to a distance of 4 mm away from the interface in the presence of the magnetic field. Conclusion: The BSDF does not depend upon the thickness of the homogeneous material above the inhomogeneity for either the B-Off or B-On cases. For all the materials investigated, the BSDF is lower at the interface for the B-On case. The exponential fall-off of the BSDF away from the interface is not valid for all the materials when the magnetic field is turned ON. Funding support for this research was provided by Elekta.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TH-AB-BRA-10: The Physics of Interface Effects for Radiation Treatments in a MRI-Linac: A Monte Carlo Study

Purpose: To investigate and explain the interface effects for clinically relevant materials being irradiated in the presence of a 1.5 T transverse magnetic field. Methods: Interface effects were investigated using Geant4.10.1 both with (B-On) and without (B-Off) a magnetic field for an Elekta MRI-Linac. A slab of thickness 8 cm, representing inhomogeneity, was placed at a depth of 4 cm in a 20×20×20 cm water phantom. Backscattered electron fluence was calculated through a 20×20 cm plane aligned with the surface of the inhomogeneity. Inhomogeneities investigated were lung, bone, aluminum, titanium, stainless steel, and dental filling. A photon beam with field size of 2×2 cm at the isocenter and SAD of 143.5 cm was generated from a point source with energy distribution sampled from a histogram representing the true Elekta MRI-Linac photon spectrum. Results: In the B-On case, if the heterogeneity is a low Z{sub eff} material, such as lung, the backscattered electron fluence is increased considerably, i.e. by 54 %, and the corresponding dose is expected to be higher near the interface compared to the B-Off case. On the contrary, if the heterogeneity is a high Z{sub eff} material then the backscattered electron fluence is reduced in the B-On electron fluence is reduced in the B-On case. This reduction leads to a lower dose deposition at the interface compared to the B-Off case. Conclusion: The reduction in dose at the interface, in the B-On case, is directly related to the reduction in backscattered electron fluence. The reduction in backscattered electron fluence occurs due to two different reasons. First, the electron energy spectrum hitting the interface is changed for the B-On case which changes the electron scattering probability. Second, some electrons that are looping under the influence of the magnetic field are captured by the higher density side of the interface and no longer contribute to the backscattered electron stream. Funding support for this study was provided by ElektaTM.

60 APPLIED LIFE SCIENCES↗