Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “regression”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Utilizing machine learning to improve the precision of fluorescence imaging of cavity-generated spin squeezed states

We present a supervised learning model to calibrate the photon collection rate during the fluorescence imaging of cold atoms. The linear regression model finds the collection rate at each location on the sensor such that the atomic population difference equals that of a highly precise optical cavity measurement. This 192 variable regression results in a measurement variance 27% smaller than our previous single variable regression calibration. The measurement variance is now in agreement with the theoretical limit due to other known noise sources. This model efficiently trains in less than a minute on a standard personal computer's CPU and requires less than 10 min of data collection. Furthermore, the model is applicable across a large change in population difference and across data collected on different days.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Empirical radius formulas for canonical neutron stars from bidirectionally selecting features of equations of state in extended Bayesian analyses of observational data

Significant advancement in Bayesian inference of nuclear equation of state (EOS) from gravitational wave and x-ray observations of neutron stars (NSs) has been made by the nuclear astrophysics community especially since GW170817. By extending the traditional Bayesian analysis which normally ends at presenting the marginalized posterior probability distribution functions (PDFs) of individual EOS parameters and their correlations (or sometimes only the Pearson correlation coefficients which are only reliably useful when the variables are linearly correlated while they are actually often not), we search for a data-driven and robust empirical formula for the radius 𝑅 1.4 of canonical NSs in terms of the characteristic EOS parameters (features). We also identify the single most important but currently poorly known EOS parameter for determining the 𝑅 1.4 . Using three regression-model-building methodologies: bidirectional stepwise feature selection, least absolute shrinkage selection operator (LASSO) regression, and neural network regression on a large set of posterior EOSs and the corresponding 𝑅 1.4 values inferred from earlier comprehensive Bayesian analyses of NS observational data, we systematically and rigorously develop the most probable 𝑅 1.4 formulas with varying statistical accuracy and technical complexity. Here, the most important EOS parameters for determining 𝑅 1.4 are found consistently in each of the feature selection processes to be (in order of decreasing importance): curvature 𝐾 sym , slope 𝐿, skewness 𝐽 sym of nuclear symmetry energy, skewness 𝐽 0 , incompressibility 𝐾 0 of symmetric nuclear matter, and the magnitude 𝐸 sym ⁡(𝜌 0 ) of symmetry energy at the saturation density 𝜌 0 of nuclear matter.

Bayesian methods↗
DOE OSTI2026

A Machine Learning based Approach of Estimating Equivalent Circuit Model Parameters at Different SoCs of Li-ion Batteries from Voltage Relaxation

Abstract: In this study, an approach of estimating the equivalent circuit model (ECM) parameters for Li-ion batteries (LIBs) is proposed based on the voltage value at different intervals while relaxing the LIB after discharge. The typical approach for estimating ECM parameters of a LIB is to conduct electrochemical impedance spectroscopy (EIS) measurements at different frequencies and fit them to a predefined circuit model, which requires additional measuring arrangements and specialized devices. The proposed methodology utilizes four different voltages at 0s, 60s, 360s, and 1800s alongside the specific state of charge (SoC) value for a specific constant discharge current value of ~1C until the relaxation stage to train and evaluate three regression-based machine learning models— Support Vector Regression (SVR), Extreme Gradient Boosting (XGBoost), and Gaussian Process Regression (GPR)—for estimating the ECM parameters of the selected model. Bayesian optimization is employed for hyperparameter tuning to achieve optimal performance for all the regressor models, among which, the GPR provided the best performance with the root-mean-squared error (RMSE) of less than 4x10-4 on average for the resistive components and less than 0.27 for capacitive components with excellent R2 scores. The simplicity of the approach enables it to eliminate the need for sophisticated measuring equipment and computation power.

Sagar, Md. Samiul [The University of Alabama (UA)]↗
DOE OSTI2018

Temporal ensemble learning of univariate methods for short term load forecasting

