Search NASA⌕ Search

SEARCH · Search NASA

Results for “statistical visualization”

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 217 records · Page 12

An affinity-structure database of helix-turn-helix: DNA complexes with a universal coordinate system

Molecular interactions between proteins and DNA molecules underlie many cellular processes, including transcriptional regulation, chromosome replication, and nucleosome positioning. Computational analyses of protein-DNA interactions rely on experimental data characterizing known protein-DNA interactions structurally and biochemically. While many databases exist that contain either structural or biochemical data, few integrate these two data sources in a unified fashion. Such integration is becoming increasingly critical with the rapid growth of structural and biochemical data, and the emergence of algorithms that rely on the synthesis of multiple data types to derive computational models of molecular interactions. We have developed an integrated affinity-structure database in which the experimental and quantitative DNA binding affinities of helix-turn-helix proteins are mapped onto the crystal structures of the corresponding protein-DNA complexes. This database provides access to: (i) protein-DNA structures, (ii) quantitative summaries of protein-DNA binding affinities using position weight matrices, and (iii) raw experimental data of protein-DNA binding instances. Critically, this database establishes a correspondence between experimental structural data and quantitative binding affinity data at the single basepair level. Furthermore, we present a novel alignment algorithm that structurally aligns the protein-DNA complexes in the database and creates a unified residue-level coordinate system for comparing the physico-chemical environments at the interface between complexes. Using this unified coordinate system, we compute the statistics of atomic interactions at the protein-DNA interface of helix-turn-helix proteins. We provide an interactive website for visualization, querying, and analyzing this database, and a downloadable version to facilitate programmatic analysis. Lastly, this database will facilitate the analysis of protein-DNA interactions and the development of programmatic computational methods that capitalize on integration of structural and biochemical datasets. The database can be accessed at http://ProteinDNA.hms.harvard.edu.

59 BASIC BIOLOGICAL SCIENCES↗

BioSig3D: High Content Screening of Three-Dimensional Cell Culture Models

BioSig3D is a computational platform for high-content screening of three-dimensional (3D) cell culture models that are imaged in full 3D volume. It provides an end-to-end solution for designing high content screening assays, based on colony organization that is derived from segmentation of nuclei in each colony. BioSig3D also enables visualization of raw and processed 3D volumetric data for quality control, and integrates advanced bioinformatics analysis. The system consists of multiple computational and annotation modules that are coupled together with a strong use of controlled vocabularies to reduce ambiguities between different users. It is a web-based system that allows users to: design an experiment by defining experimental variables, upload a large set of volumetric images into the system, analyze and visualize the dataset, and either display computed indices as a heatmap, or phenotypic subtypes for heterogeneity analysis, or download computed indices for statistical analysis or integrative biology. BioSig3D has been used to profile baseline colony formations with two experiments: (i) morphogenesis of a panel of human mammary epithelial cell lines (HMEC), and (ii) heterogeneity in colony formation using an immortalized non-transformed cell line. These experiments reveal intrinsic growth properties of well-characterized cell lines that are routinely used for biological studies. BioSig3D is being released with seed datasets and video-based documentation.

Science & Technology - Other Topics↗

Large-scale tearing-mode hazard function analysis with standard matched equilibrium reconstructions

The association between features from standard tokamak equilibrium reconstructions and the onset of n = 1 tearing modes (TMs) is analyzed at scale. The TM onset rate is directly modeled with a ‘hazard’ function which gives the expected number of onsets (per unit time spent) in a given equilibrium parameter region. In particular the different statistical modeling performance achieved for magnetics-only reconstructions and motional Stark effect (MSE) enhanced reconstructions is studied. It is observed that a better hazard model for the TM onset rate can be built with the MSE-enhanced equilibria compared to the matched magnetics-only situation. This advantage disappears if internal profile details are withheld from the matched analysis. Plausibility of the hazard function is further demonstrated with visualizations of global trends in the operational space, and time-traces from specific tokamak discharges. As a result, TMs typically degrade tokamak plasma performance and may lead to plasma termination, motivating this statistical study.

