Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Histogram”

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 91 records · Page 5

A new method for diagnosing effective radiative forcing from aerosol–cloud interactions in climate models

Aerosol–cloud interactions (ACIs) are a leading source of uncertainty in estimates of the historical effective radiative forcing (ERF). One reason for this uncertainty is the difficulty in estimating the ERF from aerosol–cloud interactions (ERFaci) in climate models, which typically requires multiple calls to the radiation code. Most commonly used methods also cannot disentangle the contributions from different processes to ERFaci. Here, we develop a new, computationally efficient method for estimating the shortwave (SW) ERFaci from liquid clouds using histograms of monthly averaged cloud fraction partitioned by cloud droplet effective radius (r e ) and liquid water path (LWP). Multiplying the histograms with SW cloud radiative kernels gives the total SW ERFaci from liquid clouds, which can be decomposed into contributions from the Twomey effect, LWP adjustments, and cloud fraction (CF) adjustments. We test the method with data from five CMIP6-era models, using the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite instrument simulator to generate the histograms. Our method gives similar total SW ERFaci estimates to other established methods in regions of prevalent liquid cloud and indicates that the Twomey effect, LWP adjustments, and CF adjustments have contributed -0.34 ± 0.23, -0.22 ± 0.13, and -0.09 ± 0.11 W m -2 , respectively, to the effective radiative forcing of the climate since 1850 in the ensemble mean (95 % confidence). These results demonstrate that widespread adoption of a MODIS r e –LWP joint histogram diagnostic would allow the SW ERFaci and its components to be quickly and accurately diagnosed from climate model outputs, a crucial step for reducing uncertainty in the historical ERF.

54 ENVIRONMENTAL SCIENCES↗

Statistical circuit design for yield improvement in CMOS circuits

This paper addresses the statistical design of CMOS integrated circuits for improved parametric yield. The work uses the Monte Carlo technique of circuit simulation to obtain an unbiased estimation of the yield. A simple graphical analysis tool, the yield factor histogram, is presented. The yield factor histograms are generated by a new computer program called SPICENTER. Using the yield factor histograms, the most sensitive circuit parameters are noted, and their nominal values are changed to improve the yield. Two basic CMOS example circuits, one analog and one digital, are chosen and their designs are 'centered' to illustrate the use of the yield factor histograms for statistical circuit design.

Kamath, H. J.↗

Statistical Properties of Line Centroid Velocity Increments in the rho Ophiuchi Cloud

We present a comparison of histograms of CO (2-1) line centroid velocity increments in the rho Ophiuchi molecular cloud with those computed for spectra synthesized from a three-dimensional, compressible, but non-starforming and non-gravitating hydrodynamic simulation. Histograms of centroid velocity increments in the rho Ophiuchi cloud show clearly non-Gaussian wings, similar to those found in histograms of velocity increments and derivatives in experimental studies of laboratory and atmospheric flows, as well as numerical simulations of turbulence. The magnitude of these wings increases monotonically with decreasing separation, down to the angular resolution of the data. This behavior is consistent with that found in the phase of the simulation which has most of the properties of incompressible turbulence. The time evolution of the magnitude of the non-Gaussian wings in the histograms of centroid velocity increments in the simulation is consistent with the evolution of the vorticity in the flow. However, we cannot exclude the possibility that the wings are associated with the shock interaction regions. Moreover, in an active starforming region like the rho Ophiuchi cloud, the effects of shocks may be more important than in the simulation. However, being able to identify shock interaction regions in the interstellar medium is also important, since numerical simulations show that vorticity is generated in shock interactions.

Lis, D. C.↗

Neutron camera employing row and column summations

