Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrices”

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,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Polarimetric Interferometric SAR Change Detection Discrimination

A coherent change detection (CCD) image, computed from a geometrically matched, temporally separated pair of complex-valued synthetic aperture radar (SAR) image sets, conveys the pixel-level equivalence between the two observations. Low-coherence values in a CCD image are typically due to either some physical change in the corresponding pixels or a low signal-to-noise observation. A CCD image does not directly convey the nature of the change that occurred to cause low coherence. In this paper, we introduce a mathematical framework for discriminating between different types of change within a CCD image. We utilize the extra degrees of freedom and information from polarimetric interferometric SAR (PolInSAR) data and PolInSAR processing techniques to define a 29-dimensional feature vector that contains information capable of discriminating between different types of change in a scene. We also propose two change-type discrimination functions that can be trained with feature vector training data and demonstrate change-type discrimination on an example image set for three different types of change. In conclusion, we also describe and characterize the performance of the two proposed change-type discrimination functions by way of receiver operating characteristic curves, confusion matrices, and pass matrices.

58 GEOSCIENCES↗
DOE OSTI2017

D–Q Impedance Based Stability Analysis and Parameter Design of Three-Phase Inverter-Based AC Power Systems

Small-signal stability is an important concern in three-phase inverter-based ac power systems. The impedance-based approach based on the generalized Nyquist stability criterion (GNC) can analyze the stability related with the medium and high-frequency modes of the systems. However,. the GNC involves the right-half-plane (RHP) pole calculation of return-ratio transfer function matrices, which cannot be avoided for stability analysis of complicated ac power systems. Therefore, it necessitates the detailed internal control information of the inverters, which is not normally available for commercial inverters. To address this issue, this paper introduces the component connection method (CCM) in the frequency domain for stability analysis in the synchronous d-q frame, by proposing a method of deriving the impedance matrix of the connection networks of inverter-based ac power systems. Demonstration on a two-area system and a microgrid shows that: The CCM-enabled approach can avoid the RHP pole calculation of return-ratio matrices and enables the stability analysis by using only the impedances of system components, which could be measured without the need for the internal information. A stability analysis method based on d-q impedances, the CCM, and the determinant-based GNC is also proposed to further simplify the analysis process. Inverter controller parameters can be designed as stability regions in parameter spaces, by repetitively applying the proposed stability analysis method. Simulation and experimental results verify the validity of the proposed stability analysis method and the parameter design approach.

42 ENGINEERING↗
DOE OSTI2020

Reduce Operations: Send Volume Balancing While Minimizing Latency

Communication hypergraph model was introduced in a two-phase setting for encapsulating multiple communication cost metrics (bandwidth and latency), which are proven to be important in parallelizing irregular applications. In the first phase, computational-task-to-processor assignment is performed with the objective of minimizing total volume while maintaining computational load balance. In the second phase, communication-task-to-processor assignment is performed with the objective of minimizing total number of messages while maintaining communication-volume balance. The reduce-communication hypergraph model suffers from failing to correctly encapsulate send-volume balancing. We propose a novel vertex weighting scheme that enables part weights to correctly encode send-volume loads of processors for send-volume balancing. The model also suffers from increasing the total communication volume during partitioning. To decrease this increase, we propose a method that utilizes the recursive bipartitioning framework and refines each bipartition by vertex swaps. For performance evaluation, we consider column-parallel SpMV, which is one of the most widely known applications in which the reduce-task assignment problem arises. Extensive experiments on 313 matrices show that, compared to the existing model, the proposed models achieve considerable improvements in all communication cost metrics. Furthermore, these improvements lead to an average decrease of 30% in parallel SpMV time on 512 processors for 70 matrices with high irregularity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

GSoFa: Scalable Sparse Symbolic LU Factorization on GPUs