equilibrium↗

A generalized machine learning workflow to visualize mechanical discontinuity

Accurate detection and mapping of mechanical discontinuity in materials has widespread industrial and research applications. Herein, we developed a generalized machine-learning framework for visualizing single mechanical discontinuity embedded in material of any composition, velocity, density, porosity, and size with limited data. The proposed visualization of discontinuity requires accurate estimations of the length, location, and orientation of the embedded discontinuity by processing multipoint wave-transmission measurements. k-Wave simulator is used to create a large dataset of elastic waveforms recorded during multi-point wave-transmission measurements through materials containing single mechanical discontinuity. k-Wave simulator considers the wave attenuation, dispersion, and mode conversion in wave motion. Discrete wavelet transform (DWT) and statistical feature extraction are essential for data preprocessing prior to the data-driven model development. DWT also minimizes the effect of noise. Using hyper-parameter tuning and cross validation, gradient boosting regression can visualize the mechanical discontinuity with an accuracy of 0.85, in terms of coefficient of determination. A double-layered neural network-based regression has better performance with an accuracy of 0.95. Use of convolutional neural network converts the predictive task from a waveform processing to an image processing problem. Convolutional neural network achieved a generalization performance of 0.91. The proposed generalized workflow requires robust simulation of wave propagation, signal processing, feature engineering, and model evaluation. Sensors closest to the source and those located opposite the source are the most significant for the desired visualization. Notably, the sensors closest to the source capture the non-linear associations, whereas the sensor on the border opposite to the source capture the linear associations between the measured waveforms and the properties of the mechanical discontinuity.

42 ENGINEERING↗

Unsupervised anomaly detection in MeV ultrafast electron diffraction

MeV ultrafast electron diffraction (MUED) is a pump-probe technique used to study the dynamic structural evolution of materials. An ultrashort laser pulse triggers structural changes, which are then probed by an ultrashort relativistic electron beam. To overcome low signal-to-noise ratios, diffraction patterns are averaged over thousands of shots. However, shot-to-shot instabilities in the electron beam can distort individual patterns, introducing uncertainty. Improving MUED accuracy requires detecting and removing these anomalous patterns from large datasets. In this work, we developed a fully unsupervised methodology for the detection of anomalous diffraction patterns. Using a convolutional autoencoder, we calculate the reconstruction mean squared error of the diffraction patterns. Based on the statistical analysis of this error, we provide the user an estimation of the probability that the pattern is normal, which also allows a posterior visual inspection of the images that are difficult to classify. This method has been trained with only 100 diffraction patterns and tested on 1521 patterns, resulting in a false positive rate between 0.2% and 0.4%, with a training time of 10 s per image and a test time of about 1 s per image. Here, the proposed methodology can also be applied to other diffraction techniques in which large datasets are collected that include faulty images due to instrumental instabilities.

43 PARTICLE ACCELERATORS↗

Gamma, electron beam and X-ray irradiation effects on polymers in an advanced bone cement mixer device

The medical device industry has been investigating ways to increase the use of alternatives to cobalt-60 gamma radiation and ethylene-oxide gas sterilization, such as electron beam (E-beam) and X-ray radiation, due to regulatory and market pressures. One impediment to switching to E-beam or X-ray technology for sterilization is the lack of data on the effects of these radiation sources on medical device polymers. To provide such data this work considers irradiation and testing of a common single-use medical bone cement mixing system, the Stryker Advanced Cement Mixer (ACM®), that is composed of seven polymer materials. The ACM® devices considered here were processed to sterilization-relevant doses (15, 25, 50, and 70 kGy) using three radiation technologies: gamma, E-beam, and X-ray. The system and its polymer components were tested for product functionality, as well as mechanical and visual properties to determine how exposure effects may be influenced by radiation technology and dose level. We found that although there were instances of statistically significant differences in effects between the gamma-irradiated products and those irradiated with E-beam and X-ray, those effects were negligible in terms of retained functionality of the product and retained mechanical properties of the polymer components. Overall, results of this study demonstrate that, for of the effects studied, E-beam and X-ray are viable alternatives to cobalt-60 gamma radiation for sterilization of the polymer-based device investigated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Optimization Studies of Radiation Shielding for the PIP-II Project at Fermilab

