SEARCH · Search NASA

Results for “Maximum Likelihood Estimation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 541 records · Page 30

All sources NASA NTRS DOE OSTI

DOE OSTI2016

SU-C-207A-01: A Novel Maximum Likelihood Method for High-Resolution Proton Radiography/proton CT

Purpose: Multiple Coulomb scattering is the largest contributor to blurring in proton imaging. Here we tested a maximum likelihood least squares estimator (MLLSE) to improve the spatial resolution of proton radiography (pRad) and proton computed tomography (pCT). Methods: The object is discretized into voxels and the average relative stopping power through voxel columns defined from the source to the detector pixels is optimized such that it maximizes the likelihood of the proton energy loss. The length spent by individual protons in each column is calculated through an optimized cubic spline estimate. pRad images were first produced using Geant4 simulations. An anthropomorphic head phantom and the Catphan line-pair module for 3-D spatial resolution were studied and resulting images were analyzed. Both parallel and conical beam have been investigated for simulated pRad acquisition. Then, experimental data of a pediatric head phantom (CIRS) were acquired using a recently completed experimental pCT scanner. Specific filters were applied on proton angle and energy loss data to remove proton histories that underwent nuclear interactions. The MTF10% (lp/mm) was used to evaluate and compare spatial resolution. Results: Numerical simulations showed improvement in the pRad spatial resolution for the parallel (2.75 to 6.71 lp/cm) and conical beam (3.08 to 5.83 lp/cm) reconstructed with the MLLSE compared to averaging detector pixel signals. For full tomographic reconstruction, the improved pRad were used as input into a simultaneous algebraic reconstruction algorithm. The Catphan pCT reconstruction based on the MLLSE-enhanced projection showed spatial resolution improvement for the parallel (2.83 to 5.86 lp/cm) and conical beam (3.03 to 5.15 lp/cm). The anthropomorphic head pCT displayed important contrast gains in high-gradient regions. Experimental results also demonstrated significant improvement in spatial resolution of the pediatric head radiography. Conclusion: The proposed MLLSE shows promising potential to increase the spatial resolution (up to 244%) in proton imaging.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2017

The Formation and Physical Origin of Highly Ionized Cooling Gas

We present a simple model that explains the origin of warm, diffuse gas seen primarily as highly ionized absorption-line systems in the spectra of background sources. We predict the observed column densities of several highly ionized transitions such as O vi, O vii, Ne viii, N v, and Mg x, and we present a unified comparison of the model predictions with absorption lines seen in the Milky Way disk, Milky Way halo, starburst galaxies, the circumgalactic medium, and the intergalactic medium at low and high redshifts. We show that diffuse gas seen in such diverse environments can be simultaneously explained by a simple model of radiatively cooling gas. We show that most such absorption-line systems are consistent with being collisionally ionized, and we estimate the maximum-likelihood temperature of the gas in each observation. This model satisfactorily explains why O vi is regularly observed around star-forming low- z L* galaxies, and why N v is rarely seen around the same galaxies. We further present some consequences of this model in quantifying the dynamics of the cooling gas around galaxies and predict the shock velocities associated with such flows. A unique strength of this model is that while it has only one free (but physically well-constrained) parameter, it nevertheless successfully reproduces the available data on O vi absorbers in the interstellar, circumgalactic, intragroup, and intergalactic media, as well as the available data on other absorption lines from highly ionized species.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1971

The analysis of incomplete data.

In this paper, we attempt to provide a simple taxonomy for incomplete-data problems and at the same time develop unified methods of analysis. The emphasis is on techniques which are natural extensions of the complete-data analysis and which will handle rather general classes of incomplete-data problems as opposed to custom-made techniques for special problems. The principle of estimation is either maximum likelihood or is at least based on maximum likelihood.

Hartley, H. O.↗

NASA NTRS1979

Large Area Crop Inventory Experiment (LACIE). Development of procedure M for multicrop inventory, with tests of a spring-wheat configuration

The author has identified the following significant results. An outgrowth of research and development activities in support of LACIE was a multicrop area estimation procedure, Procedure M. This procedure was a flexible, modular system that could be operated within the LACIE framework. Its distinctive features were refined preprocessing (including spatially varying correction for atmospheric haze), definition of field like spatial features for labeling, spectral stratification, and unbiased selection of samples to label and crop area estimation without conventional maximum likelihood classification.

Horvath, R.↗

NASA NTRS1979