Decomposing a matrix $\mathbf {A}$ into a lower matrix $\mathbf {L}$ and an upper matrix $\mathbf {U}$, which is also known as LU decomposition, is an essential operation in numerical linear algebra. For a sparse matrix, LU decomposition often introduces more nonzero entries in the $\mathbf {L}$ and $\mathbf {U}$ factors than in the original matrix. A symbolic factorization step is needed to identify the nonzero structures of $\mathbf {L}$ and $\mathbf {U}$ matrices. Attracted by the enormous potentials of the Graphics Processing Units (GPUs), an array of efforts have surged to deploy various LU factorization steps except for the symbolic factorization, to the best of our knowledge, on GPUs. This article introduces gSoFa, the first GPU-based symbolic factorization design with the following three optimizations to enable scalable LU symbolic factorization for nonsymmetric pattern sparse matrices on GPUs. First, here we introduce a novel fine-grained parallel symbolic factorization algorithm that is well suited for the Single Instruction Multiple Thread (SIMT) architecture of GPUs. Second, we tailor supernode detection into a SIMT friendly process and strive to balance the workload, minimize the communication and saturate the GPU computing resources during supernode detection. Third, we introduce a three-pronged optimization to reduce the excessive space consumption problem faced by multi-source concurrent symbolic factorization. Taken together, gSoFa achieves up to 31× speedup from 1 to 44 Summit nodes (6 to 264 GPUs) and outperforms the state-of-the-art CPU project, on average, by 5×. Notably, gSoFa also achieves up to 47 percent of the peak memory throughput of a V100 GPU in the Summit Supercomputer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

A Sparse Tableau Formulation for Node-Breaker Representations in Security-Constrained Optimal Power Flow

Realistic representations of contingencies in AC optimal power flow (ACOPF) often challenge traditional bus-branch network models. Derived from nodal analysis, such approaches represent network constraints in terms of the familiar bus admittance matrix, Ybus. A fixed Ybus is unable to represent common circuit breaker actions such as bus splitting. Work-arounds for Ybus-based analysis typically rely on topology processing, switching between different Ybus matrices depending on breaker settings. In this paper, we propose a very general node-breaker approach, employing multi-port element models and using a sparse tableau formulation (STF) for network constraints. Instead of treating breaker action as altering network topology, and hence changing the structure of Kirchhoff’s voltage law (KVL) and Kirchhoff’s current law (KCL) equations, this approach represents a breaker’s position as impacting only constraints associated with a single component, thereby maintaining fixed structure in the KVL and KCL constraints. While larger numbers of variables are required, STF proves sparser than Ybus formulations. Numerical case studies herein demonstrate that STF provides computational efficiency comparable to a Ybus-based ACOPF at the scale of several hundred buses, and lower computational cost in an example of over one thousand buses.

Park, Byungkwon (ORCID:000000033189816X)↗
DOE OSTI2023

Sensitivity Analysis of a Polyphase Wireless Power Transfer System under Off-Nominal Conditions

Misaligned and/or variable coil airgaps cause coupling coefficient variation in wireless power transfer (WPT) systems, resulting in a decrease in the system’s transmitted power and efficiency. This paper presents sensitivity analyses of a three-phase, Y-Y connected, series-tuned WPT system in the frequency domain in terms of several different electric vehicle wireless charging off-nominal conditions (misalignments in Δ x and Δ y directions, change of airgap Δ z , and the roll Δψ, pitch Δθ, and yaw ΔΦ angles) as specified in Society of Automotive Engineers (SAE) J2954 Standard. Coil inductance matrices were obtained by measuring the self- and mutual inductances of the primary and secondary coil phase windings at variable airgap classes (from 5 cm to 30 cm) for five different charging positions as identified in SAE J2954. These 6 × 6 inductance matrices were used in sensitivity and Plexim/PLECS simulation analyses. The sensitivity of the WPT system was analyzed analytically using the input impedance and phase angle, voltage gain, current gain, quality factor, and coupling coefficient parameters of the series-tuned WPT system. The results were confirmed experimentally on a 50 kW WPT system.