Short term load forecasting (STLF) is a fundamental component of demand response programs in smart grids. Recent literature has focused on complex neural networks and deep learning models for solving STLF. While these models work well for load forecasting with complex non-linear relationships, they have been shown to be less effective than simpler univariate models for STLF problems with under a 6 hours horizon. Given the lack of multivariate data (such as temperature) in many practical datasets, we need better univariate prediction models for STLF. By partitioning the dataset by temporal features, we develop a novel ensemble learning method that exploits daily seasonality in electricity consumption data to improve accuracy of popular univariate models. We train an ensemble of models from the dataset partitions. We develop a variety of methods, including Ridge Regression, to increase the robustness of the ensemble prediction. To show the effectiveness of our approach, we perform detailed evaluation using an aggregated user electricity consumption dataset collected by the Los Angeles Department of Water and Power (LADWP). We select four popular prediction algorithms in literature for our experiments, including Kernel Regression (KR), Support Vector Regression (SVR) and neural network approaches. We compare the performance of these algorithms applying our ensemble approach to training only one single model. Our approach leads to an 11.2% decrease in mean absolute percentage error (MAPE) and 21.3% decrease in root mean squared error (RMSE) over the single model approach for KR, and a 30% and 32.4% decrease in MAPE and RMSE respectively for SVR. These ensemble models also outperform the neural network approaches.

14 SOLAR ENERGY↗
DOE OSTI2016

PV Degradation Methodology Comparison – A Basis for a Standard

What is the best method to determine long-term PV system performance and degradation rates? Ideally, one universally applicable methodology would be desirable so that a single number could be derived. However, data sets vary in their attributes and evidence is presented that defining two methodologies may be preferable. Monte Carlo simulations of artificial performance data allowed investigation of different methodologies and their respective confidence intervals. Tradeoffs between different approaches were delineated, elucidating as to why two separate approaches may need to be included in a standard. Regression approaches tend to be preferable when data sets are less contaminated by seasonality, noise and occurrence of outliers although robust regression can significantly improve the accuracy when outliers are present. In the presence of outliers, marked seasonality, or strong soiling events, year-on-year approaches tend to outperform regression approaches.

14 SOLAR ENERGY↗
DOE OSTI2021

Distribution Grid Modeling Using Smart Meter Data

The knowledge of distribution grid models, including topologies and line impedances, is essential for grid monitoring, control and protection. However, such information is often unavailable, incomplete or outdated. The increasing deployment of smart meters (SMs) provides a unique opportunity to tackle this issue. This paper proposes a two-stage framework for distribution grid modeling using SM data. In the first stage, the network topology is identified by reconstructing a weighted Laplacian matrix of distribution networks. In the second stage, a least absolute deviations (LAD) regression model is developed for estimating line impedance of a single branch based on the nonlinear (inverse) power flow model, wherein a conductor library is leveraged to narrow down the solution space. The LAD regression model is originally a mixed-integer nonlinear program whose continuous relaxation is still non-convex. Furthermore, we specially address its convex relaxation and discuss the exactness. The modified regression model is then embedded within a bottom-up sweep algorithm to achieve the identification across the network in a branch-wise manner. Numerical results on the IEEE 13-bus, 37-bus and 69-bus test feeders validate the effectiveness of the proposed methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Method for estimating the density of high-level nuclear waste glass

A database of over 1100 silicate glass compositions and densities was compiled and used to evaluate the efficacy of an algorithm for estimating the density of silicate glass compositions. We sought to develop a parsimonious algorithm based on the additivity of partial molar volumes of individual oxide components weighted by their mole fraction in a glass composition. Bound molar volumes were used for oxides in which the density of the oxide bound in a glass matrix was previously determined. The bound molar volumes were known for oxides covering 97.5 mole percent of the database compositional space. The measured glass densities were plotted against the estimated glass densities and a linear regression yielded an R 2 adj. = 0.95 and a slope and intercept of approximately one and zero, respectively. This regression suggests that glass densities estimated by the algorithm, within analysis uncertainty, are equal to the measured densities of the glasses. In addition to the development of the density estimation, we corroborated many of the referenced bound molar volume data used in the parameterization of the estimation algorithm via linear regression of the individual partial molar volumes versus the inverse measured densities (specific volumes) of the glasses in the database.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Pregnancy outcome after first trimester exposure to domperidone—An observational cohort study

