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 649 records · Page 36

FREDA: A Web Application for the Processing, Analysis, and Visualization of Fourier‐Transform Mass Spectrometry Data

The high-resolution measurement capability of Fourier-transform mass spectrometry (FT-MS) has made it a necessity for exploring the molecular composition of complex organic mixtures, like soil, plant, aquatic, and petroleum samples. This demand has driven a need for informatics tools to explore and analyze FT-MS data in a robust and reproducible manner. FREDA is an interactive web application developed to enable spectrometrists to format, process, and explore their FT-MS data without the need for statistical programming expertise. FREDA was built to explore outputs from a molecular identification tool, like CoreMS, and provide a suite of methods to filter data, compute chemical properties of peaks, statistically compare samples and groups of samples, conduct exploratory data analysis, and download the results with a report detailing all steps conducted. To demonstrate the utility of FREDA, an example analysis was conducted using FT-MS data from a soil microbiology study of samples collected in two different soil depths at the Sphagnum bog forest north of Grand Rapids, Minnesota. Differences between the two depths are observed using Kendrick, Gibbs free energy, and van Krevelen plots. G-tests are used to quantify a significant difference between the groups. All analyses and plotting are conducted using only the FREDA application. FREDA is an open-source and readily available web application that allows users to explore and make statistically valid conclusions about their FT-MS data. The application is available online (https://map.emsl.pnnl.gov/app/freda) with a tutorial web series (https://youtu.be/k5HLE2kNSBY?si=yB6sGoyvzxrFf5MP) and freely accessible code on Github (https://github.com/EMSL-Computing/FREDA).

47 OTHER INSTRUMENTATION↗

Classification of Cloud Particle Imagery from Aircraft Platforms Using Convolutional Neural Networks

Abstract A vast amount of ice crystal imagery exists from a variety of field campaign initiatives that can be utilized for cloud microphysical research. Here, nine convolutional neural networks are used to classify particles into nine regimes on over 10 million images from the Cloud Particle Imager probe, including liquid and frozen states and particles with evidence of riming. A transfer learning approach proves that the Visual Geometry Group (VGG-16) network best classifies imagery with respect to multiple performance metrics. Classification accuracies on a validation dataset reach 97% and surpass traditional automated classification. Furthermore, after initial model training and preprocessing, 10 000 images can be classified in approximately 35 s using 20 central processing unit cores and two graphics processing units, which reaches real-time classification capabilities. Statistical analysis of the classified images indicates that a large portion (57%) of the dataset is unusable, meaning the images are too blurry or represent indistinguishable small fragments. In addition, 19% of the dataset is classified as liquid drops. After removal of fragments, blurry images, and cloud drops, 38% of the remaining ice particles are largely intersecting the image border (≥10% cutoff) and therefore are considered unusable because of the inability to properly classify and dimensionalize. After this filtering, an unprecedented database of 1 560 364 images across all campaigns is available for parameter extraction and bulk statistics on specific particle types in a wide variety of storm systems, which can act to improve the current state of microphysical parameterizations.

54 ENVIRONMENTAL SCIENCES↗

GOES-R AWG GLM Val Tool Development

We are developing tools needed to enable the validation of the Geostationary Lightning Mapper (GLM). In order to develop and test these tools, we have need of a robust, high-fidelity set of GLM proxy data. Many steps have been taken to ensure that the proxy data are high quality. LIS is the closest analog that exists for GLM, so it has been used extensively in developing the GLM proxy. We have verified the proxy data both statistically and algorithmically. The proxy data are pixel (event) data, called Level 1B. These data were then clustered into flashes by the Lightning Cluster-Filter Algorithm (LCFA), generating proxy Level 2 data. These were then compared with the data used to generate the proxy, and both the proxy data and the LCFA were validated. We have developed tools to allow us to visualize and compare the GLM proxy data with several other sources of lightning and other meteorological data (the so-called shallow-dive tool). The shallow-dive tool shows storm-level data and can ingest many different ground-based lightning detection networks, including: NLDN, LMA, WWLLN, and ENTLN. These are presented in a way such that it can be seen if the GLM is properly detecting the lightning in location and time comparable to the ground-based networks. Currently in development is the deep-dive tool, which will allow us to dive into the GLM data, down to flash, group and event level. This will allow us to assess performance in comparison with other data sources, and tell us if there are detection, timing, or geolocation problems. These tools will be compatible with the GLM Level-2 data format, so they can be used beginning on Day 0.

Bateman, Monte↗

Methyl-CpG island-associated genome signature tags

Disclosed is a method for analyzing the organismic complexity of a sample through analysis of the nucleic acid in the sample. In the disclosed method, through a series of steps, including digestion with a type II restriction enzyme, ligation of capture adapters and linkers and digestion with a type IIS restriction enzyme, genome signature tags are produced. The sequences of a statistically significant number of the signature tags are determined and the sequences are used to identify and quantify the organisms in the sample. Various embodiments of the invention described herein include methods for using single point genome signature tags to analyze the related families present in a sample, methods for analyzing sequences associated with hyper- and hypo-methylated CpG islands, methods for visualizing organismic complexity change in a sampling location over time and methods for generating the genome signature tag profile of a sample of fragmented DNA.

59 BASIC BIOLOGICAL SCIENCES↗

Data Integrity Challenges in NASA Giovanni

The Geospatial Interactive Online Visualization ANd aNalysis Infrastructure (Giovanni) is an online tool developed by the NASA Goddard Earth Sciences (GES) Data and Information Services Center (DISC), one of 12 NASA Science Mission Directorate Data Centers (DAACs) to analyze and visualize NASA remote sensing and model data without downloading data and software. As of this writing, over 2000 Earth satellite and model variables are available in Giovanni, including several well-known NASA satellite missions (e.g., TRMM, GPM) and projects (e.g., MERRA-2, GPCP). There are twenty-two plots provided by Giovanni that can be used to analyze, compare, and explore Earth data across disciplines. Results can be shared with colleagues and downloaded for further analysis. Giovanni has helped publish over 3000 referral papers over the years. As open science policies roll in, data integrity has become a major challenge for Giovanni and other tools. For integrity, both data and workflows must be transparent. FAIR-compliant data, including input, intermediate, and result products, as well as their associated statistics, metadata, and information, are needed. The NASA Data Product Development Guide for Data Producers provides a key resource on how to develop FAIR-compliant data products. Data quality information is also needed from data producers and analysis services like Giovanni. The workflow part is quite challenging and requires workflow management improvements, such as recording workflows and making them available to users. In this presentation, we will discuss the data integrity challenges in Giovanni.

data analysis, visualization↗

Band Excitation in Scanning Probe Microscopy: Recognition and Functional Imaging

Field confinement at the junction between a biased scanning probe microscope s (SPM) tip and solid surface enables local probing of various bias-induced transformations such as polarization switching, ionic motion, or electrochemical reactions to name a few. The nanoscale size of the biased region is smaller or comparable to features like grain boundaries and dislocations, potentially allows for the study of kinetics and thermodynamics at the level of a single defect. In contrast to classical statistically averaged approaches, this allows one to link structure to functionality and deterministically decipher associated mesoscopic and atomistic mechanisms. Furthermore, this type of information can serve as a fingerprint of local material functionality, allowing for local recognition imaging. Here, current progress in multidimensional SPM techniques based on band-excitation time and voltage spectroscopies is illustrated, including discussions on data acquisition, dimensionality reduction, and visualization along with future challenges and opportunities for the field.

Jesse, Stephen↗

The Effects of Gravity on Wrinkled Laminar Flames

The effects of gravity are significant to the dynamics of idealized unconfined open premixed flames. Moderate to low turbulence Reynolds number flames, i.e., wrinkled laminar flames, of various unconfined geometries have been used extensively for investigating fundamental processes of turbulent flame propagation and to validate theoretical models. Without the wall constraints, the flames are free to expand and interact with surrounding ambient air. The flow field in which the flame exists is determined by a coupling of burner geometry, flame orientation and the gravity field. These complex interactions raise serious questions regarding the validity of comparing the experimental data of open flames with current theoretical and numerical models that do not include the effects of gravity nor effects of the larger aerodynamic flowfield. Therefore, studies of wrinkled laminar flame in microgravity are needed for a better understanding of the role of gravity on flame characteristics such as the orientation, mean aerodynamics stretch, flame wrinkle size and burning rate. Our approach to characterize and quantify turbulent flame structures under microgravity is to exploit qualitative and quantitative flow visualization techniques coupled with video recording and computer controlled image analysis technologies. The experiments will be carried out in the 2.2 second drop tower at the NASA Lewis Research Center. The longest time scales of typical wrinkled laminar flames in the geometries considered here are in the order of 10 msec. Hence, the duration of the drop is sufficient to obtain the amount of statistical data necessary for characterize turbulent flame structures.

Kostiuk, Larry W.↗

Suite of Comprehensive Rapid Analysis Tools for Estimation (SOCRATES) - 19076

Long-term monitoring and remedy implementation require ongoing assessment and analytics, much of which relies on intensive and time-consuming data integration and engineering, data manipulation, and 'big-data' analytics. Many of these analyses, such as groundwater flow direction, contaminant rate attenuation, and non-parametric statistics, are spatial-temporal in nature and involve coalescing disparate data sets from a variety of sources with measurement frequencies from milliseconds to years, as well as total data volumes often exceeding the petabyte scale. Web-based analytical tools are advantageous for ongoing real-time or near-real-time ingestion of data from various sources for general data access and visualization, and also allow for complex, data rich calculations by leveraging expansive networks and virtually unlimited compute power. A web-based framework offers a natural architecture for such endeavors, as access to compute nodes through auto-scaling and dynamic frameworks enables speedy data manipulation and unlimited access to any data set that is made available via an http request. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Z-checker: A framework for assessing lossy compression of scientific data

Because of the vast volume of data being produced by today’s scientific simulations and experiments, lossy data compressor allowing user-controlled loss of accuracy during the compression is a relevant solution for significantly reducing the data size. However, lossy compressor developers and users are missing a tool to explore the features of scientific data sets and understand the data alteration after compression in a systematic and reliable way. To address this gap, we have designed and implemented a generic framework called Z-checker. On the one hand, Z-checker combines a battery of data analysis components for data compression. On the other hand, Z-checker is implemented as an open-source community tool to which users and developers can contribute and add new analysis components based on their additional analysis demands. In this article, we present a survey of existing lossy compressors. Then, we describe the design framework of Z-checker, in which we integrated evaluation metrics proposed in prior work as well as other analysis tools. Specifically, for lossy compressor developers, Z-checker can be used to characterize critical properties (such as entropy, distribution, power spectrum, principal component analysis, and autocorrelation) of any data set to improve compression strategies. For lossy compression users, Z-checker can detect the compression quality (compression ratio and bit rate) and provide various global distortion analysis comparing the original data with the decompressed data (peak signal-to-noise ratio, normalized mean squared error, rate–distortion, rate-compression error, spectral, distribution, and derivatives) and statistical analysis of the compression error (maximum, minimum, and average error; autocorrelation; and distribution of errors). Z-checker can perform the analysis with either coarse granularity (throughout the whole data set) or fine granularity (by user-defined blocks), such that the users and developers can select the best fit, adaptive compressors for different parts of the data set. Z-checker features a visualization interface displaying all analysis results in addition to some basic views of the data sets such as time series. To the best of our knowledge, Z-checker is the first tool designed to assess lossy compression comprehensively for scientific data sets.

97 MATHEMATICS AND COMPUTING↗

Federal Emergency Management Agency Mobile Housing Unit Running Gear Testing (Final Report)

The Federal Emergency Management Agency (FEMA) stages manufactured housing units (MHUs) at key locations to be used during disasters when other housing accommodations cannot be arranged for individuals or families. Due to the unpredictability of natural disasters and the limitations of the manufactured housing industry to produce large numbers of MHUs in very short periods of time, it is necessary for FEMA to store a large number of units for many months or years. With these long storage periods come degradation issues with frames, tires, bearing, and braking systems. These types of issues are also of interest to the U.S. Department of Transportation (DOT) Federal Motor Carrier Safety Administration (FMCSA) because regulated carriers normally transport these units and 49 CFR 393.75 Tires specifies tire loading restrictions for manufactured homes. To compare the running gear (tires, brakes, drums, bearings, axles, springs, and spring mounting hardware) currently specified in the 2014 Rugged Based Performance Requirements (RBPR) and the newly recommended running gear specified in the upcoming revision to the RBPR, the Oak Ridge National Laboratory (ORNL) outlined a Longevity Test Plan, which focused on the assessment of component longevity (failure rates by component type). The Test Plan also studied the forces imparted to the MHU from the running gear from interaction with the roadway for both the running gear in the 2014 RBPR and the running gear specified in the upcoming RBPR revision. Since the assessment of component longevity and failures are also issues of concern for FMCSA, both agencies collaborated to support the research presented in this report. The tests and data collection were conducted at the Navistar Proving Grounds (NPG) in New Carlisle (near South Bend), Indiana, during the May-November 2018 period. Four alternatives for the suspension system of a typical three-bedroom MHU were considered in this test. Alternative 0 was the MHU that was selected from the FEMA staging area in Selma, Alabama. Empirical data had shown that, on average, this alternative experiences one tire failure every 150 miles (or 0.0067 tire failures per mile traveled). Before testing, Alternative 0 was equipped with similar running gear hardware, although in new condition and installed following the manufacturer’s specifications. This became the Baseline Unit, or Alternative 1. Another unit, also selected from the FEMA staging area in Selma, Alabama, was equipped with better axles (i.e., axles that were rated for 7,000lbs as opposed to 6,000lbs for the Baseline Unit) and better tires. This became the Upgrade Unit, or Alternative 2. At the end of the Longevity Test, the Upgrade Unit was mounted with radial tires. This became Alternative 3.Alternative 1 (Baseline Unit) had, on average, one tire failure every 251 miles (or 0.0040 tire failures per mile traveled). So, by simply following the correct manufacturer’s specifications for the assembly of the transportation system, a reduction of 41% in the per mile tire failure was achieved. In the same test, Alternative 2 (Upgrade Unit) showed, on average, one tire failure every 728 miles traveled (or 0.0014 tire failure s per mile traveled). When the Upgrade Unit was mounted with radial tires in the last part of the test, no tire failures were observed for the 2,400 miles tested. That is, in the case of the radial tires (Alternative 3) the tire failure was less than one per 2,400 miles traveled (or less than 0.0004 tire failures per mile traveled). The data collected for the Longevity Test showed that the maximum accelerations were registered in all cases at the center of the MHU in the vertical direction. The overall maximum was registered by the Baseline Unit, Alternative 1. If that value is normalized to 100%, then the Upgrade Unit (Alternative 2) showed, on average, 81% at the same location, and Alternative 3 (radial-tire test) 59%. That is, the upgraded transportation system showed a reduction of 19% and 41% in the vertical accelerations transmitted to MHU at its center when regular and radial tires, respectively, were used. As part of the Test Protocol, NPG conducted a visual inspection of the MHUs and transportation system at the end of each test day. ORNL reviewed all the inspection reports in an attempt to determine if there was any difference in the damage that the Baseline and Upgrade MHUs showed during the test. The Upgrade Unit always showed lower accelerations transmitted to the MHU than the Baseline Unit. In most cases, these differences were statistically significant. The FEMA inspection report, conducted by a FEMA certified inspector at the end of the tests, included a shorter list of damaged items for the Baseline Unit than for the Upgrade Unit, although similar type of damages was described for both units. This may be an indication that, although lower than those of the Baseline Unit, the accelerations transmitted by the Upgrade Unit were above the threshold at which MHU elements start to fail.

42 ENGINEERING↗

Polarbrdf: A General Purpose Python Package for Visualization Quantitative Analysis of Multi-Angular Remote Sensing Measurements

The Bidirectional Reflectance Distribution Function (BRDF) is a fundamental concept for characterizing the reflectance property of a surface, and helps in the analysis of remote sensing data from satellite, airborne and surface platforms. Multi-angular remote sensing measurements are required for the development and evaluation of BRDF models for improved characterization of surface properties. However, multi-angular data and the associated BRDF models are typically multidimensional involving multi-angular and multi-wavelength information. Effective visualization of such complex multidimensional measurements for different wavelength combinations is presently somewhat lacking in the literature, and could serve as a potentially useful research and teaching tool in aiding both interpretation and analysis of BRDF measurements. This article describes a newly developed software package in Python (PolarBRDF) to help visualize and analyze multi-angular data in polar and False Color Composite (FCC) forms. PolarBRDF also includes functionalities for computing important multi-angular reflectance/albedo parameters including spectral albedo, principal plane reflectance and spectral reflectance slope. Application of PolarBRDF is demonstrated using various case studies obtained from airborne multi-angular remote sensing measurements using NASA's Cloud Absorption Radiometer (CAR). Our visualization program also provides functionalities for untangling complex surface/atmosphere features embedded in pixel-based remote sensing measurements, such as the FCC imagery generation of BRDF measurements of grasslands in the presence of wild fire smoke and clouds. Furthermore, PolarBRDF also provides quantitative information of the angular distribution of scattered surface/atmosphere radiation, in the form of relevant BRDF variables such as sunglint, hotspot and scattering statistics.

CAR↗

Probing and Tuning Strain‐Localized Exciton Emission in 2D Material Bubbles at Room Temperature

In monolayer transition metal dichalcogenides bubbles-nanoscale deformations typically exhibiting a dome-like shape-Excitons are confined by the strain effect, which exhibits extraordinary emission properties, such as single photon generation, enhanced light emission, and spectrally tunable excitonic states. While the strain profiles of these bubbles are extensively studied, this work provides an approach 1) to directly visualize the associated exciton properties in bubbles formed in WSe2 monolayer, revealing an intrinsic emission wavelength shift of ≈40 nm, and 2) actively modify local strain, enabling further exciton emission tuning over a range of 50 nm. These are achieved by emission mapping and nanoindentation using a dielectric near-field probe, which enables the detection of local emission spectra and emission lifetimes within individual bubbles. Statistical analysis of 67 bubbles uncovers an emission wavelength distribution centered around 780 nm. Furthermore, saturation behavior in the power-dependent studies and the associated lifetime change reveal the localized nature of the strain-induced states. These findings provide direct insights into the strain-localized emission dynamics in bubbles and establish a robust framework for non-destructive, reversible, and predictable nanoscale emission control, presenting a potential avenue for developing next-generation tunable quantum optical sources.