The PIP-II project at Fermilab, which includes an 800-MeV superconducting LINAC, demands rigorous radiation shielding optimization to meet safety requirements. We updated the MARS geometry model to reflect new magnet and collimator designs and introduced high-resolution detector planes to better capture radiation field distributions. To overcome the significant computational demands, we implemented a well-known branching technique that drastically reduced simulation runtimes while maintaining statistical integrity. This was achieved through particle splitting and the application of Russian Roulette techniques. Additionally, new graphical tools were created to streamline data visualization and MARS code usability.

Makovec, Alajos [Fermilab] (ORCID:0000000286157492↗

Alternatives to Contour Visualizations for Power Systems Data

Electrical grids are geographical and topological structures whose voltage states are challenging to represent accurately and efficiently for visual analysis. The current common practice is to use colored contour maps, yet these can misrepresent the data. We examine the suitability of four alternative visualization methods for depicting voltage data in a geographically dense distribution system - Voronoi polygons, H3 tessellations, S2 tessellations, and a network-weighted contour map. We find that Voronoi tessellations and network-weighted contour maps more accurately represent the statistical distribution of the data than regular contour maps.

MATHEMATICS AND COMPUTING,POWER TRANSMISSION AND D↗

Exploration of thermal counterflow in He II using particle tracking velocimetry

Flow visualization using particle image velocimetry (PIV) and particularly particle tracking velocimetry (PTV) has been applied to thermal counterflow in He II for nearly two decades now, but the results remain difficult to interpret because tracer particle motion can be influenced by both the normal fluid and superfluid components of He II as well as the quantized vortex tangle. For instance, in one early experiment it was observed (using PTV) that tracer particles move at the normal fluid velocity v n , while in another it was observed (using PIV) that particles move at v n /2. Besides the different visualization methods, the range of applied heat flux investigated by these experiments differed by an order of magnitude. To resolve this apparent discrepancy and explore the statistics of particle motion in thermal counterflow, we apply the PTV method to a wide range of heat flux at a number of different fluid temperatures. In our analysis, we introduce a scheme for analyzing the velocity of particles presumably moving with the normal fluid separately from those presumably influenced by the quantized vortex tangle. Our results show that for lower heat flux there are two distinct peaks in the streamwise particle velocity probability density function (PDF), with one centered at the normal fluid velocity v n (named G2 for convenience) while the other is centered near v n /2 (G1). For higher heat flux there is a single peak centered near v n /2 (G3). Using our separation scheme, we show quantitatively that there is no size difference between the particles contributing to G1 and G2. We also show that nonclassical features of the transverse particle velocity PDF arise entirely from G1, while the corresponding PDF for G2 exhibits the classical Gaussian form. The G2 transverse velocity fluctuation, backed up by second sound attenuation in decaying counterflow, suggests that large-scale turbulence in the normal fluid is absent from the two-peak region. We offer a brief discussion of the physical mechanisms that may be responsible for our observations, revealing that G1 velocity fluctuations may be linked to fluctuations of quantized vortex line velocity, and suggest a number of numerical simulations that may reveal the underlying physics in detail.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multi-Level Anomaly Detection on Time-Varying Graph Data

This work presents a novel modeling and analysis framework for graph sequences which addresses the challenge of detecting and contextualizing anomalies in labelled, streaming graph data. We introduce a generalization of the BTER model of Seshadhri et al. by adding flexibility to community structure, and use this model to perform multi-scale graph anomaly detection. Specifically, probability models describing coarse subgraphs are built by aggregating probabilities at finer levels, and these closely related hierarchical models simultaneously detect deviations from expectation. This technique provides insight into a graph's structure and internal context that may shed light on a detected event. Additionally, this multi-scale analysis facilitates intuitive visualizations by allowing users to narrow focus from an anomalous graph to particular subgraphs or nodes causing the anomaly. For evaluation, two hierarchical anomaly detectors are tested against a baseline Gaussian method on a series of sampled graphs. We demonstrate that our graph statistics-based approach outperforms both a distribution-based detector and the baseline in a labeled setting with community structure, and it accurately detects anomalies in synthetic and real-world datasets at the node, subgraph, and graph levels. To illustrate the accessibility of information made possible via this technique, the anomaly detector and an associated interactive visualization tool are tested on NCAA football data, where teams and conferences that moved within the league are identified with perfect recall, and precision greater than 0.786.

anomaly detection↗

SU-E-I-21: Deformation Mapping and Shape Prediction with 3D Tumor Volume Morphing

Purpose: To compare projection-based versus global correction that compensate for deadtime count loss in SPECT/CT images. Methods: SPECT/CT images of an IEC phantom (2.3GBq 99mTc) with ∼10% deadtime loss containing the 37mm (uptake 3), 28 and 22mm (uptake 6) spheres were acquired using a 2 detector SPECT/CT system with 64 projections/detector and 15 s/projection. The deadtime, Ti and the true count rate, Ni at each projection, i was calculated using the monitor-source method. Deadtime corrected SPECT were reconstructed twice: (1) with projections that were individually-corrected for deadtime-losses; and (2) with original projections with losses and then correcting the reconstructed SPECT images using a scaling factor equal to the inverse of the average fractional loss for 5 projections/detector. For both cases, the SPECT images were reconstructed using OSEM with attenuation and scatter corrections. The two SPECT datasets were assessed by comparing line profiles in xyplane and z-axis, evaluating the count recoveries, and comparing ROI statistics. Higher deadtime losses (up to 50%) were also simulated to the individually corrected projections by multiplying each projection i by exp(-a*Ni*Ti), where a is a scalar. Additionally, deadtime corrections in phantoms with different geometries and deadtime losses were also explored. The same two correction methods were carried for all these data sets. Results: Averaging the deadtime losses in 5 projections/detector suffices to recover >99% of the loss counts in most clinical cases. The line profiles (xyplane and z-axis) and the statistics in the ROIs drawn in the SPECT images corrected using both methods showed agreement within the statistical noise. The count-loss recoveries in the two methods also agree within >99%. Conclusion: The projection-based and the global correction yield visually indistinguishable SPECT images. The global correction based on sparse sampling of projections losses allows for accurate SPECT deadtime loss correction while keeping the study duration reasonable.

60 APPLIED LIFE SCIENCES↗

SU-F-R-25: Automatic Identification of Suspicious Recurrent/residual Disease Regions After Prostatectomy

Purpose: To automatically identify and outline suspicious regions of recurrent or residual disease in the prostate bed using Dynamic Contrast Enhanced-MRI (DCE-MRI) in patients after prostatectomy. Methods: Twenty-two patients presenting for salvage radiotherapy and with identified Gross Tumor Volume (GTV) in the prostate bed were retrospectively analyzed. The MRI data consisted of Axial T2weighted-MRI (T2w) of the pelvis: resolution 1.25×1.25×2.5 mm; Field of View (FOV): 320×320 mm; slice thickness=2.5mm; 72 slices; and Dynamic Contrast Enhanced MRI (DCE-MRI)–12 series of T1w with identical spatial resolution to T2w and at 30–34s temporal resolution. Unsupervised pattern recognition was used to decompose the 4D DCE data as the product W.H of weights W of k patterns H. A well-perfused pattern Hwp was identified and the weight map Wwp associated to Hwp was used to delineate suspicious volumes. Threshold of Wwp set at mean(Wwp)+S*std(Wwp), S=1,1.5,2 and 2.5 defined four volumes labeled as DCE1.0 to DCE2.5. These volumes were displayed on T2w and, along with GTV, were correlated with the highest pre-treatment PSA values, and with pharmacokinetic analysis constants. Results: GTV was significantly correlated with DCE2.0(ρ= 0.60, p<0.003), and DCE 2.5 (ρ=0.58, p=0.004)). Significant correlation was found between highest pre-treatment PSA and GTV(ρ=0.42, p<0.049), DCE2.0(ρ= 0.52, p<0.012), and DCE 2.5 (ρ=0.67, p<<0.01)). Kruskal-Wallis analysis showed that Ktrans median value was statistically different between non-specific prostate bed tissue NSPBT and both GTV (p<<0.001) and DCE2.5 (p<<0.001), but while median Ve was statistically different between DCE2.5 and NSPBT (p=0.002), it was not statistically different between GTV and NSPBT (p=0.054), suggesting that automatic volumes capture more accurately the area of malignancy. Conclusion: Software developed for identification and visualization of suspicions regions in DCE-MRI from post-prostatectomy patients has been validated by PSA and pharmacokinetic constants analysis showing that it generates clinically relevant volumes.

