Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical change detection”

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 109 records · Page 6

Remote Sensing Time Series Product Tool

The TSPT (Time Series Product Tool) software was custom-designed for NASA to rapidly create and display single-band and band-combination time series, such as NDVI (Normalized Difference Vegetation Index) images, for wide-area crop surveillance and for other time-critical applications. The TSPT, developed in MATLAB, allows users to create and display various MODIS (Moderate Resolution Imaging Spectroradiometer) or simulated VIIRS (Visible/Infrared Imager Radiometer Suite) products as single images, as time series plots at a selected location, or as temporally processed image videos. Manually creating these types of products is extremely labor intensive; however, the TSPT development tool makes the process simplified and efficient. MODIS is ideal for monitoring large crop areas because of its wide swath (2330 km), its relatively small ground sample distance (250 m), and its high temporal revisit time (twice daily). Furthermore, because MODIS imagery is acquired daily, rapid changes in vegetative health can potentially be detected. The new TSPT technology provides users with the ability to temporally process high-revisit-rate satellite imagery, such as that acquired from MODIS and from its successor, the VIIRS. The TSPT features the important capability of fusing data from both MODIS instruments onboard the Terra and Aqua satellites, which drastically improves cloud statistics. With the TSPT, MODIS metadata is used to find and optionally remove bad and suspect data. Noise removal and temporal processing techniques allow users to create low-noise time series plots and image videos and to select settings and thresholds that tailor particular output products. The TSPT GUI (graphical user interface) provides an interactive environment for crafting what-if scenarios by enabling a user to repeat product generation using different settings and thresholds. The TSPT Application Programming Interface provides more fine-tuned control of product generation, allowing experienced programmers to bypass the GUI and to create more user-specific output products, such as comparison time plots or images. This type of time series analysis tool for remotely sensed imagery could be the basis of a large-area vegetation surveillance system. The TSPT has been used to generate NDVI time series over growing seasons in California and Argentina and for hurricane events, such as Hurricane Katrina.

Predos, Don↗

Acoustic space occupancy: Combining ecoacoustics and lidar to model biodiversity variation and detection bias across heterogeneous landscapes

There is global interest in quantifying changing biodiversity in human-modified landscapes. Ecoacoustics may offer a promising pathway for supporting multi-taxa monitoring, but its scalability has been hampered by the sonic complexity of biodiverse ecosystems and the imperfect detectability of animal-generated sounds. The acoustic signature of a habitat, or soundscape, contains information about multiple taxa and may circumvent species identification, but robust statistical technology for characterizing community-level attributes is lacking. Here, we present the Acoustic Space Occupancy Model, a flexible hierarchical framework designed to account for detection artifacts from acoustic surveys in order to model biologically relevant variation in acoustic space use among community assemblages. We illustrate its utility in a biologically and structurally diverse Amazon frontier forest landscape, a valuable test case for modeling biodiversity variation and acoustic attenuation from vegetation density. We use complementary airborne lidar data to capture aspects of 3D forest structure hypothesized to influence community composition and acoustic signal detection. Our novel analytic framework permitted us to model both the assembly and detectability of soundscapes using lidar-derived estimates of forest structure. Our empirical predictions were consistent with physical models of frequency-dependent attenuation, and we estimated that the probability of observing animal activity in the frequency channel most vulnerable to acoustic attenuation varied by over 60%, depending on vegetation density. There were also large differences in the biotic use of acoustic space predicted for intact and degraded forest habitats, with notable differences in the soundscape channels predominantly occupied by insects. This study advances the utility of ecoacoustics by providing a robust modeling framework for addressing detection bias from remote audio surveys while preserving the rich dimensionality of soundscape data, which may be critical for inferring biological patterns pertinent to multiple taxonomic groups in the tropics. Our methodology paves the way for greater integration of remotely sensed observations with high-throughput biodiversity data to help bring routine, multi-taxa monitoring to scale in dynamic and diverse landscapes.

Airborne lidar↗

The Ocean Colour Climate Change Initiative: I. A Methodology for Assessing Atmospheric Correction Processors Based on In-Situ Measurements