2D materials bubbles↗

External Threat Risk Assessment Algorithm (ExTRAA)

Two risk assessment algorithms and philosophies have been augmented and combined to form a new algorit hm, the External Threat Risk Assessment Algorithm (ExTRAA), that allows for effective and statistically sound analysis of external threat sources in relation to individual attack methods . In addition to the attack method use probability and the attack method employment consequence, t he concept of defining threat sources is added to the risk assessment process. Sample data is tabulated and depicted in radar plots and bar graphs for algorithm demonstration purposes. The largest success of ExTRAA is its ability to visualize the kind of r isk posed in a given situation using the radar plot method.

97 MATHEMATICS AND COMPUTING↗

Content analysis to detect high stress in oral interviews and text documents

A system of interrogation to estimate whether a subject of interrogation is likely experiencing high stress, emotional volatility and/or internal conflict in the subject's responses to an interviewer's questions. The system applies one or more of four procedures, a first statistical analysis, a second statistical analysis, a third analysis and a heat map analysis, to identify one or more documents containing the subject's responses for which further examination is recommended. Words in the documents are characterized in terms of dimensions representing different classes of emotions and states of mind, in which the subject's responses that manifest high stress, emotional volatility and/or internal conflict are identified. A heat map visually displays the dimensions manifested by the subject's responses in different colors, textures, geometric shapes or other visually distinguishable indicia.