60 APPLIED LIFE SCIENCES↗

Using Stage Level Microseismic Analysis to Gain Insight into Fracture Efficiency and Completion Effectiveness

Microseismic event density heat maps at sequential time intervals during the treatment of individual completion stages have been developed to better understand fracture behavior on a stage level. Animating these heat maps through time on a single stage provides fascinating visual insights into the dynamic sequencing of populations of events, event timing and their likely genesis, providing a new way of visualizing and interpreting microseismic data sets. This technique leverages heat-map benefits from URTEC 2674376, excluding microseismic event clusters that may obscure or skew statistical analysis of the event densities, thus allowing for more accurate assessment of geologic features, stimulated rock volume, near-field concentration of events due to varying cluster spacing, fracture efficiency and stress shadowing. Far-field microseismic events prior to ramp up of treatment pressure are clearly indicated from this technique that cannot be observed with standard heat maps with certain completion designs. Unambiguous changes in microseismic event character from these far-field pre-ramp up events to those near-field events occurring at the onset of pressure are also observed. Applying a stage by stage approach on packages of wells and different locations provides robust statistical conclusions pertaining to reservoir behavior during well completion.

Wang, Iris↗

A multi-level anomaly detection algorithm for time-varying graph data with interactive visualization

This work presents a novel modeling and analysis framework for graph sequences which addresses the challenge of detecting and contextualizing anomalies in labelled, streaming graph data. We introduce a generalization of the BTER model of Seshadhri et al. by adding flexibility to community structure, and use this model to perform multi-scale graph anomaly detection. Specifically, probability models describing coarse subgraphs are built by aggregating node probabilities, and these related hierarchical models simultaneously detect deviations from expectation. This technique provides insight into a graph's structure and internal context that may shed light on a detected event. Additionally, this multi-scale analysis facilitates intuitive visualizations by allowing users to narrow focus from an anomalous graph to particular subgraphs or nodes causing the anomaly. For evaluation, two hierarchical anomaly detectors are tested against a baseline Gaussian method on a series of sampled graphs. We demonstrate that our graph statistics-based approach outperforms both a distribution-based detector and the baseline in a labeled setting with community structure, and it accurately detects anomalies in synthetic and real-world datasets at the node, subgraph, and graph levels. Furthermore, to illustrate the accessibility of information made possible via this technique, the anomaly detector and an associated interactive visualization tool are tested on NCAA football data, where teams and conferences that moved within the league are identified with perfect recall, and precision greater than 0.786.

