Search NASA⌕ Search

SEARCH · Search NASA

Results for “functional Principal Component Analysis”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Use of UV Sources for Detection and Identification of Explosives

Measurement of Raman and native fluorescence emission using ultraviolet (UV) sources (<400 nm) on targeted materials is suitable for both sensitive detection and accurate identification of explosive materials. When the UV emission data are analyzed using a combination of Principal Component Analysis (PCA) and cluster analysis, chemicals and biological samples can be differentiated based on the geometric arrangement of molecules, the number of repeating aromatic rings, associated functional groups (nitrogen, sulfur, hydroxyl, and methyl), microbial life cycles (spores vs. vegetative cells), and the number of conjugated bonds. Explosive materials can be separated from one another as well as from a range of possible background materials, which includes microbes, car doors, motor oil, and fingerprints on car doors, etc. Many explosives are comprised of similar atomic constituents found in potential background samples such as fingerprint oils/skin, motor oil, and soil. This technique is sensitive to chemical bonds between the elements that lead to the discriminating separability between backgrounds and explosive materials.

Hug, William↗
DOE OSTI2019

Function Preserving Projection

With function preserving projection (FPP), we aim to address high-dimensional data exploration challenges common in many scientific applications. Instead of preserving the global or local geometric structure of a dataset like traditional dimensionality reduction methods (e.g., principal component analysis, t-stochastic neighborhood embedding), the proposed approach finds a linear projection that reveal interpretable (yet potentially nonlinear) pattern of the properties of interests defined on a highdimensional domain (e.g., simulation yield defined on the input parameter space). The key insight driving FPP is the fact that while non-linear embeddings are difficult to interpret, humans are highly skilled in understanding complex visual patterns. FPP lets users select the complexity of a pattern they consider interpretable, allows to find a joint projection that explains multiple functions, and provides a confidence statistic (akin to a p-value) indicating how salient the resulting pattern is. We demonstrate the efficacy of the proposed method on several synthetic and scientific datasets for both continuous (regression) and discrete (classification) problems of interests and illustrate the scalability of FPP by finding an optimal projection for close to 40 gigabytes in around 10 minutes.

Liu, SHUSEN↗
DOE OSTI2022

Training and projecting: A reduced basis method emulator for many-body physics

Here, we present the reduced basis method as a tool for developing emulators for equations with tun able parameters within the context of the nuclear many-body problem. The method uses a basis expansion informed by a set of solutions for a few values of the model parameters and then projects the equations over a well-chosen low-dimensional subspace. We connect some of the results in the eigenvector continuation literature to the formalism of reduced basis methods and show how these methods can be applied to a broad set of problems. As we illustrate, the possible success of the formalism on such problems can be diagnosed beforehand by a principal component analysis. We apply the reduced basis method to the one-dimensional Gross-Pitaevskii equation with a harmonic trap ping potential and to nuclear density functional theory for 48 Ca, achieving speed-ups of more than x150 and x250, respectively, when compared to traditional solvers. The outstanding performance of the approach, together with its straightforward implementation, show promise for its application to the emulation of computationally demanding calculations, including uncertainty quantification.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

COSMIC SHEAR RESULTS FROM THE DEEP LENS SURVEY. I. JOINT CONSTRAINTS ON {Omega} {sub M} AND {sigma}{sub 8} WITH A TWO-DIMENSIONAL ANALYSIS

We present a cosmic shear study from the Deep Lens Survey (DLS), a deep BVRz multi-band imaging survey of five 4 deg{sup 2} fields with two National Optical Astronomy Observatory (NOAO) 4 m telescopes at Kitt Peak and Cerro Tololo. For both telescopes, the change of the point-spread-function (PSF) shape across the focal plane is complicated, and the exposure-to-exposure variation of this position-dependent PSF change is significant. We overcome this challenge by modeling the PSF separately for individual exposures and CCDs with principal component analysis (PCA). We find that stacking these PSFs reproduces the final PSF pattern on the mosaic image with high fidelity, and the method successfully separates PSF-induced systematics from gravitational lensing effects. We calibrate our shears and estimate the errors, utilizing an image simulator, which generates sheared ground-based galaxy images from deep Hubble Space Telescope archival data with a realistic atmospheric turbulence model. For cosmological parameter constraints, we marginalize over shear calibration error, photometric redshift uncertainty, and the Hubble constant. We use cosmology-dependent covariances for the Markov Chain Monte Carlo analysis and find that the role of this varying covariance is critical in our parameter estimation. Our current non-tomographic analysis alone constrains the {Omega} {sub M}-{sigma}{sub 8} likelihood contour tightly, providing a joint constraint of {Omega} {sub M} = 0.262 {+-} 0.051 and {sigma}{sub 8} = 0.868 {+-} 0.071. We expect that a future DLS weak-lensing tomographic study will further tighten these constraints because explicit treatment of the redshift dependence of cosmic shear more efficiently breaks the {Omega} {sub M}-{sigma}{sub 8} degeneracy. Combining the current results with the Wilkinson Microwave Anisotropy Probe 7 year (WMAP7) likelihood data, we obtain {Omega} {sub M} = 0.278 {+-} 0.018 and {sigma}{sub 8} = 0.815 {+-} 0.020.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Simplified metrics for the identification of the Madden-Julian oscillation in models

