Search NASA⌕ Search

SEARCH · Search NASA

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 145 records · Page 8

TU-F-CAMPUS-J-01: Dosimetric Effects of HU Changes During the Course of Proton Therapy for Lung Cancer

Purpose: To characterize the changes in Hounsfield unit (HU) in lung radiotherapy with proton beams during the course of treatment and to study the effect on the proton plan dose distribution. Methods: Twenty consecutive patients with non-small cell lung cancer treated with proton radiotherapy who underwent multiple CT scans including the planning CT and weekly verification CTs were studied. HU histograms were computed for irradiated lung volumes in beam paths for all scans using the same treatment plan. Histograms for un-irradiated lung volume were used as control to characterize inter-scan variations. HU statistics were calculated for both irradiated and un-irradiated lung volumes for each patient scan. Further, multiple CT scans based on the same planning CT were generated by replacing the HU of the lung based on the verification CT scans HU values. Using the same beam arrangement, we created plans for each of the altered CT scans to study the dosimetric effect using the dose volume histogram. Results: Lung HU decreased for irradiated lung volume during the course of radiotherapy. The magnitude of this change increased with total irradiation dose. On average, HU changed by −53.8 in the irradiated volume. This change resulted in less than 0.5mm of beam overshoot in tissue for every 1cm beam traversed in the irradiated lung. The dose modification is about +3% for the lung, and less than +1% for the primary tumor. Conclusion: HU of the lung decrease throughout the course of radiation therapy. This change results in a beam overshoot (e.g. 3mm for 6cm of lung traversed) and causes a small dose modification in the overall plan. However, this overshoot does not affect the quality of plans since the margins used in planning, based on proton range uncertainty, are greater. HU needs to change by 150 units before re-planning is warranted.

60 APPLIED LIFE SCIENCES↗

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↗

Langston University - High Energy Physics (LU-HEP)

This final report is presented by Langston University (LU) for the project entitled "Langston University High Energy Physics" (LUHEP) under the direction of principal investigator (PI) and project director Professor Joel Snow. The project encompassed high energy physics research performed at hadron colliders. The PI is a collaborator on the DZero experiment at Fermi National Accelerator Laboratory in Batavia, IL, USA and the ATLAS experiment at CERN in Geneva, Switzerland and was during the entire project period from April 1, 1999 until May 14, 2012. Both experiments seek to understand the fundamental constituents of the physical universe and the forces that govern their interactions. In 1999 as member of the Online Systems group for Run 2 the PI developed a cross-platform Python-based, Graphical User Interface (GUI) application for monitoring and control of EPICS based devices for control room use. This served as a model for other developers to enhance and build on for further monitoring and control tasks written in Python. Subsequently the PI created and developed a cross-platform C++ GUI utilizing a networked client-server paradigm and based on ROOT, the object oriented analysis framework from CERN. The GUI served as a user interface to the Examine tasks running in the D\O\ control room which monitored the status and integrity of data taking for Run 2. The PI developed the histogram server/control interface to the GUI client for the EXAMINE processes. The histogram server was built from the ROOT framework and was integrated into the D\O\ framework used for online monitoring programs and offline analysis. The PI developed the first implementation of displaying histograms dynamically generated by ROOT in a Web Browser. The PI's work resulted in several talks and papers at international conferences and workshops. The PI established computing software infrastructure at LU and U. Oklahoma (OU) to do analysis of DZero production data and produce simulation data for the experiment. Eventually this included the FNAL SAM data grid system, the SAMGrid (SG) infrastructure, and the Open Science Grid software stacks for computing and storage elements. At the end of 2003 Snow took on the role of global Monte Carlo production coordinator for the DØ experiment. A role which continues til this day. In January of 2004 Snow started working with the SAMGrid development team to help debug, deploy, and integrate SAMGrid with DØ Monte Carlo production. Snow installed and configured SG execution and client sites at LUHEP and OUHEP, and a SG scheduler site at LUHEP. The PI developed a python based GUI (DAJ) that acts as a front end for job submission to SAMGrid. The GUI interfaces to the DZero Mone Carlo (MC) request system that uses SAM to manage MC requests by the physics analysis groups. DAJ significantly simplified SG job submission and was deployed in DZero in an effort to increase the user base of SG. The following year was the advent of SAMGrid job submission to the Open Science Grid (OSG) and LHC Computing Grid (LCG) through a forwarding mechanism. The PI oversaw the integration of these grids into the existing production infrastructure. The PI developed an automatic MC (Automc) request processing system capable of operating without user intervention (other than getting grid credentials), and able to submit to any number of sites on various grids. The system manages production at all but 2 sites. The system was deployed at Fermilab and remains operating there today. The PI's work in distributed computing resulted in several talks at international conferences. UTA, OU, and LU were chosen as the collaborating institutions that form the Southwest Tier 2 Center (SWT2) for ATLAS. During the project period the PI contributed to the online and offline software infrastructure through his work with the Run 2 online group, and played a major role in Monte Carlo production for DZero. During the part of the project period in which the PI served as MC production coordinator MC production increased very significantly. In the first year of the PI's tenure as production coordinator production was 159M events and 6.7~TB of data. During the last year of the project period production was 2,342~M events and 262~TB of data. That is a factor of 15 increase in events and 39 in data volume. The increase occurred with improvements in computer hardware and networks, through the use of grid technology on diverse resources, and through increased automation and efficiency of the production process. LU HEP developed and deployed the automatic MC request processing system in use at FNAL. The complementary strategies of automation and grid production served DZero well. Fermilab has recognized LU HEP's contribution to DZero by allowing the PI to devote full time to research activities by appointing him a guest scientist for the last six years of the project period.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

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