Thirumalainambi, Rajkumar↗

Watching Without Seeing a Tool to Surveil Astronaut Health Outcomes While Maintaining Astronaut Medical Privacy

BACKGROUND The Privacy Act of 1974 regulates the use a nd disclosure of personally identifiable information by US Federal agencies. The Act applies to biographical, financial, a nd other identity-linked information, a s well a s personal health information (PHI). As such, the use of astronaut PHI is limited to authorized personnel for preapproved uses, with data reporting often limited to aggregated information about groups. These limitations on the use a nd reporting of astronaut PHI complicates surveillance efforts, wherein epidemiologists a t the National Aeronautics and Space Administration (NASA)monitor the incidence of targeted health conditions in the astronaut population, or to discover emerging trends of aging and disease. Stratification on one or more covariates –particularly time-period, sex, a nd mission participation –can lead to extremely small datasets such that the reporting of results is potentially attributable to individuals. An additional challenge is the small size of the astronaut population, both in terms of numbers of individuals a s well a s in terms of density of exposure time. Such small datasets yield volatile rate estimates that are difficult to interpret. To a id the epidemiological surveillance efforts, a surveillance tool is required that can (a) satisfy the need for rapid computation of condition-specific incidence and mortality rates; (b) improve the statistical estimates of these estimated rates; and (c) maintain astronaut privacy. Here we describe a nd demonstrate such a tool. METHODS We devised a system that models incidence a nd mortality rates rather than calculating them directly. This ha s the advantage of using all the available data to derive the estimates, lea ding to rates that a re not attributable to any one individual, a nd a re a s numerically stable a s they can be given the extremely limited data. The system models disease endpoints using a Poisson regression model with exposure density (measured in person-years) a s a n offset term. By doing so the model is estimating event counts per person-year, equivalent to modeling the rates directly. It uses a standard (pre-specified)set of covariates; the system does not engage in “model-building” as model parsimony is not the goa l. Instead, it is explicitly recognized that if a covariate is not statistically significant a nd not a confounder then it will likely have very little effect on the estimate of the incidence a nd mortality rates. Users are able to specify the disease endpoint of interest and the covariates over which they would like to stratify. The system then uses the resulting model to compute the estimated rates for the user-chosen configuration of variables as visualizes those either over an age range within a specified time-period, or over time for astronauts with a specified age range. RESULTS The first iteration of the tool computes incidence a nd mortality rates for cardiovascular conditions and cancers. Code ha s been developed to retrieve the appropriate data from the IMPALA analysis platform, compute the models for incidence a nd mortality, a nd then use those models to generate the corresponding rate curves. A companion graphical user interface allows the user to specify the curves and visualize the results. CONCLUSIONS It is important to note that the rapid surveillance tool described here is neither meant to be a definitive assessment of the incidence or mortality of any particular disease or condition in the astronaut population, nor is it meant to be used for research purposes. Rather, it is meant as an early indicator that in-depth investigation may be warranted. By automating a repetitive process and leveraging carefully curated astronaut health outcomes, the tool makes possible a rapid “first look” into known areas of concern, and, if used judiciously, may surface new areas of concern for long-term astronaut health. This work is supported in part by the Translational Research Institute for Space Health (TRISH) through NASA Cooperative Agreement NNX16AO69A.

