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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Texture analysis of small renal cell carcinomas at MDCT for predicting relevant histologic and protein biomarkers

Purpose: To assess CT texture features of small renal cell carcinomas (≤ 4cm) for association with key pathologic features including protein biomarkers. Methods: Quantitative CT texture analysis (CTTA) of small renal cancers (≤ 4cm) was performed on non-contrast and portal venous phase abdominal MDCT scans with an ROI drawn at the largest cross-sectional diameter of the tumor using commercially available software. Texture parameters including mean pixel attenuation, the standard deviation (SD) of the pixel distribution histogram, entropy, the mean of positive pixels, the skewness (i.e., asymmetry) of the pixel histogram, kurtosis (i.e., peakness) of the pixel histogram, and the percentage of positive pixels were correlated with pathologic data from surgical resection, including histology and nuclear grade, as well as microarray analysis in a subset (n = 40) including Ki67 index, CRP, and neovascularization (CD105/CD31). Results: Portal venous phase images were available in 249 patients (105 women, 144 men; mean age, 56.7 years) with tumors ≤ 4cm (mean, median, range, ± SD; 2.66, 2.60, 0.3–4.0 ± 0.85 cm). CT texture features of standard deviation, mean of the positive pixels, and entropy of the pixel histogram were significantly associated with histologic cell type (clear vs. non-clear; p < 0.001). Entropy and mean of the positive pixels also showed an association with nuclear grade, although not statistically significant. In the microarray analysis subset, kurtosis of the pixel histogram was associated with CD105/CD31 (p = 0.05). SD also showed some association with CD 105 positivity (p = 0.02) and CAIX expression (p = 0.01). Non-contrast CT images were available in 174 patients (72 women, 102 men; mean age, 57.5 years). Although the association with histology was not as strong as on the portal venous phase, in the subset of patients with microarray data, SD was found to correlate with CRP (p = 0.08), kurtosis with CRP (p = 0.004), CD105/CD31 (p = 0.002), and with Ki 67 index (p < 0.001). Conclusion: CT texture features were significantly associated with important histopathologic features in small renal cancers. These non-invasive measures can be performed retrospectively and may provide useful information when determining follow-up and treatment of small renal cancers.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS1992

Cosmic-ray exposure history of ordinary chondrites

The exposure age histograms for H-, L-, and LL-chondrites are discussed. None of these histograms is consistent with a continuous delivery of asteroidal material to the earth, as the observed T(e) histograms clearly disagree with expected exponential distributions for a variety of orbital lifetimes. It is concluded that T(e) histograms are dominated by stochastic events and that the continuous supply of asteroidal material can account only for a minor background of the T(e) histograms. An attempt is made to identify major collisional levels among the major classes of ordinary chondrites in order to estimate the frequency of stochastic events. Orbital maturity in the inner solar system is documented by the p.m./total fall ratio among observed meteorite falls. All chondrite classes exhibit a uniform ratio of 2/3 except type H5 chondrites, which reveal a 0.5 or lower p.m./total fall ratio. This shift in the time of fall statistics suggests a strongly evolved orbit for the H5 parent at the time of collision about 7 Ma ago.

Marti, K.↗
NASA NTRSDate not supplied

Insights into Long-Term Global Precipitation from IMERG and GPCP