Abstract Aim To assess the teratogenic risk of domperidone by comparing the incidence of major malformation with domperidone to a control. Methods Pregnancy outcome data were obtained for women at two Japanese facilities that provide counseling on drug use during pregnancy between April 1988 and December 2017. The incidence of major malformation was calculated among infants born to women taking domperidone ( n = 519), nonteratogenic drugs (control, n = 1673), or metoclopramide (reference, n = 241) during the first trimester of pregnancy. Using the control group as reference, the crude odds ratio (OR) of the incidence of major malformation in the domperidone and metoclopramide groups was calculated using univariable logistic regression analysis. Adjusted OR was also calculated using multivariable logistic regression analysis adjusted for various other factors. Results The incidence of major malformation was 2.9% (14/485, 95% confidence interval [CI]: 1.6–4.8) in the domperidone group, 1.7% (27/1554, 95%CI: 1.1–2.5) in the control group, and 3.6% (8/224, 95%CI: 1.6–6.9) in the metoclopramide group. The adjusted multivariable logistic regression analysis showed no significant difference in incidence between the control and domperidone groups (adjusted OR: 1.86 [95%CI: 0.73–4.70], p = 0.191) or between the control and metoclopramide groups (adjusted OR: 2.20 [95%CI: 0.69–6.98], p = 0.183). Conclusions This observational cohort study showed that domperidone exposure during the first trimester was not associated with increased risk of major malformation in infants. These results may help alleviate the anxiety of patients who took domperidone during pregnancy.

Hishinuma, Kayoko↗
DOE OSTI2014

Quantitative ultrasound characterization of locally advanced breast cancer by estimation of its scatterer properties

Purpose: Tumor grading is an important part of breast cancer diagnosis and currently requires biopsy as its standard. Here, the authors investigate quantitative ultrasound parameters in locally advanced breast cancers that can potentially separate tumors from normal breast tissue and differentiate tumor grades. Methods: Ultrasound images and radiofrequency data from 42 locally advanced breast cancer patients were acquired and analyzed. Parameters related to the linear regression of the power spectrum—midband fit, slope, and 0-MHz-intercept—were determined from breast tumors and normal breast tissues. Mean scatterer spacing was estimated from the spectral autocorrelation, and the effective scatterer diameter and effective acoustic concentration were estimated from the Gaussian form factor. Parametric maps of each quantitative ultrasound parameter were constructed from the gated radiofrequency segments in tumor and normal tissue regions of interest. In addition to the mean values of the parametric maps, higher order statistical features, computed from gray-level co-occurrence matrices were also determined and used for characterization. Finally, linear and quadratic discriminant analyses were performed using combinations of quantitative ultrasound parameters to classify breast tissues. Results: Quantitative ultrasound parameters were found to be statistically different between tumor and normal tissue (p < 0.05). The combination of effective acoustic concentration and mean scatterer spacing could separate tumor from normal tissue with 82% accuracy, while the addition of effective scatterer diameter to the combination did not provide significant improvement (83% accuracy). Furthermore, the two advanced parameters, including effective scatterer diameter and mean scatterer spacing, were found to be statistically differentiating among grade I, II, and III tumors (p = 0.014 for scatterer spacing, p = 0.035 for effective scatterer diameter). The separation of the tumor grades further improved when the textural features of the effective scatterer diameter parametric map were combined with the mean value of the map (p = 0.004). Conclusions: Overall, the binary classification results (tumor versus normal tissue) were more promising than tumor grade assessment. Combinations of advanced parameters can further improve the separation of tumors from normal tissue compared to the use of linear regression parameters. While the linear regression parameters were sufficient for characterizing breast tumors and normal breast tissues, advanced parameters and their textural features were required to better characterize tumor subtypes.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-E-J-03: A Comprehensive Comparison Between Alpha and Beta Emitters for Cancer Radioimmunotherapy

Purpose: The purpose of this study is to perform a comprehensive comparison of the therapeutic efficacy and cytotoxicity of alpha and beta emitters for Radioimmunotherapy (RIT). For each stage of cancer development, specific models were built for the separate objectives of RIT to be addressed:a) kill isolated cancer cells in transit in the lymphatic and vascular circulation,b) regress avascular cell clusters,c) regress tumor vasculature and tumors. Methods: Because of the nature of short range, high LET alpha and long energy beta radiation and heterogeneous antigen expression among cancer cells, the microdosimetric approach is essential for the RIT assessment. Geant4 based microdosimetric models are developed for the three different stages of cancer progression: cancer cells, cell clusters and tumors. The energy deposition, specific energy resulted from different source distribution in the three models was calculated separately for 4 alpha emitting radioisotopes ({sup 211}At, {sup 213}Bi, {sup 223}Ra and {sup 225}Ac) and 6 beta emitters ({sup 32}P, {sup 33}P, {sup 67}Cu, {sup 90}Y, {sup 131}I and {sup 177}Lu). The cell survival, therapeutic efficacy and cytotoxicity are determined and compared between alpha and beta emitters. Results: We show that internal targeted alpha radiation has advantages over beta radiation for killing isolated cancer cells, regressing small cell clusters and also solid tumors. Alpha particles have much higher dose specificity and potency than beta particles. They can deposit 3 logs more dose than beta emitters to single cells and solid tumor. Tumor control probability relies on deep penetration of radioisotopes to cancer cell clusters and solid tumors. Conclusion: The results of this study provide a quantitative understanding of the efficacy and cytotoxicity of RIT for each stage of cancer development.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-342: Use of Patient Geometry Measurements to Predict Dosimetric Gain with VMAT Over 3D for Chestwall and Regional Nodal Radiation

