Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “maximum likelihood”

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

How Well Does the Latest Anthropomorphic Test Device Mimic Human Impact Responses?

One of the goals of the NASA Occupant Protection Group is to understand the human tolerance to dynamic loading. This knowledge has to come through indirect approaches such as existing human response databases, anthropometric test devices (ATD), animal testing, post‐mortem human subjects, and models. This study investigated the biofidelity of the National Highway Traffic Safety Administration's ATD named the THOR (test device for human occupant restraint). If THOR responds comparably to humans, then it could potentially be used as a human surrogate to help validate space vehicle requirements for occupant protection. The THOR responses to frontal and spinal impacts (ranging from 8 to 12 G with rise times of 40, 70, and 100 ms) were measured and compared to human volunteer responses (95 trials in frontal and 58 in spinal) previously collected by the U. S. Air Force on the same horizontal impact accelerator. The impact acceleration profiles tested are within the expected range of multi‐purpose crew vehicle (MPCV) landing dynamics. A correlation score was calculated for each THOR to human comparison using CORA (CORrelation and Analysis) software. A two‐parameter beta distribution model fit was obtained for each dependent variable using maximum likelihood estimation. For frontal impacts, the THOR head x‐acceleration peak response correlated with the human response at 8 and 10‐G 100 ms but not 10‐G 70 ms. The phase lagged the human response. Head z‐acceleration was not correlated. Chest x‐acceleration was in phase, had a higher peak response, and was well correlated with lighter subjects (Cora = 0.8 for 46 kg vs. Cora = 0.4 for 126 kg). Head x‐displacement had a leading phase. Several subjects responded with the same peak displacement but the mean of the group was lower. The shoulder x‐displacement was in phase but had higher peaks than the human response. For spinal impacts, the THOR head x‐acceleration was not well correlated. Head and chest z‐acceleration was in phase but had a higher peak response. Chest z‐acceleration was highly correlated with heavier subjects at lower G pulses (Cora = 0.86 for 125 kg at 8 G). The human response was variable in shoulder z‐displacement but the THOR was in phase and was comparable to the mean peak response. Head xand z‐displacement was in phase but had higher peaks. Seat pan forces were well correlated, were in phase, but had a larger peak response than most subjects. The THOR does not respond to frontal and spinal impacts exactly the same way that a human does. Some responses are well matched and others are not. Understanding the strengths and weaknesses of this ATD is an important first step in determining its usefulness in occupant protection at NASA

Newby, N.↗

How Well Does the Latest Anthropomorphic Test Device Mimic Human Impact Responses?

One of the goals of the NASA Occupant Protection Group is to understand the human tolerance to dynamic loading. This knowledge has to come through indirect approaches such as existing human response databases, anthropometric test devices (ATD), animal testing, post-­‐mortem human subjects, and models. This study investigated the biofidelity of the National Highway Traffic Safety Administration's ATD named the THOR (test device for human occupant restraint). If THOR responds comparably to humans, then it could potentially be used as a human surrogate to help validate space vehicle requirements for occupant protection. The THOR responses to frontal and spinal impacts (ranging from 8 to 12 G with rise times of 40, 70, and 100 ms) were measured and compared to human volunteer responses (95 trials in frontal and 58 in spinal) previously collected by the U. S. Air Force on the same horizontal impact accelerator. The impact acceleration profiles tested are within the expected range of multi-­‐purpose crew vehicle (MPCV) landing dynamics. A correlation score was calculated for each THOR to human comparison using CORA (CORrelation and Analysis) software. A two-­‐parameter beta distribution model fit was obtained for each dependent variable using maximum likelihood estimation. For frontal impacts, the THOR head x-­‐acceleration peak response correlated with the human response at 8 and 10-­‐G 100 ms but not 10-­‐G 70 ms. The phase lagged the human response. Head z-­‐acceleration was not correlated. Chest x-­‐acceleration was in phase, had a higher peak response, and was well correlated with lighter subjects (Cora = 0.8 for 46 kg vs. Cora = 0.4 for 126 kg). Head x-­‐displacement had a leading phase. Several subjects responded with the same peak displacement but the mean of the group was lower. The shoulder x-­‐displacement was in phase but had higher peaks than the human response. For spinal impacts, the THOR head x-­‐acceleration was not well correlated. Head and chest z-­‐acceleration was in phase but had a higher peak response. Chest z-­‐acceleration was highly correlated with heavier subjects at lower G pulses (Cora = 0.86 for 125 kg at 8 G). The human response was variable in shoulder z-­‐displacement but the THOR was in phase and was comparable to the mean peak response. Head x-­‐ and z-­‐displacement was in phase but had higher peaks. Seat pan forces were well correlated, were in phase, but had a larger peak response than most subjects. The THOR does not respond to frontal and spinal impacts exactly the same way that a human does. Some responses are well matched and others are not. Understanding the strengths and weaknesses of this ATD is an important first step in determining its usefulness in occupant protection at NASA

