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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2016

MO-FG-202-06: Improving the Performance of Gamma Analysis QA with Radiomics- Based Image Analysis

Purpose: The use of gamma analysis for IMRT quality assurance has well-known limitations. Traditionally, a simple thresholding technique is used to evaluated passing criteria. However, like any image the gamma distribution is rich in information which thresholding mostly discards. We therefore propose a novel method of analyzing gamma images that uses quantitative image features borrowed from radiomics, with the goal of improving error detection. Methods: 368 gamma images were generated from 184 clinical IMRT beams. For each beam the dose to a phantom was measured with EPID dosimetry and compared to the TPS dose calculated with and without normally distributed (2mm sigma) errors in MLC positions. The magnitude of 17 intensity histogram and size-zone radiomic features were derived from each image. The features that differed most significantly between image sets were determined with ROC analysis. A linear machine-learning model was trained on these features to classify images as with or without errors on 180 gamma images.The model was then applied to an independent validation set of 188 additional gamma distributions, half with and half without errors. Results: The most significant features for detecting errors were histogram kurtosis (p=0.007) and three size-zone metrics (p<1e-6 for each). The sizezone metrics detected clusters of high gamma-value pixels under mispositioned MLCs. The model applied to the validation set had an AUC of 0.8, compared to 0.56 for traditional gamma analysis with the decision threshold restricted to 98% or less. Conclusion: A radiomics-based image analysis method was developed that is more effective in detecting error than traditional gamma analysis. Though the pilot study here considers only MLC position errors, radiomics-based methods for other error types are being developed, which may provide better error detection and useful information on the source of detected errors. This work was partially supported by a grant from the Agency for Healthcare Research and Quality, grant number R18 HS022244-01.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

WE-AB-202-06: Correlating Lung CT HU with Transformation-Based and Xe-CT Derived Ventilation

Purpose: Regional lung ventilation is useful to reduce radiation-induced function damage during lung cancer radiation therapy. Recently a new direct HU (Hounsfield unit)-based method was proposed to estimate the ventilation potential without image registration. The purpose of this study is to examine if there is a functional dependence between HU values and transformation-based or Xe-CT derived ventilation. Methods: 4DCT images acquired from 13 patients prior to radiation therapy and 4 mechanically ventilated sheep subjects which also have associated Xe-CT images were used for this analysis. Transformation-based ventilation was computed using Jacobian determinant of the transformation field between peak-exhale and peak-inhale 4DCT images. Both transformation and Xe-CT derived ventilation was computed for each HU bin. Color scatter plot and cumulative histogram were used to compare and validate the direct HU-based method. Results: There was little change of the center and shape of the HU histograms between free breathing CT and 4DCT average, with or without smoothing, and between the repeated 4DCT scans. HU of −750 and −630 were found to have the greatest transformation-based ventilation for human and sheep subjects, respectively. Maximum Xe-CT derived ventilation was found to locate at HU of −600 in sheep subjects. The curve between Xe-CT ventilation and HU was noisy for tissue above HU −400, possibly due to less intensity change of Xe gas during wash-out and wash-in phases. Conclusion: Both transformation-based and Xe-CT ventilation demonstrated that lung tissues with HU values in the range of (-750, −600) HU have the maximum ventilation potential. The correlation between HU and ventilation suggests that HU might be used to help guide the ventilation calculation and make it more robust to noise and image registration errors. Research support from NIH grants CA166703 and CA166119 and a gift from Roger Koch.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

WE-AB-202-12: Voxel-Wise Analysis of Apparent Diffusion Coefficient and Perfusion Maps in Multi-Parametric MRI of Prostate Cancer