R J Reynolds↗

SU-C-18C-03: Dual-Energy X-Ray Fluoroscopy Imaging System

Purpose: This work studies the clinical utility of dual energy (DE) subtraction fluoroscopy for fiducial-free tumor tracking in lung radiation therapy (RT). Improvement in tumor visualization and quantification of tumor shift within a breathing cycle were analyzed. Methods: Twenty subjects who were undergoing RT for lung cancer were recruited following institutional review board approval. The subjects had a range of tumor sizes, locations in the lungs, and body sizes. An x-ray imaging system was setup with the following components: (a) x-ray tube (Varian G-242, Varian Medical Systems (VMS), CA) (b) flat panel detector (4030CB, VMS, CA) and (c) x-ray generator (EPS 50RF, EMD, Canada). Firmware and software modifications were made to the generator to allow 10 x-ray pulse pairs with alternating low/high kV, 100 ms apart for ∼4s (one breathing cycle). Images were obtained at 4 angles: 0°, 45°, 90° and 135°. Weighted subtraction of a kV-pair image set was used to create a “bone-free” image of the lungs. The 2D tumor-shift in each subtracted image and the 3D shift during a breathing cycle was calculated using all views. Results: The subjects enrolled had the following statistics: average age 62.3±7.1 years, 5 female/15 male, 11 had tumors on the right and 9 on the left and the average tumor size was ∼31.4±10.8 mm. X-ray imaging conditions for the pulse pairs were: 70/120 kVp, 280/221 mA and 65/8 ms. For views where these parameters were insufficient 80/130 kVp, 280/221 mA and 60/12 ms was used. Tumor visibility improved for 0°, 45°, 90° and 135° in 100%, 55%, 75% and 80% of the cases respectively. Tumor shift during a breathing cycle was: 2.4±1.0 mm AP, 2.7±1.4 mm LR and 7.6±4.8 mm IS. Conclusion: DE subtraction fluoroscopy allowed improved visualization and quantification of movement of tumors in the lungs during a breathing cycle. This study was entirely funded by Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗

