Search NASA⌕ Search

SEARCH · Search NASA

Results for “Filtering algorithms”

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 1,117 records · Page 62

Total Variation Majorization Minimization (TV-MM) Approach to Radiometer Brightness Temperature Gridding and Reconstruction

This paper presents the implementation of an algorithm to enhance the image resolution of the Earth's surface brightness temperature (T B ) data measured by radiometers such as the one onboard of the Soil Moisture Active Passive (SMAP) mission. A key step in radiometer T B processing is the conversion of the swath-based calibrated antenna temperature (T A ) measurements to the Level 3 Earth-centered grid. The simplest algorithm to transform this data from swath to gridded format is called drop-in-the-bucket which simply averages surrounding noisy T A samples to form a T B value at the gridded location. This method reduces noise, however produces low resolution products. To obtain a higher resolution product, SMAP uses other techniques such the Backus-Gilbert (BG) algorithm, which is the conventional method used in microwave radiometry. Although this method performs the required interpolation, it is not effective in denoising and removing blurring effects due to antenna filtering of the radiometer image data. Our motivation for this development is to further improve the resolution through post-processing of the radiometer T B image, a highly cost-effective method of image enhancement. The approach adapted in this work is based on the minimization of the Total Variation (TV) regularized objective function that is used extensively in solving general ill-posed linear inverse problems in image processing. Since the TV-based objective function is convex but not everywhere differentiable, there exists many numerical algorithms that can estimate the solution and the one selected for this work is called Majorization- Minimization (MM). By applying this algorithm, simulation experiments were performed based on synthetic data from the Geophysical model as well as real SMAP data to demonstrate the effectiveness of the technique. Results were then compared against the BG method.

Wing Lee↗

Real-time optical correlator using computer-generated holographic filter on a liquid crystal light valve

Limitations associated with the binary phase-only filter often used in optical correlators are presently circumvented in the writing of complex-valued data on a gray-scale spatial light modulator through the use of a computer-generated hologram (CGH) algorithm. The CGH encodes complex-valued data into nonnegative real CGH data in such a way that it may be encoded in any of the available gray-scale spatial light modulators. A CdS liquid-crystal light valve is used for the complex-valued CGH encoding; computer simulations and experimental results are compared, and the use of such a CGH filter as the synapse hologram in a holographic optical neural net is discussed.

Chao, Tien-Hsin↗

Exploring the Use of Alfven Waves in Magnetometer Calibration at Geosynchronous Orbit

An Alfven wave is a type magnetohydrodynamicwave that travels through a conducting fluid under the influence of a magnetic field. Researchers have successfully calculated offset vectors of magnetometers in interplanetary space by optimizing the offset to maximize certain Alfvenic properties of observed waves (Leinweber, Belcher). If suitable Alfven waves can be found in the magnetosphere at geosynchronous altitude then these techniques could be used to augment the overall calibration plan for magnetometers in this region such as on the GOES spacecraft, possibly increasing the time between regular maneuvers. Calibration maneuvers may be undesirable because they disrupt the activities of other instruments. Various algorithms to calculate an offset using Alfven waves were considered. A new variation of the Davis-Smith method was derived because it can be mathematically shown that the Davis-Smith method tolerates filtered data, which expands potential applications. The variant developed was designed to find only the offset in the plane normal to the main field because the overall direction of Earth's magnetic field rarely changes, and theory suggests the Alfvenic disturbances occur transverse to the main field. Other variations of the Davis-Smith method encounter problems with data containing waves that propagate in mostly the same direction. A searching algorithm was then designed to look for periods of time with potential Alfven waves in GOES 15 data based on parameters requiring that disturbances be normal to the main field and not change field magnitude. Final waves for calculation were hand-selected. These waves produced credible two-dimensional offset vectors when input to the Davis-Smith method. Multiple two-dimensional solutions in different planes can be combined to get a measurement of the complete offset. The resulting three dimensional offset did not show sufficient precision over several years to be used as a primary calibration method, but reflected changes in the offset fairly well, suggesting that the method could be helpful in monitoring trends of the offset vector when maneuvers cannot be used.

calibration↗

Kalman Filtering USNO's GPS Observations for Improved Time Transfer Predictions

The Global Positioning System (GPS) Master Control Station (MCS) performs the Coordinated Universal Time (UTC) time transfer mission by uploading and broadcasting predictions of the GPS-UTC offset in subframe 4 of the GS navigation message. These predictions are based on only two successive daily data points obtained from the US Naval Observatory (USNO). USNO produces these daily smoothed data points by performing a least-squares fit on roughly 38 hours worth of data from roughly 160 successive 13-minute tracks of GPS satellites. Though sufficient for helping to maintain a time transfer error specification of 28 ns (1 Sigma), the MCS's prediction algorithm does not make the best use of the available data from from USNO, and produces data that can degrade quickly over extended prediction spans. This paper investigates how, by applying Kalman filtering to the same available tracking data, the MCS could improve its estimate of GPS-UTC, and in particular, the GPS-UTC A(sub 1) term. By refining the A(sub 1) (frequency) estimate for GPS-UTC predictions, error in GPS time transfer could drop significantly. Additional, the risk of future spikes in GPS's time transfer error could similarly be minimized, by employing robust Kalman filtering for GPS-UTC predictions.