Newby, Nate↗

A Parameter Estimation Scheme for Multiscale Kalman Smoother (MKS) Algorithm Used in Precipitation Data Fusion

A new approach is presented in this paper to effectively obtain parameter estimations for the Multiscale Kalman Smoother (MKS) algorithm. This new approach has demonstrated promising potentials in deriving better data products based on data of different spatial scales and precisions. Our new approach employs a multi-objective (MO) parameter estimation scheme (called MO scheme hereafter), rather than using the conventional maximum likelihood scheme (called ML scheme) to estimate the MKS parameters. Unlike the ML scheme, the MO scheme is not simply built on strict statistical assumptions related to prediction errors and observation errors, rather, it directly associates the fused data of multiple scales with multiple objective functions in searching best parameter estimations for MKS through optimization. In the MO scheme, objective functions are defined to facilitate consistency among the fused data at multiscales and the input data at their original scales in terms of spatial patterns and magnitudes. The new approach is evaluated through a Monte Carlo experiment and a series of comparison analyses using synthetic precipitation data. Our results show that the MKS fused precipitation performs better using the MO scheme than that using the ML scheme. Particularly, improvements are significant compared to that using the ML scheme for the fused precipitation associated with fine spatial resolutions. This is mainly due to having more criteria and constraints involved in the MO scheme than those included in the ML scheme. The weakness of the original ML scheme that blindly puts more weights onto the data associated with finer resolutions is overcome in our new approach.

multiscale↗

Assessing the Effect of Spaceflight on the Propensity for Astronauts to Develop Disk Herniation

BACKGROUND: A previous study [1] reported that the instantaneous risk of developing a Herniated Nucleus Pulposus (HNP) was higher in astronauts who had flown at least one mission, as compared with those in the corps who had not yet flown. However, the study only analyzed time to HNP after the first mission (if any) and did not account for the possible effects of multiple missions. While many HNP's occurred well into astronauts' careers or in some cases years after retirement, the higher incidence of HNPs relatively soon after completion of space missions appears to indicate that spaceflight may lead to an increased risk of HNP. The purpose of this study was to support the Human System Risk Board assessment of back pain, evaluate the risk of injury due to dynamic loads, and update the previous dataset which contained events up to December 31, 2006. METHODS: Data was queried from the electronic medical record and provided by the Lifetime Surveillance of Astronaut Health. The data included all 330 United States astronauts from 1959 through February 2014. Cases were confirmed by Magnetic Resonance Imaging, Computerized Tomography, Myelography, operative findings, or through clinical confirmation with a neurologist or neurosurgeon. In this analysis, astronauts who had an HNP at selection into the corps or had an HNP diagnosis prior to their first flight were excluded. The statistical challenges in using the available data to separate effects of spaceflight from those associated with general astronaut training and lifestyle on propensity to develop HNPs are many. The primary outcome is reported date of first HNP (if any), which at best is only an approximation to the actual time of occurrence. To properly analyze this data with a survival analysis model, one must also know the "exposure" time - i.e. how long each astronaut has been at risk for developing an HNP. If an HNP is reported soon after a mission, is it mission caused or general? If the former, exposure time should be counted from the time of landing (assuming the risk of HNP occurring during a mission is zero). If the latter, exposure time should be counted from the time of selection; however we can't directly know which one to use. In our analysis we take both of these possibilities into account with a competing risks model, wherein two distinct stochastic processes are going on: TG = time to HNP (general) and TS = time to HNP (spaceflight). Under this type of model, whichever of these occurs first is what we observe; in other words we don't observe TG or TS, only min(TG, TS). Here, we parameterized the model in terms of separate Weibull hazard functions for each process and estimated all parameters using maximum likelihood. In addition, we allowed for a "cured fraction" - i.e. the possibility that some astronauts may never develop an HNP. RESULTS: Results will include a depiction of the competing hazard functions as well as a probability curve for the relative likelihood that an HNP reported at a given time after a mission is actually mission caused. Other factors, such as dwell time in microgravity and vehicle landing environment will be explored. An overall assessment as to whether spaceflight truly exacerbates HNP risk will be made.