For each photomultiplier tube in an Anger camera, an R.times.S array of preamplifiers is provided to detect electrons generated within the photomultiplier tube. The outputs of the preamplifiers are digitized to measure the magnitude of the signals from each preamplifier. For each photomultiplier tube, a corresponding summation circuitry including R row summation circuits and S column summation circuits numerically add the magnitudes of the signals from preamplifiers for each row and for each column to generate histograms. For a P.times.Q array of photomultiplier tubes, P.times.Q summation circuitries generate P.times.Q row histograms including R entries and P.times.Q column histograms including S entries. The total set of histograms include P.times.Q.times.(R+S) entries, which can be analyzed by a position calculation circuit to determine the locations of events (detection of a neutron).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Cylindrical Kernel Density Estimators For Monte Carlo Neutron Transport Reactor Physics Problems

Kernel Density Estimators (KDEs) have been investigated as an alternative to histogram tallies in Monte Carlo radiation transport simulations for capturing spatially-resolved quantities such as scalar flux and reaction rate densities. KDEs are non-parametric estimators that show potential for producing smooth estimates of quantities of interest with reduced variance compared to traditional histogram tallies. Furthermore, KDEs obtain estimates of the underlying distribution at user-defined points without requiring an underlying mesh, thus making KDEs an attractive potential alternative to histogram tallies. Recent work has introduced the Mean Free Path based KDE (MFP KDE) in order to capture distributions with large curvature at material interfaces, e.g. absorption reaction rates in UO{sub 2} fuel in light water reactors. However, the multivariate MFP KDE produces inaccuracies in the form of localized spikes in problems with material cross sections that contain large resonances. Since the MFP KDE works well in one dimension, it is hypothesized that a different MFP KDE formulation where the normalization coefficient only contains one factor of the cross section would be able to capture these distributions with large curvature without producing spikes in the estimate. This can be accomplished by using a cylindrical KDE whose 2-D support area is wedge-shaped and applying the MFP KDE to the radial component of the cylindrical KDE. These cylindrical KDEs use 2-D (r,θ) cylindrical coordinates with no dependence in axial dimension z. This paper derives the cylindrical KDE and the cylindrical MFP KDE and tests them on a reactor physics pin-cell benchmark problem, producing accurate estimates of flux and reaction rate densities with differences less than 0.5% when compared to a reference histogram. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An Efficient Pipeline for Abdomen Segmentation in CT Images

Computed tomography (CT) scans usually include some disadvantages due to the nature of the imaging procedure, and these handicaps prevent accurate abdomen segmentation. Discontinuous abdomen edges, bed section of CT, patient information, closeness between the edges of the abdomen and CT, poor contrast, and a narrow histogram can be regarded as the most important handicaps that occur in abdominal CT scans. Currently, one or more handicaps can arise and prevent technicians obtaining abdomen images through simple segmentation techniques. In other words, CT scans can include the bed section of CT, a patient’s diagnostic information, low-quality abdomen edges, low-level contrast, and narrow histogram, all in one scan. These phenomena constitute a challenge, and an efficient pipeline that is unaffected by handicaps is required. In addition, analysis such as segmentation, feature selection, and classification has meaning for a real-time diagnosis system in cases where the abdomen section is directly used with a specific size. A statistical pipeline is designed in this study that is unaffected by the handicaps mentioned above. Intensity-based approaches, morphological processes, and histogram-based procedures are utilized to design an efficient structure. Performance evaluation is realized in experiments on 58 CT images (16 training, 16 test, and 26 validation) that include the abdomen and one or more disadvantage(s). The first part of the data (16 training images) is used to detect the pipeline’s optimum parameters, while the second and third parts are utilized to evaluate and to confirm the segmentation performance. The segmentation results are presented as the means of six performance metrics. Thus, the proposed method achieves remarkable average rates for training/test/validation of 98.95/99.36/99.57% (jaccard), 99.47/99.67/99.79% (dice), 100/99.91/99.91% (sensitivity), 98.47/99.23/99.85% (specificity), 99.38/99.63/99.87% (classification accuracy), and 98.98/99.45/99.66% (precision). In summary, a statistical pipeline performing the task of abdomen segmentation is achieved that is not affected by the disadvantages, and the most detailed abdomen segmentation study is performed for the use before organ and tumor segmentation, feature extraction, and classification.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