Hutsell, Steven T.↗

Hybrid eigensolvers for nuclear configuration interaction calculations

We examine and compare several iterative methods for solving large-scale eigenvalue problems arising from nuclear structure calculations. In particular, we discuss the possibility of using block Lanczos method, a Chebyshev filtering based subspace iterations and the residual minimization method accelerated by direct inversion of iterative subspace (RMM-DIIS) and describe how these algorithms compare with the standard Lanczos algorithm and the locally optimal block preconditioned conjugate gradient (LOBPCG) algorithm. Although the RMM-DIIS method does not exhibit rapid convergence when the initial approximations to the desired eigenvectors are not sufficiently accurate, it can be effectively combined with either the block Lanczos or the LOBPCG method to yield a hybrid eigensolver that has several desirable properties. We will describe a few practical issues that need to be addressed to make the hybrid solver efficient and robust.

97 MATHEMATICS AND COMPUTING↗

The Potential Benefits of Handling Mixture Statistics via a Bi-Gaussian EnKF: Tests With All-Sky Satellite Infrared Radiances

The meteorological characteristics of cloudy atmospheric columns can be very different from their clear counterparts. Thus, when a forecast ensemble is uncertain about the presence/absence of clouds at a specific atmospheric column (i.e., some members are clear while others are cloudy), that column's ensemble statistics will contain a mixture of clear and cloudy statistics. Such mixtures are inconsistent with the ensemble data assimilation algorithms currently used in numerical weather prediction. Hence, ensemble data assimilation algorithms that can handle such mixtures can potentially outperform currently used algorithms. In this study, we demonstrate the potential benefits of addressing such mixtures through a bi-Gaussian extension of the ensemble Kalman filter (BGEnKF). The BGEnKF is compared against the commonly used ensemble Kalman filter (EnKF) using perfect model observing system simulated experiments (OSSEs) with a realistic weather model (the Weather Research and Forecast model). Synthetic all-sky infrared radiance observations are assimilated in this study. In these OSSEs, the BGEnKF outperforms the EnKF in terms of the horizontal wind components, temperature, specific humidity, and simulated upper tropospheric water vapor channel infrared brightness temperatures. This study is one of the first to demonstrate the potential of a Gaussian mixture model EnKF with a realistic weather model. Our results thus motivate future research toward improving numerical Earth system predictions though explicitly handling mixture statistics.

54 ENVIRONMENTAL SCIENCES↗

An Iterative Approach for Solving the SCOPF Problem Applying LP, SOCP, and NLP Subproblems

We propose to develop efficient algorithms and software for the SCOPF problem. We will employ an iterative approach that will: a) use linear subproblems and other active set filtering techniques to identify the most important contingencies and drastically reduce the SCOPF model size; b) solve SOCP relaxations of the reduced SCOPF to converge to the neighborhood of the global optimal solution and establish a lower bound on the solution, and; c) use a non-convex, nonlinear interior-point solver, Artelys Knitro, to converge quickly to the optimal solution. To identify the most effective approach, we will experiment with several techniques to identify the tradeoffs between contingency subproblem complexity and fast solvability.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Microprogrammed digital filters for strapdown guidance application.

Discussion of an approach for implementing digital filters using microprogrammed control logic with read only memory (ROM) for strapdown guidance applications, and description of a second-order multiplexed system using multiprogrammed control instruction. The microprogramming technique for control using ROM is shown to enhance higher order digital filter realization. The high speed of digital circuits reduces propagation time and facilitates multiplexing. Programmed ROM's can be altered easily for different algorithms.

Kapadia, K.↗

A parameter-estimation subroutine package

Estimation subroutine package comprises fast, efficient, and simple least-squares data-processing algorithms for use in orbit determination and related analyses. Very reliable and general algorithms have been documented. Package contains collection of streamlined subroutines that can be used to solve large variety of parameter-estimation and filtering problems. Special routines are included for problems with colored process noise and covariance (factor) mapping.

Bierman, G.↗

Unsupervised Classification of Global Radar Units on Venus