Multiple satellite-based datasets provide estimates of the long-term record of global precipitation. Each necessarily includes both the real atmospheric behavior and a collection of artifacts driven by the input data sources and design choices in the retrievals and dataset construction. For this presentation, the Integrated Multi-satellitE Retrievals for GPM (IMERG) products from the Global Precipitation Measurement (GPM) mission’s U.S. Science Team are designed as a modern High-Resolution Precipitation Product (HRPP), while the Global Precipitation Climatology Project’s (GPCP) products are designed as a modern Climate Data Record (CDR). Although they share some common input data sources, the computational process for each is rather different. Comparing the precipitation estimates from each algorithm gives us insights into likely artifacts and the natural variations that are common to both. The goal of an HRPP is to give the best estimate of precipitation at each time step in the dataset, generally meaning that “all possible available” data are used. In contrast, the CDR is intended to provide a precipitation record that has relatively homogeneous statistics, necessary for climate analysis. Each analysis strives for both goals, of course, but this means that the HRPP uses the disparate collection of satellites whose statistics don’t quite match the ideal record, while the CDR is computed from fewer satellites that provide a relatively homogeneous set of precipitation statistics. This presentation employs the current versions of the IMERG and GPCP products, V07B and V3.2, respectively. A variety of analyses have demonstrated that IMERG has more artifacts, as expected, which tend to result from changes in the source of intercalibration for the various input satellites, namely the Combined Radar-Radiometer Algorithm using TMI and PR during the TRMM era and using GMI and DPR-Ku during the GPM era. As well, a shift occurs when the altitudes of the TRMM and GPM satellites change as a result of orbit boosts. For the most part the mean precipitation shows good continuity across these boundaries, but there are noticeable changes in the respective histograms. The latter has taken on more importance in recent years due to increased scrutiny on extremes which implicitly focus on the upper end of the precipitation histograms. The GPCP product is more homogeneous, with only one major change in calibrator, transitioning from the SSMI series of satellites to the SSMIS series of satellites in 2009. The GPCP analysis does a good job of minimizing artifacts in the means at this boundary, while the histograms show a smaller, but still noticeable shift in the histograms, a result that is similar to data boundaries in IMERG. The GPCP Daily, which is a month-by-month rescaling of the daily accumulated IMERG Final product, shows more consistency than the equivalent daily IMERG, but inherits the histogram shifts from IMERG. Given this overlay of artifacts, there is enough consistency between the products to illustrate some important long-term variations, including interannual variations, not all of which are easily attributable to ENSO events, and trends that are large regionally, but comparatively small when averaged across the globe.

precipitation↗
DOE OSTI2021

Validation of the TOPAS Monte Carlo toolkit for HDR brachytherapy simulations

The goal of this work is to validate the user-friendly Geant4-based Monte Carlo toolkit TOol for PArticle Simulation (TOPAS) for brachytherapy applications. Brachytherapy simulations performed with TOPAS were systematically compared with published TG-186 reference data. The photon emission energy spectrum, the air-kerma strength, and the dose-rate constant of the model-based dose calculation algorithm (MBDCA)-WG generic Ir-192 source were extracted. For dose calculations, a track-length estimator was implemented. Additionally, the four Joint AAPM/ESTRO/ABG MBDCA-WG test cases were evaluated through histograms of the local and global dose difference volumes. A prostate, a palliative lung, and a breast case were simulated. For each case, the dose ratio map, the histogram of the global dose difference volume, and cumulative dose-volume histograms were calculated. The air-kerma strength was (9.772 ± 0.001) × 10 -8 U Bq -1 (within 0.3% of the reference value). The dose-rate constant was 1.1107 ± 0.0005 cGy h -1 U -1 (within 0.01% of the reference value). For all cases, at least 96.9% of voxels had a local dose difference within [-1%, 1%] and at least 99.9% of voxels had a global dose difference within [-0.1%, 0.1%]. The implemented track-length estimator scorer was more efficient than the default analog dose scorer by a factor of 237. For all clinical cases, at least 97.5% of voxels had a global dose difference within [-1%, 1%]. Dose-volume histograms were consistent with the reference data. TOPAS was validated for high-dose-rate brachytherapy simulations following the TG-186 recommended approach for MBDCAs. Built on top of Geant4, TOPAS provides broad access to a state-of-the-art Monte Carlo code for brachytherapy simulations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

Modeling the dosimetry of organ-at-risk in head and neck IMRT planning: An intertechnique and interinstitutional study