42 ENGINEERING↗
DOE OSTI2020

Pore-scale view of microbial turnover: combining 14 C imaging, μCT, and zymography after adding soluble carbon to soil pores of specific sizes

The location of microorganisms and substrates within soil pore networks plays a crucial role for organic carbon (C) processing, its microbial utilization and turnover, and has direct consequences for C and nutrient cycling. An optimal approach to quantify responses to new C inputs from microorganisms residing in specific pores is the addition of new C to pores of target sizes in undisturbed soil cores. We used the matric potential approach to add 14 C labeled glucose to small (< 40 μm, root free) or large (60-180 μm, potentially inhabited by roots) pores of undisturbed soil cores. Localization of glucose-derived C via 14 C imaging was related to pore size distributions and connectivity, characterized via X-ray computed micro-tomography (μCT), and to β-glucosidase activity, characterized via zymography. After two-week incubations, 1.3 times more glucose was mineralized ( 14 CO 2 ) when it was added to the large pores, however more 14 C remained in microbial biomass when glucose was added to the small pores. Consequently, while utilizing the same amounts of easily available C, the microorganisms localized in the large pores had faster turnover compared to microorganisms in small pores. Stronger associations between β-glucosidase activity and glucose-derived C were observed when glucose was added to the large pores. We conclude that 1) the matric potential approach allows placing, albeit not exactly, soluble substrates into pores of target diameter range, and 2) microorganisms localized in large pores respond to new C inputs with faster turnover, greater growth, and more intensive enzyme production compared to those inhabiting the small pores.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Grain size variation in nanocrystalline silicon carbide irradiated at elevated temperatures

This study reports on ion irradiation-induced grain growth in nanocrystalline SiC films on Si substrates. The SiC grains with average size ranging from ~2 to 20 nm were embedded in amorphous SiC matrices. Irradiation was performed using 5 MeV Xe23+ ions to 115 ions/nm2 at 700 K. The irradiated films were characterized using x-ray diffraction, transmission electron microscopy and Raman spectroscopy. Significant grain growth is observed for smaller grains that tend to saturate at ~8 nm. In contrast, irradiation of larger grains (~20 nm in size) leads to a decrease in the grain size, which could be associated with the production of lattice disorder within the grains. Homonuclear C-C bonds in the irradiated amorphous SiC matrix are found to be graphitized. This bonding transformation could limit or inhibit grain growth and contribute to the size saturation. The results from this study may suggest nanocrystalline SiC as a promising candidate structural or cladding material for applications in advanced nuclear reactors.

Grain growth, ion irradiation, SiC↗
DOE OSTI2013

Investigation of processing advanced waste forms in Cold Crucible Induction Melter systems

The Cold Crucible Induction Melter (CCIM) technology has been identified in several independent reviews and studies as a potential alternative or supplemental vitrification technology for specific U.S. Department of Energy (DOE) high level radioactive waste (HLW) streams. These inventories have unique chemistries, such as high alumina, high iron, high chromium, etc., that are not amenable to efficient immobilization in conventional borosilicate glass (BSG) matrices. Advanced waste forms, such as iron phosphate glass (FePG), high crystalline BSG (HC-BSG), and alumino-silicate glass-ceramic (GC) have been shown to provide significantly improved immobilization matrices for these challenging waste streams, including dramatically improved waste loadings; however, they are difficult to efficiently process in the baseline joule-heated ceramic melters (JHCMs). The CCIM technology can effectively process these advanced waste forms at the desired high waste loadings. This paper describes the most recent collaborative activities between the DOE Office of Environmental Management (DOE EM) and various Russian institutions in modeling, development and testing of advanced waste forms in CCIMs, including FePG and HC-BSG formulations. Additionally, past interactions with other international partners, including France and Korea, are described. Finally, recent enhancements to CCIM processing related to modeling, draining techniques and automated control, which were collaboratively developed, will be discussed. (authors)

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