The Ocean Colour Climate Change Initiative intends to provide a long-term time series of ocean colour data and investigate the detectable climate impact. A reliable and stable atmospheric correction procedure is the basis for ocean colour products of the necessary high quality. In order to guarantee an objective selection from a set of four atmospheric correction processors, the common validation strategy of comparisons between in-situ and satellite derived water leaving reflectance spectra, is extended by a ranking system. In principle, the statistical parameters such as root mean square error, bias, etc. and measures of goodness of fit, are transformed into relative scores, which evaluate the relationship of quality dependent on the algorithms under study. The sensitivity of these scores to the selected database has been assessed by a bootstrapping exercise, which allows identification of the uncertainty in the scoring results. Although the presented methodology is intended to be used in an algorithm selection process, this paper focusses on the scope of the methodology rather than the properties of the individual processors.

Climate Change↗

Recent Results from the Goldstone Orbital Debris Radar: 2016-2017

Since 1993, the NASA Orbital Debris Program Office has used the Goldstone Orbital Debris Radar (Goldstone) to sample statistically the orbital debris environment. Due to the sensitivity of this radar, which can detect an approximately 3 mm-diameter conducting sphere at 1,000 km, it has filled an important role in the characterization of the sub-centimeter-sized orbital debris population. Through the years, the capabilities of this system have increased – recent updates include increased receiver bandwidth and a change in the bi-static observation geometry – both of which enhance the radar’s ability to estimate orbital parameters. In 2016, dual polarization capability was added, making this the first year where both right- and left-hand circularly polarized information was available from this sensor. This additional polarization information may enable better characterization of sub-centimeter-sized particles in low Earth orbit, particularly since the receiver triggers on reflected energy from both left- and right-handed circular polarizations independently. In this paper, we present measurements and results derived from data taken during the calendar years (CY) 2016-2017 by Goldstone and compare this dataset to measurements taken by the Haystack Ultra-wideband Satellite Imaging Radar (HUSIR) during a similar timeframe.

Murray, James↗

Detecting Risk and Anomalies in Airplane Dynamics Through Entropic Analysis of Time Series Data

Despite recent efforts to move away from traditional threshold exceedance detection methods for aircraft state monitoring, modern aircraft still rely on safety thresholds to communicate to pilots the identification of an anomaly in the aircraft when a threshold is surpassed. Current anomaly detection methods mainly depend on uninterpretable machine learning models to learn complex patterns and relationships contained in the time series data of aircraft. Although these methods are capable of identifying known anomalies, their deficiency in interpretability presents a challenge when translating them to different aircraft. To overcome this deficiency, entropic analysis of aircraft dynamics seeks to characterize the complexity, or lack thereof, of the aircraft dynamics prior to the development of a risk scenario. This complexity characterization provides a more straightforward summary of state changes in the dynamics of flight variables. To build a foundation for entropic analysis, we analyzed the complexity of unstable approaches, an anomalous event present in many of today’s aviation accidents. The analysis revealed a statistically significant difference in the complexity distribution of flight variables under a stable approach versus an unstable approach. These differences in complexity were especially notable minutes before an approach was identified as unstable. Moreover, the multiscale entropic analysis revealed the presence of signal complexity at multiple time scales across multiple time windows before landing. By capturing state changes and corrections in the aircraft dynamics using entropy, advanced, yet still interpretable, sensor systems based on entropic frameworks from this study can be constructed in the future using classical machine learning approaches.

Risk detection↗

Pulse profile and refined orbital elements for SMC X-1

Results of a reanalysis of SAS-3 data for SMC X-1, including greatly improved orbital elements, are presented. These elements, when combined with newly available data for the optical companion, Sk 160, yield improved estimates for the masses of SMC X-1 and Sk 160. A significant rate of change of intrinsic pulse period is detected during the course of the four-day observation which is consistent with the average value that was previously deduced from observations separated by 5 years. The X-ray pulse profile, averaged over one orbital cycle and with approximately 15-ms time resolution, is also presented. The pulse profile is found to be constant to within statistics, as a function of orbital phase. The implications of the small observed eccentricities of SMC X-1 and other binary X-ray pulsar systems are discussed.

Primini, F.↗

The Resource Prospector Neutron Spectrometer System: RP's Bloodhound