The role of nuclear medicine in assessments of cardiac dyssynchrony

Radionuclide imaging has an advantage for quantitative analyses of the tracer concentration and its temporal changes. Myocardial perfusion and function have been adapted for synchrony analyses. Extracted parameters have been demonstrated to measure ventricular synchrony and even to predict CRT outcomes. ERNA has the advantages of higher temporal resolution, greater reproducibility, and the volumetric analysis of both ventricles that can be applied for analyses of intraventricular synchrony and interventricular synchrony. Several software packages such as Quantitative Gated SPECT, the Emory Cardiac Toolbox, cardioREPO, and Heart Function View are available to assess the LV dyssynchrony parameters from GSPECT. A count-based method is applied to extract the amplitude and phase from each of the reconstructed GSPECT short-axis datasets throughout the cardiac cycle and then subjected to a Fourier analysis, the results of which are displayed on a polar map and histogram. Some of the parameters such as the bandwidth (expressed as the 95% width of the phase histogram) and the standard deviation of the phase are obtained by the phase histogram to assess the intraventricular synchrony. This review paper focuses on the application of the LV dyssynchrony parameters estimated by cardiac SPECT in patients with a heart disease.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SAMBA: Sparse Approximation of Moment-Based Arbitrary Polynomial Chaos

A new arbitrary Polynomial Chaos (aPC) method is presented for moderately high-dimensional problems characterised by limited input data availability. The proposed methodology improves the algorithm of aPC and extends the method, that was previously only introduced as tensor product expansion, to moderately high-dimensional stochastic problems. The fundamental idea of aPC is to use the statistical moments of the input random variables to develop the polynomial chaos expansion. This approach provides the possibility to propagate continuous or discrete probability density functions and also histograms (data sets) as long as their moments exist, are finite and the determinant of the moment matrix is strictly positive. For cases with limited data availability, this approach avoids bias and fitting errors caused by wrong assumptions. In this work, an alternative way to calculate the aPC is suggested, which provides the optimal polynomials, Gaussian quadrature collocation points and weights from the moments using only a handful of matrix operations on the Hankel matrix of moments. It can therefore be implemented without requiring prior knowledge about statistical data analysis or a detailed understanding of the mathematics of polynomial chaos expansions. The extension to more input variables suggested in this work, is an anisotropic and adaptive version of Smolyak's algorithm that is solely based on the moments of the input probability distributions. It is referred to as SAMBA (PC), which is short for Sparse Approximation of Moment-Based Arbitrary Polynomial Chaos. It is illustrated that for moderately high-dimensional problems (up to 20 different input variables or histograms) SAMBA can significantly simplify the calculation of sparse Gaussian quadrature rules. SAMBA's efficiency for multivariate functions with regard to data availability is further demonstrated by analysing higher order convergence and accuracy for a set of nonlinear test functions with 2, 5 and 10 different input distributions or histograms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Characterizing the Stability of NISQ Devices

In this study, we focus on the question of stability of NISQ devices. The parameters that define the device stability profile are motivated by the work of DiVincenzo in [9] where the requirements for physical implementation of quantum computing are discussed. We develop the metrics and theoretical framework to quantify the DiVincenzo requirements and study the stability of those key metrics. The basis of our assessment is histogram similarity (in time and space). For identical experiments, devices which produce reproducible histograms in time, and similar histograms in space, are considered more reliable. To investigate such reliability concerns robustly, we propose a moment-based distance (MBD) metric. We illustrate our methodology using data collected from IBM's Yorktown device. Two types of assessments are discussed: spatial stability and temporal stability.

Dasgupta, Samudra↗

A method for normalizing pathology images to improve feature extraction for quantitative pathology

