Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiresolution”

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 91 records · Page 5

Multiscale Morphological Filtering for Analysis of Noisy and Complex Images

Images acquired with passive sensing techniques suffer from illumination variations and poor local contrasts that create major difficulties in interpretation and identification tasks. On the other hand, images acquired with active sensing techniques based on monochromatic illumination are degraded with speckle noise. Mathematical morphology offers elegant techniques to handle a wide range of image degradation problems. Unlike linear filters, morphological filters do not blur the edges and hence maintain higher image resolution. Their rich mathematical framework facilitates the design and analysis of these filters as well as their hardware implementation. Morphological filters are easier to implement and are more cost effective and efficient than several conventional linear filters. Morphological filters to remove speckle noise while maintaining high resolution and preserving thin image regions that are particularly vulnerable to speckle noise were developed and applied to SAR imagery. These filters used combination of linear (one-dimensional) structuring elements in different (typically four) orientations. Although this approach preserves more details than the simple morphological filters using two-dimensional structuring elements, the limited orientations of one-dimensional elements approximate the fine details of the region boundaries. A more robust filter designed recently overcomes the limitation of the fixed orientations. This filter uses a combination of concave and convex structuring elements. Morphological operators are also useful in extracting features from visible and infrared imagery. A multiresolution image pyramid obtained with successive filtering and a subsampling process aids in the removal of the illumination variations and enhances local contrasts. A morphology-based interpolation scheme was also introduced to reduce intensity discontinuities created in any morphological filtering task. The generality of morphological filtering techniques in extracting information from a wide variety of images obtained with active and passive sensing techniques is discussed. Such techniques are particularly useful in obtaining more information from fusion of complex images by different sensors such as SAR, visible, and infrared.

Kher, A.↗

Analysis of SSME inspection imagery using AI approaches

The automated analysis of SSME injector assemblies has been investigated for the cases of LOX post surface defects and injector-baffle deterioration. Defects are isolated via 2D feature extraction from borescope and camera images; temporal-frequency transforms are then used to create a multiresolution set of feature vectors representing image contents. The potential flaws thus discriminated are then segmented and classified according to known categories. AI is applied in the form of a blackboard architecture that is controlled by a rule-based production system.

Finegan, Michael K., Jr.↗

Hierarchical multisensor analysis for robotic exploration

Robotic vehicles for lunar and Mars exploration will carry an array of complex instruments requiring real-time data interpretation and fusion. The system described here uses hierarchical multiresolution analysis of visible and multispectral images to extract information on mineral composition, texture and object shape. This information is used to characterize the site geology and choose interesting samples for acquisition. Neural networks are employed for many data analysis steps. A decision tree progressively integrates information from multiple instruments and performs goal-driven decision making. The system is designed to incorporate more instruments and data types as they become available.

Eberlein, Susan↗

Engineering a visual system for seeing through fog

We examine the requirements for on-board aircraft sensor systems that would allow pilots to 'see through' poor weather, especially fog, and land and rollout aircraft under conditions that currently cause flight cancellations and airport closures. Three visual aspects of landing and rollout are distinguished: guidance, hazard detection and hazard recognition. The visual features which support the tasks are discussed. Three broad categories of sensor technology are examined: passive millimeter wave (PMMW), imaging radar, and passive infrared (IR). PMMW and imaging radar exhibit good weather penetration, but poor spatial and temporal resolution. Imaging radar exhibits good weather penetration, but typically relies on a flat-earth assumption which can lead to interpretive errors. PMMW systems have a narrow field of view. IR has poorer weather penetration but good spatial resolution. We recommend using both millimeter-wave and infrared sensor systems, blending the images using multiresolution digital-image pyramid-processing technology, and fusing the resulting real-time images with stored database imagery of the same scene.

Larimer, J.↗

Remote sensing of multilevel clouds during FIRE IFO 2

