Search NASA⌕ Search

SEARCH · Search NASA

Results for “Function 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.

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

At least 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Compensation of Oxygen Transmittance Effects for Proximal Sensing Retrieval of Canopy-Leaving Sun-Induced Chlorophyll Fluorescence

Estimates of Sun-Induced vegetation chlorophyll Fluorescence (SIF) using remote sensing techniques are commonly determined by exploiting solar and/or telluric absorption features. When SIF is retrieved in the strong oxygen (O2) absorption features, atmospheric effects must always be compensated. Whereas correction of atmospheric effects is a standard airborne or satellite data processing step, there is no consensus regarding whether it is required for SIF proximal-sensing measurements nor what is the best strategy to be followed. Thus, by using simulated data, this work provides a comprehensive analysis about how atmospheric effects impact SIF estimations on proximal sensing, regarding: (1) the sensor height above the vegetated canopy; (2) the SIF retrieval technique used, e.g., Fraunhofer Line Discriminator (FLD) family or Spectral Fitting Methods (SFM); and (3) the instrument's spectral resolution. We demonstrate that for proximal-sensing scenarios compensating for atmospheric effects by simply introducing the O2 transmittance function into the FLD or SFM formulations improves SIF estimations. However, these simplistic corrections still lead to inaccurate SIF estimations due to the multiplication of spectrally convolved atmospheric transfer functions with absorption features. Consequently, a more rigorous oxygen compensation strategy is proposed and assessed by following a classic airborne atmospheric correction scheme adapted to proximal sensing. This approach allows compensating for the O2 absorption effects and, at the same time, convolving the high spectral resolution data according to the corresponding Instrumental Spectral Response Function (ISRF) through the use of an atmospheric radiative transfer model. Finally, due to the key role of O2 absorption on the evaluated proximal-sensing SIF retrieval strategies, its dependency on surface pressure (p) and air temperature (T) was also assessed. we combined simulated spectral data with p and T measurements obtained for a one-year period in the Hyytiala Forestry Field Station in Finland. Of importance hereby is that seasonal dynamics in terms of T and p, if not appropriately considered as part of the retrieval strategy, can result in erroneous SIF seasonal trends that mimic those of known dynamics for temperature-dependent physiological responses of vegetation.

"sun–induced chlorophyll fluorescence (SIF); pro↗
DOE OSTI2024

Increasing the power of weak lensing data with multipole-based intrinsic alignment estimators

ABSTRACT It has long been known that galaxy shapes align coherently with the large-scale density field. Characterizing this effect is essential to interpreting measurements of weak gravitational lensing, the deflection of light from distant galaxies by matter overdensities along the line of sight, as it also produces coherent galaxy alignments that we wish to interpret in terms of a cosmological model. Existing direct measurements of intrinsic alignments using galaxy samples with high-quality shape and redshift measurements typically use well-understood but sub-optimal projected estimators, which do not make good use of the information in the data when comparing those estimators to theoretical models. We demonstrate a more optimal estimator, based on a multipole expansion of the correlation functions or power spectra, for direct measurements of galaxy intrinsic alignments. We show that even using the lowest order multipole alone increases the significance of inferred model parameters using simulated and real data, without any additional modelling complexity. We apply this estimator to measurements of parameters of the non-linear alignment model using data from the Sloan Digital Sky survey, demonstrating consistent results with a factor of ∼2 greater precision in parameter fits to intrinsic alignments models. This result is functionally equivalent to quadrupling the survey area, but without the attendant costs – thereby demonstrating the value in using this new estimator in current and future intrinsic alignments measurements using spectroscopic galaxy samples.

Astronomy & Astrophysics↗
DOE OSTI2024

Learning Functions Varying along a Central Subspace

Many functions of interest are in a high-dimensional space but exhibit low-dimensional structures. This paper studies regression of an s-Hölder function in $R^D$ which varies along a central subspace of dimension $d$ while $d \ll D$. A direct approximation of $f$ in $R^D$ with an accuracy $\varepsilon$ requires the number of samples in the order of $\varepsilon^{-(2s+D)/s}$. In this paper, we analyze the generalized contour regression (GCR) algorithm for the estimation of the central subspace and use piecewise polynomials for function approximation. GCR is among the best estimators for the central subspace, but its sample complexity is an open question. In this paper, we partially answer this questions by proving that if a variance quantity is exactly known, GCR leads to a mean squared estimation error of $O(n^{-1})$ for the central subspace. The estimation error of this variance quantity is also given in this paper. The mean squared regression error of $f$ is proved to be in the order of $(n/\log n)^{-\frac{2s}{2s+d}}$, where the exponent depends on the dimension of the central subspace instead of the ambient space . This result demonstrates that GCR is effective in learning the low-dimensional central subspace. We also propose a modified GCR with improved efficiency. Here, the convergence rate is validated through several numerical experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