Purpose: With the advent of digital slide scanning technologies and the potential proliferation of large repositories of digital pathology images, many research studies can leverage these data for biomedical discovery and to develop clinical applications. However, quantitative analysis of digital pathology images is impeded by batch effects generated by varied staining protocols and staining conditions of pathological slides. Methods: To overcome this problem, this paper proposes a novel, fully automated stain normalization method to reduce batch effects and thus aid research in digital pathology applications. Their method, intensity centering and histogram equalization (ICHE), normalizes a diverse set of pathology images by first scaling the centroids of the intensity histograms to a common point and then applying a modified version of contrast-limited adaptive histogram equalization. Normalization was performed on two datasets of digitized hematoxylin and eosin (H&E) slides of different tissue slices from the same lung tumor, and one immunohistochemistry dataset of digitized slides created by restaining one of the H&E datasets. Results: The ICHE method was evaluated based on image intensity values, quantitative features, and the effect on downstream applications, such as a computer aided diagnosis. For comparison, three methods from the literature were reimplemented and evaluated using the same criteria. The authors found that ICHE not only improved performance compared with un-normalized images, but in most cases showed improvement compared with previous methods for correcting batch effects in the literature. Conclusions: ICHE may be a useful preprocessing step a digital pathology image processing pipeline.

61 RADIATION PROTECTION AND DOSIMETRY↗

SU-F-R-20: Image Texture Features Correlate with Time to Local Failure in Lung SBRT Patients

Purpose: To explore possible correlation between CT image-based texture and histogram features and time-to-local-failure in early stage non-small cell lung cancer (NSCLC) patients treated with stereotactic body radiotherapy (SBRT).Methods and Materials: From an IRB-approved lung SBRT registry for patients treated between 2009–2013 we selected 48 (20 male, 28 female) patients with local failure. Median patient age was 72.3±10.3 years. Mean time to local failure was 15 ± 7.1 months. Physician-contoured gross tumor volumes (GTV) on the planning CT images were processed and 3D gray-level co-occurrence matrix (GLCM) based texture and histogram features were calculated in Matlab. Data were exported to R and a multiple linear regression model was used to examine the relationship between texture features and time-to-local-failure. Results: Multiple linear regression revealed that entropy (p=0.0233, multiple R2=0.60) from GLCM-based texture analysis and the standard deviation (p=0.0194, multiple R2=0.60) from the histogram-based features were statistically significantly correlated with the time-to-local-failure. Conclusion: Image-based texture analysis can be used to predict certain aspects of treatment outcomes of NSCLC patients treated with SBRT. We found entropy and standard deviation calculated for the GTV on the CT images displayed a statistically significant correlation with and time-to-local-failure in lung SBRT patients.

60 APPLIED LIFE SCIENCES↗

SU-F-J-94: Development of a Plug-in Based Image Analysis Tool for Integration Into Treatment Planning

Purpose: To extend the functionality of a commercial treatment planning system (TPS) to support (i) direct use of quantitative image-based metrics within treatment plan optimization and (ii) evaluation of dose-functional volume relationships to assist in functional image adaptive radiotherapy. Methods: A script was written that interfaces with a commercial TPS via an Application Programming Interface (API). The script executes a program that performs dose-functional volume analyses. Written in C#, the script reads the dose grid and correlates it with image data on a voxel-by-voxel basis through API extensions that can access registration transforms. A user interface was designed through WinForms to input parameters and display results. To test the performance of this program, image- and dose-based metrics computed from perfusion SPECT images aligned to the treatment planning CT were generated, validated, and compared. Results: The integration of image analysis information was successfully implemented as a plug-in to a commercial TPS. Perfusion SPECT images were used to validate the calculation and display of image-based metrics as well as dose-intensity metrics and histograms for defined structures on the treatment planning CT. Various biological dose correction models, custom image-based metrics, dose-intensity computations, and dose-intensity histograms were applied to analyze the image-dose profile. Conclusion: It is possible to add image analysis features to commercial TPSs through custom scripting applications. A tool was developed to enable the evaluation of image-intensity-based metrics in the context of functional targeting and avoidance. In addition to providing dose-intensity metrics and histograms that can be easily extracted from a plan database and correlated with outcomes, the system can also be extended to a plug-in optimization system, which can directly use the computed metrics for optimization of post-treatment tumor or normal tissue response models. Supported by NIH - P01 - CA059827.