Purpose: Multi-parametric MRI (mp-MRI) is being introduced in radiotherapy (RT) of prostate cancer, including for tumour delineation in focal boosting strategies. We recently developed an image-based tumour control probability model, based on cell density distributions derived from apparent diffusion coefficient (ADC) maps. Beyond tumour volume and cell densities, tumour hypoxia is also an important determinant of RT response. Since tissue perfusion from mp-MRI has been related to hypoxia we have explored the patterns of ADC and perfusion maps, and the relations between them, inside and outside prostate index lesions. Methods: ADC and perfusion maps from 20 prostate cancer patients were used, with the prostate and index lesion delineated by a dedicated uro-radiologist. To reduce noise, the maps were averaged over a 3×3×3 voxel cube. Associations between different ADC and perfusion histogram parameters within the prostate, inside and outside the index lesion, were evaluated with the Pearson’s correlation coefficient. In the voxel-wise analysis, scatter plots of ADC vs perfusion were analysed for voxels in the prostate, inside and outside of the index lesion, again with the associations quantified with the Pearson’s correlation coefficient. Results: Overall ADC was lower inside the index lesion than in the normal prostate as opposed to ktrans that was higher inside the index lesion than outside. In the histogram analysis, the minimum ktrans was significantly correlated with the maximum ADC (Pearson=0.47; p=0.03). At the voxel level, 15 of the 20 cases had a statistically significant inverse correlation between ADC and perfusion inside the index lesion; ten of the cases had a Pearson < −0.4. Conclusion: The minimum value of ktrans across the tumour was correlated to the maximum ADC. However, on the voxel level, the ‘local’ ktrans in the index lesion is inversely (i.e. negatively) correlated to the ‘local’ ADC in most patients. Research agreement with Varian Medical Systems, not related to the work presented in this abstract.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

TH-CD-207A-10: Using the Gamma Index to Flag Changes in Anatomy During Radiation Therapy of Head and Neck Cancer

Purpose: This article presents a fast algorithm for comparing 3-D anatomy from Cone-Beam CT (CBCT) imaging using the gamma comparison index and to demonstrate how this can be used to flag patients for possible re-planning of treatment. Methods: CBCT scans acquired on a Varian linear accelerator during treatment were used as input to the gamma comparator using thresholds of 5 mm distance to agreement and 30 Hounsfield Unit CT number difference. The fraction 1 CBCT study was initially used as the reference. Should there be a re-plan during treatment, the reference resets to the CBCT study acquired on the day 1 of the re-plan. Histograms of failing pixels (γ > 1) were generated from each 3-D gamma map. An indicator of anatomy congruence, the match quality parameter (MQP), was derived from failed pixel histograms using the 90th percentile gamma value. The MQP was plotted versus fraction number and related to actual repeat computed tomography (re-CT) order dates as decided by a radiation oncologist. From this, decision criteria were derived for the algorithm to “trigger” re-CT consideration and predictive power was scored using receiver-operator characteristic (ROC) analysis. Results: The MQP plot generally showed that the on-line match from CBCT image guidance deteriorated as the treatment progressed due to weight loss and tumor regression. The optimized MQP criteria for triggering re-CT consideration demonstrated high sensitivity and specificity, consistent with actual re-CT order dates within ± 3 fractions. Out of 20 patients that were actually re-planned, the algorithm failed to trigger a re-CT recommendation only twice and this was caused by CBCT ring artifacts. Conclusion: We have demonstrated that gamma comparisons can be used to evaluate CBCT-acquired anatomy pairs and, from this, an algorithm can be “trained” to flag patients for possible re-planning in a manner consistent with local radiation oncology practice.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

TH-CD-209-05: Impact of Spot Size and Spacing On the Quality of Robustly-Optimized Intensity-Modulated Proton Therapy Plans for Lung Cancer

Purpose: To investigate how spot size and spacing affect plan quality, especially, plan robustness and the impact of interplay effect, of robustly-optimized intensity-modulated proton therapy (IMPT) plans for lung cancer. Methods: Two robustly-optimized IMPT plans were created for 10 lung cancer patients: (1) one for a proton beam with in-air energy dependent large spot size at isocenter (σ: 5–15 mm) and spacing (1.53σ); (2) the other for a proton beam with small spot size (σ: 2–6 mm) and spacing (5 mm). Both plans were generated on the average CTs with internal-gross-tumor-volume density overridden to irradiate internal target volume (ITV). The root-mean-square-dose volume histograms (RVH) measured the sensitivity of the dose to uncertainties, and the areas under RVH curves were used to evaluate plan robustness. Dose evaluation software was developed to model time-dependent spot delivery to incorporate interplay effect with randomized starting phases of each field per fraction. Patient anatomy voxels were mapped from phase to phase via deformable image registration to score doses. Dose-volume-histogram indices including ITV coverage, homogeneity, and organs-at-risk (OAR) sparing were compared using Student-t test. Results: Compared to large spots, small spots resulted in significantly better OAR sparing with comparable ITV coverage and homogeneity in the nominal plan. Plan robustness was comparable for ITV and most OARs. With interplay effect considered, significantly better OAR sparing with comparable ITV coverage and homogeneity is observed using smaller spots. Conclusion: Robust optimization with smaller spots significantly improves OAR sparing with comparable plan robustness and similar impact of interplay effect compare to larger spots. Small spot size requires the use of larger number of spots, which gives optimizer more freedom to render a plan more robust. The ratio between spot size and spacing was found to be more relevant to determine plan robustness and the impact of interplay effect than spot size alone. This research was supported by the National Cancer Institute Career Developmental Award K25CA168984, by the Fraternal Order of Eagles Cancer Research Fund Career Development Award, by The Lawrence W. and Marilyn W. Matteson Fund for Cancer Research, by Mayo Arizona State University Seed Grant, and by The Kemper Marley Foundation.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2017