The Variable Saturation Hydraulic Conductivity Method for Improving Soil Water Content Simulation in EPIC and APEX Models

Core Ideas We propose VSHC for calculating soil water percolation in continuous models. VSHC relates permeability to daily hydraulic properties to estimate hydraulic conductivity. VSHC performed well at field scale and improved estimates of peak streamflows. Soil water percolation is a key process in the life cycle of water in fields, watersheds, and river basins. The Environmental Policy Integrated Climate (EPIC) and the Agricultural Policy/Environmental eXtender (APEX) are continuous models developed for evaluating the environmental effects of agricultural management. Traditionally, these models have simulated soil water percolation processes using a tipping‐bucket approach, with the rate of flow limited by the saturated hydraulic conductivity. This simple approach often leads to inaccuracy in simulating elevated soil water conditions where soil water content (SWC) levels may remain above field capacity under prolonged wet weather periods or limited drainage. To overcome this deficiency, a new submodel, the variable saturation hydraulic conductivity (VSHC) method, was developed for simulating soil water percolation processes using a nonlinear equation to estimate the effective hydraulic conductivity as a function of the SWC and soil properties. The VSHC method was evaluated at three sites in the United States and two sites in Europe. In addition, a numerical solution of the Richards equation was used as a benchmark for SWC comparison. Results show that the VSHC method substantially improves the accuracy of the SWC simulation in long‐term simulations, particularly during wet periods. At the watershed scale, results on the Riesel Y2 watershed indicate that the VSHC method enhances model performance in the high‐flow regime of channel peak flows because of the improved estimation of SWC, which implies that the improved SWC simulation at the field scale is beneficial to hydrologic modeling at the watershed scale.

Doro, Luca↗
DOE OSTI2022

Advanced Laboratory and Field Arrays: Debris Modeling, Detection,& Mitigation (Task 1)

The statement of project objectives for this task was: develop tools, methods and models to assess, and mitigate the risk of damage to MHK infrastructure from woody debris. Develop the capability to detect woody debris using sonar and/or physical methods for purposes of characterizing debris statistics in river (at UAF’s Tanana River Test Site) and near-shore wave (at Yakutat, AK) environments and to activate debris mitigation measures. Develop debris impact risk maps and tables using statistics on debris size, geometry, type, prevalence, mobility and location. Improve and apply the COUPi discrete element method (DEM) to develop models of debris movement and impact on MHK infrastructure to evaluate risk of damage, and interference, to operations from debris. The proposed final deliverable for the task was a set of tools or techniques for providing estimates of the probability of debris impact, and resulting impact forces, on MHK infrastructure as a function of debris size, type, wave regime, and current velocity. Such estimates are required to assess damage risks to operational MHK infrastructure.

13 HYDRO ENERGY↗
DOE OSTI2016

LOW-METALLICITY YOUNG CLUSTERS IN THE OUTER GALAXY. II. SH 2-208

We obtained deep near-infrared images of Sh 2-208, one of the lowest-metallicity H ii regions in the Galaxy, [O/H] = −0.8 dex. We detected a young cluster in the center of the H ii region with a limiting magnitude of K = 18.0 mag (10 σ ), which corresponds to a mass detection limit of ∼0.2 M {sub ⊙}. This enables the comparison of star-forming properties under low metallicity with those of the solar neighborhood. We identified 89 cluster members. From the fitting of the K -band luminosity function (KLF), the age and distance of the cluster are estimated to be ∼0.5 Myr and ∼4 kpc, respectively. The estimated young age is consistent with the detection of strong CO emission in the cluster region and the estimated large extinction of cluster members ( A{sub V} ∼ 4–25 mag). The observed KLF suggests that the underlying initial mass function (IMF) of the low-metallicity cluster is not significantly different from canonical IMFs in the solar neighborhood in terms of both high-mass slope and IMF peak (characteristic mass). Despite the very young age, the disk fraction of the cluster is estimated at only 27% ± 6%, which is significantly lower than those in the solar metallicity. Those results are similar to Sh 2-207, which is another star-forming region close to Sh 2-208 with a separation of 12 pc, suggesting that their star-forming activities in low-metallicity environments are essentially identical to those in the solar neighborhood, except for the disk dispersal timescale. From large-scale mid-infrared images, we suggest that sequential star formation is taking place in Sh 2-207, Sh 2-208, and the surrounding region, triggered by an expanding bubble with a ∼30 pc radius.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Stress observations and driving force models for the South American plate