Purpose: To predict the dosimetric gain of VMAT over 3D for the treatment ofchestwall/IMN/supraclavicular nodes using geometric parameters acquired during simulation Methods: CT scans for 20 left and 20 right sided patients were retrospectively analyzed toobtain percent ipsilateral lung volume included in the PWT and supraclavicular fields, central lung depth (CLD), maximum lung depth (MLD), separation, chestwall concavity (defined here as the product of CLD and separation) and the maximum heart depth (MHD). VMAT, PWT and P/E plans were done for each case. The ipsilateral lung V20 Gy and mean, total lung V20 Gy and mean, heart V25 Gy and mean were noted for each plan. Correlation coefficients were obtained and linear regression models were built using data from the above training set of patients and then tested on 4 new patients. Results: The decrease in ipsilateral lung V20 Gy, total lung V20 Gy, ipsilateral lung mean and total lung mean with VMAT over PWT significantly (p<0.05) correlated with the percent volume of ipsilateral lung included in the PWT and supraclavicular fields with correlation coefficient values of r = 0.83, r = 0.77, r = 0.78 and r = 0.75 respectively. Significant correlations were also found between MHD and the decrease in heart V25 Gy and mean of r = 0.77 and r = 0.67 respectively. Dosimetric improvement with VMAT over P/E plans showed no correlation to any of the geometric parameters investigated in this study. The dosimetric gain predicted for the 4 test cases by the linear regression models given their respective percent ipsilateral lung volumes fell within the 95% confidence intervals around the best regression fit. Conclusion: The percent ipsilateral lung volume appears to be a strong predictor of the dosimetric gain on using VMAT over PWT apriori.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

TU-C-12A-02: Development of a Multiparametric Statistical Response Map for Quantitative Imaging

Purpose: Quantitative imaging biomarkers (QIB) are becoming increasingly utilized in early phase clinical trials as a means of non-invasively assessing treatment response and associated response heterogeneity. The aim of this study was to develop a flexible multiparametric statistical framework to predict voxel-by-voxel response of several potential MRI QIBs. Methods: Patients with histologically proven glioblastomas (n=11) were treated with chemoradiation (with/without bevacizumab) and underwent one baseline and two mid-treatment (3–4wks) MRIs. Dynamic contrast-enhanced (3D FSPGR, 6.3sec/phase, 0.1 mmol/kg Gd-DTPA), dynamic susceptibility contrast (2D GRE-EPI, 1.5sec/phase, 0.2mmol/kg Gd-DTPA), and diffusion tensor (2D DW-EPI, b=0, 1200 s/mm{sup 2}, 27 directions) imaging acquisitions were obtained during each study. Mid-treatment and pre-treatment images were rigidly aligned, and regions of partial response (PR), stable disease (SD), and progressive disease (PD) were contoured in consensus by two experienced radiation oncologists. Voxels in these categories were used to train ordinal (PR<SD<PD) and logistic (PR/SD<PD) models incorporating the apparent diffusion coefficient, fractional anisotropy, Ktrans, ve, vp, and relative cerebral blood volume imaging biomarkers, as well as treatment type. Leave-one-out cross-validation was performed at the patient level to assess model prediction accuracy. Results: Ordinal regression resulted in model prediction accuracies of 60% (PR), 0% (SD), 81% (PD), and 69% (overall), with coefficients of variation (COV) of 9.4%, 9.6%, and 23.6%, respectively. Logistic regression resulted in accuracies of 82.0% (PR/SD), 46.2% (PD), and 76.2% (overall) with COVs of 22.4%, 45.7%, and 23.8%, respectively. Conclusion: Despite limited patient numbers, this feasibility pilot study demonstrates that ordinal and logistic regression models potentially provide a flexible statistical framework for incorporating longitudinal multiparametric quantitative imaging data and, for these particular data, can potentially distinguish PR from PD. Additional training data exhibiting a full range of treatment responses is expected to allow prediction of all three treatment outcomes. Moreover, these models are extensible to multi-modality, multiparametric analyses.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

