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 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS2003

A New Objective Technique for Verifying Mesoscale Numerical Weather Prediction Models

This report presents a new objective technique to verify predictions of the sea-breeze phenomenon over east-central Florida by the Regional Atmospheric Modeling System (RAMS) mesoscale numerical weather prediction (NWP) model. The Contour Error Map (CEM) technique identifies sea-breeze transition times in objectively-analyzed grids of observed and forecast wind, verifies the forecast sea-breeze transition times against the observed times, and computes the mean post-sea breeze wind direction and speed to compare the observed and forecast winds behind the sea-breeze front. The CEM technique is superior to traditional objective verification techniques and previously-used subjective verification methodologies because: It is automated, requiring little manual intervention, It accounts for both spatial and temporal scales and variations, It accurately identifies and verifies the sea-breeze transition times, and It provides verification contour maps and simple statistical parameters for easy interpretation. The CEM uses a parallel lowpass boxcar filter and a high-order bandpass filter to identify the sea-breeze transition times in the observed and model grid points. Once the transition times are identified, CEM fits a Gaussian histogram function to the actual histogram of transition time differences between the model and observations. The fitted parameters of the Gaussian function subsequently explain the timing bias and variance of the timing differences across the valid comparison domain. Once the transition times are all identified at each grid point, the CEM computes the mean wind direction and speed during the remainder of the day for all times and grid points after the sea-breeze transition time. The CEM technique performed quite well when compared to independent meteorological assessments of the sea-breeze transition times and results from a previously published subjective evaluation. The algorithm correctly identified a forecast or observed sea-breeze occurrence or absence 93% of the time during the two- month evaluation period from July and August 2000. Nearly all failures in CEM were the result of complex precipitation features (observed or forecast) that contaminated the wind field, resulting in a false identification of a sea-breeze transition. A qualitative comparison between the CEM timing errors and the subjectively determined observed and forecast transition times indicate that the algorithm performed very well overall. Most discrepancies between the CEM results and the subjective analysis were again caused by observed or forecast areas of precipitation that led to complex wind patterns. The CEM also failed on a day when the observed sea- breeze transition affected only a very small portion of the verification domain. Based on the results of CEM, the RAMS tended to predict the onset and movement of the sea-breeze transition too early and/or quickly. The domain-wide timing biases provided by CEM indicated an early bias on 30 out of 37 days when both an observed and forecast sea breeze occurred over the same portions of the analysis domain. These results are consistent with previous subjective verifications of the RAMS sea breeze predictions. A comparison of the mean post-sea breeze winds indicate that RAMS has a positive wind-speed bias for .all days, which is also consistent with the early bias in the sea-breeze transition time since the higher wind speeds resulted in a faster inland penetration of the sea breeze compared to reality.

Case, Jonathan L.↗

NASA NTRS2003

Design of a Borescope for Extravehicular Non-Destructive Applications

Anomalies such as corrosion, structural damage, misalignment, cracking, stress fiactures, pitting, or wear can be detected and monitored by the aid of a borescope. A borescope requires a source of light for proper operation. Today s current lighting technology market consists of incandescent lamps, fluorescent lamps and other types of electric arc and electric discharge vapor lamp. Recent advances in LED technology have made LEDs viable for a number of applications, including vehicle stoplights, traffic lights, machine-vision-inspection, illumination, and street signs. LEDs promise significant reduction in power consumption compared to other sources of light. This project focused on comparing images taken by the Olympus IPLEX, using two different light sources. One of the sources is the 50-W internal metal halide lamp and the other is a 1 W LED placed at the tip of the insertion tube. Images acquired using these two light sources were quantitatively compared using their histogram, intensity profile along a line segment, and edge detection. Also, images were qualitatively compared using image registration and transformation [l]. The gray-level histogram, edge detection, image profile and image registration do not offer conclusive results. The LED light source, however, produces good images for visual inspection by an operator. Analysis using pattern recognition using Eigenfaces and Gaussian Pyramid in face recognition may be more useful.

Bachnak, Rafic↗

NASA NTRS1997

Anterior segment sparing to reduce charged particle radiotherapy complications in uveal melanoma

