Engineering PapersSearch

SEARCH · Engineering Papers

Results for “image decomposition”

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 55 records · Page 3

Stereo imaging velocimetry for microgravity applications

Stereo imaging velocimetry is the quantitative measurement of three-dimensional flow fields using two sensors recording data from different vantage points. The system described in this paper, under development at NASA Lewis Research Center in Cleveland, Ohio, uses two CCD cameras placed perpendicular to one another, laser disk recorders, an image processing substation, and a 586-based computer to record data at standard NTSC video rates (30 Hertz) and reduce it offline. The flow itself is marked with seed particles, hence the fluid must be transparent. The velocimeter tracks the motion of the particles, and from these we deduce a multipoint (500 or more), quantitative map of the flow. Conceptually, the software portion of the velocimeter can be divided into distinct modules. These modules are: camera calibration, particle finding (image segmentation) and centroid location, particle overlap decomposition, particle tracking, and stereo matching. We discuss our approach to each module, and give our currently achieved speed and accuracy for each where available.

Miller, Brian B.

Characterization of the Inner Knot of the Crab: the Site of the Gamma-ray Flares?

One of the most intriguing recent discoveries has been the detection of powerful gamma-ray flares from the Crab Nebula. Such events, with a recurrence time of about once per year, can be so dramatic to make the system the brightest source in the gamma‐ray sky as occurred, e.g. in April 2011. These flares challenge our understanding of how pulsar wind nebulae work and defy current astrophysical models for particle acceleration. We present here our study of the inner knot located within a fraction of an arcsecond from the pulsar with the aim of characterizing the feature and asking if this might be the site of the origin of the gamma-ray flares. We took data using Keck, HST, and Chandra obtained as part of our multi‐wavelength campaign to identify the source of the enigmatic flares. We set an upper limit as to the gamma-ray flux from the knot. We also find that the dimensions, surface brightness, flux, etc. of the optical and infrared knot are all correlated with distance from the pulsar. This distance, in turn, varies with time. In addition to this most thorough characterization of the inner knot's properties, we examine the hypothesis that the knot may be the site of the flares by examining the knot separation versus the Fermi/LAT gamma-ray flux. Finally, as part of this research, we make use of a new approach employing singular value decomposition (SVD) for analyzing time series of images and compare the approach to more traditional methods. Our conclusions are only refined but not impacted by using the new approach.

Gamma-ray Flares

Advanced Background Subtraction Applied to Aeroacoustic Wind Tunnel Testing

An advanced form of background subtraction is presented and applied to aeroacoustic wind tunnel data. A variant of this method has seen use in other fields such as climatology and medical imaging. The technique, based on an eigenvalue decomposition of the background noise cross-spectral matrix, is robust against situations where isolated background auto-spectral levels are measured to be higher than levels of combined source and background signals. It also provides an alternate estimate of the cross-spectrum, which previously might have poor definition for low signal-to-noise ratio measurements. Simulated results indicate similar performance to conventional background subtraction when the subtracted spectra are weaker than the true contaminating background levels. Superior performance is observed when the subtracted spectra are stronger than the true contaminating background levels. Experimental results show limited success in recovering signal behavior for data where conventional background subtraction fails. They also demonstrate the new subtraction technique's ability to maintain a proper coherence relationship in the modified cross-spectral matrix. Beam-forming and de-convolution results indicate the method can successfully separate sources. Results also show a reduced need for the use of diagonal removal in phased array processing, at least for the limited data sets considered.

Bahr, Christopher J.

Relating Lightning Flash Characteristics to Convective Structure in Tropical Cyclones