A comparison of stress observations to predictions of driving force models for the South American plate is presented. Conditions for dynamical equilibrium are used to estimate the magnitudes of forces as a function of age, with ridge push and slab pull as active forces, and plate drag and slab resistance as passive forces. The results are consistent with earlier published models, and it is found that the primary balance for an oceanic plate is between slab pull and the resisting drag forces. The model is generalized for the possible forces driving the South American plate, and it is compared with observed stress measurements and kinematics. The plate balance is found to be largely between ridge push and plate drag. The magnitude of a possible trench suction force is estimated using a model with an eddy between the Nazca slab and the overlying South American plate.

Stefanick, Michael↗
NASA NTRS1981

Evaluation of a segment-based LANDSAT full-frame approach to corp area estimation

As the registration of LANDSAT full frames enters the realm of current technology, sampling methods should be examined which utilize other than the segment data used for LACIE. The effect of separating the functions of sampling for training and sampling for area estimation. The frame selected for analysis was acquired over north central Iowa on August 9, 1978. A stratification of he full-frame was defined. Training data came from segments within the frame. Two classification and estimation procedures were compared: statistics developed on one segment were used to classify that segment, and pooled statistics from the segments were used to classify a systematic sample of pixels. Comparisons to USDA/ESCS estimates illustrate that the full-frame sampling approach can provide accurate and precise area estimates.

Bauer, M. E.↗
DOE OSTI2025

A Class of Sparse Johnson–Lindenstrauss Transforms and Analysis of their Extreme Singular Values

The Johnson–Lindenstrauss (JL) lemma is a powerful tool for dimensionality reduction in modern algorithm design. The lemma states that any set of high-dimensional points in a Euclidean space can be projected into lower dimensions while approximately preserving pairwise Euclidean distances. Random matrices satisfying this lemma are called JL transforms (JLTs). Inspired by existing $s$-hashing JLTs with exactly $s$ nonzero elements on each column, the present work introduces an ensemble of sparse matrices encompassing so-called $s$-hashing-like matrices whose expected number of nonzero elements on each column is $s$. The independence of the sub-Gaussian entries of these matrices and the knowledge of their exact distribution play an important role in their analyses. Using properties of independent sub-Gaussian random variables, these matrices are demonstrated to be JLTs, and their smallest nontrivial singular values and largest singular values are estimated nonasymptotically using a technique from geometric functional analysis. As the dimensions of the matrix grow to infinity, these singular values are proved to converge almost surely to fixed quantities (by using the universal Bai–Yin law) and in distribution to the Gaussian orthogonal ensemble Tracy–Widom law after proper rescalings. Understanding the behaviors of extreme singular values is important in general because they are often used to define a measure of stability of matrix algorithms. For example, JLTs were recently used in derivative-free optimization algorithmic frameworks to select random subspaces in which are constructed random models or poll directions to achieve scalability, and hence estimating their smallest singular value in particular helps determine the dimension of these subspaces.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Maximum likelihood tuning of a vehicle motion filter

This paper describes the use of maximum likelihood parameter estimation unknown parameters appearing in a nonlinear vehicle motion filter. The filter uses the kinematic equations of motion of a rigid body in motion over a spherical earth. The nine states of the filter represent vehicle velocity, attitude, and position. The inputs to the filter are three components of translational acceleration and three components of angular rate. Measurements used to update states include air data, altitude, position, and attitude. Expressions are derived for the elements of filter matrices needed to use air data in a body-fixed frame with filter states expressed in a geographic frame. An expression for the likelihood functions of the data is given, along with accurate approximations for the function's gradient and Hessian with respect to unknown parameters. These are used by a numerical quasi-Newton algorithm for maximizing the likelihood function of the data in order to estimate the unknown parameters. The parameter estimation algorithm is useful for processing data from aircraft flight tests or for tuning inertial navigation systems.

Trankle, Thomas L.↗
NASA NTRS2007

Broadband Phase Retrieval for Image-Based Wavefront Sensing