The primary goal of the Resource Prospector Neutron Spectrometer System (NSS) is to locate and characterize hydrogen-bearing volatile deposits, especially subsurface ice, that may exist at the lunar poles. A key objective is to detect water-equivalent hydrogen concentrations of 0.5 wt% or greater while on a moving rover. A second objective is to determine approximate burial depth of enhanced hydrogen-bearing materials up to 1 meter below otherwise dry regolith. The instrument will be carried aboard a landed mobility system at the lunar poles. The instrument operates by measuring the changes in the leakage flux of low energy neutrons out of the regolith. These neutrons are produced by galactic cosmic rays, which are so energetic that they shatter the nuclei in surface materials. The neutrons interact with other nuclei and lose energy, becoming thermalized in the process. Hydrogen is most efficient at thermalizing neutrons owing to protons' similar mass - statistically, neutrons lose half their energy per collision with protons. With hydrogen in the soil, leakage fluxes of neutrons in the 0.5 eV to 500 keV energy range are reduced. A concentration of only1-2 wt% water-equivalent hydrogen results in a decrease in epithermal leakage flux of a factor of two. The leakage flux of thermal neutrons, from 0 to 0.5 eV in energy, can either increase or decrease depending on the hydrogen abundance and stratigraphy. As with the highly successful Lunar Prospector Neutron Spectrometer, the RP NSS detects both thermal and epithermal neutrons by using two helium-3 gas proportional counters, one covered by cadmium and the other uncovered. The former detects only epithermal neutrons with energies above approximately 0.5 eV, the latter detects both thermal (less than 0.5 eV) and epithermal energies (greater than 0.5 eV). When a neutron enters the detector tube and interacts with a helium-3 nucleus, the resulting reaction produces an energetic proton and triton that ionize the gas. The resulting electrons are accelerated toward a high-voltage anode and cascade, amplifying the net charge, which is collected at the anode. The number of electrons produced is proportional to the energy that the triton and proton deposit in the gas. A charge sensitive pre-amplifier converts the total charge to a step voltage output. A shaper amplifier then shapes this step into a uni-polar waveform with peaking time appropriate for the detection depending on the event rate. The integrated shaped waveform, representing the deposited triton/proton energy, is then measured. A histogram, or pulse height analysis, is performed to record the main capture peak and wall effect pulses. A threshold for detection is also required to limit the low amplitude counting rate such as noise floor. The system electronics consists of 2 modules - the Sensor Module (SM) front-end and the Data Processing Module (DPM) back-end circuits. SM is designed as a light-weight and low power front-end housing the two helium-3 proportional counter detectors, preamp and HVPS. It is mounted external to the rover body to detect the thermalized neutron flux with a minimum of host background. The DPM is located inside the rover; it digitizes the SM signals, performs pulse height analysis and accumulates the count rate for each spectral channel. The DPM controls high voltage and thresholding, and sends the science data to the host craft via an RS422 serial asynchronous protocol. The payload host provides all thermal management and control for the SM and DPM.

hydrogen-bearing volatile deposits↗

Analysis of Pulsed Airborne Lidar Measurements of Atmospheric CO2 Column Absorption from 3-13 km Altitudes

We have developed a pulsed lidar technique for measuring the tropospheric CO2 concentrations as a candidate for NASA's ASCENDS space mission [1]. It uses two pulsed laser transmitters allowing simultaneous measurement of a CO2 absorption line in the 1575 nm band, O2 extinction in the Oxygen A-band, surface height and backscatter profile. The lasers are precisely stepped in wavelength across the CO2 line and an O2 line region during the measurement. The direct detection receiver measures the energies of the laser echoes from the surface along with the range profile of scattering in the path. The column densities for the CO2 and O2 gases are estimated from the ratio of the on- and off-line signals via the integrated path differential absorption (IPDA) technique. The time of flight of the laser pulses is used to estimate the height of the scattering surface and to reject laser photons scattered in the atmosphere. We developed an airborne lidar to demonstrate an early version of the CO2 measurement from the NASA Glenn Lear-25 aircraft. The airborne lidar stepped the pulsed laser's wavelength across the selected CO2 line with 20 wavelength steps per scan. The line scan rate is 450 Hz, the laser pulse widths are 1 usec, and laser pulse energy is 24 uJ. The time resolved laser backscatter is collected by a 20 cm telescope, detected by a NIR photomultiplier and is recorded on every other reading by a photon counting system [2]. During August 2009 we made a series of 2.5 hour long flights and measured the atmospheric CO2 absorption and line shapes using the 1572.33 nm CO2 line. Measurements were made at stepped altitudes from 3-13 km over locations in the US, including the SGP ARM site in Oklahoma, central Illinois, north-eastern North Carolina, and over the Chesapeake Bay and the eastern shore of Virginia. Although the received signal energies were weaker than expected for ASCENDS, clear CO2 line shapes were observed at all altitudes, and some measurements were made through thin clouds. The Oklahoma and east coast flights were coordinated with a LaRC/ITT CO2 lidar on the LaRC UC-12 aircraft, and in-situ measurements were made using its CO2 sensor and radiosondes. We have conducted an analysis of the ranging and IPDA lidar measurements from these four flights. Most flights had 5-6 altitude steps with 200-300 seconds of recorded measurements per step. We used a cross-correlation approach to process the laser echo records. This was used to estimate the range to the scattering surface, to define the edges of the laser pulses and to determine echo pulse energy at each wavelength. We used a minimum mean square approach to fit an instrument response function and to solve for the best-fit CO2 absorption line shape. We then calculated the differential optical depth (DOD) of the fitted CO2 line. We computed its statistics at the various altitude steps, and compare them to the DODs calculated from spectroscopy based on HITRAN 2008 and the column conditions calculated from the airborne in-situ readings. The results show the lidar and in-situ measurements have very similar DOD change with altitude and greater than 10 segments per flight where the scatter in the lidar measurements are less than or equal to 1ppm. We also present the results from subsequent CO2 column absorption measurements, which were made with stronger detected signals during three flights on the NASA DC-8 over the southwestern US in during July 2010.