For decades, lightning has been used with some success as an indicator of tropical cyclone (TC) intensity change. An increase in inner-core lightning is usually associated with TC intensification, but this is not always the case. Significant lightning outbreaks have been observed in weakening storms, and hurricanes can rapidly intensify in the total absence of inner-core lightning. These conflicting results arise because lightning production depends on convective and microphysical characteristics that can arise during both intensification and weakening phases in TCs. Recent work using flash size and energy metrics provided by the Geostationary Lightning Mapper (GLM) aboard the GOES-R series of satellites shows some promise in distinguishing between lightning outbreaks that correspond to intensification and those that correspond to weakening. If a hurricane is rapidly intensifying in the presence of inner-core lightning, the size and optical energy of the flashes generally increases as intensification proceeds. It is hypothesized that this relationship is governed by changes in the TC-scale circulation that produce enhanced upper-level outflow and large, strongly helical convective updrafts in the inner core. These changes increase the quantity and spatial extent of charged upper-level ice, which allows large intracloud flashes to occur when electric field breakdown occurs. This presentation will present objective diagnostics of convective asymmetry using azimuthal wavenumber decomposition of GOES-R Advanced Baseline Imager brightness temperatures. The degree of symmetry will be related to GLM flash characteristics to better understand how changes in the symmetry of TC inner-core convection relates to flash rate, size, and energy. It will also show how the evolution of inner-core asymmetries relates to the spatial distribution and characteristics of lightning flashes in the inner core.

Patrick Duran

Self-Aligned Focusing Schlieren and OH Planar Laser-Induced Fluorescence Flow Visualization in a Dual-Mode Scramjet

Ahigh-speed self-aligned focusing schlieren (SAFS) system was used to visualize density gradients in and around the cavity flameholder of the combustor section of the University of Virginia Supersonic Combustion Facility (UVASCF). Images with this system were acquired at a framing rate of 110 kHz with no fuel injection, with fuel injection but no flame, and for fuel injection with combustion corresponding to a global equivalence ratio of 𝜙 = 0.18. Images with an air throttle in operation to modify the shock train location with fuel injection and with flame were also acquired. Simultaneous OH planar laser-induced fluorescence (PLIF) images were also acquired at a framing rate of 20 Hz. Results obtained with both visualization techniques are compared to one another to highlight how SAFS can complement more advanced flow visualization techniques and resolve dynamic behavior that may not otherwise be captured. Both proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) analysis techniques are applied to the SAFS image sequences to identify coherent periodic structures for the runs with fuel injection and combustion.

Brett F Bathel

Vegetation Biochemistry: What Can Imaging Spectrometry Tell Us About Canopies?

Changes in ecosystem processes such as productivity and decomposition may be expressed in the canopy foliar chemistry resulting from altered carbon allocation patterns, metabolic processes and nutrient availability. Understanding carbon balance on land over large regions requires quantitative determination of leaf constituents such as lignin and total nitrogen from remote sensing imaging systems. Results from spectral reflectance measurements of stacked leaves in the laboratory show that spectrum matching techniques are applicable to the derivation of the equivalent liquid water thickness in plants as well as to the extraction of dry leaf matter reflectance spectra from spectra of green leaves. The residual spectra derived by subtracting water spectra from the spectra of green leaves shows a feature at 1.72 micrometers that can be related to the lignin content of the leaves. Oak leaves have a deeper residual absorption feature than do cotton leaves which is consistent with their relative lignin content. Similar results are achieved when deriving the residuals from images taken over areas of grass and pine trees. Imaging spectrometry provides promise in developing images of various foliar biochemical constituents.

Goetz, Alexander F. H.

Subspace-based Background Subtraction Applied to Aeroacoustic Wind Tunnel Testing

A subspace-based form of background subtraction is presented and applied to aeroacoustic wind tunnel data. A variant of this method has seen use in other fields such as climatology and medical imaging. The technique is based on an eigenvalue decomposition of the background noise cross-spectral matrix. Simulated results indicate similar performance to conventional background subtraction when the subtracted spectra are weaker than the true contaminating background levels. Superior performance is observed when the subtracted spectra are stronger than the true contaminating background levels, and when background data do not match between measurements. Experimental results show limited success in recovering signal behavior for data in which conventional background subtraction fails. The results also demonstrate the subspace subtraction technique's ability to maintain a physical coherence relationship in the modified cross-spectral matrix. Deconvolution results from microphone phased array data indicate that array integration methods are largely insensitive to subtraction type, and that background subtraction with appropriate background data is an effective alternative to diagonal removal.

deconvolution

Hierarchical image coding with diamond-shaped sub-bands