Large Area Crop Inventory Experiment (LACIE). Evaluation of three-category classification

The author has identified the following signficant results. Examination of both machine estimates and stratified areal estimates produced by clustering and classification reveal no significant differences between the proportion estimates and ground truth estimates. When testing the variances of the machine estimates, a significant reduction in the variances was found when the number of starting dots was increased from 30 to 45. The variances were again reduced, although not significantly, when the number of starting dots was increased from 45 to 60. From these results, 60 starting dots are recommended for a three-category classifier. When examining the variances of the estimates for the four estimation procedures (using 60 dots), no significant differences were found between procedures. Thus, only the machine clustering may be used to produce an estimate, and the stratified areal estimate computations and maximum likelihood classification can be deleted.

Havens, K. A.↗

NASA NTRS1980

Formulation and implementation of a practical algorithm for parameter estimation with process and measurement noise

A new formulation is proposed for the problem of parameter estimation of dynamic systems with both process and measurement noise. The formulation gives estimates that are maximum likelihood asymptotically in time. The means used to overcome the difficulties encountered by previous formulations are discussed. It is then shown how the proposed formulation can be efficiently implemented in a computer program. A computer program using the proposed formulation is available in a form suitable for routine application. Examples with simulated and real data are given to illustrate that the program works well.

Maine, R. E.↗

NASA NTRS1983

A Review of System Identification Methods Applied to Aircraft

Airplane identification, equation error method, maximum likelihood method, parameter estimation in frequency domain, extended Kalman filter, aircraft equations of motion, aerodynamic model equations, criteria for the selection of a parsimonious model, and online aircraft identification are addressed.

Klein, V.↗

NASA NTRS1987

Dynamic data acquisition, reduction, and analysis for the identification of high-speed compressor component post-stability characteristics

A compressor test was conducted in which transient data were obtained for the purpose of identifying the high-speed post-stability characteristics. The transient, surge-cycle nature of high-speed post-stability operation precludes the possibility of obtaining the characteristics in a steady-state manner, as is possible during low-speed poststability operation, which is characterized by quasi-steady rotating-stall behavior. Specialized compressor instrumentation was developed and was used to obtain the necessary surge-cycle performance data, which were then digitized, filtered, and analyzed. The high-speed post-stability characteristics were obtained through the use of a maximum likelihood-parameter estimation technique. The estimated characteristics were found to be insensitive to the presence of measurement noise and unmodelled system dynamics, but the compressor time-response constants, which were also estimated, were more sensitive to these same disturbances.

Dvorak, S. D.↗

NASA NTRS1998

Guest Investigator Studies with the Compton Gamma Ray Observatory

The cumulative all-sky survey by the Compton Gamma Ray Observatory (CGRO), composed of data acquired during the first three years of the mission, included a number of regions of very limited exposure. The most glaring deficiency in coverage was toward the region of the South Galactic Pole (SGP), which received significantly less exposure than other directions- by a factor of at least 2 to 3. Furthermore, nearly all of the SGP exposure was accumulated during the first year of the mission. Since blazars are known to be time-variable, and of unknown duty cycle, a pointing of the CCRO in that direction was considered highly desirable, and long overdue. In addition, data gathered from a pointing toward the SGP and its comparison with comprehensive data available for the North Galactic Pole would be extremely valuable to investigators studying the extragalactic diffuse emission. The reasons outlined above prompted our initiation of a Cycle 4 campaign to systematically search with EGRET and COMPTEL for gamma-ray emission from sources near the South Galactic Pole. The Cycle 4 SGP campaign consisted of tnvo 14-day observations separated in in time by approximately 10 months. The temporal separation of the observations was requested to allow a test for possible variations in the detected sources. Our primary targets were 38 FSRQs which lie within 30 degrees of the SGP, and which satisfy the basic criteria for candidate gamma-ray AGNs,flat-spectrum radio sources, many of which exhibit blazar-type properties). These targets were selected from the standard references, and from the available on-line databases (e.g., the NASA Extragalactic Database, NED), as the most promising AGN targets in the vicinity of the SGP. A 30 radius from the SGP was chosen as the boundary of our survey, since the selected targets would then fall within the most sensitive portion of the fields of view of EGRET and COMPTEL (i.e., within a 30 zenith angle), for a CGRO pointing directed exactly at the SGP. Our South Galactic Pole Survey yielded a number of exciting results. The EGRET data were analyzed using the maximum likelihood techniques to estimate the intensity, spectrum, and position of gamma-ray sources in the field of view. Our analysis revealed four sources at energies greater than 100 MeV with likelihood ratios corresponding to greater than 30 detections (Vestrand et al. 1996). One of the sources is associated with the well known gamma-ray blazar PKS 0208-512, but the other three were previously unknown. Among the new detections was PKS 2155-304 which is often considered a prototype of the x-ray selected BL Lacs. PKS 2155-304, which was also detected at hard x-ray energies by CGRO/OSSE, is one of the brightest BL Lac objects in the sky at optical through x-ray energies and has a history of rapid, strong multiwavelength variability. As such, it has been the subject of intensive, contemporaneous, multiwavelength monitoring covering radio frequencies to x-ray energies.