Purpose: To build a statistical model to quantitatively correlate the anatomic features of structures and the corresponding dose-volume histogram (DVH) of head and neck (HN) Tomotherapy (Tomo) plans. To study if the model built upon one intensity modulated radiation therapy (IMRT) technique (such as conventional Linac) can be used to predict anticipated organs-at-risk (OAR) DVH of patients treated with a different IMRT technique (such as Tomo). To study if the model built upon the clinical experience of one institution can be used to aid IMRT planning for another institution. Methods: Forty-four Tomotherapy intensity modulate radiotherapy plans of HN cases (Tomo-IMRT) from Institution A were included in the study. A different patient group of 53 HN fixed gantry IMRT (FG-IMRT) plans was selected from Institution B. The analyzed OARs included the parotid, larynx, spinal cord, brainstem, and submandibular gland. Two major groups of anatomical features were considered: the volumetric information and the spatial information. The volume information includes the volume of target, OAR, and overlapped volume between target and OAR. The spatial information of OARs relative to PTVs was represented by the distance-to-target histogram (DTH). Important anatomical and dosimetric features were extracted from DTH and DVH by principal component analysis. Two regression models, one for Tomotherapy plan and one for IMRT plan, were built independently. The accuracy of intratreatment-modality model prediction was validated by a leave one out cross-validation method. The intertechnique and interinstitution validations were performed by using the FG-IMRT model to predict the OAR dosimetry of Tomo-IMRT plans. The dosimetry of OARs, under the same and different institutional preferences, was analyzed to examine the correlation between the model prediction and planning protocol. Results: Significant patient anatomical factors contributing to OAR dose sparing in HN Tomotherapy plans have been analyzed and identified. For all the OARs, the discrepancies of dose indices between the model predicted values and the actual plan values were within 2.1%. Similar results were obtained from the modeling of FG-IMRT plans. The parotid gland was spared in a comparable fashion during the treatment planning of two institutions. The model based on FG-IMRT plans was found to predict the median dose of the parotid of Tomotherapy plans quite well, with a mean error of 2.6%. Predictions from the FG-IMRT model suggested the median dose of the larynx, median dose of the brainstem and D2 of the brainstem could be reduced by 10.5%, 12.8%, and 20.4%, respectively, in the Tomo-IMRT plans. This was found to be correlated to the institutional differences in OAR constraint settings. Re-planning of six Tomotherapy patients confirmed the potential of optimization improvement predicted by the FG-IMRT model was correct. Conclusions: The authors established a mathematical model to correlate the anatomical features and dosimetric indexes of OARs of HN patients in Tomotherapy plans. The model can be used for the setup of patient-specific OAR dose sparing goals and quality control of planning results. The institutional clinical experience was incorporated into the model which allows the model from one institution to generate a reference plan for another institution, or another IMRT technique.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

TH-E-BRF-09: Gaussian Mixture Model Analysis of Radiation-Induced Parotid-Gland Injury: An Ultrasound Study of Acute and Late Xerostomia in Head-And-Neck Radiotherapy