A GPU based high-resolution multilevel biomechanical head and neck model for validating deformable image registration

Purpose: Validating the usage of deformable image registration (DIR) for daily patient positioning is critical for adaptive radiotherapy (RT) applications pertaining to head and neck (HN) radiotherapy. The authors present a methodology for generating biomechanically realistic ground-truth data for validating DIR algorithms for HN anatomy by (a) developing a high-resolution deformable biomechanical HN model from a planning CT, (b) simulating deformations for a range of interfraction posture changes and physiological regression, and (c) generating subsequent CT images representing the deformed anatomy. Methods: The biomechanical model was developed using HN kVCT datasets and the corresponding structure contours. The voxels inside a given 3D contour boundary were clustered using a graphics processing unit (GPU) based algorithm that accounted for inconsistencies and gaps in the boundary to form a volumetric structure. While the bony anatomy was modeled as rigid body, the muscle and soft tissue structures were modeled as mass–spring-damper models with elastic material properties that corresponded to the underlying contoured anatomies. Within a given muscle structure, the voxels were classified using a uniform grid and a normalized mass was assigned to each voxel based on its Hounsfield number. The soft tissue deformation for a given skeletal actuation was performed using an implicit Euler integration with each iteration split into two substeps: one for the muscle structures and the other for the remaining soft tissues. Posture changes were simulated by articulating the skeletal structure and enabling the soft structures to deform accordingly. Physiological changes representing tumor regression were simulated by reducing the target volume and enabling the surrounding soft structures to deform accordingly. Finally, the authors also discuss a new approach to generate kVCT images representing the deformed anatomy that accounts for gaps and antialiasing artifacts that may be caused by the biomechanical deformation process. Accuracy and stability of the model response were validated using ground-truth simulations representing soft tissue behavior under local and global deformations. Numerical accuracy of the HN deformations was analyzed by applying nonrigid skeletal transformations acquired from interfraction kVCT images to the model’s skeletal structures and comparing the subsequent soft tissue deformations of the model with the clinical anatomy. Results: The GPU based framework enabled the model deformation to be performed at 60 frames/s, facilitating simulations of posture changes and physiological regressions at interactive speeds. The soft tissue response was accurate with a R{sup 2} value of >0.98 when compared to ground-truth global and local force deformation analysis. The deformation of the HN anatomy by the model agreed with the clinically observed deformations with an average correlation coefficient of 0.956. For a clinically relevant range of posture and physiological changes, the model deformations stabilized with an uncertainty of less than 0.01 mm. Conclusions: Documenting dose delivery for HN radiotherapy is essential accounting for posture and physiological changes. The biomechanical model discussed in this paper was able to deform in real-time, allowing interactive simulations and visualization of such changes. The model would allow patient specific validations of the DIR method and has the potential to be a significant aid in adaptive radiotherapy techniques.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Cyst-based measurements for assessing lymphangioleiomyomatosis in computed tomography

Purpose: To investigate the efficacy of a new family of measurements made on individual pulmonary cysts extracted from computed tomography (CT) for assessing the severity of lymphangioleiomyomatosis (LAM). Methods: CT images were analyzed using thresholding to identify a cystic region of interest from chest CT of LAM patients. Individual cysts were then extracted from the cystic region by the watershed algorithm, which separates individual cysts based on subtle edges within the cystic regions. A family of measurements were then computed, which quantify the amount, distribution, and boundary appearance of the cysts. Sequential floating feature selection was used to select a small subset of features for quantification of the severity of LAM. Adjusted R{sup 2} from multiple linear regression and R{sup 2} from linear regression against measurements from spirometry were used to compare the performance of our proposed measurements with currently used density based CT measurements in the literature, namely, the relative area measure and the D measure. Results: Volumetric CT data, performed at total lung capacity and residual volume, from a total of 49 subjects enrolled in the MILES trial were used in our study. Our proposed measures had adjusted R{sup 2} ranging from 0.42 to 0.59 when regressing against the spirometry measures, with p < 0.05. For previously used density based CT measurements in the literature, the best R{sup 2} was 0.46 (for only one instance), with the majority being lower than 0.3 or p > 0.05. Conclusions: The proposed family of CT-based cyst measurements have better correlation with spirometric measures than previously used density based CT measurements. They show potential as a sensitive tool for quantitatively assessing the severity of LAM.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Online updating of context-aware landmark detectors for prostate localization in daily treatment CT images