60 APPLIED LIFE SCIENCES↗

TU-H-CAMPUS-JeP3-02: Automated Dose Accumulation and Dose Accuracy Assessment for Online Or Offline Adaptive Replanning

Purpose: With introduction of high-quality treatment imaging during radiation therapy (RT) delivery, e.g., MR-Linac, adaptive replanning of either online or offline becomes appealing. Dose accumulation of delivered fractions, a prerequisite for the adaptive replanning, can be cumbersome and inaccurate. The purpose of this work is to develop an automated process to accumulate daily doses and to assess the dose accumulation accuracy voxel-by-voxel for adaptive replanning. Methods: The process includes the following main steps: 1) reconstructing daily dose for each delivered fraction with a treatment planning system (Monaco, Elekta) based on the daily images using machine delivery log file and considering patient repositioning if applicable, 2) overlaying the daily dose to the planning image based on deformable image registering (DIR) (ADMIRE, Elekta), 3) assessing voxel dose deformation accuracy based on deformation field using predetermined criteria, and 4) outputting accumulated dose and dose-accuracy volume histograms and parameters. Daily CTs acquired using a CT-on-rails during routine CT-guided RT for sample patients with head and neck and prostate cancers were used to test the process. Results: Daily and accumulated doses (dose-volume histograms, etc) along with their accuracies (dose-accuracy volume histogram) can be robustly generated using the proposed process. The test data for a head and neck cancer case shows that the gross tumor volume decreased by 20% towards the end of treatment course, and the parotid gland mean dose increased by 10%. Such information would trigger adaptive replanning for the subsequent fractions. The voxel-based accuracy in the accumulated dose showed that errors in accumulated dose near rigid structures were small. Conclusion: A procedure as well as necessary tools to automatically accumulate daily dose and assess dose accumulation accuracy is developed and is useful for adaptive replanning. Partially supported by Elekta, Inc.

60 APPLIED LIFE SCIENCES↗

The Poisson tensor completion parametric estimator

We introduce the Poisson tensor completion (PTC) estimator that exploits inter-sample relationships to compute a low-rank Poisson tensor decomposition of the frequency histogram for samples of a multivariate distribution. Our crucial observation is that the histogram bins are an instance of a space partitioning of counts and thus can be identified with a spatial non-homogeneous Poisson process. The Poisson tensor decomposition leads to a completion of the mean measure over all bins—including those containing few to no samples—and leads to our proposed estimator. A Poisson tensor decomposition models the underlying distribution of the count data and guarantees non-negative estimated values obviating the need for additional constraints to ensure non-negativity. Furthermore, we demonstrate that our PTC estimator is a substantial improvement over standard histogram-based estimators for sub-Gaussian probability distributions because of the concentration of norm phenomenon.

97 MATHEMATICS AND COMPUTING↗

Laser-Induced Fluorescence for Monitoring Environmental Contamination and Stress in the Moss Thuidium plicatile

The ability to detect, measure, and locate the source of contaminants, especially heavy metals and radionuclides, is of ongoing interest. A common tool for contaminant identification and bioremediation is vegetation that can accumulate and indicate recent and historic pollution. However, large-scale sampling can be costly and labor-intensive. Hence, non-invasive in-situ techniques such as laser-induced fluorescence (LIF) are becoming useful and effective ways to observe the health of plants through the excitation of organic molecules, e.g., chlorophyll. The technique presented utilizes images collected of LIF in moss to identify different metals and environmental stressors. Analysis through image processing of LIF response was key to identifying Cu, Zn, Pb, and a mixture of the metals at nmol/cm 2 levels. Specifically, the RGB values from each image were used to create density histograms of each color channel’s relative pixel abundance at each decimal code value. These histograms were then used to compare color shifts linked to the successful identification of contaminated moss samples. Photoperiod and extraneous environmental stressors had minimal impact on the histogram color shift compared to metals and presented with a response that differentiated them from metal contamination.