97 MATHEMATICS AND COMPUTING↗

Visual HPC Workflows for the Analysis of System Dynamics Models

Visual analytics supported by high performance computing (HPC) accelerates and enhances the discovery, exploration, and analysis of causal patterns in complex system dynamics (SD) models. We present a suite of visualization-assisted ensemble-based techniques for hypothesis generation and testing, and for sensitivity analysis. By employing HPC to provide parallel, on-demand simulation of SD models, one can “steer” an ensemble of simulated scenarios in real time as one first formulates and then informally tests those hypotheses: this provides rapid feedback for analysts to refine their understanding of the causal relationships emergent from a model. Such understandings can be followed and augmented by rigorous application of statistical methods, namely global variance-based sensitivity analysis, Monte-Carlo filtering, adaptive regional sensitivity analysis, and self-organized maps: here timely computation relies on HPC, while effective presentation emphasizes high-dimensional multivariate data visualization. Immersive visualization in virtual 3D environments provides an excellent adjunct to the traditional 2D graphics typically used for SD models, as it generates an embodied understanding of model behavior and facilitates an active, collaborative critique of model structure and output. Finally, we summarize prospects for HPC-enabled visual analytics applied to SD modeling.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Unsupervised Machine Learning for Exploratory Data Analysis of Exoplanet Transmission Spectra