Realistic simulation of reduced-dose CT with noise modeling and sinogram synthesis using DICOM CT images

Purpose: Reducing the patient dose while maintaining the diagnostic image quality during CT exams is the subject of a growing number of studies, in which simulations of reduced-dose CT with patient data have been used as an effective technique when exploring the potential of various dose reduction techniques. Difficulties in accessing raw sinogram data, however, have restricted the use of this technique to a limited number of institutions. Here, we present a novel reduced-dose CT simulation technique which provides realistic low-dose images without the requirement of raw sinogram data. Methods: Two key characteristics of CT systems, the noise equivalent quanta (NEQ) and the algorithmic modulation transfer function (MTF), were measured for various combinations of object attenuation and tube currents by analyzing the noise power spectrum (NPS) of CT images obtained with a set of phantoms. Those measurements were used to develop a comprehensive CT noise model covering the reduced x-ray photon flux, object attenuation, system noise, and bow-tie filter, which was then employed to generate a simulated noise sinogram for the reduced-dose condition with the use of a synthetic sinogram generated from a reference CT image. The simulated noise sinogram was filtered with the algorithmic MTF and back-projected to create a noise CT image, which was then added to the reference CT image, finally providing a simulated reduced-dose CT image. The simulation performance was evaluated in terms of the degree of NPS similarity, the noise magnitude, the bow-tie filter effect, and the streak noise pattern at photon starvation sites with the set of phantom images. Results: The simulation results showed good agreement with actual low-dose CT images in terms of their visual appearance and in a quantitative evaluation test. The magnitude and shape of the NPS curves of the simulated low-dose images agreed well with those of real low-dose images, showing discrepancies of less than +/−3.2% in terms of the noise power at the peak height and +/−1.2% in terms of the spatial frequency at the peak height. The magnitudes of the noise measured for 12 different combinations the phantom size, tube current, and reconstruction kernel for the simulated and real low-dose images were very similar, with differences of 0.1 to 4.7%. Thep value for a statistical testing of the difference in the noise magnitude ranged from 0.99 to 0.11, showing that there was no difference statistically between the noise magnitudes of the real and simulated low-dose images using our method. The strength and pattern of the streak noise in an anthropomorphic phantom was also consistent with expectations. Conclusions: A novel reduced-dose CT simulation technique was developed which uses only CT images while not requiring raw sinogram data. Our method can provide realistic simulation results under reduced-dose conditions both in terms of the noise magnitude and the textual appearance. This technique has the potential to promote clinical research for patient dose reductions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Towards the distributed burning regime in turbulent premixed flames