A focus-diverse phase-retrieval algorithm has been shown to perform adequately for the purpose of image-based wavefront sensing when (1) broadband light (typically spanning the visible spectrum) is used in forming the images by use of an optical system under test and (2) the assumption of monochromaticity is applied to the broadband image data. Heretofore, it had been assumed that in order to obtain adequate performance, it is necessary to use narrowband or monochromatic light. Some background information, including definitions of terms and a brief description of pertinent aspects of image-based phase retrieval, is prerequisite to a meaningful summary of the present development. Phase retrieval is a general term used in optics to denote estimation of optical imperfections or aberrations of an optical system under test. The term image-based wavefront sensing refers to a general class of algorithms that recover optical phase information, and phase-retrieval algorithms constitute a subset of this class. In phase retrieval, one utilizes the measured response of the optical system under test to produce a phase estimate. The optical response of the system is defined as the image of a point-source object, which could be a star or a laboratory point source. The phase-retrieval problem is characterized as image-based in the sense that a charge-coupled-device camera, preferably of scientific imaging quality, is used to collect image data where the optical system would normally form an image. In a variant of phase retrieval, denoted phase-diverse phase retrieval [which can include focus-diverse phase retrieval (in which various defocus planes are used)], an additional known aberration (or an equivalent diversity function) is superimposed as an aid in estimating unknown aberrations by use of an image-based wavefront-sensing algorithm. Image-based phase-retrieval differs from such other wavefront-sensing methods, such as interferometry, shearing interferometry, curvature wavefront sensing, and Shack-Hartmann sensing, all of which entail disadvantages in comparison with image-based methods. The main disadvantages of these non-image based methods are complexity of test equipment and the need for a wavefront reference.

Dean, Bruce H.↗
DOE OSTI2017

Optimal Bandwidth for Multitaper Spectrum Estimation

A systematic method for bandwidth parameter selection is desired for Thomson multitaper spectrum estimation. We give a method for determining the optimal bandwidth based on a mean squared error (MSE) criterion. When the true spectrum has a second-order Taylor series expansion, one can express quadratic local bias as a function of the curvature of the spectrum, which can be estimated by using a simple spline approximation. This is combined with a variance estimate, obtained by jackknifing over individual spectrum estimates, to produce an estimated MSE for the log spectrum estimate for each choice of time-bandwidth product. The bandwidth that minimizes the estimated MSE then gives the desired spectrum estimate. Additionally, the bandwidth obtained using our method is also optimal for cepstrum estimates. We give an example of a damped oscillatory (Lorentzian) process in which the approximate optimal bandwidth can be written as a function of the damping parameter. Furthermore, the true optimal bandwidth agrees well with that given by minimizing estimated the MSE in these examples.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2010

Estimation of Supersonic Stage Separation Aerodynamics of Winged-Body Launch Vehicles Using Response Surface Methods

Response surface methodology was used to estimate the longitudinal stage separation aerodynamic characteristics of a generic, bimese, winged multi-stage launch vehicle configuration at supersonic speeds in the NASA LaRC Unitary Plan Wind Tunnel. The Mach 3 staging was dominated by shock wave interactions between the orbiter and booster vehicles throughout the relative spatial locations of interest. The inference space was partitioned into several contiguous regions within which the separation aerodynamics were presumed to be well-behaved and estimable using central composite designs capable of fitting full second-order response functions. The underlying aerodynamic response surfaces of the booster vehicle in belly-to-belly proximity to the orbiter vehicle were estimated using piecewise-continuous lower-order polynomial functions. The quality of fit and prediction capabilities of the empirical models were assessed in detail, and the issue of subspace boundary discontinuities was addressed. Augmenting the central composite designs to full third-order using computer-generated D-optimality criteria was evaluated. The usefulness of central composite designs, the subspace sizing, and the practicality of fitting lower-order response functions over a partitioned inference space dominated by highly nonlinear and possibly discontinuous shock-induced aerodynamics are discussed.

Erickson, Gary E.↗
NASA NTRS2006

Estimating Temperature Retrieval Accuracy Associated With Thermal Band Spatial Resolution Requirements for Center Pivot Irrigation Monitoring and Management