Multicanonical sampling of the space of states of ℋ(2, n)-vector models

Problems of temperature behavior of specific heat are solved by the entropy simulation method for Ising models on a simple square lattice and a square spin ice (SSI) lattice with nearest neighbor interaction, models of hexagonal lattices with short-range (SR) dipole interaction, as well as with long-range (LR) dipole interaction and free boundary conditions, and models of spin quasilattices with finite interaction radius. It is established that systems of a finite number of Ising spins with LR dipole interaction can have unusual thermodynamic properties characterized by several specific-heat peaks in the absence of an external magnetic field. For a parallel multicanonical sampling method, optimal schemes are found empirically for partitioning the space of states into energy bands for Ising and SSI models, methods of concatenation and renormalization of histograms are discussed, and a flatness criterion of histograms is proposed. It is established that there is no phase transition in a model with nearest neighbor interaction on a hexagonal lattice, while the temperature behavior of specific heat exhibits singularity in the same model, in case of LR interaction. A spin quasilattice is found that exhibits a nonzero value of residual entropy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

Simulation of polarized optical radiation transport in time and frequency representations

Monte Carlo simulation of backward scattering of polarized optical radiation has been performed in the time and frequency representations based on the description of radiation transport in terms of the Bethe–Salpeter equation. The histograms of an ultrashort pulse of linearly or circularly polarized radiation have been calculated in the time representation; it is found that for distances between the source and detector on the order of the transport length, the histograms exhibit a strong dependence on the polarization type; in the case of circular polarization, the direction of the polarization plane rotation changes upon an increase in the scattering anisotropy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

Formation of Radiation Defects by Proton Braking in Lightly Doped n- and p-SiC Layers

Mathematical simulation of the cascade of displacements in SiC is used to consider the specific features of Frenkel-pair generation upon the scattering of 8- and 15-MeV protons. The distribution histograms of energies acquired not only by primary knocked-out atoms, but also by recoil atoms generated in displacement cascades, are calculated. An analysis of the histograms considers two energy ranges. In the first range of “low” energies, the spontaneous recombination of genetically related Frenkel pairs is dominant. Recoil atoms in the second range have a higher energy, which enables these atoms to leave the spontaneousrecombination zone and dissociate into isolated components. The compensation of lightly doped n- and p-4H-SiC samples grown by gas-phase epitaxy is experimentally studied under irradiation with 8- and 15-MeV protons. The carrier removal rates are measured. The calculated and experimental data are compared and estimates are obtained for the size of the spontaneous-recombination zone.

36 MATERIALS SCIENCE↗

DOE OSTI2019

pnnl/rhizo_phosph

Python package to be used when working with the DDAO method of visualizing phosphatase activity in the rhizosphere. Includes a module with functions that generate figures for analysis of phosphatase levels in images of soil and root based on pixel intensity. These functions can be used on images with a single root or multiple roots. Figures include distance maps to show how phosphatase levels change as you move away from the root and along the roots themselves, line plots to show average phosphatase levels moving down the root and the soil surrounding the root, 2D histograms to show how phosphatase levels change moving away from the root, and 1D histograms to show overall differences in phosphatase distribution depending on whether a root was grown next to a phosphorus source or not.

McGrady, Monee↗

DOE OSTI2018

pnnl/pg-histogram

PostgreSQL histograms for ActiveRecord relations. A plugin for the Ruby on Rails (RoR) web application development framework web application development framework that leverages the "Active Relation" relational algebra component to automatically formulate SQL statements for the construction of histograms.

Borkum, Mark↗

