Search NASA⌕ Search

SEARCH · Search NASA

Results for “Directional Gaussian smoothing”

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.

Improving the reconstruction of the velocity potential and primordial density fluctuations by choice of smoothing windows

Recently considerable attention has been focused on improving algorithms for restoring primordial density fluctuations in the universe or investigating large-scale velocity fields by going to higher order approximations relating nonlinear states to their initial conditions. An alternative approach is to investigate the use of alternative smoothing windows within the context of first-order approximations. I present evidence that the universally used Gaussian smoothing window is far from optimal. A sharp truncation P(k) = 0 for k greater than k(c) leads to a much more direct connection to initial conditions for either 'Gaussianization' methods or use of the Zel'dovich approximation, as discussed in the text.

Melott, Adrian L.↗

Controlled experiments in cosmological gravitational clustering

A systematic study is conducted of gravitational instability in 3D on the basis of power-law initial spectra with and without spectral cutoff, emphasizing nonlinear effects and measures of nonlinearity; effects due to short and long waves in the initial conditions are separated. The existence of second-general pancakes is confirmed, and it is noted that while these are inhomogeneous, they generate a visually strong signal of filamentarity. An explicit comparison of smoothed initial conditions with smoothed envelope models also reconfirms the need to smooth over a scale larger than any nonlinearity, in order to extrapolate directly by linear theory from Gaussian initial conditions.

Melott, Adrian L.↗

Theory for computing the field scattered from a smooth inflected surface

A theory is described for computing the reflected or scattered field from a smooth body with inflection points. These inflections occur in certain directions at each surface point for which the total (Gaussian) curvature is zero or negative. For surface illumination in one of these critical directions, the usual reflection formulas obtained by the high-frequency approximation are inapplicable, and a shadow zone exists in the reflected field. Scattering into the shadow zone is treated, as well as specular reflection. This theory should have a variety of applications such as for certain optics problems, computer graphics modeling of three-dimensional shapes, and the design and analysis of specialized microwave reflector antennas.

Barger, R. L.↗

Analysis of high-incidence separated flow past airfoils

An unsteady Navier-Stokes (NS) analysis is developed and used to carefully examine high-incidence aerodynamic separated flows past airfoils. Clustered conformal C-grids are employed for the 12 percent thick symmetric Joukowski airfoil as well as for the NACA 0012 airfoil with a sharp trailing edge. The clustering is controlled by appropriate one-dimensional stretching transformations. An attempt is made to resolve many of the dominant scales of an unsteady flow with massive separation, while maintaining the transformation metrics to be smooth and continuous in the entire flow field. A fully implicit time-marching alternating-direction implicit-block Gaussian elimination (ADI-BGE) method is employed, in which no use is made of any explicit artificial dissipation. Detailed results are obtained for massively separated, unsteady flow past symmetric Joukowski and NACA 0012 airfoils.

Chia, K. N.↗

Partial and Interaction Spline Models for the Semiparametric Estimation of Functions of Several Variables

A partial spline model is a model for a response as a function of several variables, which is the sum of a smooth function of several variables and a parametric function of the same plus possibly some other variables. Partial spline models in one and several variables, with direct and indirect data, with Gaussian errors and as an extension of GLIM to partially penalized GLIM models are described. Application to the modeling of change of regime in several variables is described. Interaction splines are introduced and described and their potential use for modeling non-linear interactions between variables by semiparametric methods is noted. Reference is made to recent work in efficient computational methods.

Wahba, Grace↗

The topology of large-scale structure. I - Topology and the random phase hypothesis

Many models for the formation of galaxies and large-scale structure assume a spectrum of random phase (Gaussian), small-amplitude density fluctuations as initial conditions. In such scenarios, the topology of the galaxy distribution on large scales relates directly to the topology of the initial density fluctuations. Here a quantitative measure of topology - the genus of contours in a smoothed density distribution - is described and applied to numerical simulations of galaxy clustering, to a variety of three-dimensional toy models, and to a volume-limited sample of the CfA redshift survey. For random phase distributions the genus of density contours exhibits a universal dependence on threshold density. The clustering simulations show that a smoothing length of 2-3 times the mass correlation length is sufficient to recover the topology of the initial fluctuations from the evolved galaxy distribution. Cold dark matter and white noise models retain a random phase topology at shorter smoothing lengths, but massive neutrino models develop a cellular topology.

Weinberg, David H.↗

International Conference on Advances in Communication and Control Systems, 1st, Washington, DC, June 18-20, 1987, Proceedings

