Search NASA⌕ Search

SEARCH · Search NASA

Results for “sampling algorithm”

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 685 records · Page 38

High resolution observations of the L1551 bipolar outflow

The nearby dark cloud Lynds 1551 contains one of the closest examples of a well-collimated bipolar molecular outflow. This source has the largest angular size of any known outflow and was the first bipolar outflow to be detected. The outflow originates from a low-luminosity young stellar object, IRS-5. Optical and radio continuum observations show the presence of a highly collimated, ionized stellar wind orginating from close to IRS-5 and aligned with the molecular outflow. However, we have little information on the actual mechanism that generates the stellar wind and collimates it into opposed jets. The Very Large Array (VLA) observations indicate that the winds originate within 10(15) cm of IRS-5, unfortunately at a size scale difficult to resolve. For these reasons, observations of the structure and dynamics of the hypersonic molecular gas may provide valuable information on the origin and evolution of these outflows. In addition, the study of the impact of the outflowing gas on the surrounding molecular material is essential to understand the consequence these outflows have on the evolution and star formation history of the entire cloud. Moriarty-Schieven et al. (1986) obtained a oversampled map of the CO emission of a portion of both the blueshifted and redshifted outflows in LI551 using Five College Radio Astronomy Observatory 14 m telescope. The oversampled maps have been reconstructed to an effective angular resolution of 20 arcsec using a maximum entropy algorithm. A continuation of the study of Moriarty-Schieven et al. is presented. The entire L1551 outflow has now been mapped at 12 arcsec sampling requiring roughly 4000 spectra. This data has been constructed to 20 arcsec resolution to provide the first high resolution picture of the entire L1551 outflow. This new data has shown that the blueshifted lobe is more extended than previously thought and has expanded downstream sufficiently to break out of the dense molecular cloud, but the redshifted outflow is still confined within the molecular cloud. Details of the structure and kinematics of the high velocity gas are used to test the various models of the origin and evolution of outflows.

Snell, R.↗

Studies in astronomical time series analysis. III - Fourier transforms, autocorrelation functions, and cross-correlation functions of unevenly spaced data

This paper develops techniques to evaluate the discrete Fourier transform (DFT), the autocorrelation function (ACF), and the cross-correlation function (CCF) of time series which are not evenly sampled. The series may consist of quantized point data (e.g., yes/no processes such as photon arrival). The DFT, which can be inverted to recover the original data and the sampling, is used to compute correlation functions by means of a procedure which is effectively, but not explicitly, an interpolation. The CCF can be computed for two time series not even sampled at the same set of times. Techniques for removing the distortion of the correlation functions caused by the sampling, determining the value of a constant component to the data, and treating unequally weighted data are also discussed. FORTRAN code for the Fourier transform algorithm and numerical examples of the techniques are given.

Scargle, Jeffrey D.↗

VLSI chip-set for data compression using the Rice algorithm

A full custom VLSI implementation of a data compression encoder and decoder which implements the lossless Rice data compression algorithm is discussed in this paper. The encoder and decoder reside on single chips. The data rates are to be 5 and 10 Mega-samples-per-second for the decoder and encoder respectively.

Venbrux, J.↗

On the optimality of a universal noiseless coder

Rice developed a universal noiseless coding structure that provides efficient performance over an extremely broad range of source entropy. This is accomplished by adaptively selecting the best of several easily implemented variable length coding algorithms. Variations of such noiseless coders have been used in many NASA applications. Custom VLSI coder and decoder modules capable of processing over 50 million samples per second have been fabricated and tested. In this study, the first of the code options used in this module development is shown to be equivalent to a class of Huffman code under the Humblet condition, for source symbol sets having a Laplacian distribution. Except for the default option, other options are shown to be equivalent to the Huffman codes of a modified Laplacian symbol set, at specified symbol entropy values. Simulation results are obtained on actual aerial imagery over a wide entropy range, and they confirm the optimality of the scheme. Comparison with other known techniques are performed on several widely used images and the results further validate the coder's optimality.

Yeh, Pen-Shu↗

Kravchuk functions for the finite oscillator approximation

Kravchuk orthogonal functions - Kravchuk polynomials multiplied by the square root of the weight function - simplify the inversion algorithm for the analysis of discrete, finite signals in harmonic oscillator components. They can be regarded as the best approximation set. As the number of sampling points increases, the Kravchuk expansion becomes the standard oscillator expansion.

Atakishiyev, Natig M.↗

Learning In networks

