Search NASA⌕ Search

SEARCH · Search NASA

Results for “statistical inference”

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 343 records · Page 19

A two step linear statistical technique using leaps and bounds procedure for retrieval of geophysical parameters from microwave radiometric data

A linear statistical technique using a 'leaps and bounds' procedure (Furnival and Wilson, 1974) is developed for retrieving geophysical parameters from remote measurements. It is used for retrieving sea surface temperatures from the Scaning Multichannel Microwave Radiometer (SMMR) on Seasat. The technique uses an efficient algorithm to select the best fixed-size subset of the 10 SMMR channels for linearly retrieving a given geophysical parameter. The 5-channel subset (6.6V, 6.6H 10H, 18V, 21H), where V and H refer to, respectively, the vertical and horizontal polarizations and the numbers are the channel frequencies in gigahertz, gives the minimum rms error in estimating the sea surface temperature. A comparison with ground truth indicates that the algorithm infers the temperature with an rms accuracy of better than 1.5 K under most environmental conditions. A quality control procedure which is seen as holding promise for further improving the accuracy is proposed.

Kakar, R. K.↗

A statistical examination of Nimbus-7 SMMR data and remote sensing of sea surface temperature, liquid water content in the atmosphere and surface wind speed

Nimbus 7 Scanning Multichannel Microwave Radiometer (SMMR) brightness temperature measurements over the global oceans have been examined with the help of statistical and empirical techniques. Such analyses show that zonal averages of brightness temperature measured by SMMR over the oceans on a large scale are primarily influenced by the water vapor in the atmosphere. Liquid water in the clouds and rain, which has a much smaller spatial and temporal scale, contributes substantially to the variability of the SMMR measurements within the latitudinal zones. The surface wind not only increases the surface emissivity, but through its interactions with the atmosphere produces correlations in the SMMR brightness temperature data that have significant meteorological implications. It is found that a simple meteorological model can explain the general characteristics of the SMMR data. With the help of this model, methods to infer over the global oceans, the surface temperature, liquid water content in the atmosphere, and surface wind speed are developed. Monthly mean estimates of the sea surface temperature and surface winds are compared with the ship measurements. Estimates of liquid water content in the atmosphere are consistent with earlier satellite measurements. Previously announced in STAR as N83-19187

Prabhakara, C.↗

Maps of Magnetic Field Strength in the OMC-1 Using HAWC+FIR Polarimetric Data

Far-infrared dust polarimetry enables the study of interstellar magnetic fields via tracing of the polarized emission from dust grains that are partially aligned with the direction of the field. The advent of high-quality polarimetric data has permitted the use of statistical methods to extract both the direction and magnitude of the magnetic field. In this work, the Davis–Chandrasekhar–Fermi technique is used to make maps of the plane-of-sky (POS) component of the magnetic field in the Orion Molecular Cloud (OMC-1) by combining polarization maps at 53, 89, 154 and 214 μm from HAWC+/SOFIA with maps of density and velocity dispersion. In addition, maps of the local dispersion of polarization angles are used in conjunction with Zeeman measurements to estimate a map of the strength of the line-of-sight (LOS) component of the field. Combining these maps, information about the threedimensional magnetic field configuration (integrated along the LOS) is inferred over the OMC-1 region. POS magnetic field strengths of up to 2 mG are observed near the BN/KL object, while the OMC-1 bar shows strengths of up to a few hundred μG. These estimates of the magnetic field components are used to produce maps of the mass-to-magnetic-flux ratio (M/Φ)—a metric for probing the conditions for star formation in molecular clouds— and determine regions of sub- and supercriticality in OMC-1. Such maps can provide invaluable input and comparison to MHD simulations of star formation processes in filamentary structures of molecular clouds. Unified Astronomy Thesaurus concepts: Molecular clouds (1072); Giant molecular clouds (653); Interstellar magnetic fields (845); Far infrared astronomy (529)

Jordan A Guerra↗

Wind-Driven Angular Dependence of Sea-Surface Reflectance Measured with an Airborne Doppler Lidar