Theoretical models of communication and control systems are discussed in reviews and reports. Topics addressed include smoothing and identification for random fields, the information and coding capacities of mismatched Gaussian channels, recursive least-squares estimation and Kalman filtering by systolic arrays, Kemp echo digital filters, a periodic test-scheduling scheme for communication and queuing processes, and receivers for direct-sequence systems. Consideration is given to a distributed-parameter model for detecting cracks in rotors, active control of aeroelastic systems governed by functional differential equations, robust multivariable control of large space structures, finite-rank relatively bounded perturbations of semigroup generators, and sensitivity analysis of convex optimal-control problems.

Declaris, Nicholas↗

Application of super smooth optics to extra-solar planet detection

The observational constraints on direct detection of extrasolar planets around nearby stars are defined. The design of an instrument capable of meeting these goals requires the efficient reduction of diffracted and scattered light by about three orders of magnitude below the airy wings of the unapodized aperture. The diffraction reduction is achieved through the high-efficiency coronagraph, which uses Gaussian occulting masks in the first focus to concentrate better diffracted light in the pupil plane. Laboratory experiments confirm this high efficiency to the limits set by the experimental setup. Scattering due to midspatial scale figure error in the mirrors must also be reduced to a comparable level.

Terrile, Richard J.↗

High spatial resolution 100 micron observations of the M83 bar

A program of high spatial resolution far-infrared observations of galaxies using the Kuiper Airborne Observatory (KAO), was conducted to better understand the role of star formation, the general interstellar radiation field, and non-thermal activity in powering the prodigious far-infrared luminosities seen in spiral and interacting galaxies. Here, researchers present observations of the central region of the well-known barred spiral M83 (NGC 5236). The resultant channel 3 scans for M83 and IRC + 10216, after co-addition and smoothing, are shown. These data show that M83 is extended at 100 microns compared to a point source. A simple Gaussian deconvolution of the M83 data with the point source profile from IRC+10216 gives a full width half maximum (FWHM) of about 19 seconds for M83. By comparison with IRC+10216, researchers obtain a flux for the unresolved component in M83 of about 110 Jy. This is about 1/6 the total flux for M83 (Rice et al. 1988) and about 1/2 the PSC flux. The M83 and IRC+10216 profiles in the cross-scan direction (SE-NW) were also compared, and show that M83 is extended in this direction as well, with a width of about 18 seconds. A comparison of the different channel profiles for M83 and IRC+10216 shows that there is an asymmetry in the M83 data, in that the maximum in the profiles shifts from southeast to northwest as channel number increases. This corresponds to the extension in the bar seen in the CO data. Thus the far-infrared emission in the central region of M83 tends to trace the CO bar. The new 100 micron data is also compared with previous H alpha observations from the literature, to determine how well the far-infrared traces the stellar structure, the star formation as measured by H alpha, and the optical colors.

Smith, B. J.↗

Gravitational clustering in the expanding universe - Controlled high-resolution studies in two dimensions

Results are presented from a series of gravitational clustering simulations in two dimensions. These simulations are a significant departure from previous work, since in two dimensions one can have large dynamic range in both length scale and mass using present computer technology. Controlled experiments were conducted by varying the slope of power-law initial density fluctuation spectra and varying cutoffs at large k, while holding constant the phases of individual Fourier components and the scale of nonlinearity. Filaments are found in many different simulations, even with pure power-law initial conditions. By direct comparison, filaments, called 'second-generation pancakes' are shown to arise as a consequence of mild nonlinearity on scales much larger than the correlation length and are not relics of an initial lattice or due to sparse sampling of the Fourier components. Bumps of low amplitude in the two-point correlation are found to be generic but usually only statistical fluctuations. Power spectra are much easier to relate to initial conditions, and seem to follow a simple triangular shape (on log-log plot) in the nonlinear regime. The rms density fluctuation with Gaussian smoothing is the most stable indicator of nonlinearity.

Beacom, John Francis↗

Pursuit Latency for Chromatic Targets