The phyto-remediation of radioactively contaminated land - a feasible approach or just bananas?

Soil is an essential component of all terrestrial ecosystems and is under increasing threat from human activity. Techniques available for removing radioactive contamination from soil and aquatic substrates are limited and often costly to implement; particularly over large areas. Frequently, bulk soil removal, with its attendant consequences, is a significant component of the majority of contamination incidents. Alternative techniques capable of removing contamination or exposure pathways without damaging or removing the soil are therefore of significant interest. An increasing number of old nuclear facilities are entering 'care and maintenance', with significant ground contamination issues. Phyto-remediation - the use of plants' natural metabolic processes to remediate contaminated sites is one possible solution. Its key mechanisms include phyto-extraction and phyto-stabilisation. These are analogues of existing remedial techniques. Further, phyto-remediation can improve soil quality and stability and restore functionality. Information on the application of phyto-remediation in the nuclear industry is widely distributed over an extended period of time and sources. It is therefore difficult to quickly and effectively identify which plants would be most suitable for phyto-remediation on a site by site basis. In response, a phyto-remediation tool has been developed to address this issue. Existing research and case studies were reviewed to understand the mechanisms of phyto-remediation, its effectiveness and the benefits and limitations of implementation. The potential for cost recovery from a phyto-remediation system is also briefly considered. An overview of this information is provided here. From this data, a set of matrices was developed to guide potential users through the plant selection process. The matrices take the user through a preliminary screening process to determine whether the contamination present at their site is amenable to phyto-remediation, and to give a rough indication as to what plants might be suitable. The second two allow the user to target specific plant species that would be most likely to successfully establish based on prevailing site conditions. The outcome of this study is a phyto-remediation tool that can facilitate the development of phyto-remediation projects, avoiding the need for in-depth research to identify optimal plant species on a case-by-case basis. (authors)

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

An algorithm to calculate a collapsed arc dose matrix in volumetric modulated arc therapy

Purpose: The delivery of volumetric modulated arc therapy (VMAT) is more complex than other conformal radiotherapy techniques. In this work, the authors present the feasibility of performing routine verification of VMAT delivery using a dose matrix measured by a gantry mounted 2D ion chamber array and corresponding dose matrix calculated by an inhouse developed algorithm.Methods: Pinnacle, v9.0, treatment planning system (TPS) was used in this study to generate VMAT plans for a 6 MV photon beam from an Elekta-Synergy linear accelerator. An algorithm was developed and implemented with inhouse computer code to calculate the dose matrix resulting from a VMAT arc in a plane perpendicular to the beam at isocenter. The algorithm was validated using measurement of standard patterns and clinical VMAT plans with a 2D ion chamber array. The clinical VMAT plans were also validated using ArcCHECK measurements. The measured and calculated dose matrices were compared using gamma ({gamma}) analysis with 3%/3 mm criteria and {gamma} tolerance of 1.Results: The dose matrix comparison of standard patterns has shown excellent agreement with the mean {gamma} pass rate 97.7 ({sigma}= 0.4)%. The validation of clinical VMAT plans using the dose matrix predicted by the algorithm and the corresponding measured dose matrices also showed good agreement with the mean {gamma} pass rate of 97.6 ({sigma}= 1.6)%. The validation of clinical VMAT plans using ArcCHECK measurements showed a mean pass rate of 95.6 ({sigma}= 1.8)%.Conclusions: The developed algorithm was shown to accurately predict the dose matrix, in a plane perpendicular to the beam, by considering all possible leaf trajectories in a VMAT delivery. This enables the verification of VMAT delivery using a 2D array detector mounted on a treatment head.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

Expediting model-based optoacoustic reconstructions with tomographic symmetries