PURPOSE: The purpose of this investigation is to delineate the risk factors in the development of neovascular glaucoma (NVG) after helium-ion irradiation of uveal melanoma patients and to propose treatment technique that may reduce this risk. METHODS AND MATERIALS: 347 uveal melanoma patients were treated with helium-ions using a single-port treatment technique. Using univariate and multivariate statistics, the NVG complication rate was analyzed according to the percent of anterior chamber in the radiation field, tumor size, tumor location, sex, age, dose, and other risk factors. Several University of California San Francisco-Lawrence Berkeley National Laboratory (LBNL) patients in each size category (medium, large, and extralarge) were retrospectively replanned using two ports instead of a single port. By using appropriate polar and azimuthal gaze angles or by treating patients with two ports, the maximum dose to the anterior segment of the eye can often be reduced. Although a larger volume of anterior chamber may receive a lower dose by using two ports than a single port treatment. We hypothesize that this could reduce the level of complications that result from the irradiation of the anterior chamber of the eye. Dose-volume histograms were calculated for the lens, and compared for the single and two-port techniques. RESULTS: NVG developed in 121 (35%) patients. The risk of NVG peaked between 1 and 2.5 years posttreatment. By univariate and multivariate analysis, the percent of lens in the field was strongly correlated with the development of NVG. Other contributing factors were tumor height, history of diabetes, and vitreous hemorrhage. Dose-volume histogram analysis of single-port vs. two-port techniques demonstrate that for some patients in the medium and large category tumor groups, a significant decrease in dose to the structures in the anterior segment of the eye could have been achieved with the use of two ports. CONCLUSION: The development of NVG after helium-ion irradiation is correlated to the amount of lens, anterior chamber in the treatment field, tumor height, proximity to the fovea, history of diabetes, and the development of vitreous hemorrhage. Although the influence of the higher LET deposition of helium-ions is unclear, this study suggests that by reducing the dose to the anterior segment of the eye may reduce the NVG complications. Based on this retrospective analysis of LBNL patients, we have implemented techniques to reduce the amount of the anterior segment receiving a high dose in our new series of patients treated with protons using the cyclotron at the UC Davis Crocker Nuclear Laboratory (CNL).

NASA Discipline Radiation Health↗

NASA NTRS2005

NSTAR Extended Life Test Discharge Chamber Flake Analysis

The Extended Life Test (ELT) of the NASA Solar Electric Propulsion Technology Readiness (NSTAR) ion thruster was concluded after 30,352 hours of operation. The ELT was conducted using the Deep Space 1 (DS1) back-up flight engine, a 30 cm diameter xenon ion thruster. Post-test inspection of the ELT engine revealed numerous contaminant flakes distributed over the bottom of the cylindrical section of the anode within the discharge chamber (DC). Extensive analyses were conducted to determine the source of the particles, which is critical to the understanding of degradation mechanisms of long life ion thruster operation. Analyses included: optical microscopy (OM) and particle length histograms, field emission scanning electron microscopy (FESEM) combined with energy dispersive spectroscopy (EDS), and atomic oxygen plasma exposure tests. Analyses of the particles indicate that the majority of the DC flakes consist of a layered structure, typically with either two or three layers. The flakes comprising two layers were typically found to have a molybdenum-rich (Mo-rich) layer on one side and a carbon-rich (C-rich) layer on the other side. The flakes comprising three layers were found to be sandwich-like structures with Mo-rich exterior layers and a C-rich interior layer. The presence of the C-rich layers indicates that these particles were produced by sputter deposition build-up on a surface external to the discharge chamber from ion sputter erosion of the graphite target in the test chamber. This contaminant layer became thick enough that particles spalled off, and then were electro-statically attracted into the ion thruster interior, where they were coated with Mo from internal sputter erosion of the screen grid and cathode components. Atomic oxygen tests provided evidence that the DC chamber flakes are composed of a significant fraction of carbon. Particle size histograms further indicated that the source of the particles was spalling of carbon flakes from downstream surfaces. Analyses of flakes taken from the downstream surface of the accelerator grid provided additional supportive information. The production of the downstream carbon flakes, and hence the potential problems associated with the flake particles in the ELT ion thruster engine is a facility induced effect and would not occur in the space environment.