Vestrand, W. T.↗

NASA NTRS2004

Historical Study of Radiation Exposures and the Incidence of Cataracts in Astronauts

For over 35 years, astronauts in low Earth orbit or on missions to the moon have been exposed to space radiation comprised of high-energy protons, heavy ions, and secondary neutrons. We reviewed the radiation exposures received by astronauts in space and on Earth, and presented results from the first epidemiological study of cataract incidence in the astronauts. Our data suggested an increased risk for cataracts from space radiation exposures. Using parametric survival analysis and the maximum likelihood method, we estimated the dose-response and age distribution for cataract incidence in astronauts by space radiation. Considering the high-LET dose contributions on specific space missions as well as data from animal studies with neutrons and heavy ions, suggested a linear response with no dose-threshold for cataracts. However, there are unanswered questions related to the importance and the definition of clinically significant cataracts commonly used in radiation protection, especially in light of epidemiological data suggesting that the probability that sub-clinical cataracts will progress is highly dependent on the age at which cataracts appear. We briefly describe a new study that will address the measurement of cataract progression-rates in astronauts and a ground-based comparison group.

Cucinotta, F. A.↗

NASA NTRS2004

Historical Study of Radiation Exposures and the Incidence of Cataracts in Astronauts

For over 35 years, astronauts in low Earth orbit or on missions to the moon have been exposed to space radiation comprised of high-energy protons, heavy ions, and secondary neutrons. We reviewed the radiation exposures received by astronauts in space and on Earth, and presented results from the first epidemiological study of cataract incidence in the astronauts. Our data suggested an increased risk for cataracts from space radiation exposures*. Using parametric survival analysis and the maximum likelihood method, we estimated the dose-response and age distribution for cataract incidence in astronauts by space radiation. Considering the high-LET dose contributions on specific space missions as well as data from animal studies with neutrons and heavy ions, suggested a linear response with no dose-threshold for cataracts. However, there are unanswered questions related to the importance and the definition of "clinically significant" cataracts commonly used in radiation protection, especially in light of epidemiological data suggesting that the probability that "sub-clinical" cataracts will progress is highly dependent on the age at which cataracts appear. We briefly describe a new study that will address the measurement of cataract progression-rates in astronauts and a ground-based comparison group.

Cucinotta, F. A.↗

NASA NTRS1976

Selecting class weights to minimize classification bias in acreage estimation

Preliminary results of experiments being performed to select optimal class weights for use with the maximum likelihood classifier in acreage estimation using remote sensor imagery are presented. These weights will be optimal in the sense that the bias will be minimized in the proportion estimate obtained from the classification results by sample counting. The procedure was tested using Landsat MSS data from an 8 by 9.6 km area of ground truth in Finney County, Kansas.

Belcher, W. M.↗

NASA NTRS1984

Bayesian estimation of normal mixture parameters

A Bayesian, or penalized maximum likelihood, approach to the problem of estimating the parameters of a mixture of multivariate normal distributions is proposed. The Bayesian formulation eliminates the problem of singularities in the likelihood function and results in an attractive EM-like procedure. Although the question of consistency is not settled, it is suggested that the proposed method has certain advantages over both the constrained and unconstrained maximum likelihood procedures.

Peters, C.↗

NASA NTRS1995

On-line estimation of error covariance parameters for atmospheric data assimilation