An unresolved difficulty in the remote sensing of clouds concerns the inability of the cloud retrieval algorithms to adequately recognize and analyze scenes containing overlapping cloud layers. Most cloud retrieval schemes, such as that used by the International Satellite Cloud Climatology Project (ISCCP) assume that each picture element (pixel) contains a single cloud layer. The current study begins to address the complexities of multilayered cloud property retrieval through the application of a modified multispectral, multiresolution (MSMR) method, first detailed in Baum et al. (1992), which merges 1.1 -km (at nadir) spectral data from the Advanced Very High Resolution Radiometer (AVHRR) with 17.4-km (at nadir) High Resolution Infrared Radiometer Sounder (HIRS/2, henceforth HIRS). Both instruments are flown aboard the National Oceanic and Atmospheric Administration (NOAA) polar-orbiting platforms. An ideal case study for this investigation is provided by the NOAA-11 overpass at 20:48 UTC on November 28, 1991. At this time, a large-scale cirrostratus veil overlaid a low-level stratus deck over much of the IFO region. There were both surface lidar and radar observations of the clouds as well as University of North Dakota (UND) Citation aircraft measurements. The presence of overlapping cloud layers within a HIRS FOV is determined from colocated AVHRR spectral data through the use of a fuzzy logic expert system. Conventional algorithms such as spatial coherence and CO2 slicing are used to retrieve cloud pressure and height for each identified cloud layer. The results from the satellite cloud retrieval analysis are compared to results from both surface- and aircraft-based measurements.

Baum, Bryan A.↗

Atmospheric parameterization schemes for satellite cloud property retrieval during FIRE IFO 2

Satellite cloud retrieval algorithms generally require atmospheric temperature and humidity profiles to determine such cloud properties as pressure and height. For instance, the CO2 slicing technique called the ratio method requires the calculation of theoretical upwelling radiances both at the surface and a prescribed number (40) of atmospheric levels. This technique has been applied to data from, for example, the High Resolution Infrared Radiometer Sounder (HIRS/2, henceforth HIRS) flown aboard the NOAA series of polar orbiting satellites and the High Resolution Interferometer Sounder (HIS). In this particular study, four NOAA-11 HIRS channels in the 15-micron region are used. The ratio method may be applied to various channel combinations to estimate cloud top heights using channels in the 15-mu m region. Presently, the multispectral, multiresolution (MSMR) scheme uses 4 HIRS channel combination estimates for mid- to high-level cloud pressure retrieval and Advanced Very High Resolution Radiometer (AVHRR) data for low-level (is greater than 700 mb) cloud level retrieval. In order to determine theoretical upwelling radiances, atmospheric temperature and water vapor profiles must be provided as well as profiles of other radiatively important gas absorber constituents such as CO2, O3, and CH4. The assumed temperature and humidity profiles have a large effect on transmittance and radiance profiles, which in turn are used with HIRS data to calculate cloud pressure, and thus cloud height and temperature. For large spatial scale satellite data analysis, atmospheric parameterization schemes for cloud retrieval algorithms are usually based on a gridded product such as that provided by the European Center for Medium Range Weather Forecasting (ECMWF) or the National Meteorological Center (NMC). These global, gridded products prescribe temperature and humidity profiles for a limited number of pressure levels (up to 14) in a vertical atmospheric column. The FIRE IFO 2 experiment provides an opportunity to investigate current atmospheric profile parameterization schemes, compare satellite cloud height results using both gridded products (ECMWF) and high vertical resolution sonde data from the National Weather Service (NWS) and Cross Chain Loran Atmospheric Sounding System (CLASS), and suggest modifications in atmospheric parameterization schemes based on these results.

Titlow, James↗

Applications of squeezed states: Bogoliubov transformations and wavelets to the statistical mechanics of water and its bubbles