Mission load dynamic tests of two undensified Space shuttle thermal protection system tiles

Two tests of undensified Space Shuttle thermal protection tiles under combined static and dynamic loads were conducted. The tiles had a density of approximately 144 Kg/cum (LI900 tiles) and were mounted on a strain isolation pad which was 0.41 cm (.160 inch) thick. A combined static and dynamic mission stress histogram representative of the W-3 area of the wing of the orbiter vehicle was applied. The stress histogram was provided by the space shuttle project. Results presented include: tabulation of measured peak and root-mean-square (RMS) accelerations in both compression and tension; peak SIP stress in compression and tension, peak and RMS amplitude response ratios; lateral to vertical response ratios; response time histories; peak stress distributions (histograms), and SIP extension measured both with and without static tension at various mission times.

Leatherwood, J. D.↗

Analyzing Multidimensional Image Data

Six computer programs perform histogram cluster analysis. Histogram Cluster Analysis Procedure (HICAP) developed to perform unsupervised classification of multidimensional image data. Clustering approach used in HICAP based on algorithm which uses multidimensional histogram to perform unsupervised classification of four-dimensional Landsat multispectral-scanner data. HICAP generalizes this procedure to process up to 32-bit data with arbitrary number of dimensions. Also incorporates efficiency improvements so classification requires less computation than original algorithm. Computational savings afforded by HICAP increase with number of dimensions in data. HICAP programs written in FORTRAN 77 for batch or interactive execution.

Wharton, S. W.↗

Nearest neighbor spacing of fair weather cumulus clouds

Histograms of nearest neighbor spacings of fair weather cumulus at 15 locations over the world's oceans are presented based on the analysis of high resolution Landsat 3 Multispectral Scanner images for amounts of cloud cover ranging from 0.6 to 37.6 percent. These histograms are found to be essentially the same at all locations analyzed, similarly to previous findings on the size distributions and the fractal dimensions of the perimeters for this cloud type. The nearest neighbor spacings are linearly dependent on the effective cloud radii, with a proportionality factor ranging from five to twenty. The histograms peak at about 0.5 km. Nearest-neighbor spacings smaller than about a kilometer, associated with cumulus clouds with an effective radius less than a few hundred meters, have a distribution of cloud centers that is almost indepedent in the horizontal plane and show a tendency for the formation of clumps. Larger spacings of up to thirty kilometers occur and are associated with the larger clouds. These latter spacings are not independent.