Intelligent systems require software incorporating probabilistic reasoning, and often times learning. Networks provide a framework and methodology for creating this kind of software. This paper introduces network models based on chain graphs with deterministic nodes. Chain graphs are defined as a hierarchical combination of Bayesian and Markov networks. To model learning, plates on chain graphs are introduced to model independent samples. The paper concludes by discussing various operations that can be performed on chain graphs with plates as a simplification process or to generate learning algorithms.

Buntine, Wray L.↗

Application of Self Nulling Eddy Current Probe Technique to the Detection of Fatigue Crack Initiation and Control of Test Procedures

A major part of fracture mechanics is concerned with studying the initiation and propagation of fatigue cracks. This typically requires constant monitoring of crack growth during fatigue cycles and the knowledge of the precise location of the crack tip at any given time. One technique currently available for measuring fatigue crack length is the Potential Drop method. The method, however, may be inaccurate if the direction of crack growth deviates considerably from what was assumed initially or the curvature of the crack becomes significant. Another popular approach is to optically view the crack using a high magnification microscope, but this entails a person constantly monitoring it. The present proposed technique uses an automated scheme, in order to eliminate the need for a person to constantly monitor the experiment. Another technique under development elsewhere is to digitize an optical image of the test specimen surface and then apply a pattern recognition algorithm to locate the crack tip. A previous publication showed that the self nulling eddy current probe successfully tracked a simulated crack in an aluminum sample. This was the impetus to develop an online real time crack monitoring system. An automated system has been developed which includes a two axis scanner mounted on the tensile testing machine, the probe and its instrumentation and a personal computer (PC) to communicate and control all the parameters. The system software controls the testing parameters as well as monitoring the fatigue crack as it propagates. This paper will discuss the experimental setup in detail and demonstrate its capabilities. A three dimensional finite element model is utilized to model the magnetic field distribution due to the probe and how the probe voltage changes as it scans the crack. Experimental data of the probe for different samples under zero load, static load and high cycle fatigue load will be discussed. The final section summarizes the major accomplishments of the present work, the elements of the future R&D needs and the advantages and disadvantages of using this system in the laboratory and field.

Namkung, M.↗

ISSM: Ice Sheet System Model

In order to have the capability to use satellite data from its own missions to inform future sea-level rise projections, JPL needed a full-fledged ice-sheet/iceshelf flow model, capable of modeling the mass balance of Antarctica and Greenland into the near future. ISSM was developed with such a goal in mind, as a massively parallelized, multi-purpose finite-element framework dedicated to ice-sheet modeling. ISSM features unstructured meshes (Tria in 2D, and Penta in 3D) along with corresponding finite elements for both types of meshes. Each finite element can carry out diagnostic, prognostic, transient, thermal 3D, surface, and bed slope simulations. Anisotropic meshing enables adaptation of meshes to a certain metric, and the 2D Shelfy-Stream, 3D Blatter/Pattyn, and 3D Full-Stokes formulations capture the bulk of the ice-flow physics. These elements can be coupled together, based on the Arlequin method, so that on a large scale model such as Antarctica, each type of finite element is used in the most efficient manner. For each finite element referenced above, ISSM implements an adjoint. This adjoint can be used to carry out model inversions of unknown model parameters, typically ice rheology and basal drag at the ice/bedrock interface, using a metric such as the observed InSAR surface velocity. This data assimilation capability is crucial to allow spinning up of ice flow models using available satellite data. ISSM relies on the PETSc library for its vectors, matrices, and solvers. This allows ISSM to run efficiently on any parallel platform, whether shared or distrib- ISSM: Ice Sheet System Model NASA's Jet Propulsion Laboratory, Pasadena, California uted. It can run on the largest clusters, and is fully scalable. This allows ISSM to tackle models the size of continents. ISSM is embedded into MATLAB and Python, both open scientific platforms. This improves its outreach within the science community. It is entirely written in C/C++, which gives it flexibility in its design, and the power/speed that C/C++ allows. ISSM is svn (subversion) hosted, on a JPL repository, to facilitate its development and maintenance. ISSM can also model propagation of rifts using contact mechanics and mesh splitting, and can interface to the Dakota software. To carry out sensitivity analysis, mesh partitioning algorithms are available, based on the Scotch, Chaco, and Metis partitioners that ensure equal area mesh partitions can be done, which are then usable for sampling and local reliability methods.

Larour, Eric↗

Hardware Implementation of Lossless Adaptive and Scalable Hyperspectral Data Compression for Space

