SEARCH · Search NASA
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.
The statistical properties of a cubic smoothing spline and its derivative are analyzed. It is shown that unless unnatural boundary conditions hold, the integrated squared bias is dominated by local effects near the boundary. Similar effects are shown to occur in the regularized solution of a translation-kernel intergral equation. These results are derived by developing a Fourier representation for a smoothing spline.
Several fields of orbital gamma ray spectroscopy data have been deconvolved in order to model the distribution of Th over the Mare Imbrium and northern Oceanus Procellarum portions of the Apollo 15 lunar ground track, which in combination with a prior study of the Apenninus region covers a continuous swath from 10 deg E to 60 deg W in the northwest quadrant. The crater of the Aristarchus region dominates the Th distribution, with a concentration of 20 ppm, and substantial enhancements are also found in the mare regions around Brayley and at the ejecta blankets of Timocharis and Lambert. The existence of enhanced Th concentrations in mare basalt regions suggests that reservoirs of some late stage mare basalts incorporated KREEP-rich material during formation or transit.
A solution is presented to the problem of deconvolving wide field of view measurements of earth-emitted radiation for the case of limb-darkening functions which are geographically dependent. The method of linear perturbations is used and the effects of limb-darkening function variations at various scales are evaluated. A simple model for the geographical variation of the limb-darkening function based on the radiant exitance is introduced. It is shown that the effect of these variations on the global average radiant exitance vanishes identically, and is insignificant for large spatial scales. For wave numbers six through twelve, the limb-darkening effects are small but possibly significant. At wave numbers greater than twelve, the noise due to meteorological variability dominates for an orbit altitude of 1100 km. The smallness of these effects is due to the smoothness of the monthly averaged radiant exitance field.
The development of vector processing computers of the streaming array architectures has made possible a dramatic decrease in the time required for the solution of problems, i.e., for those algorithms which readily lend themselves to sequential operations on long vectors. There has concurrently been rapid growth in the applications of the techniques generally known as mathematical digital filtering, also called signal analysis, digital signal processing, time series analysis, and digital image processing. The best known applications of these techniques is to seismic data and two-dimensional images, data types consisting of very large collections of numbers. A major limitation for these cases is the size and speed of the computer available. Therefore, a move to the streaming array architecture can result in marked improvement in the data analysis techniques which may be employed.
Linear filtering techniques currently used for the restoration of noisy, blurred or otherwise degraded image data are discussed and new techniques related to the iterative techniques of Morrison and van Cittert are developed and implemented. Programs written for the implementation are discussed in the appendices. It is shown that the new techniques are convergent for any system response function, and they are applied to the task of restoring a severely blurred image.
Current observational questions in asteroseismology require high spectral resolution that can be only be obtained with observations spanning many days or months. The primary constraint in the full utilization of single mid-latitude observing sites is the presence of diurnal data gaps. Several methods for removing the effect of these gaps in the spectral obtained from velocity observations have been suggested. In the case of data coverage of less than 50%, none of the methods considered has been successful at unambiguously recovering the true spectrum from test spectra with realistic complexity. The limitations of these methods and their applicability to the helioseismology problem is discussed.
Explore the source record for details and available documents.
Visible to near infrared reflectance spectra of macroscopic mixtures have been shown to be linear combinations of the reflections of the pure mineral components in the mixture. However, for microscopic mixtures the mixing systematics are in general nonlinear. The systematics may be linearized by conversion of reflectance to single scattering albedo (w), where the equations which relate reflectance to w depend on the method of data collection. Several proposed mixing models may be used to estimate mineral abundances from the reflectance spectra of intimate mixtures. These models are summarized and a revised model is presented. A noniterative (linear) least squares approach was used for curve fitting and the data, measured as bi-directional reflectance with incidence and emergence angles of 30 and 0 deg were converted to w by a simplified version of Hapke's equation for bi-directional reflectance. This model was tested with two mixture series composed of 45 to 75 micron particles: an anorthite-enstatite series and an olivine-magnetite series. The data indicate that the simplified Hapke's equation may be used to convolve reflectance spectra into mineral abundances if appropriate endmembers are known or derived from other techniques. For surfaces that contain a significant component of very low albedo material, a somewhat modified version of this technique will need to be developed. Since the abundances are calculated using a noniterative approach, the application of this method is especially efficient for large spectral data sets, such as those produced by mapping spectrometers.
Current observational questions in asteroseismology require high spectral resolution that can only be obtained with observations spanning many days or months. The primary constraint in the full utilization of single mid-latitude observing sites is the presence of diurnal data gaps. Several methods for removing the effect of these gaps in the spectral obtained from velocity observations have been suggested. In the case of data coverage of less than 50 percent, none of the methods considered has been successful at unambiguously recovering the true spectrum from test spectra with realistic complexity. The limitations of these methods and their applicability to the helioseismology problem is discussed.
An algorithm has been developed which reduces the effects of (deconvolves) instrument-induced spatial crosstalk in satellite image data by several orders of magnitude where highly precise radiometry is required. The algorithm is based upon radiance transfer ratios which are defined as the fractional bilateral exchange of energy betwen pixels A and B.
The measurement of reflected solar radiation by spacecraft wide-field-of-view radiometers is considered analytically, and the method proposed by Smith (1987) for the solution of the discretized governing equations is described and demonstrated. The theoretical basis of the singular-value decomposition method (partitioning the albedo field into observable and unobservable components) is outlined; its application to Nimbus 7 ERB data is explained in detail; and numerical results are presented in graphs and maps and briefly characterized.
The block iterative method (Young 1971) of algebraic image restoration is applied to a photographic plate of the field around SN 1987A obtained 4 yr before outburst by Chu at the CTIO 4-m telescope. By setting appropriate constraints for the solution on a pixel-by-pixel basis, the image of a starfield centered on Sk -69 202, the precursor to the supernova, is restored. It is found that this star has two companions, which others (e.g., Walborn et al., 1987 and West et al., 1987) have also noted and designated as Star 2 and Star 3. Also found is marginal evidence for Star 4, a weak source at theta = 222 deg, r = 1.9 arcsec. Correction for the contributions of these companions yields an estimate of the magnitude for Sk -69 202 alone of V = 12.37.
The technique of maximum entropy image restoration is applied to the problem of deconvolving the point spread function from a deep, high-quality V band image of the optical jet of 3C 273. The resulting maximum entropy image has an approximate spatial resolution of 0.6 arcsec and has been used to study the morphology of the optical jet. Four regularly-spaced optical knots are clearly evident in the data, together with an optical 'extension' at each end of the optical jet. The jet oscillates around its center of gravity, and the spatial scale of the oscillations is very similar to the spacing between the optical knots. The jet is marginally resolved in the transverse direction and has an asymmetric profile perpendicular to the jet axis. The distribution of V band flux along the length of the jet, and accurate astrometry of the optical knot positions are presented.
Wide-field-of-view (WFOV) radiometers have been flown as part of the Earth Radiation Budget instrument on the Nimbus 6 and 7 spacecraft and as part of the Earth Radiation Budget Experiment (ERBE) instruments aboard the ERBE spacecraft and also the NOAA 9 and 10 operational spacecraft. The measurement is the integral of the reflected solar flux distribution at the top of the earth-atmosphere system over the field-of-view of the radiometer. This paper develops the solution to this two-dimensional integral equation for the albedo distribution in terms of the measurements.
Increasingly, aircraft system designs require that aerodynamic parameters derived from pneumatic measurements be employed as control-system feedbacks. Such high frequency pressure measurements' accuracy is compromised by pressure distortion due to frictional attenuation and pneumatic resonance within the sensing system. A pneumatic distortion model is here formulated and reduced to a low-order state-variable model which retains most of the full model's dynamic characteristics. This reduced-order model is coupled with standard results from minimum variance estimation theory to develop an algorithm to compensate for pneumatic-distortion effects.
The power spectrum for a stationary random process can be defined with the Wiener-Khintchine Theorem, which says that the power spectrum and the auto correlation function are a Fourier transform pair. To implement this theorem for signals that are discrete and of finite length we can use the Blackman-Tukey method. Blackman and Tukey (1958) show that a function w(tau), called a lag window, can be applied to the auto correlation estimates to obtain power spectrum estimates that are statistically stable. The Fourier transform of w(r) is called a spectral window. Typical choices for spectral windows show a distinct trade-off between the main lobe width and side lobe strength. A new idea for designing windows by taking linear combinations of the standard windows to produce hybrid windows was introduced by Smith (1985). We implement Smith's idea to obtain spectral windows with narrow main lobes and smaller (compared with typical windows) near side lobes. One of the main contributions of this thesis is that we show that Smith's problem is equivalent to a Quadratic Programming (QP) problem with linear equality and inequality constraints. A computer program was written to produce hybrid windows by setting up and solving the QP problem. We also developed and solved two variations of the original problem. The two variations involved changing the inequality constraints in both cases from non negativity on the combination coefficients to non negativity on the hybrid lag window itself. For the second variation, the window functions used to construct the hybrid window were changed to a frequency-variable set of truncated cosinusoids. A series of tests was run with the three computer programs to investigate the behavior of the hybrid spectral and lag windows. Emphasis was put on obtaining spectral windows with both relatively narrow main lobes and the lowest possible (for these algorithms) near side lobes. Some success was achieved for this goal. A 10 dB peak side lobe reduction over the rectangular spectral window without significant main lobe broadening was achieved. Also, average side lobe levels of -117 dB were reached at a cost of doubling the main lobe width (at the -3 dB point).
This thesis reviews the technique established to clear channels in the Power Spectral Estimate by applying linear combinations of well known window functions to the autocorrelation function. The need for windowing the auto correlation function is due to the fact that the true auto correlation is not generally used to obtain the Power Spectral Estimate. When applied, the windows serve to reduce the effect that modifies the auto correlation by truncating the data and possibly the autocorrelation has on the Power Spectral Estimate. It has been shown in previous work that a single channel has been cleared, allowing for the detection of a small peak in the presence of a large peak in the Power Spectral Estimate. The utility of this method is dependent on the robustness of it on different input situations. We extend the analysis in this paper, to include clearing up to three channels. We examine the relative positions of the spikes to each other and also the effect of taking different percentages of lags of the auto correlation in the Power Spectral Estimate. This method could have application wherever the Power Spectrum is used. An example of this is beam forming for source location, where a small target can be located next to a large target. Other possibilities extend into seismic data processing. As the method becomes more automated other applications may present themselves.
Because of the interesting science which can be performed using a satellite attached by a very long tether to a mother vehicle in orbit, such as the Space Shuttle, NASA will deploy TSS-1 (Tethered Satellite System) in 1992. A very long tether (20 km in this case) has the possibility of undergoing oscillations of several different types, or modes, and higher harmonics of these modes. The purpose of this document is to describe a method for detecting the amplitude, frequency, and phase (and predicting future motion in the steady state) of these modes, in particular, the skiprope mode, using tethered satellite dynamics measurements. Specifically the rotation rate data about two orthogonal axes, calculated from output from satellite gyroscopes, are used. The data of interest are the satellite pitch and roll rate measurements. NASA has determined to use two methods to diagnose skiprope properties and predict future values. One of these, a Fourier transform domain approach, is the subject of this notebook. The main program and all subroutines are described along with the test plan for evaluating the Frequency Domain Skiprope Observer.