Feiveson, A. H.↗

Achievable Precision for Optical Ranging Systems

Achievable RMS errors in estimating the phase, frequency, and intensity of a direct-detected intensity-modulated optical pulse train are presented. For each parameter, the Cramer-Rao-Bound (CRB) is derived and the performance of the Maximum Likelihood estimator is illustrated. Approximations to the CRBs are provided, enabling an intuitive understanding of estimator behavior as a function of the signaling parameters. The results are compared to achievable RMS errors in estimating the same parameters from a sinusoidal waveform in additive white Gaussian noise. This establishes a framework for a performance comparison of radio frequency (RF) and optical science. Comparisons are made using parameters for state-of-the-art deep-space RF and optical links. Degradations to the achievable errors due to clock phase noise and detector jitter are illustrated.

estimation theory↗

Sensitivity Analysis for Atmospheric Infrared Sounder (AIRS) CO2 Retrieval

The Atmospheric Infrared Sounder (AIRS) is a thermal infrared sensor able to retrieve the daily atmospheric state globally for clear as well as partially cloudy field-of-views. The AIRS spectrometer has 2378 channels sensing from 15.4 micrometers to 3.7 micrometers, of which a small subset in the 15 micrometers region has been selected, to date, for CO2 retrieval. To improve upon the current retrieval method, we extended the retrieval calculations to include a prior estimate component and developed a channel ranking system to optimize the channels and number of channels used. The channel ranking system uses a mathematical formalism to rapidly process and assess the retrieval potential of large numbers of channels. Implementing this system, we identifed a larger optimized subset of AIRS channels that can decrease retrieval errors and minimize the overall sensitivity to other iridescent contributors, such as water vapor, ozone, and atmospheric temperature. This methodology selects channels globally by accounting for the latitudinal, longitudinal, and seasonal dependencies of the subset. The new methodology increases accuracy in AIRS CO2 as well as other retrievals and enables the extension of retrieved CO2 vertical profiles to altitudes ranging from the lower troposphere to upper stratosphere. The extended retrieval method for CO2 vertical profile estimation using a maximum-likelihood estimation method. We use model data to demonstrate the beneficial impact of the extended retrieval method using the new channel ranking system on CO2 retrieval.

sensitivity analysis↗

Studies in Astronomical Time Series Analysis. VI. Bayesian Block Representations

This paper addresses the problem of detecting and characterizing local variability in time series and other forms of sequential data. The goal is to identify and characterize statistically significant variations, at the same time suppressing the inevitable corrupting observational errors. We present a simple nonparametric modeling technique and an algorithm implementing it-an improved and generalized version of Bayesian Blocks [Scargle 1998]-that finds the optimal segmentation of the data in the observation interval. The structure of the algorithm allows it to be used in either a real-time trigger mode, or a retrospective mode. Maximum likelihood or marginal posterior functions to measure model fitness are presented for events, binned counts, and measurements at arbitrary times with known error distributions. Problems addressed include those connected with data gaps, variable exposure, extension to piece- wise linear and piecewise exponential representations, multivariate time series data, analysis of variance, data on the circle, other data modes, and dispersed data. Simulations provide evidence that the detection efficiency for weak signals is close to a theoretical asymptotic limit derived by [Arias-Castro, Donoho and Huo 2003]. In the spirit of Reproducible Research [Donoho et al. (2008)] all of the code and data necessary to reproduce all of the figures in this paper are included as auxiliary material.