This study explores the use of synthetic thermal center pivot irrigation scenes to estimate temperature retrieval accuracy for thermal remote sensed data, such as data acquired from current and proposed Landsat-like thermal systems. Center pivot irrigation is a common practice in the western United States and in other parts of the world where water resources are scarce. Wide-area ET (evapotranspiration) estimates and reliable water management decisions depend on accurate temperature information retrieval from remotely sensed data. Spatial resolution, sensor noise, and the temperature step between a field and its surrounding area impose limits on the ability to retrieve temperature information. Spatial resolution is an interrelationship between GSD (ground sample distance) and a measure of image sharpness, such as edge response or edge slope. Edge response and edge slope are intuitive, and direct measures of spatial resolution are easier to visualize and estimate than the more common Modulation Transfer Function or Point Spread Function. For these reasons, recent data specifications, such as those for the LDCM (Landsat Data Continuity Mission), have used GSD and edge response to specify spatial resolution. For this study, we have defined a 400-800 m diameter center pivot irrigation area with a large 25 K temperature step associated with a 300 K well-watered field surrounded by an infinite 325 K dry area. In this context, we defined the benchmark problem as an easily modeled, highly common stressing case. By parametrically varying GSD (30-240 m) and edge slope, we determined the number of pixels and field area fraction that meet a given temperature accuracy estimate for 400-m, 600-m, and 800-m diameter field sizes. Results of this project will help assess the utility of proposed specifications for the LDCM and other future thermal remote sensing missions and for water resource management.

Ryan, Robert E.↗
DOE OSTI2025

SymbolFit: Automatic Parametric Modeling with Symbolic Regression

We introduce SymbolFit (API: https://github.com/hftsoi/symbolfit), a framework that automates parametric modeling by using symbolic regression to perform a machine-search for functions that fit the data while simultaneously providing uncertainty estimates in a single run. Traditionally, constructing a parametric model to accurately describe binned data has been a manual and iterative process, requiring an adequate functional form to be determined before the fit can be performed. The main challenge arises when the appropriate functional forms cannot be derived from first principles, especially when there is no underlying true closed-form function for the distribution. In this work, we develop a framework that automates and streamlines the process by utilizing symbolic regression, a machine learning technique that explores a vast space of candidate functions without requiring a predefined functional form because the functional form itself is treated as a trainable parameter, making the process far more efficient and effortless than traditional regression methods. We demonstrate the framework in high-energy physics experiments at the CERN Large Hadron Collider (LHC) using five real proton-proton collision datasets from new physics searches, including background modeling in resonance searches for high-mass dijet, trijet, paired-dijet, diphoton, and dimuon events. We show that our framework can flexibly and efficiently generate a wide range of candidate functions that fit a nontrivial distribution well using a simple fit configuration that varies only by random seed, and that the same fit configuration, which defines a vast function space, can also be applied to distributions of different shapes, whereas achieving a comparable result with traditional methods would have required extensive manual effort.

Tsoi, Ho Fung [Univ. of Pennsylvania, Philadelphia↗
DOE OSTI2020

Fermi and Swift Observations of GRB 190114C: Tracing the Evolution of High-energy Emission from Prompt to Afterglow

We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory . The prompt gamma-ray emission was detected by the Fermi GRB Monitor (GBM), the Fermi Large Area Telescope (LAT), and the Swift Burst Alert Telescope (BAT) and the long-lived afterglow emission was subsequently observed by the GBM, LAT, Swift X-ray Telescope (XRT), and Swift UV Optical Telescope. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed by the XRT at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to observe the transition from internal-shock- to external-shock-dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment. We estimate the initial bulk Lorentz factor using the observed high-energy spectral cutoff. Considering the onset of the afterglow component, we constrain the deceleration radius at which this forward shock begins to radiate in order to estimate the maximum synchrotron energy as a function of time. We find that even in the LAT energy range, there exist high-energy photons that are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high-energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy-loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2020

Fermi and Swift Observations of GRB 190114C: Tracing the Evolution of High-energy Emission from Prompt to Afterglow

We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The prompt gamma-ray emission was detected by the Fermi GRB Monitor (GBM), the Fermi Large Area Telescope (LAT), and the Swift Burst Alert Telescope (BAT) and the long-lived afterglow emission was subsequently observed by the GBM, LAT, Swift X-ray Telescope (XRT), and Swift UV Optical Telescope. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed by the XRT at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to observe the transition from internal-shock- to external-shock-dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment. We estimate the initial bulk Lorentz factor using the observed high-energy spectral cutoff. Considering the onset of the afterglow component, we constrain the deceleration radius at which this forward shock begins to radiate in order to estimate the maximum synchrotron energy as a function of time. We find that even in the LAT energy range, there exist high-energy photons that are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high-energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy-loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process.

M. Ajello↗