DOE OSTI2022

areaDetector/ADTimePix3

ADTimePix3 is an open‑source EPICS areaDetector driver for TimePix3 pixel detectors from Advanced Scientific Instruments (ASI). It provides a complete control and data‑acquisition interface between EPICS and the Serval server, supporting both real detectors (single‑chip and 2×2 quad configurations) and the TimePix3 emulator. The driver configures measurements via an HTTP/JSON REST API and acquires data through high‑rate TCP streams for raw events, images, and histograms. It offers real‑time preview and accumulated images, Time‑of‑Flight histogram processing, pixel masking (including Binary Pixel Configuration files and hot‑pixel tools), extensive health and status monitoring, and integration with standard areaDetector plugins for file saving and downstream analysis. ADTimePix3 is designed for Linux 64‑bit systems (tested on Ubuntu and RHEL) and is distributed under an MIT license with additional institutional disclaimers, enabling reuse and extension in photon science beamlines and other TimePix3‑based experimental facilities.

Gofron, Kaz↗

DOE OSTI2022

Estimating SHmax azimuth with P sources and vertical geophones: Use P-P reflection amplitudes or use SV-P reflection times?

We compared two methods for extracting the azimuth of maximum horizontal stress (SHmax) from 3D land-based seismic data generated by a P source and recorded with vertical geophones. In the first method, we used the direct-SV mode that is produced by all land-based P sources. P sources generate SV illumination that radiates in all azimuth directions from a source station and creates SV-P reflections that are recorded by vertical geophones. Unless stratigraphy has steep dip, SV-P raypaths recorded by vertical geophones are the reverse of P-SV raypaths recorded by horizontal geophones. Thus, SV-P data provide the same S-wave sensitivity to stress fields as popular P-SV data do. In the second method, we retrieved P-P reflections and then performed an amplitude-variation-with-azimuth (AVA) analysis of the amplitude-gradient behavior of P-P reflection wavelets. We did this analysis in narrow azimuth corridors to determine the gradient of reflection-wavelet amplitudes as a function of azimuth. This P-P AVA amplitude-gradient method has been of great interest in the reflection seismology community since it was introduced in the late 1990s. Each of these methods, AVA analysis of the gradient of P-P reflection amplitudes and azimuth-dependent arrival times of SV-P reflections, can be used to determine the azimuth of SHmax stress. We compare the results of the two methods with ground truth measurements of SHmax azimuth at a CO 2 sequestration site in the Michigan Basin. SHmax azimuths were determined from P-P and SV-P data at three major boundaries at depths of approximately 3500 ft (1067 m), 5500 ft (1676 m), and 7500 ft (2286 m). Two estimates of SHmax azimuth (one using SV-P data and one using P-P data) were made at each stacking bin inside a 24 mi 2 (62 km 2 ) image space. The result was approximately 98,000 estimates of SHmax azimuth across each of these three boundaries for each of these two prediction strategies. Histogram displays of PP AVA gradient estimates had peaks at correct azimuths of SHmax at all three depths, but the spread of the distributions widened with depth and split into two peaks at the deepest boundary. In contrast, each histogram of SHmax azimuth predicted by azimuth-dependent SV-P traveltimes had a single, definitive peak that was positioned at the correct SHmax azimuth at all three boundary depths.

Geochemistry & Geophysics↗

DOE OSTI2020

Rossi-alpha Uncertainty Quantification by Analytic, Bootstrap, and Sample Methods to Inform Fitting Best Practices

The prompt neutron period (the negative reciprocal of the prompt neutron decay constant) can be estimated using the Rossi-alpha technique that is predicated on fitting Rossi alpha histograms and of interest in nuclear criticality safety and nonproliferation. The histograms are traditionally fit with a one-exponential model; however, recent work has proposed a two-exponential model to account for reflector induced phenomenon. Until recently, the uncertainty quantification for either model was inadequate (inaccurate and demanded large measurement times). Measurement uncertainty quantification by sample and analytic methods was developed and validated in Ref. The purpose of this transaction is to (i) validate a new bootstrap method by comparing bin-by-bin error bar estimates and (ii) demonstrate how to choose bin widths and reset times to optimize precision and accuracy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2019

Two-Exponential Rossi-Alpha Analysis of Copper- and Polyethylene-Reflected Plutonium