Abstract Transit spectroscopy is a powerful tool for decoding the chemical compositions of the atmospheres of extrasolar planets. In this paper, we focus on unsupervised techniques for analyzing spectral data from transiting exoplanets. After cleaning and validating the data, we demonstrate methods for: (i) initial exploratory data analysis, based on summary statistics (estimates of location and variability); (ii) exploring and quantifying the existing correlations in the data; (iii) preprocessing and linearly transforming the data to its principal components; (iv) dimensionality reduction and manifold learning; (v) clustering and anomaly detection; and (vi) visualization and interpretation of the data. To illustrate the proposed unsupervised methodology, we use a well-known public benchmark data set of synthetic transit spectra. We show that there is a high degree of correlation in the spectral data, which calls for appropriate low-dimensional representations. We explore a number of different techniques for such dimensionality reduction and identify several suitable options in terms of summary statistics, principal components, etc. We uncover interesting structures in the principal component basis, namely well-defined branches corresponding to different chemical regimes of the underlying atmospheres. We demonstrate that those branches can be successfully recovered with a K-means clustering algorithm in a fully unsupervised fashion. We advocate for lower-dimensional representations of the spectroscopic data in terms of the main principal components, in order to reveal the existing structure in the data and quickly characterize the chemical class of a planet.

Matchev, Konstantin T. (ORCID:0000000341829096)↗

Determining Patterns from Radiation Portal Monitor Data: Enabling Data Insight with Visual and Interactive Exploratory Data Analysis

Smuggling of nuclear and radiological material is a complex problem that will need to be analyzed from a variety of angles in order to effectively combat something so rare, yet so hazardous. One crucial tool in the effort to detect and deter smuggling of rad/nuc materials is the radiation portal monitor (RPM)—machines that scan passengers at airports, vehicles and trains at border crossings, and shipping containers at ports throughout the world. While data from RPMs are regularly analyzed at the individual lane level to ensure monitor health and good operation, aggregate analysis of the data has not gone far beyond summary statistics such as number of occupancies and alarm rates. In particular, there are opportunities to carry out interactive exploratory data analysis with the help of tools like the R package Shiny and methods that lend themselves to the visual display of spatiotemporal correlations such as spatial PCA and detrended cross-correlation analysis. By doing so, we discover interesting spatiotemporal patterns in the data that point to specific events or anomalous moments in time that warrant further analysis.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Conductivity Features of Nanoislet Metal Films

Structural percolation transitions in ultrathin metal films at thermally activated granularity was studied. The degree of granularity of an Au/GGG film is shown to increase with decreasing its effective thickness. The film structure peculiarities are visualized upon the percolation transition via scanning electron microscopy. The feasible causes, mechanisms and models of activation conductivity in islet metal nanofilms are discussed and the statistical width distribution of inter-islet gaps is considered, as well.

77 NANOSCIENCE AND NANOTECHNOLOGY↗