Application of optimal linear estimation and control theory to the design of optical tracking systems
Optimal linear estimation and control theory application to high precision optical spacecraft tracking systems design
SEARCH · Search NASA
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.
Optimal linear estimation and control theory application to high precision optical spacecraft tracking systems design
Linear-estimator subalgorithm optimized for local terrain. Algorithm processes some of information in synthetic-aperture-radar echoes from quasi-homogeneous targets, yielding optimal estimates of Doppler frequency shifts. Differs from one described in article, "Estimating SAR Doppler Shifts From Homogeneous Targets", NPO-17869, though offers similar advantages for measurements of wind shears and speeds of ocean currents. Derivation of algorithm involves some of same terminology and simplifying assumptions, however, target area considered quasi-homogeneous in sense it consists of patchwork of subareas, each of which is homogeneous target differing from its neighbors. This assumption more realistic than assumption of homogeneity.
Extensive correlations of computer code results with experimental data are employed to illustrate the use of linearized theory attached flow methods for the estimation and optimization of the aerodynamic performance of simple hinged flap systems. Use of attached flow methods is based on the premise that high levels of aerodynamic efficiency require a flow that is as nearly attached as circumstances permit. A variety of swept wing configurations are considered ranging from fighters to supersonic transports, all with leading- and trailing-edge flaps for enhancement of subsonic aerodynamic efficiency. The results indicate that linearized theory attached flow computer code methods provide a rational basis for the estimation and optimization of flap system aerodynamic performance at subsonic speeds. The analysis also indicates that vortex flap design is not an opposing approach but is closely related to attached flow design concepts. The successful vortex flap design actually suppresses the formation of detached vortices to produce a small vortex which is restricted almost entirely to the leading edge flap itself.
Two nonparametric probability density estimators are considered. The first is the kernel estimator. The problem of choosing the kernel scaling factor based solely on a random sample is addressed. An interactive mode is discussed and an algorithm proposed to choose the scaling factor automatically. The second nonparametric probability estimate uses penalty function techniques with the maximum likelihood criterion. A discrete maximum penalized likelihood estimator is proposed and is shown to be consistent in the mean square error. A numerical implementation technique for the discrete solution is discussed and examples displayed. An extensive simulation study compares the integrated mean square error of the discrete and kernel estimators. The robustness of the discrete estimator is demonstrated graphically.
Two representative sample images of Band 4 of the Landsat Thematic Mapper are compressed with the JPEG algorithm at 8:1, 16:1 and 24:1 Compression Ratios for experimental browsing purposes. We then apply the Optimal PSNR Estimated Spectra Adaptive Postfiltering (ESAP) algorithm to reduce the DCT blocking distortion. ESAP reduces the blocking distortion while preserving most of the image's edge information by adaptively postfiltering the decoded image using the block's spectral information already obtainable from each block's DCT coefficients. The algorithm iteratively applied a one dimensional log-sigmoid weighting function to the separable interpolated local block estimated spectra of the decoded image until it converges to the optimal PSNR with respect to the original using a 2-D steepest ascent search. Convergence is obtained in a few iterations for integer parameters. The optimal logsig parameters are transmitted to the decoder as a negligible byte of overhead data. A unique maxima is guaranteed due to the 2-D asymptotic exponential overshoot shape of the surface generated by the algorithm. ESAP is based on a DFT analysis of the DCT basis functions. It is implemented with pixel-by-pixel spatially adaptive separable FIR postfilters. PSNR objective improvements between 0.4 to 0.8 dB are shown together with their corresponding optimal PSNR adaptive postfiltered images.
Lagrange multipliers in two-point boundary condition problems for trajectory optimization estimated by direct method
Absorbing aerosols play an important, but uncertain, role in the global climate. Much of this uncertainty is due to a lack of adequate aerosol measurements. While great strides have been made in observational capability in the previous years and decades, it has become increasingly apparent that this development must continue. Scanning polarimeters have been designed to help resolve this issue by making accurate, multi-spectral, multi-angle polarized observations. This work involves the use of the Research Scanning Polarimeter (RSP). The RSP was designed as the airborne prototype for the Aerosol Polarimetery Sensor (APS), which was due to be launched as part of the (ultimately failed) NASA Glory mission. Field observations with the RSP, however, have established that simultaneous retrievals of aerosol absorption and vertical distribution over bright land surfaces are quite uncertain. We test a merger of RSP and High Spectral Resolution Lidar (HSRL) data with observations of boreal forest fire smoke, collected during the Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS). During ARCTAS, the RSP and HSRL instruments were mounted on the same aircraft, and validation data were provided by instruments on an aircraft flying a coordinated flight pattern. We found that the lidar data did indeed improve aerosol retrievals using an optimal estimation method, although not primarily because of the constraints imposed on the aerosol vertical distribution. The more useful piece of information from the HSRL was the total column aerosol optical depth, which was used to select the initial value (optimization starting point) of the aerosol number concentration. When ground based sun photometer network climatologies of number concentration were used as an initial value, we found that roughly half of the retrievals had unrealistic sizes and imaginary indices, even though the retrieved spectral optical depths agreed within uncertainties to independent observations. The convergence to an unrealistic local minimum by the optimal estimator is related to the relatively low sensitivity to particles smaller than 0.1 ( m) at large optical thicknesses. Thus, optimization algorithms used for operational aerosol retrievals of the fine mode size distribution, when the total optical depth is large, will require initial values generated from table look-ups that exclude unrealistic size/complex index mixtures. External constraints from lidar on initial values used in the optimal estimation methods will also be valuable in reducing the likelihood of obtaining spurious retrievals.
We describe a variation of the Optimal Estimation (OE) method for greenhouse gas remote sensing retrievals using a singular value decomposition (SVD) and an uninformative prior. The SVD method is capable of discerning vertical information in column integrated absorption measurements. While traditional Bayesian optimal estimation (OE) assumes a prior distribution in order to regularize the inversion problem, the SVD approach identifies principal components that can be retrieved from the measurement without explicitly specifying a prior mean and prior covariance matrix. We discuss the method, illustrate its use on an integrated path differential absorption CO2 lidar measurement model, and compare it to traditional optimal estimation using numerical simulations. In the absence of forward model error, the mathematics behind the SVD method guarantee it to be bias-free, which is confirmed by the numerical simulations. In contrast, traditional OE retrievals exhibit bias when the prior mean used in the retrieval differs from the true mean. While the SVD approach can be used for most trace gas retrievals, it is particularly useful for situations where prior knowledge of the trace gas profile is poor. The SVD analysis is also useful in illustrating how vertical information is treated by the traditional OE approach.
An Optimal Estimation (OE) inversion method has been developed to retrieve vertical profiles of aerosol size distributions, volume/surface concentrations, complex refractive indices, and single scattering albedos using the High Spectral Resolution Lidar (HSRL) and polarimeter measurements. By combining the active and passive remote sensing data, we can take advantages of the high vertical resolution of the lidar measurements and the high information content of the combined observations. There are three components to the new retrieval system: 1) simultaneous aerosol vertical profile retrieval using HSRL data, 2) standalone aerosol retrieval using Research Scanning Polarimeter (RSP), and 3) combined HSRL and RSP optimal estimation retrieval. For the LIDAR retrieval, we solve for number concentrations, mode fractions, aerosol size distributions, and real and imaginary part of refractive indices simultaneously. Based on this information, other parameters such as effective particle radius, surface and volume concentration, single scattering albedo can be derived. The combined HSRL-polarimeter algorithm will provide improved single scattering albedo and refractive index retrievals. The newly developed retrieval algorithm is very fast, and it provides an optimal estimation solution with associated error estimates. By solving for the aerosol vertical profiles using HSRL extinctions and backscattering coefficients from different layers, the new OE method can effectively reduce the impact of measurement noise on the retrieved aerosol microphysical properties. It also enables us to effectively combing HSRL and RSP data in the single OE retrieval system since the passive polarimeter instruments measure path-integrated intensity and polarizations. We have performed sensitivity studies using simulated data and applied the OE algorithms to NASA airborne field data measured by the HSRL-2 and the RSP instruments.
This paper presents a sequential filtering strategy using observations corrupted with uniform measurement noise. While the Kalman filter remains the best linear estimator of the state, other filtering techniques provide minimum variance optimal estimates, a trait only enjoyed by the Kalman filter when the underlying noises are, in fact, Gaussian. This work develops a new approximate optimal estimator for uniform measurement noises. The resulting recursion requires just slightly more computational time to complete a measurement update than the Kalman filter, which generally cannot be claimed by other optimal strategies such as the particle or Gaussian mixture filters.
This paper presents a sequential filtering strategy using observations corrupted with uniform measurement noise. While the Kalman filter remains the best linear estimator of the state, other filtering techniques provide minimum variance optimal estimates, a trait only enjoyed by the Kalman filter when the underlying noises are, in fact, Gaussian. This work develops a new approximate optimal estimator for uniform measurement noises. The resulting recursion requires just slightly more computational time to complete a measurement update than the Kalman filter, which generally cannot be claimed by other optimal strategies such as the particle or Gaussian mixture filters.
This study presents and evaluates an updated algorithmfor quantification of absorbing aerosols above clouds(AACs) from passive satellite measurements. The focus isbiomass burning in the south-eastern Atlantic Ocean duringthe 2016 and 2017 ObseRvations of Aerosols above CLoudsand their intEractionS (ORACLES) field campaign deployments.The algorithm retrieves the above-cloud aerosoloptical depth (AOD) and underlying liquid cloud opticaldepth and is applied to measurements from the Sea-viewingWide Field-of-view Sensor (SeaWiFS), Moderate ResolutionImaging Spectroradiometer (MODIS), and Visible InfraredImaging Radiometer Suite (VIIRS) from 1997 to 2017. AirborneNASA Ames Spectrometers for Sky-Scanning, Sun-Tracking Atmospheric Research (4STAR) and NASA LangleyHigh Spectral Resolution Lidar 2 (HSRL2) data collectedduring ORACLES provide important validation for spectralAOD for MODIS and VIIRS; as the SeaWiFS missionended in 2010, it cannot be evaluated directly. The 4STARand HSRL2 comparisons are complementary and reveal performancegenerally in line with uncertainty estimates providedby the optimal estimation retrieval framework used. Atpresent the two MODIS-based data records seem the mostreliable, although there are differences between the deployments,which may indicate that the available data are not yetsufficient to provide a robust regional validation. Spatiotemporalpatterns in the data sets are similar, and the time seriesare very strongly correlated with each other (correlationcoefficients from 0.95 to 0.99). Offsets between the satellitedata sets are thought to be chiefly due to differences in absolutecalibration between the sensors. The available validationdata for this type of algorithm are limited to a small numberof field campaigns, and it is strongly recommended that suchairborne measurements continue to be made, both over thesouthern Atlantic Ocean and elsewhere.
Control of air contaminants is a crucial factor in the safety considerations of crewed space flight. Indoor air quality needs to be closely monitored during long range missions such as a Mars mission, and also on large complex space structures such as the International Space Station. This work mainly pertains to the detection and simulation of air contaminants in the space station, though much of the work is easily extended to buildings, and issues of ventilation systems. Here we propose a method with which to track the presence of contaminants using an accurate physical model, and also develop a robust procedure that would raise alarms when certain tolerance levels are exceeded. A part of this research concerns the modeling of air flow inside a spacecraft, and the consequent dispersal pattern of contaminants. Our objective is to also monitor the contaminants on-line, so we develop a state estimation procedure that makes use of the measurements from a sensor system and determines an optimal estimate of the contamination in the system as a function of time and space. The real-time optimal estimates in turn are used to detect faults in the system and also offer diagnoses as to their sources. This work is concerned with the monitoring of air contaminants aboard future generation spacecraft and seeks to satisfy NASA's requirements as outlined in their Strategic Plan document (Technology Development Requirements, 1996).
The sampling error of an arbitrary linear estimate of a time-averaged quantity constructed from a time series of irregularly spaced observations at a fixed located is quantified through a formalism. The method is applied to satellite observations of chlorophyll from the coastal zone color scanner. The two specific linear estimates under consideration are the composite average formed from the simple average of all observations within the averaging period and the optimal estimate formed by minimizing the mean squared error of the temporal average based on all the observations in the time series. The resulting suboptimal estimates are shown to be more accurate than composite averages. Suboptimal estimates are also found to be nearly as accurate as optimal estimates using the correct signal and measurement error variances and correlation functions for realistic ranges of these parameters, which makes it a viable practical alternative to the composite average method generally employed at present.
Given a valid system model and adequate observability, a Kalman filter will converge toward the true system state with error statistics given by the estimated error covariance matrix. The errors generally do not continue to decrease. Rather, a balance is reached between the gain of information from new measurements and the loss of information during propagation. The errors can be further reduced, however, by a second pass through the data with an optimal smoother. This algorithm obtains the optimally weighted average of forward and backward propagating Kalman filters. It roughly halves the error covariance by including future as well as past measurements in each estimate. This paper investigates whether such benefits actually accrue in the application of an optimal smoother to spacecraft attitude determination. Tests are performed both with actual spacecraft data from the Extreme Ultraviolet Explorer (EUVE) and with simulated data for which the true state vector and noise statistics are exactly known.
The most recent results in an ongoing research effort at Princeton in the area of flight dynamics in wind shear are described. The first undertaking in this project was a trajectory optimization study. The flight path of a medium-haul twin-jet transport aircraft was optimized during microburst encounters on final approach. The assumed goal was to track a reference climb rate during an aborted landing, subject to a minimum airspeed constraint. The results demonstrated that the energy loss through the microburst significantly affected the qualitative nature of the optimal flight path. In microbursts of light to moderate strength, the aircraft was able to track the reference climb rate successfully. In severe microbursts, the minimum airspeed constraint in the optimization forced the aircraft to settle on a climb rate smaller than the target. A tradeoff was forced between the objectives of flight path tracking and stall prevention.
The estimation and control of large flexible space structures pose significant new control technology problems. One major problem area is the closed-loop stability of a structure that has been modeled with truncated modal dynamics. Stability problems arise due to the control design process beginning with a deficient model. This work provides the necessary conditions for the optimal estimation of infinite dimensional (partial-differential equation) systems. This approach can be particularly useful for initial control studies and for gaining considerable insight into what the optimal estimator for a truly infinite dimensional system should be. A detailed example of the estimation of the continuous shape of a string in tension is presented.
A meaningful interpretation of airborne flux estimates over nonuniform terrain must consider local advection from nonhomogeneous source distributions (footprint correction). An empirical procedure for footprint correction of airborne flux data obtained in the First ISLSCP Field Experiment (FIFE) 1989 is presented, based on simple analytical solutions of the diffusion equation which have been adjusted to approximate more realistic solutions. Optimized estimates of surface flux 'maps' derived from airborne observations are then constructed for maximum correlation between flux estimates and independently observed surface characteristics. This paper also includes an analysis of resolution for the given procedure, in terms of amplitude recovery of hypothetical square-wave distributions of surface flux density, and discusses its implications for the interpretation of the optimized estimates of FIFE surface flux maps. Results show that corrected spatial sink distributions are obtained by the current empirical procedure at the cost of considerable reduction in amplitude recovery of small-scale variations in surface flux density. The high spatial correlation observed between corrected flux maps and surface characteristics such as greenness or surface temperature excess is ascribed to the approximately correct positioning of the main flux gradients across the site, which can be expected to have been reproduced with amplitude recovery equal to or greater than 60 percent.