The effects of wind-stress on the optical properties of the ocean surface have been studied for several decades. In particular, the classic study by Cox and Munk (1954) linking sea-surface wind field to wave slope statistics provides a phenomenology by which the sea-surface wind velocity can be estimated from direct measurement of the wave-modulated surface reflectance. A limited number of studies along these lines have been conducted using airborne or spaceborne lidar systems. In these instances, truthing was provided by in situ ship reports or satellite microwave remote sensing instruments (e.g., ERS scatterometer, SSM/I). During the second deployment of the MACAWS Doppler wind lidar in the summer of 1996 measurements of sea-surface reflectance as a function of azimuth- and nadir-viewing angles were acquired off the California coast. MACAWS data products include directly measured winds, as well as calibrated backscatter/reflectance profiles, thus enabling comparison of the winds inferred from sea-surface reflectance measurements with those deriving from the Doppler-processed direct line-of-sight (LOS) estimates. Additional validation data was extracted from the ERS and SSM/I satellite microwave sensor archives maintained by the JPL Physical Oceanography Distributed Active Archive Center (PO- DAAC).

Tratt, David M.↗

Topics in inference and decision-making with partial knowledge

Two essential elements needed in the process of inference and decision-making are prior probabilities and likelihood functions. When both of these components are known accurately and precisely, the Bayesian approach provides a consistent and coherent solution to the problems of inference and decision-making. In many situations, however, either one or both of the above components may not be known, or at least may not be known precisely. This problem of partial knowledge about prior probabilities and likelihood functions is addressed. There are at least two ways to cope with this lack of precise knowledge: robust methods, and interval-valued methods. First, ways of modeling imprecision and indeterminacies in prior probabilities and likelihood functions are examined; then how imprecision in the above components carries over to the posterior probabilities is examined. Finally, the problem of decision making with imprecise posterior probabilities and the consequences of such actions are addressed. Application areas where the above problems may occur are in statistical pattern recognition problems, for example, the problem of classification of high-dimensional multispectral remote sensing image data.

Safavian, S. Rasoul↗

On the Universality of the Kolmogorov Constant in Numerical Simulations of Turbulence

Motivated by a recent survey of experimental data, we examine data on the Kolmogorov spectrum constant in numerical simulations of isotropic turbulence, using results both from previous studies and from new direct numerical simulations over a range of Reynolds numbers (up to 240 on the Taylor scale) at grid resolutions up to 512(exp 3). It is noted that in addition to k(exp -5/3) scaling, identification of a true inertial range requires spectral isotropy in the same wavenumber range. We found that a plateau in the compensated three-dimensional energy spectrum at k(eta) approx. = 0.1 - -0.2, commonly used to infer the Kolmogorov constant from the compensated three-dimensional energy spectrum, actually does not represent proper inertial range behavior. Rather, a proper, if still approximate, inertial range emerges at k(eta) approx. = 0.02 - 0.05 when R(sub lambda) increases beyond 140. The new simulations indicate proportionality constants C(sub 1) and C in the one- and three-dimensional energy spectra respectively about 0.60 and 1.62. If the turbulence were perfectly isotropic then use of isotropy relations in wavenumber space (C(sub 1) = 18/55 C) would imply that C(sub 1) approx. = 0.53 for C = 1.62, in excellent agreement with experiments. However the one- and three-dimensional estimates are not fully consistent, because of departures (due to numerical and statistical limitations) from isotropy of the computed spectra at low wavenumbers. The inertial scaling of structure functions in physical space is briefly addressed. Since DNS is still restricted to moderate Reynolds numbers, an accurate evaluation of the Kolmogorov constant is very difficult. We focus on providing new insights on the interpretation of Kolmogorov 1941 similarity in the DNS literature and do not consider issues pertaining to the refined similarity hypotheses of Kolmogorov (K62).

Yeung, P. K.↗

Autocorrelation and regularization in digital images. I - Basic theory