We present a sub-band image coding/decoding system using a diamond-shaped pyramid frequency decomposition to more closely match visual sensitivities than conventional rectangular bands. Filter banks are composed of simple, low order IIR components. The coder is especially designed to function in a multiple resolution reconstruction setting, in situations such as variable capacity channels or receivers, where images must be reconstructed without the entire pyramid of sub-bands. We use a nonlinear interpolation technique for lost subbands to compensate for loss of aliasing cancellation.

Li, Xiaohui

Stereo Imaging Velocimetry Technique Using Standard Off-the-Shelf CCD Cameras

Stereo imaging velocimetry is a fluid physics technique for measuring three-dimensional (3D) velocities at a plurality of points. This technique provides full-field 3D analysis of any optically clear fluid or gas experiment seeded with tracer particles. Unlike current 3D particle imaging velocimetry systems that rely primarily on laser-based systems, stereo imaging velocimetry uses standard off-the-shelf charge-coupled device (CCD) cameras to provide accurate and reproducible 3D velocity profiles for experiments that require 3D analysis. Using two cameras aligned orthogonally, we present a closed mathematical solution resulting in an accurate 3D approximation of the observation volume. The stereo imaging velocimetry technique is divided into four phases: 3D camera calibration, particle overlap decomposition, particle tracking, and stereo matching. Each phase is explained in detail. In addition to being utilized for space shuttle experiments, stereo imaging velocimetry has been applied to the fields of fluid physics, bioscience, and colloidal microscopy.

McDowell, Mark

ICER-3D Hyperspectral Image Compression Software

Software has been developed to implement the ICER-3D algorithm. ICER-3D effects progressive, three-dimensional (3D), wavelet-based compression of hyperspectral images. If a compressed data stream is truncated, the progressive nature of the algorithm enables reconstruction of hyperspectral data at fidelity commensurate with the given data volume. The ICER-3D software is capable of providing either lossless or lossy compression, and incorporates an error-containment scheme to limit the effects of data loss during transmission. The compression algorithm, which was derived from the ICER image compression algorithm, includes wavelet-transform, context-modeling, and entropy coding subalgorithms. The 3D wavelet decomposition structure used by ICER-3D exploits correlations in all three dimensions of sets of hyperspectral image data, while facilitating elimination of spectral ringing artifacts, using a technique summarized in "Improving 3D Wavelet-Based Compression of Spectral Images" (NPO-41381), NASA Tech Briefs, Vol. 33, No. 3 (March 2009), page 7a. Correlation is further exploited by a context-modeling subalgorithm, which exploits spectral dependencies in the wavelet-transformed hyperspectral data, using an algorithm that is summarized in "Context Modeler for Wavelet Compression of Hyperspectral Images" (NPO-43239), which follows this article. An important feature of ICER-3D is a scheme for limiting the adverse effects of loss of data during transmission. In this scheme, as in the similar scheme used by ICER, the spatial-frequency domain is partitioned into rectangular error-containment regions. In ICER-3D, the partitions extend through all the wavelength bands. The data in each partition are compressed independently of those in the other partitions, so that loss or corruption of data from any partition does not affect the other partitions. Furthermore, because compression is progressive within each partition, when data are lost, any data from that partition received prior to the loss can be used to reconstruct that partition at lower fidelity. By virtue of the compression improvement it achieves relative to previous means of onboard data compression, this software enables (1) increased return of hyperspectral scientific data in the presence of limits on the rates of transmission of data from spacecraft to Earth via radio communication links and/or (2) reduction in spacecraft radio-communication power and/or cost through reduction in the amounts of data required to be downlinked and stored onboard prior to downlink. The software is also suitable for compressing hyperspectral images for ground storage or archival purposes.

Xie, Hua

Estimating key forest ecosystem parameters through remote sensing

Forest canopy chemistry and biomass indicators of ecosystem photosynthesis and decomposition processes are presently studied in view of Airborne Imaging Spectrometer data, which generated spectra from averaged 3 x 3 pixel areas for each of 20 sites for mutual qualitative comparison. Vegetation spectra were strongly differentiated from other cover types by an apparent absorption feature at 1500-1700 nm. Preliminary work with stepwise regression suggests that lignin may play a role in canopy reflectance, and that there is potential for remote detection of forest canopy lignin.