signal detection↗

Predicting the Redshift 2 H-Alpha Luminosity Function Using [OIII] Emission Line Galaxies

Upcoming space-based surveys such as Euclid and WFIRST-AFTA plan to measure Baryonic Acoustic Oscillations (BAOs) in order to study dark energy. These surveys will use IR slitless grism spectroscopy to measure redshifts of a large number of galaxies over a significant redshift range. In this paper, we use the WFC3 Infrared Spectroscopic Parallel Survey (WISP) to estimate the expected number of H-alpha emitters observable by these future surveys. WISP is an ongoing Hubble Space Telescope slitless spectroscopic survey, covering the 0.8 - 1.65 micrometers wavelength range and allowing the detection of H-alpha emitters up to z approximately equal to 1.5 and [OIII] emitters to z approximately equal to 2.3. We derive the H-alpha-[OIII] bivariate line luminosity function for WISP galaxies at z approximately equal to 1 using a maximum likelihood estimator that properly accounts for uncertainties in line luminosity measurement, and demonstrate how it can be used to derive the H-alpha luminosity function from exclusively fitting [OIII] data. Using the z approximately equal to 2 [OIII] line luminosity function, and assuming that the relation between H-alpha and [OIII] luminosity does not change significantly over the redshift range, we predict the H-alpha number counts at z approximately equal to 2 - the upper end of the redshift range of interest for the future surveys. For the redshift range 0.7 less than z less than 2, we expect approximately 3000 galaxies per sq deg for a flux limit of 3 x 10(exp −16) ergs per sec per sq cm (the proposed depth of Euclid galaxy redshift survey) and approximately 20,000 galaxies per sq deg for a flux limit of approximately 10(exp −16) ergs per sec per sq cm (the baseline depth of WFIRST galaxy redshift survey).

Redshift↗

Validation of Rain Rate Retrievals for the Airborne Hurricane Imaging Radiometer (HIRAD)

The NASA Hurricane and Severe Storm Sentinel (HS3) mission is an aircraft field measurements program using NASA's unmanned Global Hawk aircraft system for remote sensing and in situ observations of Atlantic and Caribbean Sea hurricanes. One of the principal microwave instruments is the Hurricane Imaging Radiometer (HIRAD), which measures surface wind speeds and rain rates. For validation of the HIRAD wind speed measurement in hurricanes, there exists a comprehensive set of comparisons with the Stepped Frequency Microwave Radiometer (SFMR) with in situ GPS dropwindsondes [1]. However, for rain rate measurements, there are only indirect correlations with rain imagery from other HS3 remote sensors (e.g., the dual-frequency Ka- & Ku-band doppler radar, HIWRAP), which is only qualitative in nature. However, this paper presents results from an unplanned rain rate measurement validation opportunity that occurred in 2013, when HIRAD flew over an intense tropical squall line that was simultaneously observed by the Tampa NEXRAD meteorological radar (Fig. 1). During this experiment, Global Hawk flying at an altitude of 18 km made 3 passes over the rapidly propagating thunderstorm, while the TAMPA NEXRAD perform volume scans on a 5-minute interval. Using the well-documented NEXRAD Z-R relationship, 2D images of rain rate (mm/hr) were obtained at two altitudes (3 km & 6 km), which serve as surface truth for the HIRAD rain rate retrievals. A preliminary comparison of HIRAD rain rate retrievals (image) for the first pass and the corresponding closest NEXRAD rain image is presented in Fig. 2 & 3. This paper describes the HIRAD instrument, which 1D synthetic-aperture thinned array radiometer (STAR) developed by NASA Marshall Space Flight Center [2]. The rain rate retrieval algorithm, developed by Amarin et al. [3], is based on the maximum likelihood estimation (MLE) technique, which compares the observed Tb's at the HIRAD operating frequencies of 4, 5, 6 and 6.6 GHz with corresponding theoretical Tb values from a forward radiative transfer model (RTM). The optimum solution is the integrated rain rate that minimizes the difference between RTM and observed values. Because the excess Tb from rain comes from the direct upwelling and the indirect reflected downwelling paths through the atmosphere, there are several assumptions made for the 2D rain distribution in the antenna incident plane (crosstrack to flight direction). The opportunity to knowing 2D rain surface truth from NEXRAD at two different altitudes will enable a comprehensive evaluation to be preformed and reported in this paper.