Spatial structure occurs in remotely sensed images when the imaged scenes contain discrete objects that are identifiable in that their spectral properties are more homogeneous within than between them and other scene elements. The spatial structure introduced is manifest in statistical measures such as the autocovariance function and variogram associated with the scene, and it is possible to formulate these measures explicitly for scenes composed of simple objects of regular shapes. Digital images result from sensing scenes by an instrument with an associated point spread function (PSF). Since there is averaging over the PSF, the effect, termed regularization, induced in the image data by the instrument will influence the observable autocovariance and variogram functions of the image data. It is shown how the autocovariance or variogram of an image is a composition of the underlying scene covariance convolved with an overlap function, which is itself a convolution of the PSF. The functional form of this relationship provides an analytic basis for scene inference and eventual inversion of scene model parameters from image data.

Jupp, David L. B.↗

Non-Gaussian Distributions of TOA SW Flux as Observed by MISR and CERES

The Top of Atmosphere (TOA) shortwave (SW) flux, converted from Terra Multi angle Imaging‐SpectroRadiometer (MISR) narrow band albedos, is compared with that measured from Clouds and the Earth’s Radiant Energy System (CERES). We describe the probability density function (PDF) of the monthly TOA SW flux and how the statistical third moment, skewness, can impact the quantification of the flux. The PDF of the SW flux is not normally distributed but positively skewed. In both sets of observations, the near-global (80S-80N) median value of the SW flux is ≈3 W/m 2 less than the mean value, due to the positive skewness of the distribution. The near-global mean TOA SW flux converted from MISR is about 7 W/m 2 (≈7%) less than CERES measured flux during the last two decades. Surprisingly, hemispheric asymmetry exists with TOA SW observations from Terra platform. SH reflects 3.92 W/m 2 and 1.15 W/m 2 more mean SW flux than NH, from MISR and CERES Single Scanner Footprint products, respectively. We can infer that the offsetting by morning clouds in the SH is greater than the effect of hemispheric imbalance of SW flux caused by different land masses in two hemispheres. While the characteristics of the two SW fluxes broadly agree with each other, differences in the regional PDF from two different SW fluxes are substantially different over high cloud regions and high altitude regions. Our analysis shows that some parts of the different skewness from two measurements may be attributed to the different calibration of the radiance anisotropy over high cloud scenes.

satellite measurements↗

X-ray studies of quasars and active galaxies with the Einstein Observatory

Preliminary results of the first six months of an X-ray survey of active galaxies and quasars performed with the imaging proportional counter on board the Einstein Observatory are presented. Consideration is given to detections of Seyfert galaxies, N galaxies, BL Lac objects and quasars, and properties of selected representative objects are discussed in detail. The statistical properties, including class averages and correlation coefficients, of the various X-ray active galaxies are calculated in an attempt to determine the relation of X-ray properties to the similarities and differences between the various classes of galaxies. Evidence is obtained for a continuity of high-energy behavior among all active galaxies, and it is noted that they all exhibit flux variability, often on a short time scale. The soft X-ray emission can be explained in terms of the synchrotron self-Compton model, or alternatively by thermal emission from hot gas. A significant correlation is obtained between the X-ray and optical fluxes for all objects considered, indicating that optical quasar counts may be used to infer the contribution of quasars to the diffuse X-ray background.

Ku, W. H.-M.↗

Multiple strike slip faults sets: A case study from the Dead Sea transform

In many strike slip tectonic settings, large rotations of crust blocks about vertical axes have been inferred from paleomagnetic data. These blocks are bounded by sets of parallel faults which presumably accommodate the relative motion between the blocks as regional deformation progress. A mechanical model by Nur et al., (1986) suggests that rotations greater than phi sub c equals 25 to 45 degrees must be accommodated by more than one set of faults, with angle phi sub c between their direction; consequently the sum of the angles between sets must be roughly equal to the total tectonic material rotation. To test this model, the authors investigated the fault geometry and field relation of fault sets in the Mt. Hermon area in northern Israel, where paleomagnetic declination implies data 69 degrees plus or minus 13 degrees counter-clockwise block rotation. The statistical and field relation analysis of over 315 faults shows that the faulting is predominantly right lateral strike slip consisting of three distinct sets. The oldest set strikes 253 degrees, the second oldest set strikes 293 degrees and the youngest strikes 339 degrees. This last direction is consistent also with the current north-south direction of the maximum principle stress axis. The angle phi sub c between the first and second sets is 39 degrees and between the second and third sets 46 degrees, in good agreement with the phi sub c angle predicted from mechanical considerations. The sum of the two angles is 85 degrees, in good agreement with the 69 degrees plus or minus 13 degrees CCW paleomagnetically derived rotation. The results suggest specifically that the sequential development of multiple intersecting fault sets is responsible for the faulting in the Mt. Hermon area; and generally that the model of block rotation with multiple faults provides very good simple rules for analyzing very complex fault patterns.