On-board lossless hyperspectral data compression reduces data volume in order to meet NASA and DoD limited downlink capabilities. The technique also improves signature extraction, object recognition and feature classification capabilities by providing exact reconstructed data on constrained downlink resources. At JPL a novel, adaptive and predictive technique for lossless compression of hyperspectral data was recently developed. This technique uses an adaptive filtering method and achieves a combination of low complexity and compression effectiveness that far exceeds state-of-the-art techniques currently in use. The JPL-developed 'Fast Lossless' algorithm requires no training data or other specific information about the nature of the spectral bands for a fixed instrument dynamic range. It is of low computational complexity and thus well-suited for implementation in hardware. A modified form of the algorithm that is better suited for data from pushbroom instruments is generally appropriate for flight implementation. A scalable field programmable gate array (FPGA) hardware implementation was developed. The FPGA implementation achieves a throughput performance of 58 Msamples/sec, which can be increased to over 100 Msamples/sec in a parallel implementation that uses twice the hardware resources This paper describes the hardware implementation of the 'Modified Fast Lossless' compression algorithm on an FPGA. The FPGA implementation targets the current state-of-the-art FPGAs (Xilinx Virtex IV and V families) and compresses one sample every clock cycle to provide a fast and practical real-time solution for space applications.

FPGA implementation↗

Efficient Autonomous Learning for Statistical Pattern Recognition

We describe a neural network learning algorithm that implements differential learning in a generalized backpropagation framework. The algorithm regulates model complexity during the learning procedure, generating the best low-complexity approximation to the Bayer-optimal classifier allowed by the training sample.

Pattern↗

An Ensemble of Bayesian Neural Networks for Exoplanetary Atmospheric Retrieval

Machine learning (ML) is now used in many areas of astrophysics, from detecting exoplanets in Kepler transit signals to removing telescope systematics. Recent work demonstrated the potential of using ML algorithms for atmospheric retrieval by implementing a random forest (RF) to perform retrievals in seconds that are consistent with the traditional, computationally expensive nested-sampling retrieval method. We expand upon their approach by presenting a new ML model, plan-net, based on an ensemble of Bayesian neural networks (BNNs) that yields more accurate inferences than the RF for the same data set of synthetic transmission spectra. We demonstrate that an ensemble provides greater accuracy and more robust uncertainties than a single model. In addition to being the first to use BNNs for atmospheric retrieval, we also introduce a new loss function for BNNs that learns correlations between the model outputs. Importantly, we show that designing ML models to explicitly incorporate domain-specific knowledge both improves performance and provides additional insight by inferring the covariance of the retrieved atmospheric parameters. We apply plan-net to the Hubble Space Telescope Wide Field Camera 3 transmission spectrum for WASP-12b and retrieve an isothermal temperature and water abundance consistent with the literature. We highlight that our method is flexible and can be expanded to higher resolution spectra and a larger number of atmospheric parameters.

Adam D. Cobb↗

Powder River Basin Water Resources: Mapping Russian Olive in the Powder River Basin to Inform Invasive Species Management

Since its introduction in the late 1800s, Elaeagnus augustifolia (Russian olive) has become a widespread invasive shrub that poses a threat to native riparian species in the United States by competing with native riparian plants for space and resources. To date, limited information on the distribution of Russian olive in the Powder River Basin of Montana and Wyoming have hampered management efforts and decision making. Here, we detect and model the distribution of Russian Olive using field surveys, ocular sampling, and variables from Landsat 8 Operational Land Imager (OLI), Sentinel-2 MultiSpectral Instrument (MSI), and Shuttle Radar Topography Mission (SRTM) using the Random Forest algorithm. We derived topographic, spectral, and hydrological variables from Landsat 8 OLI, Sentinel-2 MSI, and SRTM to utilize as model inputs. The team was able to successfully create a spectral Russian olive detection map for the Powder River Basin (RMSE =15.44%, R2 = 0.6482). The team also examined change in stream channel geomorphology from 1984-2020 in a time-series analysis using Landsat visible imagery and the RivMap MATLAB package and found little change. Our results will help our partners at the Powder River County Weed Board, Gay Ranch, United States Geological Survey, and University of Northern Colorado to locate and prioritize areas for riparian habitat restoration and to understand the region’s hydrology and geomorphology.

Catherine Buczek↗

Influence of atmospherically induced random wave fronts on diffraction imagery - A computer simulation model for testing image reconstruction algorithms

This paper is devoted to the development of a two-dimensional computer-simulation model that is based on the rigid constraints of optical diffraction theory with careful attention paid to the generation of sample realizations of Gaussian-distributed, spatially random, isotropic wave fronts that have zero-mean and prescribed-covariance functions. Given a sample realization of the wave front, the corresponding centered point-spread function and optical-transfer function are evaluated. A detailed study is made of the statistics of random wave-front tilt, point-spread function, modulus squared of transfer function, and phase of transfer function.