deGroh, Kim K.↗

NASA NTRS2004

National Airspace System Delay Estimation Using Weather Weighted Traffic Counts

Assessment of National Airspace System performance, which is usually measured in terms of delays resulting from the application of traffic flow management initiatives in response to weather conditions, volume, equipment outages and runway conditions, is needed both for guiding flow control decisions during the day of operations and for post operations analysis. Comparison of the actual delay, resulting from the traffic flow management initiatives, with the expected delay, based on traffic demand and other conditions, provides the assessment of the National Airspace System performance. This paper provides a method for estimating delay using the expected traffic demand and weather. In order to identify the cause of delays, 517 days of National Airspace System delay data reported by the Federal Aviation Administration s Operations Network were analyzed. This analysis shows that weather is the most important causal factor for delays followed by equipment and runway delays. Guided by these results, the concept of weather weighted traffic counts as a measure of system delay is described. Examples are given to show the variation of these counts as a function of time of the day. The various datasets, consisting of aircraft position data, enroute severe weather data, surface wind speed and visibility data, reported delay data and number of aircraft handled by the Centers data, and their sources are described. The procedure for selecting reference days on which traffic was minimally impacted by weather is described. Different traffic demand on each reference day of the week, determined by analysis of 42 days of traffic and delay data, was used as the expected traffic demand for each day of the week. Next, the method for computing the weather weighted traffic counts using the expected traffic demand, derived from reference days, and the expanded regions around severe weather cells is discussed. It is shown via a numerical example that this approach improves the dynamic range of the weather weighted traffic counts considerably. Time histories of these new weather weighted traffic counts are used for synthesizing two statistical features, six histogram features and six time domain features. In addition to these enroute weather features, two surface weather features of number of major airports in the United States with high mean winds and low mean visibility are also described. A least squares procedure for establishing a functional relation between the features, using combinations of these features, and system delays is explored using 36 days of data. Best correlations between the estimated delays using the functional relation and the actual delays provided by the Operations Network are obtained with two different combinations of features: 1) six time domain features of weather weighted traffic counts plus two surface weather features, and 2) six histogram features and mean of weather weighted traffic counts along with the two surface weather features. Correlation coefficient values of 0.73 and 0.83 were found in these two instances.

Chatterji, Gano B.↗

NASA NTRS2004

Carrier Modulation Via Waveform Probability Density Function

Beyond the classic modes of carrier modulation by varying amplitude (AM), phase (PM), or frequency (FM), we extend the modulation domain of an analog carrier signal to include a class of general modulations which are distinguished by their probability density function histogram. Separate waveform states are easily created by varying the pdf of the transmitted waveform. Individual waveform states are assignable as proxies for digital ONEs or ZEROs. At the receiver, these states are easily detected by accumulating sampled waveform statistics and performing periodic pattern matching, correlation, or statistical filtering. No fundamental natural laws are broken in the detection process. We show how a typical modulation scheme would work in the digital domain and suggest how to build an analog version. We propose that clever variations of the modulating waveform (and thus the histogram) can provide simple steganographic encoding.

Williams, Glenn L.↗

NASA NTRS2006

Carrier Modulation Via Waveform Probability Density Function

Beyond the classic modes of carrier modulation by varying amplitude (AM), phase (PM), or frequency (FM), we extend the modulation domain of an analog carrier signal to include a class of general modulations which are distinguished by their probability density function histogram. Separate waveform states are easily created by varying the pdf of the transmitted waveform. Individual waveform states are assignable as proxies for digital one's or zero's. At the receiver, these states are easily detected by accumulating sampled waveform statistics and performing periodic pattern matching, correlation, or statistical filtering. No fundamental physical laws are broken in the detection process. We show how a typical modulation scheme would work in the digital domain and suggest how to build an analog version. We propose that clever variations of the modulating waveform (and thus the histogram) can provide simple steganographic encoding.

Williams, Glenn L.↗

NASA NTRS2007

Statistical Analyses of Satellite Cloud Object Data from CERES: Comparison with Cloud-Resolving Model Simulations of Tropical Convective Clouds - Part III