The squeezed states or Bogoliubov transformations and wavelets are applied to two problems in nonrelativistic statistical mechanics: the dielectric response of liquid water, epsilon(q-vector,w), and the bubble formation in water during insonnification. The wavelets are special phase-space windows which cover the domain and range of L(exp 1) intersection of L(exp 2) of classical causal, finite energy solutions. The multiresolution of discrete wavelets in phase space gives a decomposition into regions of time and scales of frequency thereby allowing the renormalization group to be applied to new systems in addition to the tired 'usual suspects' of the Ising models and lattice gasses. The Bogoliubov transformation: squeeze transformation is applied to the dipolaron collective mode in water and to the gas produced by the explosive cavitation process in bubble formation.

Defacio, Brian↗

On the wavelet optimized finite difference method

When one considers the effect in the physical space, Daubechies-based wavelet methods are equivalent to finite difference methods with grid refinement in regions of the domain where small scale structure exists. Adding a wavelet basis function at a given scale and location where one has a correspondingly large wavelet coefficient is, essentially, equivalent to adding a grid point, or two, at the same location and at a grid density which corresponds to the wavelet scale. This paper introduces a wavelet optimized finite difference method which is equivalent to a wavelet method in its multiresolution approach but which does not suffer from difficulties with nonlinear terms and boundary conditions, since all calculations are done in the physical space. With this method one can obtain an arbitrarily good approximation to a conservative difference method for solving nonlinear conservation laws.

Jameson, Leland↗

Multilevel cloud retrieval using multispectral HIRS and AVHRR data: Nighttime oceanic analysis

A multispectral, multiresolution (MSMR) method is developed for analyzing scenes of overlapping cloud layers. The MSMR method is applied to data from the NOAA 11 advanced very high resolution radiometer (AVHRR) and the high-resolution infrared radiometer sounder (HIRS-2). The data are from a nighttime oceanic scene in which a semitransparent cirrus veil overlays a large-scale stratus cloud. Low-cloud and clear-sky radiances are determined using a spatial coherence technique. Middle to upper level cloud pressures and radiances are estimated from HIRS-2 15 micrometer CO2 band radiometric data. The MSMR method improves the interpretation of a nighttime, oceanic scene containing thin cirrus over a large-scale stratiform cloud. If, for example, the same scene is analyzed using only the AVHRR 10.8 micrometer channel, the accompanying retrieved cloud heights are found to be between the cirrus and stratus cloud heights and are incorrectly identified as midlevel altostratus clouds. Theoretical radiative transfer model results for both water droplet spheres and randomly oriented hexagonal ice crystals are compared to observed AVHRR brightness temperature differences (BTD) between the 3.7- and 10.8 micrometer channels (BTD(sup 34)) and between the 10.8- and 12- micrometer channels (BTD(sup 45)) to distinguish among the effects of cloud optical depth, particle size, and phase for both single-layer clouds and overlapping two-layer clouds. Theoretical BTD calculations are used to estimate the range of effective particle sizes for eac h cloud layer. The data for the cirrus in the case study region near Bermuda are consistent with theoretical results for relatively small randomly oriented hexagonal ice crystals. The observed BTD(sup 34) and BTD(sup 45) values are lower for the cirrus above a lower-level cloud than for single-level cirrus with no underlying cloud. In certain cases the BTD analysis provides a way to distinguish between clouds composed of supercooled water droplets rather than ice particles. Analysis of nighttime data permits determination of stratus infrared optical depths smaller than 4.

Baum, Bryan A.↗

Deep HST imaging of distant weak radio and field galaxies