Barakat, Richard↗

Noise and drift analysis of non-equally spaced timing data

Generally, it is possible to obtain equally spaced timing data from oscillators. The measurement of the drifts and noises affecting oscillators is then performed by using a variance (Allan variance, modified Allan variance, or time variance) or a system of several variances (multivariance method). However, in some cases, several samples, or even several sets of samples, are missing. In the case of millisecond pulsar timing data, for instance, observations are quite irregularly spaced in time. Nevertheless, since some observations are very close together (one minute) and since the timing data sequence is very long (more than ten years), information on both short-term and long-term stability is available. Unfortunately, a direct variance analysis is not possible without interpolating missing data. Different interpolation algorithms (linear interpolation, cubic spline) are used to calculate variances in order to verify that they neither lose information nor add erroneous information. A comparison of the results of the different algorithms is given. Finally, the multivariance method was adapted to the measurement sequence of the millisecond pulsar timing data: the responses of each variance of the system are calculated for each type of noise and drift, with the same missing samples as in the pulsar timing sequence. An estimation of precision, dynamics, and separability of this method is given.

Vernotte, F.↗

On the effect of timing errors in run length codes

Many redundancy removal algorithms employ some sort of run length code. Blocks of timing words are coded with synchronization words inserted between blocks. The probability of incorrectly reconstructing a sample because of a channel error in the timing data is a monotonically nondecreasing function of time since the last synchronization word. In this paper we compute the 'probability that the accumulated magnitude of timing errors equal zero' as a function of time since the last synchronization word for a zero-order predictor (ZOP). The result is valid for any data source that can be modeled by a first-order Markov chain and any digital channel that can be modeled by a channel transition matrix. An example is presented.

Wilkins, L. C.↗

LACIE performance predictor FOC users manual

The LACIE Performance Predictor (LPP) is a computer simulation of the LACIE process for predicting worldwide wheat production. The simulation provides for the introduction of various errors into the system and provides estimates based on these errors, thus allowing the user to determine the impact of selected error sources. The FOC LPP simulates the acquisition of the sample segment data by the LANDSAT Satellite (DAPTS), the classification of the agricultural area within the sample segment (CAMS), the estimation of the wheat yield (YES), and the production estimation and aggregation (CAS). These elements include data acquisition characteristics, environmental conditions, classification algorithms, the LACIE aggregation and data adjustment procedures. The operational structure for simulating these elements consists of the following key programs: (1) LACIE Utility Maintenance Process, (2) System Error Executive, (3) Ephemeris Generator, (4) Access Generator, (5) Acquisition Selector, (6) LACIE Error Model (LEM), and (7) Post Processor.

Source record↗

Landsat-4 thematic mapper and thematic mapper simulator data for a porphyry copper deposit

Aircraft thematic mapper (TM) data were analyzed to evaluate the potential utility of the Landsat-4 thematic mapper for geologic mapping and detection of hydrothermal alteration zones in the Silver Bell porphyry copper deposit in southern Arizona. The data allow a comparison between aircraft TV simulator data and the Landsat-4 TM satellite data which possess similar spectral bands. A color rationcomposite of 30-m pixels was resampled, in order to clearly define a number of hydroxyl bearing minerals, (kaolinite, sericite, white mica), pyrite and iron oxide/hydroxide minerals. The iron oxide minerals have diagnostic absorption bands in the 0.45 and 0.85 micron regions of the spectrum, and the hydrous minerals are characterized by an absorption in the 2.2 micron region. The position of the spectral bands allow the TM to identify regions of hydrothermal alteration without resorting to a data processing algorithm. The comparison of the aircraft and Landsat-4 TM data showed considerable agreement, and confirmed the utility of TM data for identifying hydrothermal alteration zones. Samples of some color TM images are provided.

Abrams, M. J.↗

Development and testing of a simple assimilation technique to derive average wind fields from simulated scatterometer data

A simple algorithm is developed and tested to derive a regularly spaced wind field in a limited area from simulated multiorbit scatterometer data. The data are generated by sampling a time-varying known wind field, the 1000-mb FGGE data, from a simulated scatterometer. A simple assimilation technique is used to derive a regularly spaced (100-km grid) wind field representation of two-day averages from the simulated data. Several test cases are considered, and it is noted that the assimilation technique might be applicable to large-scale ocean or atmospheric models.

Legler, D. M.↗