Purpose: In image guided radiation therapy, it is crucial to fast and accurately localize the prostate in the daily treatment images. To this end, the authors propose an online update scheme for landmark-guided prostate segmentation, which can fully exploit valuable patient-specific information contained in the previous treatment images and can achieve improved performance in landmark detection and prostate segmentation. Methods: To localize the prostate in the daily treatment images, the authors first automatically detect six anatomical landmarks on the prostate boundary by adopting a context-aware landmark detection method. Specifically, in this method, a two-layer regression forest is trained as a detector for each target landmark. Once all the newly detected landmarks from new treatment images are reviewed or adjusted (if necessary) by clinicians, they are further included into the training pool as new patient-specific information to update all the two-layer regression forests for the next treatment day. As more and more treatment images of the current patient are acquired, the two-layer regression forests can be continually updated by incorporating the patient-specific information into the training procedure. After all target landmarks are detected, a multiatlas random sample consensus (multiatlas RANSAC) method is used to segment the entire prostate by fusing multiple previously segmented prostates of the current patient after they are aligned to the current treatment image. Subsequently, the segmented prostate of the current treatment image is again reviewed (or even adjusted if needed) by clinicians before including it as a new shape example into the prostate shape dataset for helping localize the entire prostate in the next treatment image. Results: The experimental results on 330 images of 24 patients show the effectiveness of the authors’ proposed online update scheme in improving the accuracies of both landmark detection and prostate segmentation. Besides, compared to the other state-of-the-art prostate segmentation methods, the authors’ method achieves the best performance. Conclusions: By appropriate use of valuable patient-specific information contained in the previous treatment images, the authors’ proposed online update scheme can obtain satisfactory results for both landmark detection and prostate segmentation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-251: Incorporation of Pre-Therapy 18F-FDG Uptake with CT Texture Features in a Predictive Model for Radiation Pneumonitis Development

Purpose: To determine whether the addition of standardized uptake value (SUV) statistical variables to CT lung texture features can improve a predictive model of radiation pneumonitis (RP) development in patients undergoing radiation therapy. Methods: Anonymized data from 96 esophageal cancer patients (18 RP-positive cases of Grade ≥ 2) were retrospectively collected including pre-therapy PET/CT scans, pre-/posttherapy diagnostic CT scans and RP status. Twenty texture features (firstorder, fractal, Laws’ filter and gray-level co-occurrence matrix) were calculated from diagnostic CT scans and compared in anatomically matched regions of the lung. The mean, maximum, standard deviation, and 50th–95th percentiles of the SUV values for all lung voxels in the corresponding PET scans were acquired. For each texture feature, a logistic regression-based classifier consisting of (1) the average change in that texture feature value between the pre- and post-therapy CT scans and (2) the pre-therapy SUV standard deviation (SUV{sub SD}) was created. The RP-classification performance of each logistic regression model was compared to the performance of its texture feature alone by computing areas under the receiver operating characteristic curves (AUCs). T-tests were performed to determine whether the mean AUC across texture features changed significantly when SUV{sub SD} was added to the classifier. Results: The AUC for single-texturefeature classifiers ranged from 0.58–0.81 in high-dose (≥ 30 Gy) regions of the lungs and from 0.53–0.71 in low-dose (< 10 Gy) regions. Adding SUVSD in a logistic regression model using a 50/50 data partition for training and testing significantly increased the mean AUC by 0.08, 0.06 and 0.04 in the low-, medium- and high-dose regions, respectively. Conclusion: Addition of SUVSD from a pre-therapy PET scan to a single CT-based texture feature improves RP-classification performance on average. These findings demonstrate the potential for more accurate prediction of RP using information from multiple imaging modalities. Supported, in part, by the National Institute of Biomedical Imaging and Bioengineering of the National Institutes of Health under grant number T32 EB002103; SGA receives royalties and licensing fees through the University of Chicago for computer-aided diagnosis technology. HA receives royalties through the University of Chicago for computer-aided diagnosis technology.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