The present study evaluates the ability of a cloud-resolving model (CRM) to simulate the physical properties of tropical deep convective cloud objects identified from a Clouds and the Earth s Radiant Energy System (CERES) data product. The emphasis of this study is the comparisons among the small-, medium- and large-size categories of cloud objects observed during March 1998 and between the large-size categories of cloud objects observed during March 1998 (strong El Ni o) and March 2000 (weak La Ni a). Results from the CRM simulations are analyzed in a way that is consistent with the CERES retrieval algorithm and they are averaged to match the scale of the CERES satellite footprints. Cloud physical properties are analyzed in terms of their summary histograms for each category. It is found that there is a general agreement in the overall shapes of all cloud physical properties between the simulated and observed distributions. Each cloud physical property produced by the CRM also exhibits different degrees of disagreement with observations over different ranges of the property. The simulated cloud tops are generally too high and cloud top temperatures are too low except for the large-size category of March 1998. The probability densities of the simulated top-of-the-atmosphere (TOA) albedos for all four categories are underestimated for high albedos, while those of cloud optical depth are overestimated at its lowest bin. These disagreements are mainly related to uncertainties in the cloud microphysics parameterization and inputs such as cloud ice effective size to the radiation calculation. Summary histograms of cloud optical depth and TOA albedo from the CRM simulations of the large-size category of cloud objects do not differ significantly between the March 1998 and 2000 periods, consistent with the CERES observations. However, the CRM is unable to reproduce the significant differences in the observed cloud top height while it overestimates the differences in the observed outgoing longwave radiation and cloud top temperature between the two periods. Comparisons between the CRM results and the observations for most parameters in March 1998 consistently show that both the simulations and observations have larger differences between the large- and small-size categories than between the large- and medium-size, or between the medium- and small-size categories. However, the simulated cloud properties do not change as much with size as observed. These disagreements are likely related to the spatial averaging of the forcing data and the mismatch in time and in space between the numerical weather prediction model from which the forcing data are produced and the CERES observed cloud systems.

Luo, Yali↗

NASA NTRS2009

Modulation Based on Probability Density Functions

A proposed method of modulating a sinusoidal carrier signal to convey digital information involves the use of histograms representing probability density functions (PDFs) that characterize samples of the signal waveform. The method is based partly on the observation that when a waveform is sampled (whether by analog or digital means) over a time interval at least as long as one half cycle of the waveform, the samples can be sorted by frequency of occurrence, thereby constructing a histogram representing a PDF of the waveform during that time interval.

Williams, Glenn L.↗

NASA NTRS2010

Qualification and Selection of Flight Diode Lasers for Space Applications

The reliability and lifetime of laser diodes is critical to space missions. The Nuclear Spectroscopic Telescope Array (NuSTAR) mission includes a metrology system that is based upon laser diodes. An operational test facility has been developed to qualify and select, by mission standards, laser diodes that will survive the intended space environment and mission lifetime. The facility is situated in an electrostatic discharge (ESD) certified clean-room and consist of an enclosed temperature-controlled stage that can accommodate up to 20 laser diodes. The facility is designed to characterize a single laser diode, in addition to conducting laser lifetime testing on up to 20 laser diodes simultaneously. A standard laser current driver is used to drive a single laser diode. Laser diode current, voltage, power, and wavelength are measured for each laser diode, and a method of selecting the most adequate laser diodes for space deployment is implemented. The method consists of creating histograms of laser threshold currents, powers at a designated current, and wavelengths at designated power. From these histograms, the laser diodes that illustrate a performance that is outside the normal are rejected and the remaining lasers are considered spaceborne candidates. To perform laser lifetime testing, the facility is equipped with 20 custom laser drivers that were designed and built by California Institute of Technology specifically to drive NuSTAR metrology lasers. The laser drivers can be operated in constant-current mode or alternating-current mode. Situated inside the enclosure, in front of the laser diodes, are 20 power-meter heads to record laser power throughout the duration of lifetime testing. Prior to connecting a laser diode to the current source for characterization and lifetime testing, a background program is initiated to collect current, voltage, and resistance. This backstage data collection enables the operational test facility to have full laser diode traceablity.