Wessman, Carol A.

Sub-band/transform compression of video sequences

The progress on compression of video sequences is discussed. The overall goal of the research was the development of data compression algorithms for high-definition television (HDTV) sequences, but most of our research is general enough to be applicable to much more general problems. We have concentrated on coding algorithms based on both sub-band and transform approaches. Two very fundamental issues arise in designing a sub-band coder. First, the form of the signal decomposition must be chosen to yield band-pass images with characteristics favorable to efficient coding. A second basic consideration, whether coding is to be done in two or three dimensions, is the form of the coders to be applied to each sub-band. Computational simplicity is of essence. We review the first portion of the year, during which we improved and extended some of the previous grant period's results. The pyramid nonrectangular sub-band coder limited to intra-frame application is discussed. Perhaps the most critical component of the sub-band structure is the design of bandsplitting filters. We apply very simple recursive filters, which operate at alternating levels on rectangularly sampled, and quincunx sampled images. We will also cover the techniques we have studied for the coding of the resulting bandpass signals. We discuss adaptive three-dimensional coding which takes advantage of the detection algorithm developed last year. To this point, all the work on this project has been done without the benefit of motion compensation (MC). Motion compensation is included in many proposed codecs, but adds significant computational burden and hardware expense. We have sought to find a lower-cost alternative featuring a simple adaptation to motion in the form of the codec. In sequences of high spatial detail and zooming or panning, it appears that MC will likely be necessary for the proposed quality and bit rates.

Sauer, Ken

Ionospheric Imaging From a Low Earth Orbiter Tracking GPS

Tomographic imaging of the ionosphere is examined using singular value decomposition analysis. The interdependency of the obtainable resolution, the accuracy of the solution and the noise data is explained. A simulation is illustrated where a true 2-D ionospheric structure is generated, and a tomographic inversion of the structure is carried out. Inclusion of data taken in an occultation geometry reveals the strength that is added by these data. The effect of the ionosphere on a GPS signal as viewed by a user in space is examined. Under somewhat strong solar conditions, the absolute bending of the signal is of the order of 0.01 degrees for the L1 signal; the phase advance can be as large as 90 meters.

ionosphere

Analysis and visualization of single-trial event-related potentials

In this study, a linear decomposition technique, independent component analysis (ICA), is applied to single-trial multichannel EEG data from event-related potential (ERP) experiments. Spatial filters derived by ICA blindly separate the input data into a sum of temporally independent and spatially fixed components arising from distinct or overlapping brain or extra-brain sources. Both the data and their decomposition are displayed using a new visualization tool, the "ERP image," that can clearly characterize single-trial variations in the amplitudes and latencies of evoked responses, particularly when sorted by a relevant behavioral or physiological variable. These tools were used to analyze data from a visual selective attention experiment on 28 control subjects plus 22 neurological patients whose EEG records were heavily contaminated with blink and other eye-movement artifacts. Results show that ICA can separate artifactual, stimulus-locked, response-locked, and non-event-related background EEG activities into separate components, a taxonomy not obtained from conventional signal averaging approaches. This method allows: (1) removal of pervasive artifacts of all types from single-trial EEG records, (2) identification and segregation of stimulus- and response-locked EEG components, (3) examination of differences in single-trial responses, and (4) separation of temporally distinct but spatially overlapping EEG oscillatory activities with distinct relationships to task events. The proposed methods also allow the interaction between ERPs and the ongoing EEG to be investigated directly. We studied the between-subject component stability of ICA decomposition of single-trial EEG epochs by clustering components with similar scalp maps and activation power spectra. Components accounting for blinks, eye movements, temporal muscle activity, event-related potentials, and event-modulated alpha activities were largely replicated across subjects. Applying ICA and ERP image visualization to the analysis of sets of single trials from event-related EEG (or MEG) experiments can increase the information available from ERP (or ERF) data. Copyright 2001 Wiley-Liss, Inc.

Non-NASA Center

Spatially-varying IIR filter banks for image coding