We present deep Hubble Space Telescope (HST) Wide-Field Camera (WFC) V- and I-band images of three distant weak radio galaxies with z = 0.311-2.390 and seven field galaxies with z = 0.131-0.58. The images were deconvolved with both the Lucy and multiresolution CLEAN methods, which yield a restoring Full Width at Half Maximum (FWHM) of less than or equal to 0.2 sec, (nearly) preserve photons and signal-to-noise ratio at low spatial frequencies, and produce consistent light profiles down to our 2 sigma surface brightness sensitivity limit of V approximately 27.2 and I approximately 25.9 mag/sq arcsec. Multi-component image modeling was used to provide deconvolution-independent estimates of structural parameters for symmetric galaxies. We present 12-band (m(sub 2750) UBVRIgriJHK) photometry for a subset of the galaxies and bootstrap the unknown FOC/48 zero point at 2750 A in three independent ways (yielding m(sub 2750) = 21.34 +/- 0.09 mag for 1.0 e(-)/s). Two radio galaxies with z = 0.311 and 0.528, as well as one field galaxy with z = 0.58, have the colors and spectra of early-type galaxies, and a(exp 1/4)-like light profiles in the HST images. The two at z greater than 0.5 have little or no color gradients in V - I and are likely giant ellipticals, while the z = 0.311 radio galaxy has a dim exponential disk and is likely an S0. Six of the seven field galaxies have light profiles that indicate (small) inner bulges following a(exp 1/4) laws and outer exponential disks, both with little or no color gradients. These are (early-type) spiral galaxies with z = 0.131-0.528. About half have faint companions or bars. One shows lumpy structure, possibly a merger. The compact narrow-line galaxy 53W002 at z = 2.390 has less than or = 30% +/- 10% of its HST V and I flux in the central kiloparsec (due to its weak Active Galactic Nucleus (AGN)). Most of its light (V approximately equal to 23.3) occurs in a symmetric envelope with a regular a(exp 1/4)-like profile of effective radius a approximately equal to 1.1 sec (approximately equal to 12 kpc for H(sub 0) = 50, q(sub 0) = 0.1. Its (HST) V - I color varies at most from approximately 0.3 mag at a approximately equal to 0.2 sec to approximately 1.2 mag at a approximately greater than 0.4 sec, and possibly to approximately greater than 2.2 mag at a approximately greater than 1.2 sec. Together with its I - K color (approximately equal to 2.5 mag for a approximately greater than 1.0 sec-2.0 sec), this is consistent with an aging stellar population approximately 0.3-0.5 Gyr old in the galaxy center (a approx. less than 2 kpc radius), and possibly approximately 0.5-1.0 Gyr old at a approximately greater than 10 kpc radius. While its outer part may thus have started to collapse at z = 2.5-4, its inner part still is aligned with its redshifted Ly(alpha) cloud and its radio axis, possibly caused by star formation associated with the radio jet, or by reflection from its AGN cone.

Windhorst, R. A.↗

On the Effective Construction of Compactly Supported Wavelets Satisfying Homogenous Boundary Conditions on the Interval

We construct compactly supported wavelet bases satisfying homogeneous boundary conditions on the interval (0,1). The maximum features of multiresolution analysis on the line are retained, including polynomial approximation and tree algorithms. The case of H(sub 0)(sup 1)(0, 1)is detailed, and numerical values, required for the implementation, are provided for the Neumann and Dirichlet boundary conditions.

Chiavassa, G.↗

Multilevel Algorithm for Atmospheric Data Assimilation

A multiscale algorithm for the problem of optimal statistical interpolation of observed data has been developed. This problem includes the calculation of the vector of the 'analyzed' (best estimated) atmosphere flow field w(sup a) by the formula: w(sup a) = w(sup f) + P(sup f) H(sup T) y, where the quantity y is defined by the equation (H P(sup f) H(sup T) + R)y = w(sup o) - H w(sup f), using the given model forecast first guess w(sup f) and the vector of observations w(sup o); H is an interpolation operator from the regular grid to the observation network, P(sup f) is the forecast error covariance matrix, and R is the observation error covariance matrix. At this initial stage the case of univariate analysis of single level radiosonde height data is considered. The matrix R is assumed to be diagonal, and the matrix P(sup f) is assumed to be given by the formula P(sub ij)(sup f) = sigma(sub i)(sup f) mu(sub ij) sigma(sub j)(sub f), where mu(sub ij) is a smooth, decreasing function of the distance between the i-th and the j-th points. In this paper we describe a multiscale iterative process based on a multiresolution, simultaneous displacement technique and a localized variational calculation of iteration parameters.

Brandt, Achi↗

A Subband Coding Method for HDTV

This paper introduces a new HDTV coder based on motion compensation, subband coding, and high order conditional entropy coding. The proposed coder exploits the temporal and spatial statistical dependencies inherent in the HDTV signal by using intra- and inter-subband conditioning for coding both the motion coordinates and the residual signal. The new framework provides an easy way to control the system complexity and performance, and inherently supports multiresolution transmission. Experimental results show that the coder outperforms MPEG-2, while still maintaining relatively low complexity.

Chung, Wilson↗

Advanced Signal Processing Methods Applied to Digital Mammography

The work reported here is on the extension of the earlier proposal of the same title, August 1994-June 1996. The report for that work is also being submitted. The work reported there forms the foundation for this work from January 1997 to September 1997. After the earlier work was completed there were a few items that needed to be completed prior to submission of a new and more comprehensive proposal for further research. Those tasks have been completed and two new proposals have been submitted, one to NASA, and one to Health & Human Services WS). The main purpose of this extension was to refine some of the techniques that lead to automatic large scale evaluation of full mammograms. Progress on each of the proposed tasks follows. Task 1: A multiresolution segmentation of background from breast has been developed and tested. The method is based on the different noise characteristics of the two different fields. The breast field has more power in the lower octaves and the off-breast field behaves similar to a wideband process, where more power is in the high frequency octaves. After the two fields are separated by lowpass filtering, a region labeling routine is used to find the largest contiguous region, the breast. Task 2: A wavelet expansion that can decompose the image without zero padding has been developed. The method preserves all properties of the power-of-two wavelet transform and does not add appreciably to computation time or storage. This work is essential for analysis of the full mammogram, as opposed to selecting sections from the full mammogram. Task 3: A clustering method has been developed based on a simple counting mechanism. No ROC analysis has been performed (and was not proposed), so we cannot finally evaluate this work without further support. Task 4: Further testing of the filter reveals that different wavelet bases do yield slightly different qualitative results. We cannot provide quantitative conclusions about this for all possible bases without further support. Task 5: Better modeling does indeed make an improvement in the detection output. After the proposal ended, we came up with some new theoretical explanations that helps in understanding when the D4 filter should be better. This work is currently in the review process. Task 6: N/A. This no longer applies in view of Tasks 4-5. Task 7: Comprehensive plans for further work have been completed. These plans are the subject of two proposals, one to NASA and one to HHS. These proposals represent plans for a complete evaluation of the methods for identifying normal mammograms, augmented with significant further theoretical work.