hurricane↗

Experimental Validation of Pulse Phase Tracking for X-Ray Pulsar Based

Pulsars are a form of variable celestial source that have shown to be usable as aids for autonomous, deep space navigation. Particularly those sources emitting in the X-ray band are ideal for navigation due to smaller detector sizes. In this paper X-ray photons arriving from a pulsar are modeled as a non-homogeneous Poisson process. The method of pulse phase tracking is then investigated as a technique to measure the radial distance traveled by a spacecraft over an observation interval. A maximum-likelihood phase estimator (MLE) is used for the case where the observed frequency signal is constant. For the varying signal frequency case, an algorithm is used in which the observation window is broken up into smaller blocks over which an MLE is used. The outputs of this phase estimation process were then looped through a digital phase-locked loop (DPLL) in order to reduce the errors and produce estimates of the doppler frequency. These phase tracking algorithms were tested both in a computer simulation environment and using the NASA Goddard Space flight Center X-ray Navigation Laboratory Testbed (GXLT). This provided an experimental validation with photons being emitted by a modulated X-ray source and detected by a silicon-drift detector. Models of the Crab pulsar and the pulsar B1821-24 were used in order to generate test scenarios. Three different simulated detector trajectories were used to be tracked by the phase tracking algorithm: a stationary case, one with constant velocity, and one with constant acceleration. All three were performed in one-dimension along the line of sight to the pulsar. The first two had a constant signal frequency and the third had a time varying frequency. All of the constant frequency cases were processed using the MLE, and it was shown that they tracked the initial phase within 0.15% for the simulations and 2.5% in the experiments, based on an average of ten runs. The MLE-DPLL cascade version of the phase tracking algorithm was used in the varying frequency case. This resulted in tracking of the phase and frequency by the DPLL outputs in both the simulation and experimental environments. The crab pulsar was experimentally tested with a trajectory with a higher acceleration. In this case the phase error tended toward zero as the observation extended to 250 seconds and the doppler frequency error tended to zero in under 100 seconds.

Experimental↗

Data Analysis & Statistical Methods for Command File Errors

This paper explains current work on modeling for managing the risk of command file errors. It is focused on analyzing actual data from a JPL spaceflight mission to build models for evaluating and predicting error rates as a function of several key variables. We constructed a rich dataset by considering the number of errors, the number of files radiated, including the number commands and blocks in each file, as well as subjective estimates of workload and operational novelty. We have assessed these data using different curve fitting and distribution fitting techniques, such as multiple regression analysis, and maximum likelihood estimation to see how much of the variability in the error rates can be explained with these. We have also used goodness of fit testing strategies and principal component analysis to further assess our data. Finally, we constructed a model of expected error rates based on the what these statistics bore out as critical drivers to the error rate. This model allows project management to evaluate the error rate against a theoretically expected rate as well as anticipate future error rates.

Correlation Analysis↗

ANTARES Constrains a Blazar Origin of Two IceCube PeV Neutrino Events