Abshire, James B.↗

Effects of orbital spaceflight on human osteoblastic cell physiology and gene expression

During long-term spaceflight, astronauts lose bone, in part due to a reduction in bone formation. It is not clear, however, whether the force imparted by gravity has direct effects on bone cells. To examine the response of bone forming cells to weightlessness, human fetal osteoblastic (hFOB) cells were cultured during the 17 day STS-80 space shuttle mission. Fractions of conditioned media were collected during flight and shortly after landing for analyses of glucose utilization and accumulation of type I collagen and prostaglandin E(2) (PGE(2)). Total cellular RNA was isolated from flight and ground control cultures after landing. Measurement of glucose levels in conditioned media indicated that glucose utilization occurred at a similar rate in flight and ground control cultures. Furthermore, the levels of type I collagen and PGE(2) accumulation in the flight and control conditioned media were indistinguishable. The steady-state levels of osteonectin, alkaline phosphatase, and osteocalcin messenger RNA (mRNA) were not significantly changed following spaceflight. Gene-specific reductions in mRNA levels for cytokines and skeletal growth factors were detected in the flight cultures using RNase protection assays. Steady-state mRNA levels for interleukin (IL)-1alpha and IL-6 were decreased 8 h following the flight and returned to control levels at 24 h postflight. Also, transforming growth factor (TGF)-beta(2) and TGF-beta(1) message levels were modestly reduced at 8 h and 24 h postflight, although the change was not statistically significant at 8 h. These data suggest that spaceflight did not significantly affect hFOB cell proliferation, expression of type I collagen, or PGE(2) production, further suggesting that the removal of osteoblastic cells from the context of the bone tissue results in a reduced ability to respond to weightlessness. However, spaceflight followed by return to earth significantly impacted the expression of cytokines and skeletal growth factors, which have been implicated as mediators of the bone remodeling cycle. It is not yet clear whether these latter changes were due to weightlessness or to the transient increase in loading resulting from reentry.

Non-NASA Center↗

Fatigue Crack and Porosity Measurement in Composite Materials by Thermographic and Ultrasonic Methods

Many nondestructive methods exist for the detection of localized material anomalies in an otherwise good composite structure. The problem arises when the material system as a whole has degraded during service or was improperly manufactured. Porosity and intra-ply microcracking are two such conditions that in unlined composite pressure vessels can be very troublesome to detect and when linked through the thickness can be critical to mission success. These leak paths may lead to loss of pressure/propellant, increased risk of explosion and possible cryo-pumping. Research sought nondestructive methods for quantifying porosity and microcracking in composite tankage. Both thermographic and resonance ultrasound methods have been utilized with artificial neural network and statistical approaches to analyze the data. Resonant ultrasound spectroscopy provides measurements, which are sensitive to fine details in the materials character, such as micro-cracking and porosity. Here, the higher frequency (shorter wavelength) components of the signal train provide more significant interaction with the defects causing the spectral characteristics to shift toward lower amplitudes at the higher frequencies. As the density of the defects increases more interactions occur and more drastic amplitude changes are observed. From a thermal perspective, the higher the defect density the lower the through thickness thermal diffusivity will be. Utilizing a point heat source, and thermographically recording the heat profile with time, diffusivity calculations can be made which in turn can be related to the relative quality of the material. Preliminary experiments to verify the measurable effect on the resonance spectrum of the ultrasonic data to detect microcracking and for porosity detection thermographically are presented. Methods involving supervised and unsupervised artificial neural networks as well as other clustering algorithms are developed for signal identification.