Stauduhar, Richard P.↗

The Wavelet Element Method: Realization and Additional Features in 2D and 3D - Part 2

The Wavelet Element Method (WEM) provides a construction of multiresolution systems and biorthogonal wavelets on fairly general domains. These are split into subdomains that are mapped to a single reference hypercube. Tensor products of scaling functions and wavelets defined on the unit interval are used on the reference domain. By introducing appropriate matching conditions across the interelement boundaries, a globally continuous biorthogonal wavelet basis on the general domain is obtained. This construction does not uniquely define the basis functions but rather leaves some freedom for fulfilling additional features. In this paper we detail the general construction principle of the WEM to the 1D, 2D and 3D cases. We address additional features such as symmetry, vanishing moments and minimal support of the wavelet functions in each particular dimension. The construction is illustrated by using biorthogonal spline wavelets on the interval.

Canuto, Claudio↗

Automatic registration of satellite imagery

Image registration is one of the basic image processing operations in remote sensing. With the increase in the number of images collected every day from different sensors, automated registration of multi-sensor/multi-spectral images has become an important issue. A wide range of registration techniques has been developed for many different types of applications and data. The objective of this paper is to present an automatic registration algorithm which uses a multiresolution analysis procedure based upon the wavelet transform. The procedure is completely automatic and relies on the grey level information content of the images and their local wavelet transform modulus maxima. The registration algorithm is very simple and easy to apply because it needs basically one parameter. We have obtained very encouraging results on test data sets from the TM and SPOT sensor images of forest, urban and agricultural areas.