Context. The source(s) of the neutrino excess reported by the IceCube Collaboration is unknown. The TANAMI Collaboration recently reported on the multiwavelength emission of six bright, variable blazars which are positionally coincident with two of the most energetic IceCube events. Objects like these are prime candidates to be the source of the highest-energy cosmic rays, and thus of associated neutrino emission. Aims. We present an analysis of neutrino emission from the six blazars using observations with the ANTARES neutrino telescope. Methods. The standard methods of the ANTARES candidate list search are applied to six years of data to search for an excess of muons - and hence their neutrino progenitors - from the directions of the six blazars described by the TANAMI Collaboration, and which are possibly associated with two IceCube events. Monte Carlo simulations of the detector response to both signal and background particle fluxes are used to estimate the sensitivity of this analysis for di erent possible source neutrino spectra. A maximum-likelihood approach, using the reconstructed energies and arrival directions of through-going muons, is used to identify events with properties consistent with a blazar origin. Results. Both blazars predicted to be the most neutrino-bright in the TANAMI sample (1653-329 and 1714-336) have a signal flux fitted by the likelihood analysis corresponding to approximately one event. This observation is consistent with the blazar-origin hypothesis of the IceCube event IC 14 for a broad range of blazar spectra, although an atmospheric origin cannot be excluded. No ANTARES events are observed from any of the other four blazars, including the three associated with IceCube event IC20. This excludes at a 90% confidence level the possibility that this event was produced by these blazars unless the neutrino spectrum is flatter than -2.4.

general↗

Determination of the Risk of Radiation-Associated Circulatory and Cancer Disease Mortality in a NASA Early Astronaut Cohort

Of the many possible health challenges posed during extended exploratory missions to space, the effects of space radiation on cardiovascular disease and cancer are of particular concern. There are unique challenges to estimating those radiation risks; care and appropriate and rigorous methodology should be applied when considering small cohorts such as the NASA astronaut population. The objective of this work was to determine if there was sufficient evidence for excess risk of cardiovascular disease and cancer in early NASA astronaut cohorts. NASA astronauts in selection groups 1-7 were chosen; this relatively homogeneous cohort consists of 73 white males, who unlike today's astronauts, maintained similar smoking and drinking habits to the general US population, and have published radiation doses. The participants flew in space on missions Mercury through Shuttle and received space radiation doses between 0-74.1 milligrays. Cause of death information was obtained from the Lifetime Surveillance of Astronaut Health (LSAH) program at NASA Johnson Space Center. Mortality was compared with the US male population. Trends of mortality with dose were assessed using a logistic model, fitted by maximum likelihood. Only 32 (43.84 percent) of the 73 early astronauts have died. Standard mortality ratios (SMRs) for cancer (n=7, SMR=43.4, 95 percent CI 17.8, 84.9), all circulatory disease (n=7, SMR=33.2, 95 percent CI 13.7, 65.0), and ischemic heart disease (IHD) (n=5, SMR=40.1, 95 percent CI 13.2, 89.4) were significantly lower than for the US white male population. For cerebrovascular disease, the upper confidence interval for SMR included 100, indicating it was not significantly different from the US population (n=2, SMR = 77.0, 95 percent CI 9.4, 268.2). The power of the study is low and remains below 10 percent even when risks 10 times those reported in the literature are assumed. Due to small sample size, there is currently insufficient statistical power to evaluate space radiation exposure effects on mortality in NASA astronauts. In addition to a comprehensive longitudinal study of NASA astronauts, a research strategy of low dose epidemiology data integration with cell and animal studies should be utilized for space radiation risk assessment in the astronauts.

Elgart, S. R.↗

Trajectory Dispersed Vehicle Process for Space Launch System

The Space Launch System (SLS) vehicle is part of NASA's deep space exploration plans that includes manned missions to Mars. Manufacturing uncertainties in design parameters are key considerations throughout SLS development as they have significant effects on focus parameters such as lift-off-thrust-to-weight, vehicle payload, maximum dynamic pressure, and compression loads. This presentation discusses how the SLS program captures these uncertainties by utilizing a 3 degree of freedom (DOF) process called Trajectory Dispersed (TD) analysis. This analysis biases nominal trajectories to identify extremes in the design parameters for various potential SLS configurations and missions. This process utilizes a Design of Experiments (DOE) and response surface methodologies (RSM) to statistically sample uncertainties, and develop resulting vehicles using a Maximum Likelihood Estimate (MLE) process for targeting uncertainties bias. These vehicles represent various missions and configurations which are used as key inputs into a variety of analyses in the SLS design process, including 6 DOF dispersions, separation clearances, and engine out failure studies.

