Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sequence tracking method”

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 37 records · Page 2

Application of an automatic cloud tracking technique to Meteosat water vapor and infrared observations

The automatic cloud tracking system was applied to METEOSAT 6.7 micrometers water vapor measurements to learn whether the system can track the motions of water vapor patterns. Data for the midlatitudes, subtropics, and tropics were selected from a sequence of METEOSAT pictures for 25 April 1978. Trackable features in the water vapor patterns were identified using a clustering technique and the features were tracked by two different methods. In flat (low contrast) water vapor fields, the automatic motion computations were not reliable, but in areas where the water vapor fields contained small scale structure (such as in the vicinity of active weather phenomena) the computations were successful. Cloud motions were computed using METEOSAT infrared observations (including tropical convective systems and midlatitude jet stream cirrus).

Endlich, R. M.↗

Dense Feature Tracking of Atmospheric Winds with Deep Optical Flow

Atmospheric winds are a key physical phenomenon impacting natural hazards, energy transport, ocean currents, large-scale circulation, and ecosystem fluxes. Observing winds is a complex process and presents a large gap in NASA’s Earth Observation System. Atmospheric motion vectors (AMVs) aim to fill this gap by making numerical estimates of cloud movement between sequences of multi-spectral satellite images, tracking clouds and water vapor. Recent imaging hardware and software advancements have enabled the use of numerical optical flow techniques to produce accurate and dense vector fields outperforming traditional methods. This work presents WindFlow as the first machine learning based system for feature tracking atmospheric motion using optical flow. Due to the lack of large-scale satellite-based observations, we leverage high-resolution numerical simulations from NASA's GEOS-5 Nature Run to perform supervised learning and transfer to satellite images. We demonstrate that our approach using deep learning based optical flow scales to ultra-high-resolution images of size 2881x5760 with less than 1 m/s bias and 2.5 m/s average error. Four network and learning architectures are compared and it is found that recurrent all-pairs field transforms (RAFT) produces the lowest errors on all metrics for wind speed and direction. Results on held out numerical outputs shows RAFT's good performance in each of the spatial, temporal, and physical dimensions. A comparison between WindFlow and an operational AMV product against rawinsonde observations show that RAFT transfers across simulations and thermal infrared satellite observations. This work shows that machine learning based optical flow is an efficient approach to generating robust feature tracking for AMVs consistently over large regions.

Atmospheric winds↗

Operation of Direct Drive Systems: Experiments in Peak Power Tracking and Multi-Thruster Control

Direct-drive power and propulsion systems have the potential to significantly reduce the mass of high-power solar electric propulsion spacecraft, among other advantages. Recent experimental direct-drive work has significantly mitigated or retired the technical risks associated with single-thruster operation, so attention is now moving toward systems-level areas of interest. One of those areas is the use of a Hall thruster system as a peak power tracker to fully use the available power from a solar array. A simple and elegant control based on the incremental conductance method, enhanced by combining it with the unique properties of Hall thruster systems, is derived here and it is shown to track peak solar array power very well. Another area of interest is multi-thruster operation and control. Dualthruster operation was investigated in a parallel electrical configuration, with both thrusters operating from discharge power provided by a single solar array. Startup and shutdown sequences are discussed, and it is shown that multi-thruster operation and control is as simple as for a single thruster. Some system architectures require operation of multiple cathodes while they are electrically connected together. Four different methods to control the discharge current emitted by individual cathodes in this configuration are investigated, with cathode flow rate control appearing to be advantageous. Dual-parallel thruster operation with equal cathode current sharing at total powers up to 10 kW is presented.

solar electric propulsion (SEP) vehicles↗

Joint Carrier-Phase Synchronization and LDPC Decoding