Ron, Hagai↗

Validation of TRMM Satellite Rainfall Products Over Oklahoma for a Three Year Period (1998-2000)

This study evaluates the level 2 and level 3 rainfall products from the Tropical Rainfall Measuring Mission gridded distribution of 118 high resolution tipping buckets from the Oklahoma Mesonet. The Tropical Rainfall Measuring Mission, a joint satellite mission between NASA and the National Space Development Agency (NASDA) of Japan, was designed to estimate global precipitation between 40 S and 40 N latitude. The TRMM satellite, consisting of two main precipitation sensors, a passive microwave (TMI) and precipitation radar (PR) sensors, was launched in 1997. Although the great advantage of TRMM has been its ability to sample precipitation globally over the tropical oceans in places where ground sensors do not exist, regional points over land still offer the best opportunity for validating these rain estimates. For this study, the gauge data was gridded at I x I degree resolution between 34 and 36 N and -95 and -100 W. The location of Oklahoma is somewhat unique, in that, it is located near the turning point for the satellite. This study investigates effects of temporal sampling on the satellite measurements and the resulting rainfall bias observed between space and ground based sensors. The first part of the study analyzes the monthly rainfall statistics of the 3AI2 (TMI), the 3A25 (PR) and the 31142 (combined sensors) with the rain gauges. It uses regression techniques and a computation of the probability density function (PDF) and cumulative probability density function (CDF) for each sensor to evaluate the effects of temporal sampling on the satellite products and the sufficiency of using two overpasses a day to estimate rainfall on a monthly scale. In further probing this issue, the study also looked at the direct coincidence of ground and satellite measurements by comparing instantaneous comparisons obtained from the level 2 orbital track data, e.g. 2A12, 2A25. Lastly, the effects of temporal sampling were further studied by sub-sampling the gauge data only at overpass times. This methodology provides a quantitative way of separately inferring the fraction of bias due to temporal sampling and the rainfall algorithm. In addition to showing the salient effects of temporal sampling, the study also revisits the historical problem of correlating point gauge estimates with area estimates, though, the problem is ameliorated by averaging and gridding the data to the area resolution scale of the satellite.

Fisher, Brad↗

Evaluation of Long-Term Cloud-Resolving Model Simulations Using Satellite Radiance Observations and Multi-Frequency Satellite Simulators

This paper proposes a methodology known as the Tropical Rainfall Measuring Mission (TRMM) Triple-Sensor Three-step Evaluation Framework (T3EF) for the systematic evaluation of precipitating cloud types and microphysics in a cloud-resolving model (CRM). T3EF utilizes multi-frequency satellite simulators and novel statistics of multi-frequency radiance and backscattering signals observed from the TRMM satellite. Specifically, T3EF compares CRM and satellite observations in the form of combined probability distributions of precipitation radar (PR) reflectivity, polarization-corrected microwave brightness temperature (Tb), and infrared Tb to evaluate the candidate CRM. T3EF is used to evaluate the Goddard Cumulus Ensemble (GCE) model for cases involving the South China Sea Monsoon Experiment (SCSMEX) and Kwajalein Experiment (KWAJEX). This evaluation reveals that the GCE properly captures the satellite-measured frequencies of different precipitating cloud types in the SCSMEX case but underestimates the frequencies of deep convective and deep stratiform types in the KWAJEX case. Moreover, the GCE tends to simulate excessively large and abundant frozen condensates in deep convective clouds as inferred from the overestimated GCE-simulated radar reflectivities and microwave Tb depressions. Unveiling the detailed errors in the GCE s performance provides the best direction for model improvements.