Statham, Tamara↗

Modeling of Time-Dependent Strength Degradation of SiC/SiC Ceramic Matrix Composites via MAC/GMC Computer Code

Silicon carbide fiber reinforced silicon carbide (SiC/SiC) ceramic matrix composites (CMCs) display time-dependent strength degradation at intermediate temperatures (600 to 900 degrees Centigrade). This is generally believed to be an oxidation induced phenomenon. The understanding of the effect of temperature with environment (oxidation) is key towards development of SiC/SiC CMCs with a reliable load carrying capacity. Various theories have been proposed to explain the strength degradation. One suggests that the boron nitride (BN) coating deposited on the fibers oxidizes causing fusion of fibers. Another theory proposes that the SiC fibers are oxidized forming a silica scale leading to premature fiber failure. A more recent theory suggests that SiC fiber strength is intrinsically time-dependent due to slow crack growth of flaws in the fibers. An empirical model, termed as a "fiber classic model," which is based on a standard slow crack growth type power-law, has been implemented within NASA's micromechanics-based MAC/GMC computer code as a user routine. Model parameters for this "classic model" were calibrated from stress-rupture data of Hi-Nicalon (TM) monofilaments using the maximum likelihood estimation (MLE) technique. This new capability in the MAC/GMC computer code was then used to predict the stress-rupture behavior of Hi-Nicalon (TM) tows as well as 2-D SiC/SiC composites reinforced with Hi-Nicalon (TM) fibers. Results demonstrate that the MAC/GMC with this new capability successfully predicts the time-to-failure vs. applied stress within the intermediate temperature range at various scales as well as laminated composites in an oxidizing environment.

Environmental degradation↗

Advances in Hyperspectral Image Classification Methods for Vegetation and Agricultural Cropland Studies

Hyperspectral data are becoming more widely available via sensors on airborne and unmanned aerial vehicle (UAV) platforms, as well as proximal platforms. While space-based hyperspectral data continue to be limited in availability, multiple spaceborne Earth-observing missions on traditional platforms are scheduled for launch, and companies are experimenting with small satellites for constellations to observe the Earth, as well as for planetary missions. Land cover mapping via classification is one of the most important applications of hyperspectral remote sensing and will increase in significance as time series of imagery are more readily available. However, while the narrow bands of hyperspectral data provide new opportunities for chemistry-based modeling and mapping, challenges remain. Hyperspectral data are high dimensional, and many bands are highly correlated or irrelevant for a given classification problem. For supervised classification methods, the quantity of training data is typically limited relative to the dimension of the input space. The resulting Hughes phenomenon, often referred to as the curse of dimensionality, increases potential for unstable parameter estimates, overfitting, and poor generalization of classifiers. This is particularly problematic for parametric approaches such as Gaussian maximum likelihood–based classifiers that have been the backbone of pixel-based multispectral classification methods. This issue has motivated investigation of alternatives, including regularization of the class covariance matrices, ensembles of weak classifiers, development of feature selection and extraction methods, adoption of nonparametric classifiers, and exploration of methods to exploit unlabeled samples via semi-supervised and active learning. Data sets are also quite large, motivating computationally efficient algorithms and implementations. This chapter provides an overview of the recent advances in classification methods for mapping vegetation using hyperspectral data. Three data sets that are used in the hyperspectral classification literature (e.g., Botswana Hyperion satellite data and AVIRIS airborne data over both Kennedy Space Center and Indian Pines) are described in Section 3.2 and used to illustrate methods described in the chapter. An additional high-resolution hyperspectral data set acquired by a SpecTIR sensor on an airborne platform over the Indian Pines area is included to exemplify the use of new deep learning approaches, and a multiplatform example of airborne hyperspectral data is provided to demonstrate transfer learning in hyperspectral image classification. Classical approaches for supervised and unsupervised feature selection and extraction are reviewed in Section 3.3. In particular, nonlinearities exhibited in hyperspectral imagery have motivated development of nonlinear feature extraction methods in manifold learning, which are outlined in Section 3.3.1.4. Spatial context is also important in classification of both natural vegetation with complex textural patterns and large agricultural fields with significant local variability within fields. Approaches to exploit spatial features at both the pixel level (e.g., co-occurrence–based texture and extended morphological attribute profiles [EMAPs]) and integration of segmentation approaches (e.g., HSeg) are discussed in this context in Section 3.3.2. Recently, classification methods that leverage nonparametric methods originating in the machine learning community have grown in popularity. An overview of both widely used and newly emerging approaches, including support vector machines (SVMs), Gaussian mixture models, and deep learning based on convolutional neural networks is provided in Section 3.4. Strategies to exploit unlabeled samples, including active learning and metric learning, which combine feature extraction and augmentation of the pool of training samples in an active learning framework, are outlined in Section 3.5. Integration of image segmentation with classification to accommodate spatial coherence typically observed in vegetation is also explored, including as an integrated active learning system. Exploitation of multisensor strategies for augmenting the pool of training samples is investigated via a transfer learning framework in Section 3.5.1.2. Finally, we look to the future, considering opportunities soon to be provided by new paradigms, as hyperspectral sensing is becoming common at multiple scales from ground-based and airborne autonomous vehicles to manned aircraft and space-based platforms.