Purpose: Xerostomia (dry mouth), secondary to parotid-gland injury, is a distressing side-effect in head-and-neck radiotherapy (RT). This study's purpose is to develop a novel ultrasound technique to quantitatively evaluate post-RT parotid-gland injury. Methods: Recent ultrasound studies have shown that healthy parotid glands exhibit homogeneous echotexture, whereas post-RT parotid glands are often heterogeneous, with multiple hypoechoic (inflammation) or hyperechoic (fibrosis) regions. We propose to use a Gaussian mixture model to analyze the ultrasonic echo-histogram of the parotid glands. An IRB-approved clinical study was conducted: (1) control-group: 13 healthy-volunteers, served as the control; (2) acutetoxicity group − 20 patients (mean age: 62.5 ± 8.9 years, follow-up: 2.0±0.8 months); and (3) late-toxicity group − 18 patients (mean age: 60.7 ± 7.3 years, follow-up: 20.1±10.4 months). All patients experienced RTOG grade 1 or 2 salivary-gland toxicity. Each participant underwent an ultrasound scan (10 MHz) of the bilateral parotid glands. An echo-intensity histogram was derived for each parotid and a Gaussian mixture model was used to fit the histogram using expectation maximization (EM) algorithm. The quality of the fitting was evaluated with the R-squared value. Results: (1) Controlgroup: all parotid glands fitted well with one Gaussian component, with a mean intensity of 79.8±4.9 (R-squared>0.96). (2) Acute-toxicity group: 37 of the 40 post-RT parotid glands fitted well with two Gaussian components, with a mean intensity of 42.9±7.4, 73.3±12.2 (R-squared>0.95). (3) Latetoxicity group: 32 of the 36 post-RT parotid fitted well with 3 Gaussian components, with mean intensities of 49.7±7.6, 77.2±8.7, and 118.6±11.8 (R-squared>0.98). Conclusion: RT-associated parotid-gland injury is common in head-and-neck RT, but challenging to assess. This work has demonstrated that the Gaussian mixture model of the echo-histogram could quantify acute and late toxicity of the parotid glands. This study provides meaningful preliminary data from future observational and interventional clinical research.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Postmortem validation of breast density using dual-energy mammography

Purpose: Mammographic density has been shown to be an indicator of breast cancer risk and also reduces the sensitivity of screening mammography. Currently, there is no accepted standard for measuring breast density. Dual energy mammography has been proposed as a technique for accurate measurement of breast density. The purpose of this study is to validate its accuracy in postmortem breasts and compare it with other existing techniques. Methods: Forty postmortem breasts were imaged using a dual energy mammography system. Glandular and adipose equivalent phantoms of uniform thickness were used to calibrate a dual energy basis decomposition algorithm. Dual energy decomposition was applied after scatter correction to calculate breast density. Breast density was also estimated using radiologist reader assessment, standard histogram thresholding and a fuzzy C-mean algorithm. Chemical analysis was used as the reference standard to assess the accuracy of different techniques to measure breast composition. Results: Breast density measurements using radiologist reader assessment, standard histogram thresholding, fuzzy C-mean algorithm, and dual energy were in good agreement with the measured fibroglandular volume fraction using chemical analysis. The standard error estimates using radiologist reader assessment, standard histogram thresholding, fuzzy C-mean, and dual energy were 9.9%, 8.6%, 7.2%, and 4.7%, respectively. Conclusions: The results indicate that dual energy mammography can be used to accurately measure breast density. The variability in breast density estimation using dual energy mammography was lower than reader assessment rankings, standard histogram thresholding, and fuzzy C-mean algorithm. Improved quantification of breast density is expected to further enhance its utility as a risk factor for breast cancer.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Evaluation of Matrix9 silicon photomultiplier array for small-animal PET