A simple scheme is presented for on-line estimation of covariance parameters in statistical data assimilation systems. The scheme is based on a maximum-likelihood approach in which estimates are produced on the basis of a single batch of simultaneous observations. Simple-sample covariance estimation is reasonable as long as the number of available observations exceeds the number of tunable parameters by two or three orders of magnitude. Not much is known at present about model error associated with actual forecast systems. Our scheme can be used to estimate some important statistical model error parameters such as regionally averaged variances or characteristic correlation length scales. The advantage of the single-sample approach is that it does not rely on any assumptions about the temporal behavior of the covariance parameters: time-dependent parameter estimates can be continuously adjusted on the basis of current observations. This is of practical importance since it is likely to be the case that both model error and observation error strongly depend on the actual state of the atmosphere. The single-sample estimation scheme can be incorporated into any four-dimensional statistical data assimilation system that involves explicit calculation of forecast error covariances, including optimal interpolation (OI) and the simplified Kalman filter (SKF). The computational cost of the scheme is high but not prohibitive; on-line estimation of one or two covariance parameters in each analysis box of an operational bozed-OI system is currently feasible. A number of numerical experiments performed with an adaptive SKF and an adaptive version of OI, using a linear two-dimensional shallow-water model and artificially generated model error are described. The performance of the nonadaptive versions of these methods turns out to depend rather strongly on correct specification of model error parameters. These parameters are estimated under a variety of conditions, including uniformly distributed model error and time-dependent model error statistics.

Dee, Dick P.↗

NASA NTRS1981

Characterization, parameter estimation, and aircraft response statistics of atmospheric turbulence

A nonGaussian three component model of atmospheric turbulence is postulated that accounts for readily observable features of turbulence velocity records, their autocorrelation functions, and their spectra. Methods for computing probability density functions and mean exceedance rates of a generic aircraft response variable are developed using nonGaussian turbulence characterizations readily extracted from velocity recordings. A maximum likelihood method is developed for optimal estimation of the integral scale and intensity of records possessing von Karman transverse of longitudinal spectra. Formulas for the variances of such parameter estimates are developed. The maximum likelihood and least-square approaches are combined to yield a method for estimating the autocorrelation function parameters of a two component model for turbulence.

Mark, W. D.↗

NASA NTRS1978

On-line parameter estimation using a high sensitivity estimator

An on-line parameter identification method is presented. The method is based on a recursive formulation of the maximum likelihood method, with a significant modification on the gains of the state estimator. In the conventional maximum likelihood method, the Kalman gains are used in the state estimator. This produces unbiased, minimum variance parameter estimates in the presence of process noise and measurement noise, but it also slows the convergence rate when the identification is done on-line. Here we suggest choosing the gains to maximize a measure of the sensitivities of the state estimates to parameter variations. One such criterion is to minimize the trace of the inverse information matrix. This increases the convergence rate significantly. After one or two time constants, the gains are switched to the Kalman values to assure unbiased, minimum-variance estimates. The state estimate will initially be nonoptimal, and may not be adequate for control purposes. In this case, a parallel Kalman filter which uses the identifier's parameter estimates can be used. This method is applied here for the identification of a simple first-order system, and for the identification of short-period stability derivatives of an F-8 aircraft from simulated data.

Mishne, D.↗

NASA NTRS1981

Identifiability of unknown noise covariance matrices for some special cases of a linear, time-invariant, discrete-time dynamic system

Two off-line schemes are proposed for the identification of unknown noise covariance matrices Q and R of a discrete-time dynamic system. The first scheme is based on a maximum a posteriori cost function utilizing smoothed state estimates, while the second is based on a maximum likelihood cost function utilizing filtered state estimates. Sensitivity of the cost functions to Q and R is analyzed for the following cases: (1) single-input single-output systems; (2) multiinput single-output systems; and (3) single-input multioutput systems with a diagonal R. Identifiability criteria are presented for the cases considered and demonstrated by examples.

Tsang, W. L.↗

NASA NTRS1986

Effects of wing modification on an aircraft's aerodynamic parameters as determined from flight data

A study of the effects of four wing-leading-edge modifications on a general aviation aircraft's stability and control parameters is presented. Flight data from the basic aircraft configuration and configurations with wing modifications are analyzed to determine each wing geometry's stability and control parameters. The parameter estimates and aerodynamic model forms are obtained using the stepwise regression and maximum likelihood techniques. The resulting parameter estimates and aerodynamic models are verified using vortex-lattice theory and by analysis of each model's ability to predict aircraft behavior. Comparisons of the stability and control derivative estimates from the basic wing and the four leading-edge modifications are accomplished so that the effects of each modification on aircraft stability and control derivatives can be determined.

Hess, R. A.↗