Liebe, Carl C.↗

NASA NTRS2011

Vertical Structures of Anvil Clouds of Tropical Mesoscale Convective Systems Observed by CloudSat

A global study of the vertical structures of the clouds of tropical mesoscale convective systems (MCSs) has been carried out with data from the CloudSat Cloud Profiling Radar. Tropical MCSs are found to be dominated by cloud-top heights greater than 10 km. Secondary cloud layers sometimes occur in MCSs, but outside their primary raining cores. The secondary layers have tops at 6 8 and 1 3 km. High-topped clouds extend outward from raining cores of MCSs to form anvil clouds. Closest to the raining cores, the anvils tend to have broader distributions of reflectivity at all levels, with the modal values at higher reflectivity in their lower levels. Portions of anvil clouds far away from the raining core are thin and have narrow frequency distributions of reflectivity at all levels with overall weaker values. This difference likely reflects ice particle fallout and therefore cloud age. Reflectivity histograms of MCS anvil clouds vary little across the tropics, except that (i) in continental MCS anvils, broader distributions of reflectivity occur at the uppermost levels in the portions closest to active raining areas; (ii) the frequency of occurrence of stronger reflectivity in the upper part of anvils decreases faster with increasing distance in continental MCSs; and (iii) narrower-peaked ridges are prominent in reflectivity histograms of thick anvil clouds close to the raining areas of connected MCSs (superclusters). These global results are consistent with observations at ground sites and aircraft data. They present a comprehensive test dataset for models aiming to simulate process-based upper-level cloud structure around the tropics.

Hence, Deanna A.↗

NASA NTRS2011

Vertical Structures of Anvil Clouds of Tropical Mesoscale Convective Systems Observed by CloudSat

A global study of the vertical structures of the clouds of tropical mesoscale convective systems (MCSs) has been carried out with data from the CloudSat Cloud Profiling Radar. Tropical MCSs are found to be dominated by cloud-top heights greater than 10 km. Secondary cloud layers sometimes occur in MCSs, but outside their primary raining cores. The secondary layers have tops at 6--8 and 1--3 km. High-topped clouds extend outward from raining cores of MCSs to form anvil clouds. Closest to the raining cores, the anvils tend to have broader distributions of reflectivity at all levels, with the modal values at higher reflectivity in their lower levels. Portions of anvil clouds far away from the raining core are thin and have narrow frequency distributions of reflectivity at all levels with overall weaker values. This difference likely reflects ice particle fallout and therefore cloud age. Reflectivity histograms of MCS anvil clouds vary little across the tropics, except that (i) in continental MCS anvils, broader distributions of reflectivity occur at the uppermost levels in the portions closest to active raining areas; (ii) the frequency of occurrence of stronger reflectivity in the upper part of anvils decreases faster with increasing distance in continental MCSs; and (iii) narrower-peaked ridges are prominent in reflectivity histograms of thick anvil clouds close to the raining areas of connected MCSs (superclusters). These global results are consistent with observations at ground sites and aircraft data. They present a comprehensive test dataset for models aiming to simulate process-based upper-level cloud structure around the tropics.

Yuan, J.↗

NASA NTRS2012

Histogrammatic Method for Determining Relative Abundance of Input Gas Pulse

To satisfy the Major Constituents Analysis (MCA) requirements for the Vehicle Cabin Atmosphere Monitor (VCAM), this software analyzes the relative abundance ratios for N2, O2, Ar, and CO2 as a function of time and constructs their best-estimate mean. A histogram is first built of all abundance ratios for each of the species vs time. The abundance peaks corresponding to the intended measurement and any obfuscating background are then separated via standard peak-finding techniques in histogram space. A voting scheme is then used to include/exclude this particular time sample in the final average based on its membership to the intended measurement or the background population. This results in a robust and reasonable estimate of the abundance of trace components such as CO2 and Ar even in the presence of obfuscating backgrounds internal to the VCAM device. VCAM can provide a means for monitoring the air within the enclosed environments, such as the ISS (International Space Station), Crew Exploration Vehicle (CEV), a Lunar Habitat, or another vehicle traveling to Mars. Its miniature pre-concentrator, gas chromatograph (GC), and mass spectrometer can provide unbiased detection of a large number of organic species as well as MCA analysis. VCAM s software can identify the concentration of trace chemicals and whether the chemicals are on a targeted list of hazardous compounds. This innovation s performance and reliability on orbit, along with the ground team s assessment of its raw data and analysis results, will validate its technology for future use and development.