Walker, James L.↗

Identification of interplanetary tangential and rotational discontinuities.

Discontinuities in the interplanetary magnetic field detected by the Mariner 5 magnetometer have been analyzed to distinguish tangential and rotational discontinuities. The method used was based on finding the field component normal to the plane of the discontinuities that vanishes for tangential discontinuities and is nonzero for rotational discontinuities. Slightly less than one half of the 118 discontinuities analyzed could not be classified because neither the normal component nor the changes in field magnitude were statistically significant owing to errors in the field measurements and in finding the normal direction. The 44 rotational discontinuities that were finally identified did not exhibit changes in field magnitude in excess of 20%.-

Smith, E. J.↗

A practical obstacle detection system for the Mars Rover

Discussion of an obstacle detection system which uses a laser range finder as the principal sensing device. The laser scans the scene ahead and stores the range data for the terrain in a matrix form. This matrix can be thought of as a 'range image'. This range varies from point to point, but sharp changes are caused by the presence of obstacles. Some existing algorithms and one developed by the authors were used to obtain the outlines of obstacles from the 'range image'. All these algorithms end with a thresholding operation. A theoretical analysis shows how a proper value of the threshold may be chosen, given the type of irregularity of the terrain to be detected. The analysis also shows how large the statistics of measurement noise can be allowed without leading to false alarms.

Reed, M.↗

A study of the feasibility of statistical analysis of airport performance simulation

The feasibility of conducting a statistical analysis of simulation experiments to study airport capacity is investigated. First, the form of the distribution of airport capacity is studied. Since the distribution is non-Gaussian, it is important to determine the effect of this distribution on standard analysis of variance techniques and power calculations. Next, power computations are made in order to determine how economic simulation experiments would be if they are designed to detect capacity changes from condition to condition. Many of the conclusions drawn are results of Monte-Carlo techniques.

Myers, R. H.↗

Volcanic eruptions and solar activity

The historical record of large volcanic eruptions from 1500 to 1980 is subjected to detailed time series analysis. In two weak but probably statistically significant periodicities of about 11 and 80 yr, the frequency of volcanic eruptions increases (decreases) slightly around the times of solar minimum (maximum). Time series analysis of the volcanogenic acidities in a deep ice core from Greenland reveals several very long periods ranging from about 80 to about 350 yr which are similar to the very slow solar cycles previously detected in auroral and C-14 records. Solar flares may cause changes in atmospheric circulation patterns that abruptly alter the earth's spin. The resulting jolt probably triggers small earthquakes which affect volcanism.

Stothers, Richard B.↗

Searching for hard X-ray directivity during the rise, peak, and decay phases of solar flares

We have identified 72 large solar flares (peak counting rates more than 1000 counts/s) observed by Hard X-ray Burst Spectroscopy (HXRBS) on-board the Solar Maximum Mission (SMM). Using a database of these flares, we have studied hard X-ray (50-850 keV) spectral center-to-limb variation and its evolution with time. The major results are the following: (1) During the rise phase, the center-to-limb spectral variation is small, with a hardness of delta delta = 0.02 +/- 0.25, and a statistical significance of 0.1 sigma. (2) During the peak phase, the center-to-limb variation is delta delta = 0.13 +/- 0.13, with a statistical significance of 1 sigma. (3) During the decay phase, the center-to-limb variation changes to softening. The softness is relatively large with delta delta = -0.25 +/- 0.21, and a statistical significance of 1.2 sigma. (4) The linear least-squares fits to the spectral center-to-limb variations do not have slopes significantly different from zero during all those three phases. (5) The center events and limb events spectral distributions are shown to be not different by using Kolmogorov-Smirnov two-samples test. (6) The fraction of events detected near the limb is marginally consistent with that expected from isotropically emitting flares. (7) On average, flares evolve as soft-hard-soft. These results suggest that there is no statistically significant evidence for hard X-ray directivity during the rise, peak, and decay phases of solar flares. The hard X-ray radiation pattern at those energies is almost isotropic during all those phases. This lack of directivity (or anisotropy) found in this study is not in agreement with the results discovered by Vestrand et al. (1987) in which they found energetic photon source is anisotropic, using SMM Gamma-Ray Spectrometer (GRS) data at a much higher energy band of 0.3-1 MeV. If we want to interpret the results of Vestrand et al. (1987) and our present results in a self-consistent way, we must conclude that at energies of greater than 50 keV, the Compton backscattering suppresses the hard X-ray anisotropy sufficiently so that no significant directivity can be measured during the rise, peak, and decay phases of solar flares. during the rise, peak, and decay phases of solar flares.