Fonseca, Leila M. G.↗

On the Spatial Distribution of High Velocity Al-26 Near the Galactic Center

We present results of simulations of the distribution of 1809 keV radiation from the decay of Al-26 in the Galaxy. Recent observations of this emission line using the Gamma Ray Imaging Spectrometer (GRIS) have indicated that the bulk of the AL-26 must have a velocity of approx. 500 km/ s. We have previously shown that a velocity this large could be maintained over the 10(exp 6) year lifetime of the Al-26 if it is trapped in dust grains that are reaccelerated periodically in the ISM. Here we investigate whether a dust grain velocity of approx. 500 km/ s will produce a distribution of 1809 keV emission in latitude that is consistent with the narrow distribution seen by COMPTEL. We find that dust grain velocities in the range 275 - 1000 km/ s are able to reproduce the COMPTEL 1809 keV emission maps reconstructed using the Richardson-Lucy and Maximum Entropy image reconstruction methods while the emission map reconstructed using the Multiresolution Regularized Expectation Maximization algorithm is not well fit by any of our models. The Al-26 production rate that is needed to reproduce the observed 1809 keV intensity yields in a Galactic mass of Al-26 of approx. 1.5 - 2 solar mass which is in good agreement with both other observations and theoretical production rates.

Sturner, Steven J.↗

BOREAS AFM-12 1-km AVHRR Seasonal Land Cover Classification

The Boreal Ecosystem-Atmosphere Study (BOREAS) Airborne Fluxes and Meteorology (AFM)-12 team's efforts focused on regional scale Surface Vegetation and Atmosphere (SVAT) modeling to improve parameterization of the heterogeneous BOREAS landscape for use in larger scale Global Circulation Models (GCMs). This regional land cover data set was developed as part of a multitemporal one-kilometer Advanced Very High Resolution Radiometer (AVHRR) land cover analysis approach that was used as the basis for regional land cover mapping, fire disturbance-regeneration, and multiresolution land cover scaling studies in the boreal forest ecosystem of central Canada. This land cover classification was derived by using regional field observations from ground and low-level aircraft transits to analyze spectral-temporal clusters that were derived from an unsupervised cluster analysis of monthly Normalized Difference Vegetation Index (NDVI) image composites (April-September 1992). This regional data set was developed for use by BOREAS investigators, especially those involved in simulation modeling, remote sensing algorithm development, and aircraft flux studies. Based on regional field data verification, this multitemporal one-kilometer AVHRR land cover mapping approach was effective in characterizing the biome-level land cover structure, embedded spatially heterogeneous landscape patterns, and other types of key land cover information of interest to BOREAS modelers.The land cover mosaics in this classification include: (1) wet conifer mosaic (low, medium, and high tree stand density), (2) mixed coniferous-deciduous forest (80% coniferous, codominant, and 80% deciduous), (3) recent visible bum, vegetation regeneration, or rock outcrops-bare ground-sparsely vegetated slow regeneration bum (four classes), (4) open water and grassland marshes, and (5) general agricultural land use/ grasslands (three classes). This land cover mapping approach did not detect small subpixel-scale landscape features such as fens, bogs, and small water bodies. Field observations and comparisons with Landsat Thematic Mapper (TM) suggest a minimum effective resolution of these land cover classes in the range of three to four kilometers, in part, because of the daily to monthly compositing process. In general, potential accuracy limitations are mitigated by the use of conservative parameterization rules such as aggregation of predominant land cover classes within minimum horizontal grid cell sizes of ten kilometers. The AFM-12 one-kilometer AVHRR seasonal land cover classification data are available from the Earth Observing System Data and Information System (EOSDIS) Oak Ridge National Laboratory (ORNL) Distributed Active Archive Center (DAAC). The data files are available on a CD-ROM (see document number 20010000884).

Steyaert, Lou↗