Matsui, Toshihisa↗

The analysis of a rocket tomography measurement of the N2+3914A emission and N2 ionization rates in an auroral arc

Techniques were developed for recovering two-dimensional distributions of auroral volume emission rates from rocket photometer measurements made in a tomographic spin scan mode. These tomographic inversion procedures are based upon an algebraic reconstruction technique (ART) and utilize two different iterative relaxation techniques for solving the problems associated with noise in the observational data. One of the inversion algorithms is based upon a least squares method and the other on a maximum probability approach. The performance of the inversion algorithms, and the limitations of the rocket tomography technique, were critically assessed using various factors such as (1) statistical and non-statistical noise in the observational data, (2) rocket penetration of the auroral form, (3) background sources of emission, (4) smearing due to the photometer field of view, and (5) temporal variations in the auroral form. These tests show that the inversion procedures may be successfully applied to rocket observations made in medium intensity aurora with standard rocket photometer instruments. The inversion procedures have been used to recover two-dimensional distributions of auroral emission rates and ionization rates from an existing set of N2+3914A rocket photometer measurements which were made in a tomographic spin scan mode during the ARIES auroral campaign. The two-dimensional distributions of the 3914A volume emission rates recoverd from the inversion of the rocket data compare very well with the distributions that were inferred from ground-based measurements using triangulation-tomography techniques and the N2 ionization rates derived from the rocket tomography results are in very good agreement with the in situ particle measurements that were made during the flight. Three pre-prints describing the tomographic inversion techniques and the tomographic analysis of the ARIES rocket data are included as appendices.

Mcdade, Ian C.↗

New Approaches to Quantifying Transport Model Error in Atmospheric CO2 Simulations

In recent years, much progress has been made in observing CO2 distributions from space. However, the use of these observations to infer source/sink distributions in inversion studies continues to be complicated by difficulty in quantifying atmospheric transport model errors. We will present results from several different experiments designed to quantify different aspects of transport error using the Goddard Earth Observing System, Version 5 (GEOS-5) Atmospheric General Circulation Model (AGCM). In the first set of experiments, an ensemble of simulations is constructed using perturbations to parameters in the model s moist physics and turbulence parameterizations that control sub-grid scale transport of trace gases. Analysis of the ensemble spread and scales of temporal and spatial variability among the simulations allows insight into how parameterized, small-scale transport processes influence simulated CO2 distributions. In the second set of experiments, atmospheric tracers representing model error are constructed using observation minus analysis statistics from NASA's Modern-Era Retrospective Analysis for Research and Applications (MERRA). The goal of these simulations is to understand how errors in large scale dynamics are distributed, and how they propagate in space and time, affecting trace gas distributions. These simulations will also be compared to results from NASA's Carbon Monitoring System Flux Pilot Project that quantified the impact of uncertainty in satellite constrained CO2 flux estimates on atmospheric mixing ratios to assess the major factors governing uncertainty in global and regional trace gas distributions.

Ott, L.↗

The Cosmic Shoreline

Volatile escape is the classic existential problem of planetary atmospheres. The problem has gained new currency now that we can study the cumulative effects of escape from extrasolar planets. Escape itself is likely to be a rapid process, relatively unlikely to be caught in the act, but the cumulative effects of escape in particular, the distinction between planets with and without atmospheres should show up in the statistics of the new planets. The new planets make a moving target. It can be difficult to keep up, and every day the paper boy brings more. Of course most of these will be giant planets loosely resembling Saturn or Neptune albeit hotter and nearer their stars, as big hot fast-orbiting exoplanets are the least exceedingly difficult to discover. But they are still planets, all in all, and although twenty years ago experts could prove on general principles that they did not exist, we have come round rather quickly, and they should be welcome now at LPSC. Here we will discuss the empirical division between planets with and without atmospheres. For most exoplanets the question of whether a planet has or has not an atmosphere is a fuzzy inference based on the planet's bulk density. A probably safe presumption is that a low density planet is one with abundant volatiles, in the general mold of Saturn or Neptune. On the other hand a high density low mass planet could be volatile-poor, in the general mold of Earth or Mercury. We will focus on planets, mostly seen in transit, for which both radius and mass are measured, as these are the planets with measured densities. More could be said: a lot of subtle recent work has been devoted to determining the composition of planets from equations of state or directly observing atmospheres in transit, but we will not go there. What interests us here is that, from the first, the transiting extrasolar planets appear to have fit into a pattern already seen in our own Solar System, as shown in Fig. 1. We first noticed this in 2004 when there were just two transiting exoplanets to consider. The trend was well-defined by late 2007. Figure 1 shows how matters stood in Dec 2012 with approx.240 exoplanets. The figure shows that the boundary between planets with and without active volatiles - the cosmic shoreline, as it were - is both well-defined and follows a power law.