Li, Peng↗

Optimizing Science Return from Plato's Observations of the Kepler Field

The Plato mission could substantially enhance its scientific yield by focusing a significant fraction of its observations with all 24 cameras viewing the Kepler Field of View (FOV). Kepler, while exceptionally productive, did not accomplish its goal of providing a direct measurement of the abundance of Earth-analog planet around sunlike stars. On this issue, Plato observations centered on the Kepler FOV would enable planet searches using the combined data sets to achieve more precise estimates of planetary occurrence rates than obtainable using data from either mission alone. The ecliptic latitude of the Kepler FOV is too low for Plato to observe with 24 cameras year-round, so current plan’s for Plato’s Northern Stare places the Kepler FOV in the periphery, reducing the number of photons collected from Kepler targets by factors of 2 or 4. We propose instead four ~6 - 8 month Plato stares with 24 cameras observing the Kepler FOV during the Plato primary mission, instead of a single 2 year stare with the Kepler FOV located to the periphery. This change would greatly increase the amount of 'Plato glass’ dedicated to the continued study of the Kepler FOV, while still enabling Plato to observe other parts of the sky during the majority of its prime mission. In addition to discovering planets below Kepler's detection threshold, this change would yield transit follow-up of thousands of known exoplanets with an effective throughput of a 0.6 m space telescope. It would also extend the observing baseline of the Kepler FOV to timescales of decades, providing key data for stellar cycles, the analysis of transit timing variations, and other decade-timescale variations. Such a campaign enhances the likelihood of the detection of true Earth-analog transiting planet candidates around sunlike stars, provide the best statistical study of exoplanet systems with transit timing variations and duration variations, and open up the study of stellar variability to new time-domains with precision on large numbers of faint targets that only 24 Plato cameras can provide.

Jack Lissauer↗

Automated thematic mapping and change detection of ERTS-A images

The author has identified the following significant results. For the recognition of terrain types, spatial signatures are developed from the diffraction patterns of small areas of ERTS-1 images. This knowledge is exploited for the measurements of a small number of meaningful spatial features from the digital Fourier transforms of ERTS-1 image cells containing 32 x 32 picture elements. Using these spatial features and a heuristic algorithm, the terrain types in the vicinity of Phoenix, Arizona were recognized by the computer with a high accuracy. Then, the spatial features were combined with spectral features and using the maximum likelihood criterion the recognition accuracy of terrain types increased substantially. It was determined that the recognition accuracy with the maximum likelihood criterion depends on the statistics of the feature vectors. Nonlinear transformations of the feature vectors are required so that the terrain class statistics become approximately Gaussian. It was also determined that for a given geographic area the statistics of the classes remain invariable for a period of a month but vary substantially between seasons.

Gramenopoulos, N.↗

Star Tracker Accuracy Improvement and Optimization for Attitude Measurement in Three-Axis

High precision attitude measurement systems obviate the need for the beacon from the receiver making it possible for the spacecraft to beam a laser communications signal to a ground station without the ground station advertising its location. The research presented targets new detection and estimation methods to improve the accuracy in locating stars on a focal plane detector, and an understanding of the effects of changes in the optics design parameters and aberration, including defocus, on the navigation solution itself. This understanding can lead to an optimization of the attitude solution with respect to those optics realm parameter changes. The methodology discussed includes the development of a model of a current star tracker system. Using this model, multiple algorithms are implemented, including a multi-hypothesis method (MHT), to detect and estimate the position of the stars on the focal plane detector. It will be shown that using the MHT for detection and estimation, a greater accuracy can be found for each star estimation from more traditional detection and estimation algorithms. The approach then uses the model to develop statistics of the star tracker and the attitude estimation outputs to understand the accuracy, or variance, of the system's attitude solution. This solution is repeated for a range of defocus aberration, and a lower limit to the variance of the attitude solution is shown. A Cramer Rao lower bound solution is derived for the star tracker system and the results are compared to the Monte Carlo analysis from the model and shown to correlate very well. The approach uses a star image not as a Gaussian spot on the focal plane as done in previous work, and use of an image that includes the effects of aberrations of the optic system, and the effects of under-sampling and noise from the focal plane detector as well. Analysis includes exploring a star tracker's accuracy improvement through the combination of focus error and under-sampling effects alone, possibly contradicting conventional wisdom and approaches.

Lichter, Michael J.↗