Joseph, Joachim H.↗

Estimation of area-averaged rainfall over tropical oceans from microwave radiometry - A single channel approach

This paper presents a new simple retrieval algorithm for estimating area-time averaged rain rates over tropical oceans by using single channel microwave measurements from satellites. The algorithm was tested by using the Nimbus-5 Electrically Scanning Microwave Radiometer and a simple microwave radiative transfer model to retrieve seasonal 5-deg x 5-deg area averaged rainrate over the tropical Atlantic and Pacific from December 1973 to November 1974. The brightness temperatures were collected and analyzed into histograms for each season and in each grid box from December 1973 to November 1974. The histograms suggest a normal distribution of background noise plus a skewed rain distribution at the higher brightness temperatures. By using a statistical estimation procedure based upon normally distributed background noise, the rain distribution was separated from the raw histogram. The radiative transfer model was applied to the rain-only distribution to retrieve area-time averaged rainrates throughout the tropics. Despite limitations of single channel information, the retrieved seasonal rain rates agree well in the open ocean with expectations based upon previous estimates of tropical rainfall over the oceans.

Shin, Kyung-Sup↗

Principal modes of variation of rain-rate probability distributions

Radar or satellite observations of an area generate sequences of rain-rate maps. From a gridded map a histogram of rain rates can be obtained representing the relative areas occupied by rain rates of various strengths. The histograms vary with time as precipitating systems in the area evolve and decay and amounts of convective and stratiform rain in the area change. A method of decomposing the histograms into linear combinations of a few empirical distributions with time-dependent coefficients is developed, using principal component analysis as a starting point. When applied to a tropical Atlantic dataset (GATE), two distributions emerge naturally from the analysis, resembling stratiform and convective rain-rate distributions in that they peak at low and high rain rates, respectively. The two 'modes' have different timescales and only the high-rain-rate mode has a statistically significant diurnal cycle. The ability of just two modes to describe rain variabiltiy over an area can explain why methods of estimating area-averaged rain rate from the area covered by rain rates above a certain threshold are so successful.

Bell, Thomas L.↗

The effects of cloud inhomogeneities upon radiative fluxes, and the supply of a cloud truth validation dataset

The ASTER polar cloud mask algorithm is currently under development. Several classification techniques have been developed and implemented. The merits and accuracy of each are being examined. The classification techniques under investigation include fuzzy logic, hierarchical neural network, and a pairwise histogram comparison scheme based on sample histograms called the Paired Histogram Method. Scene adaptive methods also are being investigated as a means to improve classifier performance. The feature, arctan of Band 4 and Band 5, and the Band 2 vs. Band 4 feature space are key to separating frozen water (e.g., ice/snow, slush/wet ice, etc.) from cloud over frozen water, and land from cloud over land, respectively. A total of 82 Landsat TM circumpolar scenes are being used as a basis for algorithm development and testing. Numerous spectral features are being tested and include the 7 basic Landsat TM bands, in addition to ratios, differences, arctans, and normalized differences of each combination of bands. A technique for deriving cloud base and top height is developed. It uses 2-D cross correlation between a cloud edge and its corresponding shadow to determine the displacement of the cloud from its shadow. The height is then determined from this displacement, the solar zenith angle, and the sensor viewing angle.

Welch, Ronald M.↗

Information-Adaptive Image Encoding and Restoration

The multiscale retinex with color restoration (MSRCR) has shown itself to be a very versatile automatic image enhancement algorithm that simultaneously provides dynamic range compression, color constancy, and color rendition. A number of algorithms exist that provide one or more of these features, but not all. In this paper we compare the performance of the MSRCR with techniques that are widely used for image enhancement. Specifically, we compare the MSRCR with color adjustment methods such as gamma correction and gain/offset application, histogram modification techniques such as histogram equalization and manual histogram adjustment, and other more powerful techniques such as homomorphic filtering and 'burning and dodging'. The comparison is carried out by testing the suite of image enhancement methods on a set of diverse images. We find that though some of these techniques work well for some of these images, only the MSRCR performs universally well oil the test set.

Park, Stephen K.↗