A method has been proposed to increase the degree of synchronization of a radio receiver with the phase of a suppressed carrier signal modulated with a binary- phase-shift-keying (BPSK) or quaternary- phase-shift-keying (QPSK) signal representing a low-density parity-check (LDPC) code. This method is an extended version of the method described in Using LDPC Code Constraints to Aid Recovery of Symbol Timing (NPO-43112), NASA Tech Briefs, Vol. 32, No. 10 (October 2008), page 54. Both methods and the receiver architectures in which they would be implemented belong to a class of timing- recovery methods and corresponding receiver architectures characterized as pilotless in that they do not require transmission and reception of pilot signals. The proposed method calls for the use of what is known in the art as soft decision feedback to remove the modulation from a replica of the incoming signal prior to feeding this replica to a phase-locked loop (PLL) or other carrier-tracking stage in the receiver. Soft decision feedback refers to suitably processed versions of intermediate results of iterative computations involved in the LDPC decoding process. Unlike a related prior method in which hard decision feedback (the final sequence of decoded symbols) is used to remove the modulation, the proposed method does not require estimation of the decoder error probability. In a basic digital implementation of the proposed method, the incoming signal (having carrier phase theta theta (sub c) plus noise would first be converted to inphase (I) and quadrature (Q) baseband signals by mixing it with I and Q signals at the carrier frequency [wc/(2 pi)] generated by a local oscillator. The resulting demodulated signals would be processed through one-symbol-period integrate and- dump filters, the outputs of which would be sampled and held, then multiplied by a soft-decision version of the baseband modulated signal. The resulting I and Q products consist of terms proportional to the cosine and sine of the carrier phase cc as well as correlated noise components. These products would be fed as inputs to a digital PLL that would include a number-controlled oscillator (NCO), which provides an estimate of the carrier phase, theta(sub c).

Simon, Marvin↗

Method of and means for testing a tape record/playback system

A tape record/playback system was tested by first deriving an analog test signal and a band-limited digital reference signal from a pseudo-noise sequence generator driven by a clock signal. It recorded the signals on respective tracks of the system during operation in a record mode. During the playback mode of operation of the system, a delayed analog reference signal without time base variations was reconstructed from the played back reference signal. It was compared with the played back test signal in order to obtain an error signal that was a measure of the performance of the system.

Wallace, G. R.↗

Slip Validation and Prediction for Mars Exploration Rovers

This paper presents a novel technique to validate and predict the rover slips on Martian surface for NASA's Mars Exploration Rovers. Different from the traditional approach, the proposed method uses the actual velocity profile of the wheels and the digital elevation map (DEM) from the stereo images of the terrain to formulate simplified equations of motion of the rovers. A weighted factor to the wheel-ground speed from the empirical data comprises the velocity equations of the simplified differential-algebraic system of the rover motion. Applying the discretization operator to these equations, the full kinematics state of the rover is then resolved by the configuration kinematics solution in the Robot Sequencing and Visualization Program (RSVP). This method produced accurate simulation of the rover movements compared with these of the earth testing vehicle. Using the telemetry from the onboard Visual Odometry, the simulated rover path also compares well with the actual track of the vehicle. Preliminary results indicated that the proposed computational method is very effective in planning the path of the rovers on the high-slope areas.

Mars Rovers↗

Observing rotation and deformation of sea ice with synthetic aperture radar

The ESA's ERS-1 satellite will carry SARs over the polar regions; an important component in the use of these data is an automated scheme for the extraction of sea ice velocity fields from a sequence of SAR images of the same geographical region. The image pyramid area-correlation hierarchical method is noted to be vulnerable to uncertainties for sea ice rotations greater than 10-15 deg between SAR observations. Rotation-invariant methods can successfully track isolated floes in the marginal ice zone. Hu's (1962) invariant moments are also worth considering as a possible basis for rotation-invariant tracking methods. Feature tracking is inherently robust for tracking rotating sea ice, but is limited when features are floe-lead boundaries. A variety of techniques appears neccessary.

Vesecky, J. F.↗

Kalman filter based range estimation for autonomous navigation using imaging sensors

Rotorcraft operating in high-threat environments fly close to the surface of the earth to utilize surrounding terrain, vegetation, or man-made objects to minimize the risk of being detected by the enemy. Two basic requirements for obstacle avoidance are detection and range estimation of the object from the current rotorcraft position. There are many approaches to the estimation of range using a sequence of images. The approach used in this analysis differes from previous methods in two significant ways: an attempt is not made to estimate the rotorcraft's motion from the images; and the interest lies in recursive algorithms. The rotorcraft parameters are assumed to be computed using an onboard inertial navigation system. Given a sequence of images, using image-object differential equations, a Kalman filter (Sridhar and Phatak, 1988) can be used to estimate both the relative coordinates and the earth coordinates of the objects on the ground. The Kalman filter can also be used in a predictive mode to track features in the images, leading to a significant reduction of search effort in the feature extraction step of the algorithm. The purpose is to summarize early results obtained in extending the Kalman filter for use with actual image sequences. The experience gained from the application of this algorithm to real images is very valuable and is a necessary step before proceeding to the estimation of range during low-altitude curvilinear flight. A simple recursive method is presented to estimate range to objects using a sequence of images. The method produces good range estimates using real images in a laboratory set up and needs to be evaluated further using several different image sequences to test its robustness. The feature generation part of the algorithm requires further refinement on the strategies to limit the number of features (Sridhar and Phatak, 1989). The extension of the work reported here to curvilinear flight may require the use of the extended Kalman filter.

Sridhar, Banavar↗

Real-Time Feature Tracking Using Homography

This software finds feature point correspondences in sequences of images. It is designed for feature matching in aerial imagery. Feature matching is a fundamental step in a number of important image processing operations: calibrating the cameras in a camera array, stabilizing images in aerial movies, geo-registration of images, and generating high-fidelity surface maps from aerial movies. The method uses a Shi-Tomasi corner detector and normalized cross-correlation. This process is likely to result in the production of some mismatches. The feature set is cleaned up using the assumption that there is a large planar patch visible in both images. At high altitude, this assumption is often reasonable. A mathematical transformation, called an homography, is developed that allows us to predict the position in image 2 of any point on the plane in image 1. Any feature pair that is inconsistent with the homography is thrown out. The output of the process is a set of feature pairs, and the homography. The algorithms in this innovation are well known, but the new implementation improves the process in several ways. It runs in real-time at 2 Hz on 64-megapixel imagery. The new Shi-Tomasi corner detector tries to produce the requested number of features by automatically adjusting the minimum distance between found features. The homography-finding code now uses an implementation of the RANSAC algorithm that adjusts the number of iterations automatically to achieve a pre-set probability of missing a set of inliers. The new interface allows the caller to pass in a set of predetermined points in one of the images. This allows the ability to track the same set of points through multiple frames.

Clouse, Daniel S.↗

Lunar orbit determination in the presence of unmodeled accelerations

A technique for estimating the state of an artificial satellite in the presence of unmodeled accelerations is presented. The unmodeled acceleration is approximated by a first-order Gauss-Markov sequence which can be separated into a timewise-correlated component and a purely random component. Using this approximation, a sequential procedure for estimating the position, velocity, and the unmodeled acceleration is developed. The method is evaluated by reducing range-rate observations obtained by tracking the Apollo 10 and 11 spacecraft during the lunar-orbit phase of the mission. Numerical results are presented which show that the observation residual pattern lies within the observation noise standard deviation. The values of the estimated components of the unmodeled acceleration are repeatable from orbit to orbit within a given mission and from mission to mission when the same ground track is covered. Finally, the variation in the radial component of the unmodeled acceleration shows a high correlation with the reported location of the lunar surface mascons.

Ingram, D. S.↗

Computational fluid dynamics combustion analysis evaluation

This study involves the development of numerical modelling in spray combustion. These modelling efforts are mainly motivated to improve the computational efficiency in the stochastic particle tracking method as well as to incorporate the physical submodels of turbulence, combustion, vaporization, and dense spray effects. The present mathematical formulation and numerical methodologies can be casted in any time-marching pressure correction methodologies (PCM) such as FDNS code and MAST code. A sequence of validation cases involving steady burning sprays and transient evaporating sprays will be included.

Kim, Y. M.↗

Enhanced Pseudo-Waypoint Guidance for Spacecraft Maneuvers

An enhanced version of the scheme reported in "Pseudo-Waypoint Guidance for Proximity Spacecraft Maneuvers" (NPO-42753), NASA Tech Briefs, Vol. 31, No. 6 (June 2007), page 73 was developed. To recapitulate: the scheme provides algorithms for guidance and control (G&C) of a spacecraft maneuvering near a small astronomical body. The open-loop guidance problem is solved in advance or in real time by use of the pseudo-waypoint generation (PWG) method. Feedback control is implemented to track PWG trajectories, in a manner that enables updating of G&C in a model-predictive manner. The scheme includes silent periods following each thruster firing. The original version of the scheme provides for a fire-first, followed-by-silence sequence, which is disadvantageous in that the silence after final firing precludes reduction of any remaining velocity error -- an unacceptable result in the case of a maneuver for which a specific final velocity is required. In the enhanced version, the scheme is augmented with a fire-second technique, so that the final velocity can be established with a much higher precision because both the guidance and feedback firing can be performed and ceased at the final maneuver time.

Carson, John↗

The luminosity functions of embedded stellar clusters. 1: Method of solution and analytic results

We describe a method for computing the history of the luminosity function in a young cluster still forming within a molecular cloud complex. Our method, which utilizes detailed results from stellar evolution theory, assumes that clusters arise from the continuous collapse of dense cloud cores over a protracted period of time. It is also assumed that stars reaching the main sequence are distributed in mass according to a prescribed initial mass function (IMF). We keep track separately of the contributions to the luminosity function from the populations of protostars, pre-main-sequence stars, and main-sequence stars. We derive expressions for the fractional contribution of these populations to both the total number of stars produced and the total cluster luminosity. In our model, the number of protostars rises quickly at first, but then levels off to a nearly constant value, which it maintains until the dispersal of the cloud complex. The number fraction of protostars always decreases with time. Averaged over the life of the parent cloud, this fraction is typically a few percent. The protostar mass distribution can be expressed as an integral over the IMF.

Fletcher, Andre B.↗

Measuring and Modeling Shared Visual Attention

Multi-person teams are sometimes responsible for critical tasks, such as flying an airliner. Here we present a method using gaze tracking data to assess shared visual attention, a term we use to describe the situation where team members are attending to a common set of elements in the environment. Gaze data are quantized with respect to a set of N areas of interest (AOIs); these are then used to construct a time series of N dimensional vectors, with each vector component representing one of the AOIs, all set to 0 except for the component corresponding to the currently fixated AOI, which is set to 1. The resulting sequence of vectors can be averaged in time, with the result that each vector component represents the proportion of time that the corresponding AOI was fixated within the given time interval.We present two methods for comparing sequences of this sort, one based on computing the time varying correlation of the averaged vectors, and another based on a chi-square test testing the hypothesis that the observed gaze proportions are drawn from identical probability distributions.We have evaluated the method using synthetic data sets, in which the behavior was modeled as a series of activities, each of which was modeled as a first-order Markov process. By tabulating distributions for pairs of identical and disparate activities, we are able to perform a receiver operating characteristic (ROC) analysis, allowing us to choose appropriate criteria and estimate error rates. Using these criteria, we have applied the methods to data from airline crews, collected in a high-fidelity flight simulator (Gontar Hoermann, 2014). We conclude by considering the problem of automatic (blind) discovery of activities, using methods developed for text analysis.

Mulligan, Jeffrey B.↗

Using Fractional Clock-Period Delays in Telemetry Arraying

A set of special digital all-pass finite-impulse- response (FIR) filters produces phase shifts equivalent to delays that equal fractions of the sampling or clock period of a telemetry-data-processing system. These filters have been used to enhance the arraying of telemetry signals that have been received at multiple ground stations from spacecraft (see figure). Somewhat more specifically, these filters have been used to align, in the time domain, the telemetry-data sequences received by the various antennas, in order to maximize the signal-to-noise ratio of the composite telemetric signal obtained by summing the signals received by the antennas. The term arraying in this context denotes a method of enhanced reception of telemetry signals in which several antennas are used to track a single spacecraft. Each antenna receives a signal that comprises a sum of telemetry data plus noise, and these sum data are sent to an arraying combiner for processing. Correlation is the means used to align the set of data from one antenna with that from another antenna. After the data from all the antennas have been aligned in the time domain, they are all added together, sample by sample.

Fort, David↗

Measuring and Modeling Shared Visual Attention

Multi-person teams are sometimes responsible for critical tasks, such as flying an airliner. Here we present a method using gaze tracking data to assess shared visual attention, a term we use to describe the situation where team members are attending to a common set of elements in the environment. Gaze data are quantized with respect to a set of N areas of interest (AOIs); these are then used to construct a time series of N dimensional vectors, with each vector component representing one of the AOIs, all set to 0 except for the component corresponding to the currently fixated AOI, which is set to 1. The resulting sequence of vectors can be averaged in time, with the result that each vector component represents the proportion of time that the corresponding AOI was fixated within the given time interval. We present two methods for comparing sequences of this sort, one based on computing the time-varying correlation of the averaged vectors, and another based on a chi-square test testing the hypothesis that the observed gaze proportions are drawn from identical probability distributions.We have evaluated the method using synthetic data sets, in which the behavior was modeled as a series of activities, each of which was modeled as a first-order Markov process. By tabulating distributions for pairs of identical and disparate activities, we are able to perform a receiver operating characteristic (ROC) analysis, allowing us to choose appropriate criteria and estimate error rates.We have applied the methods to data from airline crews, collected in a high-fidelity flight simulator (Haslbeck, Gontar Schubert, 2014). We conclude by considering the problem of automatic (blind) discovery of activities, using methods developed for text analysis.

attention↗

Measuring and Modeling Shared Visual Attention

Multi-person teams are sometimes responsible for critical tasks, such as flying an airliner. Here we present a method using gaze tracking data to assess shared visual attention, a term we use to describe the situation where team members are attending to a common set of elements in the environment. Gaze data are quantized with respect to a set of N areas of interest (AOIs); these are then used to construct a time series of N dimensional vectors, with each vector component representing one of the AOIs, all set to 0 except for the component corresponding to the currently fixated AOI, which is set to 1. The resulting sequence of vectors can be averaged in time, with the result that each vector component represents the proportion of time that the corresponding AOI was fixated within the given time interval. We present two methods for comparing sequences of this sort, one based on computing the time-varying correlation of the averaged vectors, and another based on a chi-square test testing the hypothesis that the observed gaze proportions are drawn from identical probability distributions. We have evaluated the method using synthetic data sets, in which the behavior was modeled as a series of "activities," each of which was modeled as a first-order Markov process. By tabulating distributions for pairs of identical and disparate activities, we are able to perform a receiver operating characteristic (ROC) analysis, allowing us to choose appropriate criteria and estimate error rates. We have applied the methods to data from airline crews, collected in a high-fidelity flight simulator (Haslbeck, Gontar & Schubert, 2014). We conclude by considering the problem of automatic (blind) discovery of activities, using methods developed for text analysis.

attention↗

Advanced data assimilation in strongly nonlinear dynamical systems

Advanced data assimilation methods are applied to simple but highly nonlinear problems. The dynamical systems studied here are the stochastically forced double well and the Lorenz model. In both systems, linear approximation of the dynamics about the critical points near which regime transitions occur is not always sufficient to track their occurrence or nonoccurrence. Straightforward application of the extended Kalman filter yields mixed results. The ability of the extended Kalman filter to track transitions of the double-well system from one stable critical point to the other depends on the frequency and accuracy of the observations relative to the mean-square amplitude of the stochastic forcing. The ability of the filter to track the chaotic trajectories of the Lorenz model is limited to short times, as is the ability of strong-constraint variational methods. Examples are given to illustrate the difficulties involved, and qualitative explanations for these difficulties are provided. Three generalizations of the extended Kalman filter are described. The first is based on inspection of the innovation sequence, that is, the successive differences between observations and forecasts; it works very well for the double-well problem. The second, an extension to fourth-order moments, yields excellent results for the Lorenz model but will be unwieldy when applied to models with high-dimensional state spaces. A third, more practical method--based on an empirical statistical model derived from a Monte Carlo simulation--is formulated, and shown to work very well. Weak-constraint methods can be made to perform satisfactorily in the context of these simple models, but such methods do not seem to generalize easily to practical models of the atmosphere and ocean. In particular, it is shown that the equations derived in the weak variational formulation are difficult to solve conveniently for large systems.

Miller, Robert N.↗