The prompt neutron decay constant $\mathcal{α}$, which is used to infer reactivity, is estimated using the Rossi-alpha technique that is predicated on fitting Rossi-alpha histograms. The histograms are traditionally fit with a one-exponential model; however, recent work has proposed a two-exponential model to account for reflector-induced phenomenon. Recent work also validated the two-exponential model for shields, a subset of reflectors, but not for reflectors in general. Ultimately, it is desirable to validate the two-exponential model for the general shield. This work will investigate if previous 3 He measurements of plutonium reflected by copper and polyethylene, and replacement measurements where the plutonium is replaced by 252 Cf are suitable for validating the new model. Since the 3 He detectors ( 3 He gas proportional counters embedded in a polyethylene matrix) inherently contain moderating/reflector material (to increase detection efficiency), this work also investigates the use of organic scintillation detectors that can directly detect fast, fission neutrons.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2022

Correcting for digital image correlation speckle inversion at high temperature using color cameras

Digital image correlation (DIC) is a popular, noncontacting technique to measure full-field deformation by using cameras to track the motion of an applied surface pattern. Because it is noncontacting, DIC can be performed for extreme temperature applications (e.g., hot-fire rocket testing of carbon composite rocket nozzles) under harsh conditions during which bonded gauges are damaged. Speckle pattern inversion is a phenomenon that sometimes occurs while performing high-temperature DIC. During speckle pattern inversion, portions of the surface pattern that were initially darker at room temperature (e.g., graphite) may emit more light due to blackbody radiation than the portions that were initially paler, thereby producing images in which the pattern appears inverted at high temperature relative to the initial pattern at room temperature. This phenomenon can prevent the correlation algorithm from being able to resolve the displacements between images. This work compares three methods to mitigate speckle pattern inversion: (A) the subtraction method, a recently-published technique in which two high-temperature images are subtracted to remove unwanted light; (B) the filtering method, a popular technique in which optical bandpass filters screen out unwanted light; and (C) the histogram rescaling method, a proposed new method that pairs a color camera with a blue light source and uses information from the green sensor of the camera to correct against inversion in the blue sensor through postprocessing. The histogram rescaling method is shown to successfully eliminate speckle pattern inversion and has the added advantages that it does not require quasi-static loading to be able to compensate for speckle pattern inversion, nor does it impose thick-glass distortions caused by the optical filter.

Rowley, Lindsey J.↗

DOE OSTI2022

NOvA Inclusive Nue CC Cross Section Data Release

NOvA inclusive electron neutrino cross section results presented at Neutrino2020Exposure:8.09E20 protons-on-target, neutrino-enhanced beamData release contains 3 ROOT files containing the double-differential measurement with respect to electron angle (cos theta) and electron energy, a single-differential measurement with respect to squared four-momentum transfer (Q2), and the total cross section as a function of neutrino energy. The 1D cross-section measurements include the same kinematic phase space restrictions present in the double-differential analysis.The files, xsec_{dThdE/Enu/dQ2}.root, contain the extracted cross section and associated covariance matrices for the double-differential cross section with respect to the observed electron kinematics, dThdE, neutrino energy, Enu, and squared four-momentum transfer, Q^2.In each file there will be the following histograms:— TH2D/TH1D xsec: unfolded measured cross section— TH2D xsec_cov: cross section covariance matrix— TH2D stat_cov: statistical covariance matrixIncludes text files, xsec_{dThdE/Enu/dQ2}.txt show the corresponding histograms in plain text format.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2016

Texture and phase analysis of deformed SUS304 by using HIPPO

These slides represent the author's research activity at Los Alamos National Laboratory (LANL), which is about texture and phase analysis of deformed SUS304 by using HIPPO. The following topics are covered: diffraction histogram at each sample position, diffraction histogram (all bank data averaged), possiblity of ε-phase, MAUD analysis with including ε-phase.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Computing a Non-trivial Lower Bound on the Joint Entropy between Two Images

In this report, a non-trivial lower bound on the joint entropy of two non-identical images is developed, which is greater than the individual entropies of the images. The lower bound is the least joint entropy possible among all pairs of images that have the same histograms as those of the given images. New algorithms are presented to compute the joint entropy lower bound with a computation time proportional to S log S where S is the number of histogram bins of the images. This is faster than the traditional methods of computing the exact joint entropy with a computation time that is quadratic in S .

97 MATHEMATICS AND COMPUTING↗