The temporal dynamics of eye movement response to a change in direction of stimulus motion has been used to compare the processing speeds of different types of stimuli (Mulligan, ARVO '97). In this study, the pursuit response to colored targets was measured to test the hypothesis that the slow response of the chromatic system (as measured using traditional temporal sensitivity measures such as contrast sensitivity) results in increased eye movement latencies. Subjects viewed a small (0.4 deg) Gaussian spot which moved downward at a speed of 6.6 deg/sec. At a variable time during the trajectory, the dot's direction of motion changed by 30 degrees, either to the right or left. Subjects were instructed to pursue the spot. Eye movements were measured using a video ophthalmoscope with an angular resolution of approximately 1 arc min and a temporal sampling rate of 60 Hz. Stimuli were modulated in chrominance for a variety of hue directions, combined with a range of small luminance increments and decrements, to insure that some of the stimuli fell in the subjects' equiluminance planes. The smooth portions of the resulting eye movement traces were fit by convolving the stimulus velocity with an exponential having variable onset latency, time constant and amplitude. Smooth eye movements with few saccades were observed for all stimuli. Pursuit responses to stimuli having a significant luminance component are well-fit by exponentials having latencies and time constants on the order of 100 msec. Increases in pursuit response latency on the order of 100-200 msec are observed in response to certain stimuli, which occur in pairs of complementary hues, corresponding to the intersection of the stimulus section with the subjects' equiluminant plane. Smooth eye movements can be made in response to purely chromatic stimuli, but are slower than responses to stimuli with a luminance component.

Mulligan, Jeffrey B.↗

Downscaling Satellite Precipitation with Emphasis on Extremes: A Variational 1-Norm Regularization in the Derivative Domain

The increasing availability of precipitation observations from space, e.g., from the Tropical Rainfall Measuring Mission (TRMM) and the forthcoming Global Precipitation Measuring (GPM) Mission, has fueled renewed interest in developing frameworks for downscaling and multi-sensor data fusion that can handle large data sets in computationally efficient ways while optimally reproducing desired properties of the underlying rainfall fields. Of special interest is the reproduction of extreme precipitation intensities and gradients, as these are directly relevant to hazard prediction. In this paper, we present a new formalism for downscaling satellite precipitation observations, which explicitly allows for the preservation of some key geometrical and statistical properties of spatial precipitation. These include sharp intensity gradients (due to high-intensity regions embedded within lower-intensity areas), coherent spatial structures (due to regions of slowly varying rainfall),and thicker-than-Gaussian tails of precipitation gradients and intensities. Specifically, we pose the downscaling problem as a discrete inverse problem and solve it via a regularized variational approach (variational downscaling) where the regularization term is selected to impose the desired smoothness in the solution while allowing for some steep gradients(called 1-norm or total variation regularization). We demonstrate the duality between this geometrically inspired solution and its Bayesian statistical interpretation, which is equivalent to assuming a Laplace prior distribution for the precipitation intensities in the derivative (wavelet) space. When the observation operator is not known, we discuss the effect of its misspecification and explore a previously proposed dictionary-based sparse inverse downscaling methodology to indirectly learn the observation operator from a database of coincidental high- and low-resolution observations. The proposed method and ideas are illustrated in case studies featuring the downscaling of a hurricane precipitation field.

Hurricanes↗

NGC 1058: Gas motions in an extended, quiescent spiral disk

Researchers investigate in detail the motion of gas in the galaxy NGC 1058 using the very large array (VLA) to map the emission in the 21-cm line. This galaxy is so nearly face-on that the contribution to the line width due to the variation of the rotational velocity across the D-array beam is small compared with the random z-motion of the gas. Researchers confirm results of earlier studies (Lewis 1987, A. and A. Suppl., 63, 515; van der Kruit and Shostak 1984, A. and A., 134, 258) of the galaxy's total neutral hydrogen (HI) and kinematics, including the fact that the rotation curve drops faster than Keplerian at the outer edge of the disk, which is interpreted as a fortuitous twist of the plane of rotation in the outer disk. However, their very high velocity resolution (2.58 km s(exp -1) after Hanning smoothing) coupled with good spatial resolution, allows researchers to measure more accurately the line width, and even to some extent its shape, throughout the disk. One of the most interesting results of this study is the remarkable constancy of the line width in the outer disk. From radius 90 to 210 seconds the Gaussian velocity dispersion (sigma sub nu) of the 21-cm line has a mean value of 5.7 km s(exp -1) (after correcting for the spectral resolution) with a dispersion of less than 0.9 km s(exp -1) (after correcting for the spectral resolution) with a dispersion of less than 0.9 km s(exp -1). Translating this directly into a kinetic temperature (Doppler temperature): T sub Dopp equals 121K (sigma sub mu exp 2/(km s(exp -1) (exp 2) gives 4000 K, with a dispersion of less than 1500 K over the outer disk. This constancy is observed even when comparing the spiral arms versus inter-arm regions, which in the radius range from 100 to 150 seconds the surface density modulates (defined as the ratio N sub peak -N sub trough/N sub peak + N sub trough) from 0.5 to 0.25 in the range 150 to 200 seconds.

Hanson, Margaret Murray↗

Transition in a Supersonic Boundary-Layer Due to Roughness and Acoustic Disturbances

The transition process induced by the interaction of an isolated roughness with acoustic disturbances in the free stream is numerically investigated for a boundary layer over a flat plate with a blunted leading edge at a free stream Mach number of 3.5. The roughness is assumed to be of Gaussian shape and the acoustic disturbances are introduced as boundary condition at the outer field. The governing equations are solved using the 5'h-~rder accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third- order total-variation-diminishing (TVD) Runge- Kutta scheme for time integration. The steady field induced by the two and three-dimensional roughness is also computed. The flow field induced by two-dimensional roughness exhibits different characteristics depending on the roughness heights. At small roughness heights the flow passes smoothly over the roughness, at moderate heights the flow separates downstream of the roughness and at larger roughness heights the flow separates upstream and downstream of the roughness. Computations also show that disturbances inside the boundary layer is due to the direct interaction of the acoustic waves and isolated roughness plays a minor role in generating instability waves.

Balakumar, P.↗

Transition Marshall Space Flight Center Wind Profiler Splicing Algorithm to Launch Services Program Upper Winds Tool

NASAs LSP customers and the future SLS program rely on observations of upper-level winds for steering, loads, and trajectory calculations for the launch vehicles flight. On the day of launch, the 45th Weather Squadron (45 WS) Launch Weather Officers (LWOs) monitor the upper-level winds and provide forecasts to the launch team via the AMU-developed LSP Upper Winds tool for launches at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station. This tool displays wind speed and direction profiles from rawinsondes released during launch operations, the 45th Space Wing 915-MHz Doppler Radar Wind Profilers (DRWPs) and KSC 50-MHz DRWP, and output from numerical weather prediction models.The goal of this task was to splice the wind speed and direction profiles from the 45th Space Wing (45 SW) 915-MHz Doppler radar Wind Profilers (DRWPs) and KSC 50-MHz DRWP at altitudes where the wind profiles overlap to create a smooth profile. In the first version of the LSP Upper Winds tool, the top of the 915-MHz DRWP wind profile and the bottom of the 50-MHz DRWP were not spliced, sometimes creating a discontinuity in the profile. The Marshall Space Flight Center (MSFC) Natural Environments Branch (NE) created algorithms to splice the wind profiles from the two sensors to generate an archive of vertically complete wind profiles for the SLS program. The AMU worked with MSFC NE personnel to implement these algorithms in the LSP Upper Winds tool to provide a continuous spliced wind profile.The AMU transitioned the MSFC NE algorithms to interpolate and fill data gaps in the data, implement a Gaussian weighting function to produce 50-m altitude intervals in each sensor, and splice the data together from both DRWPs. They did so by porting the MSFC NE code written with MATLAB software into Microsoft Excel Visual Basic for Applications (VBA). After testing the new algorithms in stand-alone VBA modules, the AMU replaced the existing VBA code in the LSP Upper Winds tool with the new algorithms. They then tested the code in the LSP Upper Winds tool with archived data. The tool will be delivered to the 45 WS after the 50-MHz DRWP upgrade is complete and the tool is tested with real-time data. The 50-MHz DRWP upgrade is expected to be finished in October 2014.

Space launch↗

On three dimensional object recognition and pose-determination: An abstraction based approach

A method of computing reliable Gaussian and mean curvature sign-map descriptors from the polynomial approximation of surfaces was demonstrated. Such descriptors which are invariant under perspective variation are suitable for hypothesis generation. A means for determining the pose of constructed geometric forms whose algebraic surface descriptors are nonlinear in terms of their orienting parameters was developed. This was done by means of linear functions which are capable of approximating nonlinear forms and determining their parameters. It was shown that biquadratic surfaces are suitable companion linear forms for cylindrical approximation and parameter estimation. The estimates provided the initial parametric approximations necessary for a nonlinear regression stage to fine tune the estimates by fitting the actual nonlinear form to the data. A hypothesis-based split-merge algorithm for extraction and pose determination of cylinders and planes which merge smoothly into other surfaces was developed. It was shown that all split-merge algorithms are hypothesis-based. A finite-state algorithm for the extraction of the boundaries of run-length regions was developed. The computation takes advantage of the run list topology and boundary direction constraints implicit in the run-length encoding.

Quek, Kok How Francis↗

HEASARC Software Archive

(1) Chandra Archive: SAO has maintained the interfaces through which HEASARC gains access to the Chandra Data Archive. At HEASARC's request, we have implemented an anonymous ftp copy of a major part of the public archive and we keep that archive up-to- date. SAO has participated in the ADEC interoperability working group, establishing guidelines or interoperability standards and prototyping such interfaces. We have provided an NVO-based prototype interface, intending to serve the HEASARC-led NVO demo project. HEASARC's Astrobrowse interface was maintained and updated. In addition, we have participated in design discussions surrounding HEASARC's Caldb project. We have attended the HEASARC Users Group meeting and presented CDA status and developments. (2) Chandra CALDB: SA0 has maintained and expanded the Chandra CALDB by including four new data file types, defining the corresponding CALDB keyword/identification structures. We have provided CALDB upgrades for the public (CIAO) and for Standard Data Processing. Approximately 40 new files have been added to the CALDB in these version releases. There have been in the past year ten of these CALDB upgrades, each with unique index configurations. In addition, with the inputs from software, archive, and calibration scientists, as well as CIAO/SDP software developers, we have defined a generalized expansion of the existing CALDB interface and indexing structure. The purpose of this is to make the CALDB more generally applicable and useful in new and future missions that will be supported archivally by HEASARC. The generalized interface will identify additional configurational keywords and permit more extensive calibration parameter and boundary condition specifications for unique file selection. HEASARC scientists and developers from SAO and GSFC have become involved in this work, which is expected to produce a new interface for general use within the current year. (3) DS9: One of the decisions that came from last year's HEADCC meeting was to make the ds9 image display program the primary vehicle for displaying line graphics (as well as images). The first step required to make this possible was to enhance the line graphics capabilities of ds9. SAO therefore spent considerable effort upgrading ds9 to use Tcl 8.4 so that the BLT line graphics package could be built and imported into ds9 from source code, rather than from a pre-built (and generally outdated) shared library. This task, which is nearly complete, allows us to extend BLT as needed for the HEAD community. Following HEADCC discussion concerning archiving and the display of archived data, we extended ds9 to support full access to many astronomical Web-based archives sites, including HEASARC, MAST, CHANDRA, SKYVIEW, ADS, NED, SIMBAD, IRAS, NVRO, SAO TDC, and FIRST. Using ds9's new internal Web access capabilities, these archives can be accessed via their Web page. FITS images, plots, spectra, and journal abstracts can be referenced, down-loaded, and displayed directly and easily in ds9. For more information, see: http://hea-www.harvard.edu/saord/ds9. Also after the HEADCC discussion concerning region filtering, we extended the Funtools sample implementation of region filtering as described in: http://hea-www.harvard.edu/saord/funtools/regions.html. In particular, we added several new composite regions for event and image filtering, including elliptical and box annuli. We also extended the panda (Pie AND Annulus) region support to include box pandas and elliptical pandas. These new composite regions are especially useful in programs that need to count photons in each separate region using only a single pass through the data. Support for these new regions was added to ds9. In the same vein, we developed new region support for filtering images using simple FITS image masks, i.e. 8-bit or 16-bit FITS images where the value of a pixel is the region id number for that pixel. Other important enhancements to DS9 this year, include supporor multiple world coordinate systems, three dimensional event file binning, image smoothing, region groups and tags, the ability to save images in a number of image formats (such as JPEG, TIFF, PNG, FITS), improvements in support for integrating external analysis tools, and support for the virtual observatory. In particular, a full-featured web browser has been implemented within D S 9 . This provides support for full access to HEASARC archive sites such as SKYVIEW and W3BROWSE, in addition to other astronomical archives sites such as MAST, CHANDRA, ADS, NED, SIMBAD, IRAS, NVRO, SA0 TDC, and FIRST. From within DS9, the archives can be searched, and FITS images, plots, spectra, and journal abstracts can be referenced, downloaded and displayed The web browser provides the basis for the built-in help facility. All DS9 documentation, including the reference manual, FAQ, Know Features, and contact information is now available to the user without the need for external display applications. New versions of DS9 maybe downloaded and installed using this facility. Two important features used in the analysis of high energy astronomical data have been implemented in the past year. The first is support for binning photon event data in three dimensions. By binning the third dimension in time or energy, users are easily able to detect variable x-ray sources and identify other physical properties of their data. Second, a number of fast smoothing algorithms have been implemented in DS9, which allow users to smooth their data in real time. Algorithms for boxcar, tophat, and gaussian smoothing are supported.

White, Nicholas↗