Three-dimensional numerical simulations of canonical statistically steady, statistically planar turbulent flames have been used in an attempt to produce distributed burning in lean methane and hydrogen flames. Dilatation across the flame means that extremely large Karlovitz numbers are required; even at the extreme levels of turbulence studied (up to a Karlovitz number of 8767) distributed burning was only achieved in the hydrogen case. In this case, turbulence was found to broaden the reaction zone visually by around an order of magnitude, and thermodiffusive effects (typically present for lean hydrogen flames) were not observed. In the preheat zone, the species compositions differ considerably from those of one-dimensional flames based a number of different transport models (mixture averaged, unity Lewis number and a turbulent eddy viscosity model). The behaviour is a characteristic of turbulence dominating non-unity Lewis number species transport, and the distinct limit is again attributed to dilatation and its effect on the turbulence. Peak local reaction rates are found to be lower in the distributed case than in the lower Karlovitz cases but higher than in the laminar flame, which is attributed to effects that arise from the modified fuel-temperature distribution that results from turbulent mixing dominating low Lewis number thermodiffusive effects. Finally, approaches to achieve distributed burning at realisable conditions are discussed; factors that increase the likelihood of realising distributed burning are higher pressure, lower equivalence ratio, higher Lewis number and lower reactant temperature.

42 ENGINEERING↗