Characterization of the Venusian surface in terms of its radar properties was accomplished by application of an unsupervised, linear discriminant algorithm to two Pioneer-Venus (PV) Orbiter radar data sets: the RMS-slope (surface roughness) and reflectivity. Both databases were spatially filtered to the same effective resolution of 100 km prior to classification. A recent supervised classification study using these data was based on presupposed morphologic significance of selected data ranges. The knowledge of both Venusian geology and the geologic significance of the radar data is so limited that the data warrant a more unsupervised approach; for this study a linear discriminant classifier was chosen. This approach is purely statistical, thereby removing any observer bias. Statistical significance of the resulting clusters was evaluated by an ancillary program in which an F test utilizing the Mahalanobis' distance.

Kozak, R. C.↗

A decentralized square root information filter/smoother

A number of developments has recently led to a considerable interest in the decentralization of linear least squares estimators. The developments are partly related to the impending emergence of VLSI technology, the realization of parallel processing, and the need for algorithmic ways to speed the solution of dynamically decoupled, high dimensional estimation problems. A new method is presented for combining Square Root Information Filters (SRIF) estimates obtained from independent data sets. The new method involves an orthogonal transformation, and an information matrix filter 'homework' problem discussed by Schweppe (1973) is generalized. The employed SRIF orthogonal transformation methodology has been described by Bierman (1977).

Bierman, G. J.↗

Effect of model uncertainty on failure detection - The threshold selector

The performance of all failure detection, isolation, and accomodation (DIA) algorithms is influenced by the presence of model uncertainty. A unique framework is presented to incorporate a knowledge of modeling error in the analysis and design of failure detection systems. The tools being used are very similar to those in robust control theory. A concept is introduced called the threshold selector, which is a nonlinear inequality whose solution defines the set of detectable sensor failure signals. The threshold selector represents an innovative tool for analysis and synthesis of DIA algorithms. It identifies the optimal threshold to be used in innovations-based DIA algorithms. The optimal threshold is shown to be a function of the bound on modeling errors, the noise properties, the speed of DIA filters, and the classes of reference and failure signals. The size of the smallest detectable failure is also determined. The results are applied to a multivariable turbofan jet engine example, which demonstrates improvements compared to previous studies.

Emami-Naeini, Abbas↗

Analysis and improved design considerations for airborne pulse Doppler radar signal processing in the detection of hazardous windshear

High resolution windspeed profile measurements are needed to provide reliable detection of hazardous low altitude windshear with an airborne pulse Doppler radar. The system phase noise in a Doppler weather radar may degrade the spectrum moment estimation quality and the clutter cancellation capability which are important in windshear detection. Also the bias due to weather return Doppler spectrum skewness may cause large errors in pulse pair spectral parameter estimates. These effects are analyzed for the improvement of an airborne Doppler weather radar signal processing design. A method is presented for the direct measurement of windspeed gradient using low pulse repetition frequency (PRF) radar. This spatial gradient is essential in obtaining the windshear hazard index. As an alternative, the modified Prony method is suggested as a spectrum mode estimator for both the clutter and weather signal. Estimation of Doppler spectrum modes may provide the desired windshear hazard information without the need of any preliminary processing requirement such as clutter filtering. The results obtained by processing a NASA simulation model output support consideration of mode identification as one component of a windshear detection algorithm.

Lee, Jonggil↗

Filtering Image Records Using Wavelets and the Zakai Equation