We propose simplified metrics to evaluate the fidelity with which the Madden–Julian oscillation (MJO) is simulated in climate models. These metrics are based on lag correlation analysis of principal component time series (PCs). The PCs are obtained by projecting simulated 20–100 day bandpass filtered daily outgoing longwave radiation onto the two leading empirical orthogonal functions of observed MJO variability. The simplified MJO metrics, the maximum positive correlation and time lag at which it occurs, provide consistent information relative to more complex diagnostics developed by the Madden–Julian Oscillation Working Group (CLIVAR MJOWG) and by Kim et al.

58 GEOSCIENCES↗
DOE OSTI2021

Material Identification Using Dual Energy X-ray Absorptiometry

Two implementations of dual energy X-ray absorptiometry were studied to identify materials using X-ray attenuation data taken with the Digital Radiography and Computed Tomography (DRCT) systems that were developed for the Recovered Chemical Materiel Directorate (RCMD). Maitrejean et al.’s approach utilizes eigen effects through Principal Component Analysis, while Osipov et al.’s approach proposed a physics-based method. Both approaches approximate mass attenuation coefficients of materials as a linear combination of basis functions (eigen effects) or physics-based equations. A set of coefficients {a 1 , a 2 , a 3 } or {B, D} were found by parameter optimization in EXCEL Solver. The identification parameters, {$\frac{a_{2}}{a_{1}}$, $\frac{a_{3}}{a_{1}}$} or estimated effective atomic number $\hat{Z}$ from {B, D}, were calculated to identify material of an aluminum 8 step wedge and a steel 8 step wedge in X-ray radiography images taken by a DRCT system. Maitrejean et al.’s approach was unable to provide reliable $\frac{a_{3}}{a_{1}}$ ratio values for identification of materials. Osipov et al.’s approach was found to be more robust in identify materials with a semi-empirical formula derived from test results in this study.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dimensional Reduction for Sampled Priors and Application to Photometric Redshift Distributions

A typical Bayesian inference on the values of some parameters of interest q from some data D involves running a Markov Chain (MC) to sample from the posterior $p$($q$,$n$|$D$) $\propto$ $\mathcal{L}$($D$|$q$,$n$)$p$(q)$p$($n$), where n are some nuisance parameters with a separable prior. In some cases, the nuisance parameters are high-dimensional, and their prior p(n) is itself defined only by a set of samples that have been drawn from some other MC. The MC for the posterior will typically require evaluation of p(n) at arbitrary values of n, i.e., one needs to provide a density estimator over the full n space from the provided samples. But the high dimensionality of n hinders both the density estimation and the efficiency of the MC for the posterior. We describe a solution to this problem: a linear compression of the n space into a much lower-dimensional space u, which projects away directions in n space that cannot appreciably alter $\mathcal{L}$. The algorithm for doing so is a slight modification to principal components analysis, and is less restrictive on p(n) than other proposed solutions to this issue. We demonstrate this “mode projection” technique using the analysis of 2-point correlation functions of weak lensing fields and galaxy density in the Dark Energy Survey, where n is a binned representation of the redshift distribution n(z) of the galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1985

Materials requirements for optical processing and computing devices

Devices for optical processing and computing systems are discussed, with emphasis on the materials requirements imposed by functional constraints. Generalized optical processing and computing systems are described in order to identify principal categories of requisite components for complete system implementation. Three principal device categories are selected for analysis in some detail: spatial light modulators, volume holographic optical elements, and bistable optical devices. The implications for optical processing and computing systems of the materials requirements identified for these device categories are described, and directions for future research are proposed.