59 BASIC BIOLOGICAL SCIENCES↗

The Herschel Orion Protostar Survey: Luminosity and Envelope Evolution

The Herschel Orion Protostar Survey obtained well-sampled 1.2–870 μm spectral energy distributions (SEDs) of over 300 protostars in the Orion molecular clouds, home to most of the young stellar objects (YSOs) in the nearest 500 pc. We plot the bolometric luminosities and temperatures for 330 Orion YSOs, 315 of which have bolometric temperatures characteristic of protostars. The histogram of the bolometric temperature is roughly flat; 29% of the protostars are in Class 0. The median luminosity decreases by a factor of four with increasing bolometric temperature; consequently, the Class 0 protostars are systematically brighter than the Class I protostars, with a median luminosity of 2.3 L{sub ⊙} as opposed to 0.87 L{sub ⊙}. At a given bolometric temperature, the scatter in luminosities is three orders of magnitude. Using fits to the SEDs, we analyze how the luminosities corrected for inclination and foreground reddening relate to the mass in the inner 2500 au of the best-fit model envelopes. The histogram of the envelope mass is roughly flat, while the median-corrected luminosity peaks at 15 L{sub ⊙} for young envelopes and falls to 1.7 L{sub ⊙} for late-stage protostars with remnant envelopes. The spread in luminosity at each envelope mass is three orders of magnitude. Envelope masses that decline exponentially with time explain the flat mass histogram and the decrease in luminosity, while the formation of a range of stellar masses explains the dispersion in luminosity.

79 ASTRONOMY AND ASTROPHYSICS↗

The Gaia-ESO survey: metal-rich bananas in the bulge

We analyze the kinematics of ∼2000 giant stars in the direction of the Galactic bulge, extracted from the Gaia-ESO survey in the region −10{sup ∘}≲ℓ≲10{sup ∘} and −11{sup ∘}≲b≲−3{sup ∘}. We find distinct kinematic trends in the metal-rich ([M/H]>0) and metal-poor ([M/H]<0) stars in the data. The velocity dispersion of the metal-rich stars drops steeply with latitude, compared to a flat profile in the metal-poor stars, as has been seen previously. We argue that the metal-rich stars in this region are mostly on orbits that support the boxy–peanut shape of the bulge, which naturally explains the drop in their velocity dispersion profile with latitude. The metal-rich stars also exhibit peaky features in their line of sight velocity histograms, particularly along the minor axis of the bulge. We propose that these features are due to stars on resonant orbits supporting the boxy–peanut bulge. This conjecture is strengthened through the comparison of the minor axis data with the velocity histograms of resonant orbits generated in simulations of buckled bars. The “banana” or 2:1:2 orbits provide strongly bimodal histograms with narrow velocity peaks that resemble the Gaia-ESO metal-rich data.

79 ASTRONOMY AND ASTROPHYSICS↗

Frequency and time dependence of the Jovian decametric radio emissions - A nineteen-year high-resolution study

A nineteen-year catalog of observations has permitted detailed comparison of Jovian source characteristics at 10 frequencies from 5.6 to 30 MHz. The spectral morphology of Jupiter's decameter-wavelength radio sources was studied by means of high-resolution histograms of occurrence probability for intervals as small as 2 deg in Jovian system III longitude and the departure of Io from superior geocentric conjunction. Longitude histograms clearly exhibit a substantial shift in the source A and source B positions between 15 and 10 MHz. Similar diagrams isolate the dependence of each source on the position of Io, illustrating changes with frequency in the structures and sizes of the main probability maxima and in the positions of surrounding secondary peaks. Individual frequency probability histograms for both longitude and Io phase reveal new phenomena, such as a bridge of radiation through source B longitudes connecting sources A and C at 15 MHz and spectral modification of portions of the Io-related source B peak.

Thieman, J. R.↗