TU-G-BRA-04: Changes in Regional Lung Function Measured by 4D-CT Ventilation Imaging for Thoracic Radiotherapy

Purpose: To test the hypothesis: 4D-CT ventilation imaging can show the known effects of radiotherapy on lung function: (1) radiation-induced ventilation reductions, and (2) ventilation increases caused by tumor regression. Methods: Repeat 4D-CT scans (pre-, mid- and/or post-treatment) were acquired prospectively for 11 thoracic cancer patients in an IRB-approved clinical trial. A ventilation image for each time point was created using deformable image registration and the Hounsfield unit (HU)-based or Jacobian-based metric. The 11 patients were divided into two subgroups based on tumor volume reduction using a threshold of 5 cm{sup 3}. To quantify radiation-induced ventilation reduction, six patients who showed a small tumor volume reduction (<5 cm{sup 3}) were analyzed for dose-response relationships. To investigate ventilation increase caused by tumor regression, two of the other five patients were analyzed to compare ventilation changes in the lung lobes affected and unaffected by the tumor. The remaining three patients were excluded because there were no unaffected lobes. Results: Dose-dependent reductions of HU-based ventilation were observed in a majority of the patient-specific dose-response curves and in the population-based dose-response curve, whereas no clear relationship was seen for Jacobian-based ventilation. The post-treatment population-based dose-response curve of HU-based ventilation demonstrated the average ventilation reductions of 20.9±7.0% at 35–40 Gy (equivalent dose in 2-Gy fractions, EQD2), and 40.6±22.9% at 75–80 Gy EQD2. Remarkable ventilation increases in the affected lobes were observed for the two patients who showed an average tumor volume reduction of 37.1 cm{sup 3} and re-opening airways. The mid-treatment increase in HU-based ventilation of patient 3 was 100.4% in the affected lobes, which was considerably greater than 7.8% in the unaffected lobes. Conclusion: This study has demonstrated that 4D-CT ventilation imaging shows the known effects of radiotherapy on lung function: radiation-induced ventilation reduction and ventilation increase caused by tumor regression, providing validation for 4D-CT ventilation imaging. This study was supported in part by a National Lung Cancer Partnership Young Investigator Research grant.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-C-BRA-05: Delineating High-Dose Clinical Target Volumes for Head and Neck Tumors Using Machine Learning Algorithms

Purpose: To develop and test population-based machine learning algorithms for delineating high-dose clinical target volumes (CTVs) in H&N tumors. Automating and standardizing the contouring of CTVs can reduce both physician contouring time and inter-physician variability, which is one of the largest sources of uncertainty in H&N radiotherapy. Methods: Twenty-five node-negative patients treated with definitive radiotherapy were selected (6 right base of tongue, 11 left and 9 right tonsil). All patients had GTV and CTVs manually contoured by an experienced radiation oncologist prior to treatment. This contouring process, which is driven by anatomical, pathological, and patient specific information, typically results in non-uniform margin expansions about the GTV. Therefore, we tested two methods to delineate high-dose CTV given a manually-contoured GTV: (1) regression-support vector machines(SVM) and (2) classification-SVM. These models were trained and tested on each patient group using leave-one-out cross-validation. The volume difference(VD) and Dice similarity coefficient(DSC) between the manual and auto-contoured CTV were calculated to evaluate the results. Distances from GTV-to-CTV were computed about each patient’s GTV and these distances, in addition to distances from GTV to surrounding anatomy in the expansion direction, were utilized in the regression-SVM method. The classification-SVM method used categorical voxel-information (GTV, selected anatomical structures, else) from a 3×3×3cm3 ROI centered about the voxel to classify voxels as CTV. Results: Volumes for the auto-contoured CTVs ranged from 17.1 to 149.1cc and 17.4 to 151.9cc; the average(range) VD between manual and auto-contoured CTV were 0.93 (0.48–1.59) and 1.16(0.48–1.97); while average(range) DSC values were 0.75(0.59–0.88) and 0.74(0.59–0.81) for the regression-SVM and classification-SVM methods, respectively. Conclusion: We developed two novel machine learning methods to delineate high-dose CTV for H&N patients. Both methods showed promising results that hint to a solution to the standardization of the contouring process of clinical target volumes. Varian Medical Systems grant.

60 APPLIED LIFE SCIENCES↗