Tanguay, A. R., Jr.↗
DOE OSTI2023

Exploring the dependence of gas cooling and heating functions on the incident radiation field with machine learning

ABSTRACT Gas cooling and heating functions play a crucial role in galaxy formation. But, it is computationally expensive to exactly compute these functions in the presence of an incident radiation field. These computations can be greatly sped up by using interpolation tables of pre-computed values, at the expense of making significant and sometimes even unjustified approximations. Here, we explore the capacity of machine learning to approximate cooling and heating functions with a generalized radiation field. Specifically, we use the machine learning algorithm XGBoost to predict cooling and heating functions calculated with the photoionization code cloudy at fixed metallicity, using different combinations of photoionization rates as features. We perform a constrained quadratic fit in metallicity to enable a fair comparison with traditional interpolation methods at arbitrary metallicity. We consider the relative importance of various photoionization rates through both a principal component analysis (PCA) and calculation of SHapley Additive exPlanation (shap) values for our XGBoost models. We use feature importance information to select different subsets of rates to use in model training. Our XGBoost models outperform a traditional interpolation approach at each fixed metallicity, regardless of feature selection. At arbitrary metallicity, we are able to reduce the frequency of the largest cooling and heating function errors compared to an interpolation table. We find that the primary bottleneck to increasing accuracy lies in accurately capturing the metallicity dependence. This study demonstrates the potential of machine learning methods such as XGBoost to capture the non-linear behaviour of cooling and heating functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The relationship between below average cognitive ability at age 5 years and the child’s experience of school at age 9

Background At age 5, while only embarking on their educational journey, substantial differences in children’s cognitive ability will already exist. The aim of this study was to examine the causal association between below average cognitive ability at age 5 years and child-reported experience of school and self-concept, and teacher-reported class engagement and emotional-behavioural function at age 9 years. Methods This longitudinal cohort study used data from 7,392 children in the Growing Up in Ireland Infant Cohort, who had completed the Picture Similarities and Naming Vocabulary subtests of the British Abilities Scales at age 5. Principal components analysis was used to produce a composite general cognitive ability score for each child. Children with a general cognitive ability score more than 1 standard deviation (SD) below the mean at age 5 were categorised as ‘Below Average Cognitive Ability’ (BACA), and those scoring above this as ‘Typical Cognitive Development’ (TCD). The outcomes of interest, measured at age 9, were child-reported experience of school, child’s self-concept, teacher-reported class engagement, and teacher-reported emotional behavioural function. Binary and multinomial logistic regression models were used to examine the association between BACA and these outcomes. Results Compared to those with TCD, those with BACA had significantly higher odds of never liking school [Adjusted odds ratio (AOR) 1.82, 95% CI 1.37–2.43, p < 0.001], of being picked on (AOR 1.27, 95% CI 1.09–1.48) and of picking on others (AOR 1.53, 95% CI 1.27–1.84). They had significantly higher odds of experiencing low self-concept (AOR 1.20, 95% CI 1.02–1.42) and emotional-behavioural difficulties (AOR 1.34, 95% CI 1.10–1.63, p = 0.003). Compared to those with TCD, children with BACA had significantly higher odds of hardly ever or never being interested, motivated and excited to learn (AOR 2.29, 95% CI 1.70–3.10). Conclusion Children with BACA at school-entry had significantly higher odds of reporting a negative school experience and low self-concept at age 9. They had significantly higher odds of having teacher-reported poor class engagement and problematic emotional-behavioural function at age 9. The findings of this study suggest BACA has a causal role in these adverse outcomes. Early childhood policy and intervention design should be cognisant of the important role of cognitive ability in school and childhood outcomes.

Bowe, Andrea K.↗
NASA NTRS1980

Eigenvector analysis of some observed and model-generated Climatological fields

An empirical orthogonal function analysis of monthly mean synoptic fields computed with a general circulation model and the corresponding observed climatological patterns is presented. A computer program was used to compute the principal components of the data matrices representing global maps of various observed and model-generated climatological fields. The variables analyzed were: 500mb geopotential height, 850 mb temperature, and sea level pressure. It is concluded that eigenvector analysis of observed and model-generated synoptic climatological fields is a useful diagnostic aid in the evaluation of global model simulation. The maps of any climatological variables are readily decomposed into a few orthogonal components which account for a large fraction of the total spatial variance, and the observed and model-generated patterns may then be compared in terms of these principal components of the original data matrices.