Zahnle, Kevin J.↗

Sensitivity of the Assimilated Ozone in the UTLS to Model and Data Selection Changes

This presentation will discuss the sensitivity of assimilated ozone fields in the upper troposphere and lower stratosphere (UTLS) to a number of factors, focusing mainly on aspects of data selection and the prediction model. This is important, because assimilation represents an attempt to construct our best estimates of the true ozone field; however, inaccuracies in the UTLS ozone distribution translate into an uncertainty in factors such as the calculated radiative forcing of climate or the inferred stratosphere-troposphere exchange (STE) of ozone. The 3D ozone data assimilation system, from NASA's Global Modeling and Assimilation Office (GMAO), combines observations of total ozone column and stratospheric profiles with predictions from an off-line, parameterized chemistry and transport model (pCTM) to produce six-hourly, global analyses. The first experiments discussed assimilate ozone retrievals from the Earth-Probe Total Ozone Mapping Spectrometer (EPTOMS) and stratospheric profiles from the Solar Backscatter UltraViolet/2 (SBUV/2) instrument. The SBUV/2 ozone data have a coarse vertical resolution, with increased uncertainty below the ozone maximum, and TOMS provides only total ozone columns. Thus, the assimilated ozone profiles in the UTLS region are only weakly constrained by the incoming SBUV and TOMS data. Consequently, the assimilated ozone distribution should be sensitive to changes in inputs to the statistical analysis scheme. Sensitivity studies have been conducted to examine the responses to TOMS and SBUV/2 data selection, modifications of the forecast and observation error covariance models, and the model formulation (turning off chemistry or using different wind analyses in the pCTM). The second set of experiments includes an additional data type: ozone retrieved from infrared limb-emission by MIPAS on Envisat. These data offer not only improved vertical resolution in the stratosphere, but also give measurements in the polar night. Comparisons of the assimilated ozone fields from both sets of experiments with independent observations, primarily ozone sondes, are used to determine the impact of each of these changes. It is shown that many of the changes have a significant impact on the UTLS ozone estimates. Implications for interpretation of STE and radiative forcing of climate are discussed.

Pawson, Steven↗

Using Machine Learning to Infer Material Properties of Debris Fragments from X-ray Images in the DebriSat Project