Purpose: The MatrixSL-9-30035-OEM (Matrix9) from SensL is a large-area silicon photomultiplier (SiPM) photodetector module consisting of a 3 × 3 array of 4 × 4 element SiPM arrays (total of 144 SiPM pixels) and incorporates SensL’s front-end electronics board and coincidence board. Each SiPM pixel measures 3.16 × 3.16 mm{sup 2} and the total size of the detector head is 47.8 × 46.3 mm{sup 2}. Using 8 × 8 polished LSO/LYSO arrays (pitch 1.5 mm) the performance of this detector system (SiPM array and readout electronics) was evaluated with a view for its eventual use in small-animal positron emission tomography (PET). Methods: Measurements of noise, signal, signal-to-noise ratio, energy resolution, flood histogram quality, timing resolution, and array trigger error were obtained at different bias voltages (28.0–32.5 V in 0.5 V intervals) and at different temperatures (5 °C–25 °C in 5 °C degree steps) to find the optimal operating conditions. Results: The best measured signal-to-noise ratio and flood histogram quality for 511 keV gamma photons were obtained at a bias voltage of 30.0 V and a temperature of 5 °C. The energy resolution and timing resolution under these conditions were 14.2% ± 0.1% and 4.2 ± 0.1 ns, respectively. The flood histograms show that all the crystals in the 1.5 mm pitch LSO array can be clearly identified and that smaller crystal pitches can also be resolved. Flood histogram quality was also calculated using different center of gravity based positioning algorithms. Improved and more robust results were achieved using the local 9 pixels for positioning along with an energy offset calibration. To evaluate the front-end detector readout, and multiplexing efficiency, an array trigger error metric is introduced and measured at different lower energy thresholds. Using a lower energy threshold greater than 150 keV effectively eliminates any mispositioning between SiPM arrays. Conclusions: In summary, the Matrix9 detector system can resolve high-resolution scintillator arrays common in small-animal PET with adequate energy resolution and timing resolution over a large detector area. The modular design of the Matrix9 detector allows it to be used as a building block for simple, low channel-count, yet high performance, small animal PET or PET/MRI systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

SU-F-I-45: An Automated Technique to Measure Image Contrast in Clinical CT Images

Purpose: To develop and validate an automated technique for measuring image contrast in chest computed tomography (CT) exams. Methods: An automated computer algorithm was developed to measure the distribution of Hounsfield units (HUs) inside four major organs: the lungs, liver, aorta, and bones. These organs were first segmented or identified using computer vision and image processing techniques. Regions of interest (ROIs) were automatically placed inside the lungs, liver, and aorta and histograms of the HUs inside the ROIs were constructed. The mean and standard deviation of each histogram were computed for each CT dataset. Comparison of the mean and standard deviation of the HUs in the different organs provides different contrast values. The ROI for the bones is simply the segmentation mask of the bones. Since the histogram for bones does not follow a Gaussian distribution, the 25th and 75th percentile were computed instead of the mean. The sensitivity and accuracy of the algorithm was investigated by comparing the automated measurements with manual measurements. Fifteen contrast enhanced and fifteen non-contrast enhanced chest CT clinical datasets were examined in the validation procedure. Results: The algorithm successfully measured the histograms of the four organs in both contrast and non-contrast enhanced chest CT exams. The automated measurements were in agreement with manual measurements. The algorithm has sufficient sensitivity as indicated by the near unity slope of the automated versus manual measurement plots. Furthermore, the algorithm has sufficient accuracy as indicated by the high coefficient of determination, R2, values ranging from 0.879 to 0.998. Conclusion: Patient-specific image contrast can be measured from clinical datasets. The algorithm can be run on both contrast enhanced and non-enhanced clinical datasets. The method can be applied to automatically assess the contrast characteristics of clinical chest CT images and quantify dependencies that may not be captured in phantom data.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

The Poisson tensor completion non-parametric differential entropy estimator

We introduce the Poisson tensor completion (PTC) estimator, a non-parametric differential entropy estimator. The PTC estimator leverages inter-sample relationships to compute a low-rank Poisson tensor decomposition of the frequency histogram. 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 Poisson process. The Poisson tensor decomposition leads to a completion of the intensity measure over all bins—including those containing few to no samples—and leads to our proposed PTC differential entropy estimator. A Poisson tensor decomposition models the underlying distribution of the count data and guarantees non-negative estimated values and so can be safely used directly in entropy estimation. Our estimator is the first tensor-based estimator that exploits the underlying spatial Poisson process related to the histogram explicitly when estimating the probability density with low-rank tensor decompositions for the purpose of tensor completion. 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.

42 ENGINEERING↗
DOE OSTI2025

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↗
NASA NTRS1990

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.↗
NASA NTRS1998

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.↗
DOE OSTI2016

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 ↗
DOE OSTI2016

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↗
DOE OSTI2018

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↗
DOE OSTI2018

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↗
DOE OSTI2016

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↗