Consider the problem of detecting and localizing a faint object moving In an "essentially stationary" background, using a sequence of two-dimensional low-SNR images of the scene. A natural approach consists of "digitizing" each snapshot into a discrete set of observations, sufficiently (perhaps not exactly) matched to the object In question, then tracking the object using an appropriate stochastic filter. The tracking would be expected to make up for the low signal-to-noise ratio, this allowing one to "coherently" process successive images in order to beat down the noise and localize the object. Thus, "tracking" here does not refer to the usual notion of detecting then tracking: rather, we track in order to detect The problem then becomes one of choosing the appropriate image representation as well as the optimal (and necessarily nonlinear filter. We propose exact and approximate solutions using wavelets and the Zakai equation. The smoothness of the wavelets used is required in the derivation of the evolution equation for the conditional density giving the filter, and their orthogonality makes it possible to carry out actual computations of the Ito- and change-of-gauge-terms in the algorithm effectively.

Haddad, Ziad S.↗

Design of an L-band Microwave Radiometer with Active Mitigation of Interference

Radio frequency interference (RFI) impairs L-band radiometry outside the protected 20 MHz frequency band around 1413 MHz. However, bandwidths of 100 MHz or more are desired for certain remote sensing applications as well as certain astronomy applications. Because much of the RFI in this band is from radars with pulse lengths on the order of microseconds, traditional radiometers (i.e., those which directly measure total power or power spectral density integrated over time scales of milliseconds or greater) are poorly-suited to this task. Simply reducing integration time and discarding contaminated outputs may not be a practical answer due to the wide variety of modulations and pulse lengths observed in L-band RFI signals, the dynamic and complex nature of the associated propagation channels, and the logistical effort associated with post-measurement data editing. This motivates the design and development of radiometers capable of coherent sampling and adaptive, real-time mitigation of interference. Such a radiometer will be described in this presentation. This design is capable of coherently-sampling up to 100 MHz bandwidth at L-band. RFI mitigation is implemented in FPGA components so that real-time suppression is achieved. The system currently uses a cascade of basic time- and frequency- domain detection and blanking techniques; more advanced algorithms are un- der consideration. The modular FPGA-based architecture provides other benefits, such as the ability to implement extremely stable digital filters and the ability to reconfigure the system "on the fly". An overview of the basic design along with on-the-air results from an initial implementation will be provided in the presentation. Related L-band RFI surveys will be described to illustrate the relevance of this approach in a variety of operating conditions.

Ellingson, S. W.↗

Next-Generation Aura/OMI NO2 and SO2 Products

The measurement of both SO2 and NO2 gases are recognized as an essential component of atmospheric composition missions. We describe current capabilities and limitations of the operational Aura/OMI NO2 and SO2 data that have been used by a large number of researchers. Analyses of the data and validation studies have brought to light a number of areas in which these products can be expanded and improved. Major improvements for new NASA standard (SP) NO2 product include more accurate tropospheric and stratospheric column amounts, along with much improved error estimates and diagnostics. Our approach uses a monthly NO2 climatology based on the NASA Global Modeling Initiative (GMI) chemistry-transport model and takes advantage of OMI data from cloudy scenes to find clean areas where the contribution from the trap NO2 column is relatively small. We then use a new filtering, interpolation and smoothing techniques for separating the stratospheric and tropospheric components of NO2, minimizing the influence of a priori information. The new algorithm greatly improves the structure of stratospheric features relative to the original SP. For the next-generation OMI SO2 product we plan to implement operationally the offline iterative spectral fitting (ISF) algorithm and re-process the OMI Level-2 SO2 dataset using a priori SO2 and aerosol profiles, clouds, and surface reflectivity appropriate for observation conditions. This will improve the ability to detect and quantify weak tropospheric SO2 loadings. The new algorithm is validated using aircraft in-situ data during field campaigns in China (2005 and 2008) and in Maryland (Frostburg, 2010 and DISCOVER-AQ in July 2011). The height of the SO2 plumes will also be estimated for high SO2 loading cases (e.g., volcanic eruptions). The same SO2 algorithm will be applied to the data from OMPS sensor to be launched on NPP satellite later this year. The next-generation NO2 and SO2 products will provide critical information (e.g., averaging kernels) for evaluation of chemistry-transport models, for data assimilation, and to impose top-down constraints on the SO2 and NO2 emission sources.

Krotkov, Nickolay↗

Robust, Practical Adaptive Control for Launch Vehicles

A modern mechanization of a classical adaptive control concept is presented with an application to launch vehicle attitude control systems. Due to a rigorous flight certification environment, many adaptive control concepts are infeasible when applied to high-risk aerospace systems; methods of stability analysis are either intractable for high complexity models or cannot be reconciled in light of classical requirements. Furthermore, many adaptive techniques appearing in the literature are not suitable for application to conditionally stable systems with complex flexible-body dynamics, as is often the case with launch vehicles. The present technique is a multiplicative forward loop gain adaptive law similar to that used for the NASA X-15 flight research vehicle. In digital implementation with several novel features, it is well-suited to application on aerodynamically unstable launch vehicles with thrust vector control via augmentation of the baseline attitude/attitude-rate feedback control scheme. The approach is compatible with standard design features of autopilots for launch vehicles, including phase stabilization of lateral bending and slosh via linear filters. In addition, the method of assessing flight control stability via classical gain and phase margins is not affected under reasonable assumptions. The algorithm s ability to recover from certain unstable operating regimes can in fact be understood in terms of frequency-domain criteria. Finally, simulation results are presented that confirm the ability of the algorithm to improve performance and robustness in realistic failure scenarios.

Orr, Jeb. S.↗

A Battery Health Monitoring Framework for Planetary Rovers

Batteries have seen an increased use in electric ground and air vehicles for commercial, military, and space applications as the primary energy source. An important aspect of using batteries in such contexts is battery health monitoring. Batteries must be carefully monitored such that the battery health can be determined, and end of discharge and end of usable life events may be accurately predicted. For planetary rovers, battery health estimation and prediction is critical to mission planning and decision-making. We develop a model-based approach utilizing computaitonally efficient and accurate electrochemistry models of batteries. An unscented Kalman filter yields state estimates, which are then used to predict the future behavior of the batteries and, specifically, end of discharge. The prediction algorithm accounts for possible future power demands on the rover batteries in order to provide meaningful results and an accurate representation of prediction uncertainty. The framework is demonstrated on a set of lithium-ion batteries powering a rover at NASA.

prediction↗