Purpose: Image quantification in optoacoustic tomography implies the use of accurate forward models of excitation, propagation, and detection of optoacoustic signals while inversions with high spatial resolution usually involve very large matrices, leading to unreasonably long computation times. The development of fast and memory efficient model-based approaches represents then an important challenge to advance on the quantitative and dynamic imaging capabilities of tomographic optoacoustic imaging. Methods: Herein, a method for simplification and acceleration of model-based inversions, relying on inherent symmetries present in common tomographic acquisition geometries, has been introduced. The method is showcased for the case of cylindrical symmetries by using polar image discretization of the time-domain optoacoustic forward model combined with efficient storage and inversion strategies. Results: The suggested methodology is shown to render fast and accurate model-based inversions in both numerical simulations andpost mortem small animal experiments. In case of a full-view detection scheme, the memory requirements are reduced by one order of magnitude while high-resolution reconstructions are achieved at video rate. Conclusions: By considering the rotational symmetry present in many tomographic optoacoustic imaging systems, the proposed methodology allows exploiting the advantages of model-based algorithms with feasible computational requirements and fast reconstruction times, so that its convenience and general applicability in optoacoustic imaging systems with tomographic symmetries is anticipated.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-E-T-465: Dose Calculation Method for Dynamic Tumor Tracking Using a Gimbal-Mounted Linac

Purpose: Dynamic tumor tracking using the gimbal-mounted linac (Vero4DRT, Mitsubishi Heavy Industries, Ltd., Japan) has been available when respiratory motion is significant. The irradiation accuracy of the dynamic tumor tracking has been reported to be excellent. In addition to the irradiation accuracy, a fast and accurate dose calculation algorithm is needed to validate the dose distribution in the presence of respiratory motion because the multiple phases of it have to be considered. A modification of dose calculation algorithm is necessary for the gimbal-mounted linac due to the degrees of freedom of gimbal swing. The dose calculation algorithm for the gimbal motion was implemented using the linear transformation between coordinate systems. Methods: The linear transformation matrices between the coordinate systems with and without gimbal swings were constructed using the combination of translation and rotation matrices. The coordinate system where the radiation source is at the origin and the beam axis along the z axis was adopted. The transformation can be divided into the translation from the radiation source to the gimbal rotation center, the two rotations around the center relating to the gimbal swings, and the translation from the gimbal center to the radiation source. After operating the transformation matrix to the phantom or patient image, the dose calculation can be performed as the no gimbal swing. The algorithm was implemented in the treatment planning system, PlanUNC (University of North Carolina, NC). The convolution/superposition algorithm was used. The dose calculations with and without gimbal swings were performed for the 3 × 3 cm{sup 2} field with the grid size of 5 mm. Results: The calculation time was about 3 minutes per beam. No significant additional time due to the gimbal swing was observed. Conclusions: The dose calculation algorithm for the finite gimbal swing was implemented. The calculation time was moderate.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

SU-E-QI-17: Dependence of 3D/4D PET Quantitative Image Features On Noise

Purpose: Quantitative imaging is a fast evolving discipline where a large number of features are extracted from images; i.e., radiomics. Some features have been shown to have diagnostic, prognostic and predictive value. However, they are sensitive to acquisition and processing factors; e.g., noise. In this study noise was added to positron emission tomography (PET) images to determine how features were affected by noise. Methods: Three levels of Gaussian noise were added to 8 lung cancer patients PET images acquired in 3D mode (static) and using respiratory tracking (4D); for the latter images from one of 10 phases were used. A total of 62 features: 14 shape, 19 intensity (1stO), 18 GLCM textures (2ndO; from grey level co-occurrence matrices) and 11 RLM textures (2ndO; from run-length matrices) features were extracted from segmented tumors. Dimensions of GLCM were 256×256, calculated using 3D images with a step size of 1 voxel in 13 directions. Grey levels were binned into 256 levels for RLM and features were calculated in all 13 directions. Results: Feature variation generally increased with noise. Shape features were the most stable while RLM were the most unstable. Intensity and GLCM features performed well; the latter being more robust. The most stable 1stO features were compactness, maximum and minimum length, standard deviation, root-mean-squared, I30, V10-V90, and entropy. The most stable 2ndO features were entropy, sum-average, sum-entropy, difference-average, difference-variance, difference-entropy, information-correlation-2, short-run-emphasis, long-run-emphasis, and run-percentage. In general, features computed from images from one of the phases of 4D scans were more stable than from 3D scans. Conclusion: This study shows the need to characterize image features carefully before they are used in research and medical applications. It also shows that the performance of features, and thereby feature selection, may be assessed in part by noise analysis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-157: CARMEN: A MatLab-Based Research Platform for Monte Carlo Treatment Planning (MCTP) and Customized System for Planning Evaluation