Pasolli, Edoardo↗

Super Resolution of Remote Sensing Images Using Edge-Directed Radial Basis Functions

Edge-Directed Radial Basis Functions (EDRBF) are used to compute super resolution(SR) image from a given set of low resolution (LR) images differing in subpixel shifts. The algorithm is tested on remote sensing images and compared for accuracy with other well-known algorithms such as Iterative Back Projection (IBP), Maximum Likelihood (ML) algorithm, interpolation of scattered points using Nearest Neighbor (NN) and Inversed Distance Weighted (IDW) interpolation, and Radial Basis Functin(RBF) . The accuracy of SR depends on various factors besides the algorithm (i) number of subpixel shifted LR images (ii) accuracy with which the LR shifts are estimated by registration algorithms (iii) and the targeted spatial resolution of SR. In our studies, the accuracy of EDRBF is compared with other algorithms keeping these factors constant. The algorithm has two steps: i) registration of low resolution images and (ii) estimating the pixels in High Resolution (HR) grid using EDRBF. Experiments are conducted by simulating LR images from a input HR image with different sub-pixel shifts. The reconstructed SR image is compared with input HR image to measure the accuracy of the algorithm using sum of squared errors (SSE). The algorithm has outperformed all of the algorithms mentioned above. The algorithm is robust and is not overly sensitive to the registration inaccuracies.

Mareboyana, Manohar M.↗

Radar Autofocus Algorithm Incorporating a priori Terrain Knowledge for Correction of Mars’ Ionospheric Distortion in MARSIS Observations

Low-frequency subsurface radar observations of Mars’ polar ice deposits by MARSIS (Mars Advanced Radar for Subsurface and Ionosphere Sounding) are heavily impacted by the electron content of Mars’ ionosphere. The resulting ionospheric distortion can be represented as attenuation and bulk delay, in addition to higher-order frequency dispersion effects. Baseline, uncorrected images are often unusable when the solar zenith angle is less than 90◦ (day side). In this work, a radar autofocus algorithm is developed that estimates and inverts ionospheric distortion, producing a focused radargram of the ice deposit subsurface. Previously published autofocus algorithms have sought to maximize peak-to-noise contrast, which may yield sub-optimal results for complex terrain. Instead, a maximum likelihood approach is developed that incorporates simulated surface clutter returns for the current spacecraft position, based on the Mars Orbiter Laser Altimeter (MOLA) elevation model of the Martian surface. An ancillary product is a surface-only clutter simulation for each orbit, which is necessary to identify true subsurface features.

McMichael, Joseph G.↗