Dennis, M.↗
DOE OSTI2014

TU-F-18C-09: Mammogram Surveillance Using Texture Analysis for Breast Cancer Patients After Radiation Therapy

Purpose: To study the feasibility of applying cancer risk model established from treated patients to predict the risk of recurrence on follow-up mammography after radiation therapy for both ipsilateral and contralateral breast. Methods: An extensive set of textural feature functions was applied to a set of 196 Mammograms from 50 patients. 56 Mammograms from 28 patients were used as training set, 44 mammograms from 22 patients were used as test set and the rest were used for prediction. Feature functions include Histogram, Gradient, Co-Occurrence Matrix, Run-Length Matrix and Wavelet Energy. An optimum subset of the feature functions was selected by Fisher Coefficient (FO) or Mutual Information (MI) (up to top 10 features) or a method combined FO, MI and Principal Component (FMP) (up to top 30 features). One-Nearest Neighbor (1-NN), Linear Discriminant Analysis (LDA) and Nonlinear Discriminant Analysis (NDA) were utilized to build a risk model of breast cancer from the training set of mammograms at the time of diagnosis. The risk model was then used to predict the risk of recurrence from mammogram taken one year and three years after RT. Results: FPM with NDA has the best classification power in classifying the training set of the mammogram with lesions versus those without lesions. The model of FPM with NDA achieved a true positive (TP) rate of 82% compared to 45.5% of using FO with 1-NN. The best false positive (FP) rates were 0% and 3.6% in contra-lateral breast of 1-year and 3-years after RT, and 10.9% in ipsi-lateral breast of 3-years after RT. Conclusion: Texture analysis offers high dimension to differentiate breast tissue in mammogram. Using NDA to classify mammogram with lesion from mammogram without lesion, it can achieve rather high TP and low FP in the surveillance of mammogram for patient with conservative surgery combined RT.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

A variance-based decomposition and global sensitivity index method for uncertainty quantification: Application to retrieved ice cloud properties

This study develops a novel uncertainty quantification (UQ) method for cloud microphysical property retrievals using variance-based decomposition and global sensitivity index. In this UQ framework, empirical orthogonal function (EOF) analysis is applied to the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) ground-based observations, which are the inputs for the cloud retrieval studied here. The principal components (PCs) in the EOF expansion are parameterized as random input variables, and hence, the input dimension is greatly reduced (up to a factor of 50), allowing large ensemble of random samplings. The EOF expansion improves the accuracy of the uncertainty estimation by taking into account the cross correlations in the input data profiles. Here, this method enables a probabilistic representation of a retrieval process by adding normally distributed perturbations into PCs of sample means of input data profiles within a time window. Therefore, it effectively facilitates objective validation of climate models against cloud retrievals under a probabilistic framework for rigorous statistical inferences. Moreover, the variance-based global sensitivity index analysis, part of this method, attributes the output uncertainties to each individual source, thus providing directions for improving retrieval algorithms and observation strategies. For demonstration, we apply this method to quantify the uncertainties of the ARM program's baseline cloud retrieval algorithm for an ice cloud case observed at the Southern Great Plains site on 9 March 2000.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Indirect dark matter signatures in the cosmic dark ages. I. Generalizing the bound on s -wave dark matter annihilation from Planck results

Recent measurements of the cosmic microwave background (CMB) anisotropies by Planck provide a sensitive probe of dark matter annihilation during the cosmic dark ages, and specifically constrain the annihilation parameter f eff <σv>=m χ . Using new results (paper II) for the ionization produced by particles injected at arbitrary energies, we calculate and provide f eff values for photons and e + e - pairs injected at keV-TeV energies; the f eff value for any dark matter model can be obtained straightforwardly by weighting these results by the spectrum of annihilation products. This result allows the sensitive and robust constraints on dark matter annihilation presented by the Planck collaboration to be applied to arbitrary dark matter models with s-wave annihilation. Here, we demonstrate the validity of this approach using principal component analysis. As an example, we integrate over the spectrum of annihilation products for a range of Standard Model final states to determine the CMB bounds on these models as a function of dark matter mass, and demonstrate that the new limits generically exclude models proposed to explain the observed high-energy rise in the cosmic ray positron fraction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Statistical Shape Modeling of Cam Femoroacetabular Impingement