This paper reports on the application of spatially variant infinite impulse response (IIR) filter banks to subband image coding. The new filter bank is based on computationally efficient recursive polyphase decompositions that dynamically change in response to the input signal. In the absence of quantization, reconstruction can be made exact. However, by proper choice of an adaptation scheme, we show that subband image coding based on time varying filter banks can yield improvement over the use of conventional filter banks.

Chung, Wilson C.

Three-Dimensional Instantaneous Flow-field Reconstruction Using Planar Spectral Proper Orthogonal Decomposition

A novel method for three-dimensional flowfield reconstruction is presented that utilizes spectral proper orthogonal decomposition (SPOD). The time-resolved PIV data are acquired in multiple two-dimensional spanwise planes, with adequate streamwise spacing to resolve the out-of-plane wavelength. Knowledge of the out-of-plane phase speed is also required, and for the present case this value was obtained using instantaneous hot-wire measurements. The method is validated by comparing phase-averaged results acquired with a forcing signal to the SPOD reconstructed results using the same data. The method is also applied to data with no forcing applied, and the results are shown to be similar to the forced case. This method is particularly useful for cases in which there are multiple frequencies of interest, or cases in which the phenomena occur over a broad range of frequencies. This method allows reconstruction of the frequencies of interest without the need for a forcing signal, so the flow can be studied in its natural (undisturbed) state.

boundary-layer transition

Robot Acting on Moving Bodies (RAMBO): Interaction with tumbling objects

Interaction with tumbling objects will become more common as human activities in space expand. Attempting to interact with a large complex object translating and rotating in space, a human operator using only his visual and mental capacities may not be able to estimate the object motion, plan actions or control those actions. A robot system (RAMBO) equipped with a camera, which, given a sequence of simple tasks, can perform these tasks on a tumbling object, is being developed. RAMBO is given a complete geometric model of the object. A low level vision module extracts and groups characteristic features in images of the object. The positions of the object are determined in a sequence of images, and a motion estimate of the object is obtained. This motion estimate is used to plan trajectories of the robot tool to relative locations rearby the object sufficient for achieving the tasks. More specifically, low level vision uses parallel algorithms for image enhancement by symmetric nearest neighbor filtering, edge detection by local gradient operators, and corner extraction by sector filtering. The object pose estimation is a Hough transform method accumulating position hypotheses obtained by matching triples of image features (corners) to triples of model features. To maximize computing speed, the estimate of the position in space of a triple of features is obtained by decomposing its perspective view into a product of rotations and a scaled orthographic projection. This allows use of 2-D lookup tables at each stage of the decomposition. The position hypotheses for each possible match of model feature triples and image feature triples are calculated in parallel. Trajectory planning combines heuristic and dynamic programming techniques. Then trajectories are created using dynamic interpolations between initial and goal trajectories. All the parallel algorithms run on a Connection Machine CM-2 with 16K processors.

Davis, Larry S.

Robot acting on moving bodies (RAMBO): Preliminary results

A robot system called RAMBO is being developed. It is equipped with a camera, which, given a sequence of simple tasks, can perform these tasks on a moving object. RAMBO is given a complete geometric model of the object. A low level vision module extracts and groups characteristic features in images of the object. The positions of the object are determined in a sequence of images, and a motion estimate of the object is obtained. This motion estimate is used to plan trajectories of the robot tool to relative locations nearby the object sufficient for achieving the tasks. More specifically, low level vision uses parallel algorithms for image enchancement by symmetric nearest neighbor filtering, edge detection by local gradient operators, and corner extraction by sector filtering. The object pose estimation is a Hough transform method accumulating position hypotheses obtained by matching triples of image features (corners) to triples of model features. To maximize computing speed, the estimate of the position in space of a triple of features is obtained by decomposing its perspective view into a product of rotations and a scaled orthographic projection. This allows the use of 2-D lookup tables at each stage of the decomposition. The position hypotheses for each possible match of model feature triples and image feature triples are calculated in parallel. Trajectory planning combines heuristic and dynamic programming techniques. Then trajectories are created using parametric cubic splines between initial and goal trajectories. All the parallel algorithms run on a Connection Machine CM-2 with 16K processors.

Davis, Larry S.