Mandrake, Lukas↗

NASA NTRS2013

The Impact of Model and Rainfall Forcing Errors on Characterizing Soil Moisture Uncertainty in Land Surface Modeling

The contribution of rainfall forcing errors relative to model (structural and parameter) uncertainty in the prediction of soil moisture is investigated by integrating the NASA Catchment Land Surface Model (CLSM), forced with hydro-meteorological data, in the Oklahoma region. Rainfall-forcing uncertainty is introduced using a stochastic error model that generates ensemble rainfall fields from satellite rainfall products. The ensemble satellite rain fields are propagated through CLSM to produce soil moisture ensembles. Errors in CLSM are modeled with two different approaches: either by perturbing model parameters (representing model parameter uncertainty) or by adding randomly generated noise (representing model structure and parameter uncertainty) to the model prognostic variables. Our findings highlight that the method currently used in the NASA GEOS-5 Land Data Assimilation System to perturb CLSM variables poorly describes the uncertainty in the predicted soil moisture, even when combined with rainfall model perturbations. On the other hand, by adding model parameter perturbations to rainfall forcing perturbations, a better characterization of uncertainty in soil moisture simulations is observed. Specifically, an analysis of the rank histograms shows that the most consistent ensemble of soil moisture is obtained by combining rainfall and model parameter perturbations. When rainfall forcing and model prognostic perturbations are added, the rank histogram shows a U-shape at the domain average scale, which corresponds to a lack of variability in the forecast ensemble. The more accurate estimation of the soil moisture prediction uncertainty obtained by combining rainfall and parameter perturbations is encouraging for the application of this approach in ensemble data assimilation systems.

Maggioni, V.↗

NASA NTRS2015

Simulations of Cloud-Radiation Interaction Using Large-Scale Forcing Derived from the CINDY/DYNAMO Northern Sounding Array

The recently completed CINDY/DYNAMO field campaign observed two Madden-Julian oscillation (MJO) events in the equatorial Indian Ocean from October to December 2011. Prior work has indicated that the moist static energy anomalies in these events grew and were sustained to a significant extent by radiative feedbacks. We present here a study of radiative fluxes and clouds in a set of cloud-resolving simulations of these MJO events. The simulations are driven by the large-scale forcing data set derived from the DYNAMO northern sounding array observations, and carried out in a doubly periodic domain using the Weather Research and Forecasting (WRF) model. Simulated cloud properties and radiative fluxes are compared to those derived from the S-PolKa radar and satellite observations. To accommodate the uncertainty in simulated cloud microphysics, a number of single-moment (1M) and double-moment (2M) microphysical schemes in the WRF model are tested. The 1M schemes tend to underestimate radiative flux anomalies in the active phases of the MJO events, while the 2M schemes perform better, but can overestimate radiative flux anomalies. All the tested microphysics schemes exhibit biases in the shapes of the histograms of radiative fluxes and radar reflectivity. Histograms of radiative fluxes and brightness temperature indicate that radiative biases are not evenly distributed; the most significant bias occurs in rainy areas with OLR less than 150 W/ cu sq in the 2M schemes. Analysis of simulated radar reflectivities indicates that this radiative flux uncertainty is closely related to the simulated stratiform cloud coverage. Single-moment schemes underestimate stratiform cloudiness by a factor of 2, whereas 2M schemes simulate much more stratiform cloud.

Indian Ocean↗

NASA NTRS2018

Fragment Size Distribution for Ice Particle Impacts on a Glass Plate