Purpose: Although there exist several radiotherapy research platforms, such as: CERR, the most widely used and referenced; SlicerRT, which allows treatment plan comparison from various sources; and MMCTP, a full MCTP system; it is still needed a full MCTP toolset that provides users complete control of calculation grids, interpolation methods and filters in order to “fairly” compare results from different TPSs, supporting verification with experimental measurements. Methods: This work presents CARMEN, a MatLab-based platform including multicore and GPGPU accelerated functions for loading RT data; designing treatment plans; and evaluating dose matrices and experimental data.CARMEN supports anatomic and functional imaging in DICOM format, as well as RTSTRUCT, RTPLAN and RTDOSE. Besides, it contains numerous tools to accomplish the MCTP process, managing egs4phant and phase space files.CARMEN planning mode assist in designing IMRT, VMAT and MERT treatments via both inverse and direct optimization. The evaluation mode contains a comprehensive toolset (e.g. 2D/3D gamma evaluation, difference matrices, profiles, DVH, etc.) to compare datasets from commercial TPS, MC simulations (i.e. 3ddose) and radiochromic film in a user-controlled manner. Results: CARMEN has been validated against commercial RTPs and well-established evaluation tools, showing coherent behavior of its multiple algorithms. Furthermore, CARMEN platform has been used to generate competitive complex treatment that has been published in comparative studies. Conclusion: A new research oriented MCTP platform with a customized validation toolset has been presented. Despite of being coded with a high-level programming language, CARMEN is agile due to the use of parallel algorithms. The wide-spread use of MatLab provides straightforward access to CARMEN’s algorithms to most researchers. Similarly, our platform can benefit from the MatLab community scientific developments as filters, registration algorithms etc. Finally, CARMEN arises the importance of grid and filtering control in treatment plan comparison.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

TU-AB-BRA-06: Texture Feature Reproducibility Between PET/CT and PET/MR Imaging Modalities

Purpose: Radiomics involves the extraction of texture features from different imaging modalities with the purpose of developing models to predict patient treatment outcomes. The purpose of this study is to investigate texture feature reproducibility across [18F]FDG PET/CT and [18F]FDG PET/MR imaging in patients with primary malignancies. Methods: Twenty five prospective patients with solid tumors underwent clinical [18F]FDG PET/CT scan followed by [18F]FDG PET/MR scans. In all patients the lesions were identified using nuclear medicine reports. The images were co-registered and segmented using an in-house auto-segmentation method. Fifty features, based on the intensity histogram, second and high order matrices, were extracted from the segmented regions from both image data sets. One-way random-effects ANOVA model of the intra-class correlation coefficient (ICC) was used to establish texture feature correlations between both data sets. Results: Fifty features were classified based on their ICC values, which were found in the range from 0.1 to 0.86, in three categories: high, intermediate, and low. Ten features extracted from second and high-order matrices showed large ICC ≥ 0.70. Seventeen features presented intermediate 0.5 ≤ ICC ≤ 0.65 and the remaining twenty three presented low ICC ≤ 0.45. Conclusion: Features with large ICC values could be reliable candidates for quantification as they lead to similar results from both imaging modalities. Features with small ICC indicates a lack of correlation. Therefore, the use of these features as a quantitative measure will lead to different assessments of the same lesion depending on the imaging modality from where they are extracted. This study shows the importance of the need for further investigation and standardization of features across multiple imaging modalities.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