In this study, statistical shape modeling (SSM) was used to quantify three-dimensional (3D) variation and morphologic differences between femurs with and without cam femoroacetabular impingement (FAI). 3D surfaces were generated from CT scans of femurs from 41 controls and 30 cam FAI patients. SSM correspondence particles were optimally positioned on each surface using a gradient descent energy function. Mean shapes for control and patient groups were defined from the resulting particle configurations. Morphological differences between group mean shapes and between the control mean and individual patients were calculated. Principal component analysis was used to describe anatomical variation present in both groups. The first 6 modes (or principal components) captured statistically significant shape variations, which comprised 84% of cumulative variation among the femurs. Shape variation was greatest in femoral offset, greater trochanter height, and the head-neck junction. The mean cam femur shape protruded above the control mean by a maximum of 3.3 mm with sustained protrusions of 2.5-3.0 mm along the anterolateral head-neck junction and distally along the anterior neck, corresponding well with reported cam lesion locations and soft-tissue damage. This study provides initial evidence that SSM can describe variations in femoral morphology in both controls and cam FAI patients and may be useful for developing new measurements of pathological anatomy. SSM may also be applied to characterize cam FAI severity and provide templates to guide patient-specific surgical resection of bone.

statistical shape modeling↗
NASA NTRS1998

Using Neural Networks for Sensor Validation

This paper presents the results of applying two different types of neural networks in two different approaches to the sensor validation problem. The first approach uses a functional approximation neural network as part of a nonlinear observer in a model-based approach to analytical redundancy. The second approach uses an auto-associative neural network to perform nonlinear principal component analysis on a set of redundant sensors to provide an estimate for a single failed sensor. The approaches are demonstrated using a nonlinear simulation of a turbofan engine. The fault detection and sensor estimation results are presented and the training of the auto-associative neural network to provide sensor estimates is discussed.

Mattern, Duane L.↗
DOE OSTI2016

Simulation methodology for fuel assembly drop accident during handling

Fresh and irradiated fuel assembly handling is a routine part of Nuclear Power Plant activities. Fuel Assembly handling occurs during fresh fuel receipt, core refueling, spent fuel assemblies transfer from spent fuel pool to the storage containers, etc. Westinghouse, ENUSA and KEPCO-NF have collaborated in the technical area of Fuel Assembly Drop Accident (FADA) Simulation to update the approach for estimating the number of failed fuel rods as a consequence of a fresh or irradiated fuel assembly drop, which might occur during a handling accident in a PWR-type reactor. The proposed FADA simulation methodology starts with an assessment of fuel assembly impact velocity at the final stage of the drop in water or air. The vertical, horizontal or inclined drop directions are considered. This velocity is an input to determine a fuel assembly structural response during impact which defines the loading function for the fuel assembly components including spacer grids. The dynamic load on the fuel rod is calculated from the detailed spacer grid model and the corresponding loading function. Finally, the load from the spacer grid spring and dimples are applied to the fuel rod model to determine the cladding parameters for use in compliance to fuel rod cladding failure criteria. A combination of first principal, finite element analysis (FEA) and computational fluid dynamic (CFD) methods has been utilized to create a technically sound methodology. A number of specific and traditional tests for the fuel assembly and its components has been used to confirm that the proposed approaches which describe the physical phenomena related to the fuel assembly drop accident are reasonably conservative.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2013

Star/galaxy separation at faint magnitudes: Application to a simulated Dark Energy Survey

We address the problem of separating stars from galaxies in future large photometric surveys. We focus our analysis on simulations of the Dark Energy Survey (DES). In the first part of the paper, we derive the science requirements on star/galaxy separation, for measurement of the cosmological parameters with the Gravitational Weak Lensing and Large Scale Structure probes. These requirements are dictated by the need to control both the statistical and systematic errors on the cosmological parameters, and by Point Spread Function calibration. We formulate the requirements in terms of the completeness and purity provided by a given star/galaxy classifier. In order to achieve these requirements at faint magnitudes, we propose a new method for star/galaxy separation in the second part of the paper. We first use Principal Component Analysis to outline the correlations between the objects parameters and extract from it the most relevant information. We then use the reduced set of parameters as input to an Artificial Neural Network. This multi-parameter approach improves upon purely morphometric classifiers (such as the classifier implemented in SExtractor), especially at faint magnitudes: it increases the purity by up to 20% for stars and by up to 12% for galaxies, at i-magnitude fainter than 23.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