This work presents the results of an experimental study of ice particle impacts on a flat glass plate. The experiment was conducted at the Ballistics Impact Laboratory of NASA Glenn Research Center. The main objective of the experiment was to gain understanding about the modifications needed to the experimental configuration for a future parametric study at a larger range of values for particle diameters and other parameters. This was achieved by studying the effect of the velocity of an impacting ice particle on the post-impact fragment size and distribution for a reduced range of impacting particle diameters. Pre-impact particle diameter and velocity data were captured with a high-speed side camera. Post-impact fragment data were captured in a single frame with a 29-megapixel camera located above and normal to the target. Repeat runs were conducted for ice particles with diameters ranging from 1.7 to 2.9 millimeters, impacting at velocities between 39 and 98 meters per second. The fragment areas were measured, and the corresponding equivalent diameters and histogram distributions were calculated. Analysis of the data showed that the average equivalent diameter for the fragments in a run was an order of magnitude smaller than the diameter of the impacting ice particle. The histograms for equivalent diameter distribution were nonnormal with long tails, with most of the fragments having equivalent diameters concentrated toward the minimum value of the fragment size that could be resolved. Factors affecting the accuracy of the data during the digital imaging analysis were identified. Needed modifications to the setup to handle small size ice particles and other testing conditions were also identified.

engine icing↗

NASA NTRS2018

ICESat-2 Photon Classification: Finding Signal Photons in the ATL03 Geolocated Photon Data Product

ICESat-2 carries NASA's next-generation laser altimeter, ATLAS, (Advanced Topographic Laser Altimeter System), designed to measure changes in ice sheet height, sea ice freeboard, and vegetation canopy height. ATLAS contains a photon-counting lidar which transmits green (532-nm) pulses at 10kHz. Each pulse is split into 3 pairs of beams (one strong and one weak). Approximately 1014 photons per pulse travel from ATLAS through the atmosphere to reflect off the Earth's surface. Some return back into the ATLAS telescope where they are recorded. Photons from sunlight and instrument noise at the same wavelength are also recorded. The flight software time tags all photons within a 500m to 6 km range window and generates histograms. Using the histograms, it selects a telemetry window which varies from 20m over flat surfaces to hundreds of meters over rougher terrain. ATL03 contains the time, height (relative to the WGS-84 ellipsoid), latitude and longitude of every photon within the telemetry window. The basic challenge is to determine which of these photons were reflected off the surface. We have developed an algorithm that identifies these signal photons and assigns a confidence level (low, medium, or high) to each signal photon based on the signal to noise ratio. We present an overview of the signal identification algorithm and show the results on actual ICESat-2 data over ice sheet, sea ice, vegetated, and water surfaces. Higher level ATLAS products work with aggregations of the photons in order to determine the ellipsoidal height of the Earth, canopy height and structure, and other quantities of geophysical interest.

Remote Sensing↗

NASA NTRS2019

Characterizing the System Impulse Response Function from Photon-Counting LiDAR Data

NASA's Multiple Altimeter Beam Experimental LiDAR (MABEL) is an aircraft-based photon-counting laser altimeter designed as a simulator to test measurement techniques and algorithms for Advanced Topographic Laser Altimeter System (ATLAS), the sole instrument on NASA's Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) mission. By measuring the time of flight, pointing angle, and absolute position for individual photons, ICESat-2 provides detailed elevation measurements of earth's surface. Calculating accurate and precise elevations requires an understanding of how photons interact with surfaces, and characterization of the photon distribution after returning from surfaces. Neither MABEL nor ATLAS records the transmitted laser pulse shape, relying instead on aggregating several pulses worth of photons, often using histograms, to characterize the pulse shape. In this paper, we assess the limitations of using histograms and propose a more robust method to describe MABEL's system impulse-response function using an exponentially modified Gaussian distribution. We also provide standard error estimates for the arithmetic mean and standard deviation calculations, and for exponentially modified Gaussian parameters using a Monte Carlo sensitivity analysis. We apply this method to photon returns from a sea ice lead and from a dry salt lake bed as case studies for estimating the standard error associated with sample size for the arithmetic mean and standard deviation, and for the exponentially modified Gaussian parameters. We use these standard errors to calculate the minimum number of photons required to find both Gaussian and exponentially modified Gaussian distribution parameters within 3 cm of their parent population values.

photoncounting↗