TU-F-CAMPUS-J-05: Effect of Uncorrelated Noise Texture On Computed Tomography Quantitative Image Features

Purpose: To investigate the relationship between quantitative image features (i.e. radiomics) and statistical fluctuations (i.e. electronic noise) in clinical Computed Tomography (CT) using the standardized American College of Radiology (ACR) CT accreditation phantom and patient images. Methods: Three levels of uncorrelated Gaussian noise were added to CT images of phantom and patients (20) acquired in static mode and respiratory tracking mode. We calculated the noise-power spectrum (NPS) of the original CT images of the phantom, and of the phantom images with added Gaussian noise with means of 50, 80, and 120 HU. Concurrently, on patient images (original and noise-added images), image features were calculated: 14 shape, 19 intensity (1st order statistics from intensity volume histograms), 18 GLCM features (2nd order statistics from grey level co-occurrence matrices) and 11 RLM features (2nd order statistics from run-length matrices). These features provide the underlying structural information of the images. GLCM (size 128x128) was calculated with a step size of 1 voxel in 13 directions and averaged. RLM feature calculation was performed in 13 directions with grey levels binning into 128 levels. Results: Adding the electronic noise to the images modified the quality of the NPS, shifting the noise from mostly correlated to mostly uncorrelated voxels. The dramatic increase in noise texture did not affect image structure/contours significantly for patient images. However, it did affect the image features and textures significantly as demonstrated by GLCM differences. Conclusion: Image features are sensitive to acquisition factors (simulated by adding uncorrelated Gaussian noise). We speculate that image features will be more difficult to detect in the presence of electronic noise (an uncorrelated noise contributor) or, for that matter, any other highly correlated image noise. This work focuses on the effect of electronic, uncorrelated, noise and future work shall examine the influence of changes in quantum noise on the features. J. Oliver was supported by NSF FGLSAMP BD award HRD #1139850 and the McKnight Doctoral Fellowship.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-D-207B-01: Radiomics Feature Reproducibility From Repeat CT Scans of Patients with Rectal Cancer

Purpose: To evaluate the reproducibility of radiomics features by repeating computed tomographic (CT) scans in rectal cancer. To choose stable radiomics features for rectal cancer. Methods: 40 rectal cancer patients were enrolled in this study, each of whom underwent two CT scans within average 8.7 days (5 days to 17 days), before any treatment was delivered. The rectal gross tumor volume (GTV) was distinguished and segmented by an experienced oncologist in both CTs. Totally, more than 2000 radiomics features were defined in this study, which were divided into four groups (I: GLCM, II: GLRLM III: Wavelet GLCM and IV: Wavelet GLRLM). For each group, five types of features were extracted (Max slice: features from the largest slice of target images, Max value: features from all slices of target images and choose the maximum value, Min value: minimum value of features for all slices, Average value: average value of features for all slices, Matrix sum: all slices of target images translate into GLCM and GLRLM matrices and superpose all matrices, then extract features from the superposed matrix). Meanwhile a LOG (Laplace of Gauss) filter with different parameters was applied to these images. Concordance correlation coefficients (CCC) and inter-class correlation coefficients (ICC) were calculated to assess the reproducibility. Results: 403 radiomics features were extracted from each type of patients’ medical images. Features of average type are the most reproducible. Different filters have little effect for radiomics features. For the average type features, 253 out of 403 features (62.8%) showed high reproducibility (ICC≥0.8), 133 out of 403 features (33.0%) showed medium reproducibility (0.8≥ICC≥0.5) and 17 out of 403 features (4.2%) showed low reproducibility (ICC≥0.5). Conclusion: The average type radiomics features are the most stable features in rectal cancer. Further analysis of these features of rectal cancer can be warranted for treatment monitoring and prognosis prediction.

60 APPLIED LIFE SCIENCES↗