The DebriSat project is a collaboration effort with the NASA Orbital Debris Program Office, the U.S. Space Force Space Systems Command Center, The Aerospace Corporation, and the University of Florida. To date, over 200,000 fragments from this ground-based, hypervelocity impact experiment have been collected, and processing is underway to determine their physical characteristics, such as material, shape, color, characteristic length, and average cross-sectional area. The x-ray process is primarily used to identify the location of the fragments and estimated size for extraction, so that these physical characteristics can be assessed. This paper proposes a machine learning-based approach to characterize materials from x-ray images of debris fragments embedded in soft-catch foam used in the DebriSat project. The novel methodology discussed in this paper will highlight the use of x-ray imagery data to characterize these fragments without extraction or a human-in-the-loop. Both supervised and unsupervised machine learning techniques are utilized with this approach to infer the physical parameters of the fragments embedded in the soft-catch foam panels used in the impact experiment based on x-ray images of the foam panels. Additionally, 3D reconstructions of the extracted fragments are created with images taken from two different angles using the structure from motion (SfM) method. The characteristic lengths and shape from the 3D reconstruction, alongside the physical characteristics of the debris, are used in the inference of the material type. To develop and test the approach, a dataset of x-ray images of debris fragments of varying sizes and materials is collected. Supervised learning methods such as convolutional neural networks (CNNs), support vector machines (SVM), decision trees, and random forest classifiers are used due to the high-dimensional feature spaces of the debris and nonlinear decision boundaries for material categorization. Given the limited pre-labeled data of embedded debris materials smaller than 10 mm, unsupervised machine learning techniques such as clustering algorithms and autoencoders are used, in addition to supervised learning methods. The clustering algorithms group similar fragments together based on their physical properties, and autoencoders reduce the dimensionality of the x ray images and extract relevant features. The performance of the proposed approach's is analyzed using a range of statistical methods, including confusion matrices, receiver operating characteristic curves, and precision-recall curves. The results are compared with those obtained using a baseline approach that relies on manual identification and classification of debris fragments. To evaluate the effectiveness of different machine learning methods, statistical tests such as t-tests, ANOVA, and cross-validation are performed, comparing the performance of CNNs, SVMs, clustering algorithms, and autoencoders. Additional analysis needs to be conducted to identify any sources of bias or variability that may affect the results, such as variations in imaging conditions or fragmentation patterns. Other topics explored are limitations, refinements, and the potential use of semi-supervised learning techniques, such as self-training to label unlabeled datasets and co-training using x-ray images taken from two different angles as two different models.

Saik Anam Siam↗

Using Machine Learning to Infer Material Properties of Debris Fragments from X-ray Images in the DebriSat Project

The DebriSat project is a collaboration effort with the NASA Orbital Debris Program Office, the U.S. Space Force Space Systems Command Center, The Aerospace Corporation, and the University of Florida. To date, over 200,000 fragments from this ground-based, hypervelocity impact experiment have been collected, and processing is underway to determine their physical characteristics, such as material, shape, color, characteristic length, and average cross-sectional area. The x-ray process is primarily used to identify the location of the fragments and estimated size for extraction, so that these physical characteristics can be assessed. This paper proposes a machine learning-based approach to characterize materials from x-ray images of debris fragments embedded in soft-catch foam used in the DebriSat project. The novel methodology discussed in this paper will highlight the use of x-ray imagery data to characterize these fragments without extraction or a human-in-the-loop. Both supervised and unsupervised machine learning techniques are utilized with this approach to infer the physical parameters of the fragments embedded in the soft-catch foam panels used in the impact experiment based on x-ray images of the foam panels. Additionally, 3D reconstructions of the extracted fragments are created with images taken from two different angles using the structure from motion (SfM) method. The characteristic lengths and shape from the 3D reconstruction, alongside the physical characteristics of the debris, are used in the inference of the material type. To develop and test the approach, a dataset of x-ray images of debris fragments of varying sizes and materials is collected. Supervised learning methods such as convolutional neural networks (CNNs), support vector machines (SVM), decision trees, and random forest classifiers are used due to the high-dimensional feature spaces of the debris and nonlinear decision boundaries for material categorization. Given the limited pre-labeled data of embedded debris materials smaller than 10 mm, unsupervised machine learning techniques such as clustering algorithms and autoencoders are used, in addition to supervised learning methods. The clustering algorithms group similar fragments together based on their physical properties, and autoencoders reduce the dimensionality of the x ray images and extract relevant features. The performance of the proposed approach's is analyzed using a range of statistical methods, including confusion matrices, receiver operating characteristic curves, and precision-recall curves. The results are compared with those obtained using a baseline approach that relies on manual identification and classification of debris fragments. To evaluate the effectiveness of different machine learning methods, statistical tests such as t-tests, ANOVA, and cross-validation are performed, comparing the performance of CNNs, SVMs, clustering algorithms, and autoencoders. Additional analysis needs to be conducted to identify any sources of bias or variability that may affect the results, such as variations in imaging conditions or fragmentation patterns. Other topics explored are limitations, refinements, and the potential use of semi-supervised learning techniques, such as self-training to label unlabeled datasets and co-training using x-ray images taken from two different angles as two different models.

Saik Anam Siam↗