Engineering Papers⌕ Search

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

Unsupervised learning-enabled pulsed infrared thermographic microscopy of subsurface defects in stainless steel

Metallic structures produced with laser powder bed fusion (LPBF) additive manufacturing method (AM) frequently contain microscopic porosity defects, with typical approximate size distribution from one to 100 microns. Presence of such defects could lead to premature failure of the structure. In principle, structural integrity assessment of LPBF metals can be accomplished with nondestructive evaluation (NDE). Pulsed infrared thermography (PIT) is a non-contact, one-sided NDE method that allows for imaging of internal defects in arbitrary size and shape metallic structures using heat transfer. PIT imaging is performed using compact instrumentation consisting of a flash lamp for deposition of a heat pulse, and a fast frame infrared (IR) camera for measuring surface temperature transients. However, limitations of imaging resolution with PIT include blurring due to heat diffusion, sensitivity limit of the IR camera. We demonstrate enhancement of PIT imaging capability with unsupervised learning (UL), which enables PIT microscopy of subsurface defects in high strength corrosion resistant stainless steel 316 alloy. PIT images were processed with UL spatial–temporal separation-based clustering segmentation (STSCS) algorithm, refined by morphology image processing methods to enhance visibility of defects. The STSCS algorithm starts with wavelet decomposition to spatially de-noise thermograms, followed by UL principal component analysis (PCA), fine-tuning optimization, and neural learning-based independent component analysis (ICA) algorithms to temporally compress de-noised thermograms. The compressed thermograms were further processed with UL-based graph thresholding K-means clustering algorithm for defects segmentation. The STSCS algorithm also includes online learning feature for efficient re-training of the model with new data. For this study, metallic specimens with calibrated microscopic flat bottom hole defects, with diameters in the range from 203 to 76 µm, were produced using electro discharge machining (EDM) drilling. While the raw thermograms do not show any material defects, using STSCS algorithm to process PIT images reveals defects as small as 101 µm in diameter. To the best of our knowledge, this is the smallest reported size of a sub-surface defect in a metal imaged with PIT, which demonstrates the PIT capability of detecting defects in the size range relevant to quality control requirements of LPBF-printed high-strength metals.

36 MATERIALS SCIENCE↗

Subband/transform functions for image processing

Functions for image data processing written for use with the MATLAB(TM) software package are presented. These functions provide the capability to transform image data with block transformations (such as the Walsh Hadamard) and to produce spatial frequency subbands of the transformed data. Block transforms are equivalent to simple subband systems. The transform coefficients are reordered using a simple permutation to give subbands. The low frequency subband is a low resolution version of the original image, while the higher frequency subbands contain edge information. The transform functions can be cascaded to provide further decomposition into more subbands. If the cascade is applied to all four of the first stage subbands (in the case of a four band decomposition), then a uniform structure of sixteen bands is obtained. If the cascade is applied only to the low frequency subband, an octave structure of seven bands results. Functions for the inverse transforms are also given. These functions can be used for image data compression systems. The transforms do not in themselves produce data compression, but prepare the data for quantization and compression. Sample quantization functions for subbands are also given. A typical compression approach is to subband the image data, quantize it, then use statistical coding (e.g., run-length coding followed by Huffman coding) for compression. Contour plots of image data and subbanded data are shown.

Glover, Daniel↗

Geometric Interpretation of the Cluster Location Problem Part II: Application to the Pahala, Hawaii, Earthquake Sequence

In the companion “Theory” article, we presented a new framing of the seismic location problem in terms of differential geometry (Harris et al., 2025). From that viewpoint, we developed a “project and correct” approach for estimating the relative locations of earthquakes. Here, in this study, we use project and correct to estimate high-precision relative locations of events from an earthquake sequence beneath the town of Pahala, Hawaii, using high-precision correlation-derived picks. The sequence was active from 2020 through 2022 and produced many highly correlated signals at Hawaii Volcano Observatory (HVO) stations on the island of Hawaii. The data we inverted consisted of 2882 events with observations at 5 HVO stations. For comparison with the travel-time image, we also produced conventional hypocenter solutions using both the Bayesloc program (Myers et al., 2007, 2009) and a purpose-built double-difference code. There were obvious structural elements in the resulting image, the resolution of which we used to test the performance of the project and the correct algorithm. For the projection step, we first produced a 3D local basis using an singular value decomposition (SVD) of the 2882 groups of times. Projection of the travel-time vectors into this basis resulted in an image with structures similar to those produced by our conventional locators, but with distortion as predicted by theory. Removing the distortion requires an inverse operator generated from the metric tensor at the geometric centroid of the events. We compared two approaches to obtaining such an inverse operator. The first uses an estimate of the geographic centroid of the event cloud from the centroid of the travel-time data. The second approach uses the centroid of the conventionally produced locations. The first approach produces a corrected image very similar to the conventional results, but with a rotation. The corrected image produced using the conventionally derived centroid is a near-exact match to the conventional locations.

Dodge, Douglas A. [Lawrence Livermore National Lab↗

Does manganese influence grass litter decomposition on a Hawaiian rainfall gradient?

Plant litter is a well-defined pool of organic matter (OM) in which the influence of manganese (Mn) on decomposition (both decomposition rate and the mix of compounds ultimately transferred to soil OM) has been clearly demonstrated in temperate forests. However, no similar study exists on grasslands and the effect of foliar Mn versus soil-derived Mn on litter decomposition is poorly known. We used a 5-month and 12-month field, and 10-month laboratory experiments to evaluate litter decomposition on the Kohala rainfall gradient (Island of Hawai‘i) in areas with different foliar and soil Mn abundances, and on which a single plant species (Pennisetum clandestinum) dominates primary production and the litter pool. The chemical imaging analyses of decomposed litter revealed that Mn 2+ oxidized to Mn 3+ and Mn 4+ on grass litter during decompositions—hallmarks of Mn-driven litter oxidation. However, these transformations and Mn abundance did not predict greater litter mass loss through decomposition. These observations demonstrate that the importance of Mn to an ecosystem’s C cycle does not rely solely on the metal’s abundance and availability.

54 ENVIRONMENTAL SCIENCES↗

Progressive Vector Quantization on a massively parallel SIMD machine with application to multispectral image data

A progressive vector quantization (VQ) compression approach is discussed which decomposes image data into a number of levels using full search VQ. The final level is losslessly compressed, enabling lossless reconstruction. The computational difficulties are addressed by implementation on a massively parallel SIMD machine. We demonstrate progressive VQ on multispectral imagery obtained from the Advanced Very High Resolution Radiometer instrument and other Earth observation image data, and investigate the trade-offs in selecting the number of decomposition levels and codebook training method.

Manohar, Mareboyana↗

Using AI Planning Techniques to Automatically Generate Image Processing Procedures: A Preliminary Report

This paper describes work on the Multimission VICAR Planner (MVP) system to automatically construct executable image processing procedures for custom image processing requests for the JPL Multimission Image Processing Lab (MIPL). This paper focuses on two issues. First, large search spaces caused by complex plans required the use of hand encoded control information. In order to address this in a manner similar to that used by human experts, MVP uses a decomposition-based planner to implement hierarchical/skeletal planning at the higher level and then uses a classical operator based planner to solve subproblems in contexts defined by the high-level decomposition.

Laboratory MIPL VICAR artificial intelligence AI↗

Model-Based Obstacle Detention From Image Sequences

This paper proposes a method for detecting obstacles on a runway from monocular image sequences. The surface of the runway is modeled as a plane and the model flow field corresponding to the runway is described by 8 coefficients. The set of 8 coefficients describing the initial model flow field is given by the data from Inertial Navigation System (INS). The uncertainties of this initial model flow field are estimated and used to obtain the accurate model flow field. The residual flow field after warping images by the accurate model flow field is computed by solving over determined gradient-based optical flow equations using a singular value decomposition (SVD). This SVD gives us a new insight of the uncertainties inherent in the optical flow computation. Those pixels with large residual flow vectors are considered as obstacles. Experimental results for two real image sequences are presented.

Sull, Sanghoon↗

Spectral decomposition of AVIRIS data

A set of techniques is presented that uses only information contained within a raw Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) scene to estimate and to remove additive components such as multiple scattering and instrument dark current. Multiplicative components (instrument gain, topographic modulation of brightness, atmospheric transmission) can then be normalized, permitting enhancement, extraction, and identification of relative reflectance information related to surface composition and mineralogy. The technique for derivation of additive-component spectra from a raw AVIRIS scene is an adaption of the 'regression intersection method' of Crippen. This method uses two surface units that are spatially extensive, and located in rugged terrain. For a given wavelength pair, subtraction of the derived additive component from individual band values will remove topography in both regions in a band/band ratio image. Normalization of all spectra in the scene to the average scene spectrum then results in cancellation of multiplicative components and production of a relative-reflectance scene. The resulting AVIRIS product contains relative-reflectance features due to mineral absorption that depart from the average spectrum. These features commonly are extremely weak and difficult to recognize, but they can be enhanced by using two simple 3-D image-processing tools. The validity of these techniques will be demonstrated by comparisons between relative-reflectance AVIRIS spectra and those derived by using JPL standard calibrations. The AVIRIS data used in this analysis were acquired over the Kelso Dunes area (34 deg 55' N, 115 deg 43' W) of the eastern Mojave Desert, CA (in 1987) and the Upheaval Dome area (38 deg 27' N, 109 deg 55' W) of the Canyonlands National Park, UT (in 1991).

Gaddis, Lisa↗

Quantitative Comparison of the Variability in Observed and Simulated Shortwave Reflectance

The Climate Absolute Radiance and Refractivity Observatory (CLARREO) is a climate observation system that has been designed to monitor the Earth's climate with unprecedented absolute radiometric accuracy and SI traceability. Climate Observation System Simulation Experiments (OSSEs) have been generated to simulate CLARREO hyperspectral shortwave imager measurements to help define the measurement characteristics needed for CLARREO to achieve its objectives. To evaluate how well the OSSE-simulated reflectance spectra reproduce the Earth s climate variability at the beginning of the 21st century, we compared the variability of the OSSE reflectance spectra to that of the reflectance spectra measured by the Scanning Imaging Absorption Spectrometer for Atmospheric Cartography (SCIAMACHY). Principal component analysis (PCA) is a multivariate decomposition technique used to represent and study the variability of hyperspectral radiation measurements. Using PCA, between 99.7%and 99.9%of the total variance the OSSE and SCIAMACHY data sets can be explained by subspaces defined by six principal components (PCs). To quantify how much information is shared between the simulated and observed data sets, we spectrally decomposed the intersection of the two data set subspaces. The results from four cases in 2004 showed that the two data sets share eight (January and October) and seven (April and July) dimensions, which correspond to about 99.9% of the total SCIAMACHY variance for each month. The spectral nature of these shared spaces, understood by examining the transformed eigenvectors calculated from the subspace intersections, exhibit similar physical characteristics to the original PCs calculated from each data set, such as water vapor absorption, vegetation reflectance, and cloud reflectance.

Roberts, Yolanda, L.↗

Fluid-Thermal-Structural Interactions Induced by an Asymmetric Shock-Wave/Boundary-Layer Interaction in a Mach-6 Compression Corner

An experimental study is conducted of the fluid-thermal-structural interaction of a clamped compliant panel exposed to a three dimensional shock-wave/boundary-layer interaction (SWBLI) induced by a Mach-6 compression ramp with a spanwise nonuniform incoming boundary layer. The nonuniform boundary layer was produced by placing trips on one side of the upstream flat plate, resulting in largely turbulent flow on the tripped side and transitional flow on the untripped side. Measurements of the flowfield confirmed that the tripped boundary layer contained elevated levels of unsteadiness, and the SWBLI was observed to vary from attached to fully separated as the ramp angle was increased from 10◦ to 38◦; the separation region on the tripped side of the panel was noticeably smaller, showing the elevated turbulence levels of the tripped-side flow to remain relatively localized rather than diffusing across the whole model. Full-field, time-resolved panel deformations were measured using high-speed photogrammetry and the vibrational response at each compression angle was characterized. Although the measured modes conformed largely to those from classical clamped-plate theory, some skewing of the mode shapes was observed. IR thermography highlighted regions of the compliant region where elevated temperatures were likely to promote thermal softening effects to the transient panel response. The quasi-static deformation and stress field was used to characterize the internal stress factor of each mode and showed a meaningful relationship between transient panel response and stress contained within each mode: modes with antinodes lying in high-stress areas of the plate tended to exhibit increases in vibrational frequency and decreases in vibrational power, whereas the opposite was true for modes with antinodes in low-stress areas.

Spectral Proper Orthogonal Decomposition↗

Polarimetric radar data decomposition and interpretation

Significant efforts have been made to decompose polarimetric radar data into several simple scattering components. The components which are selected because of their physical significance can be used to classify SAR (Synthetic Aperture Radar) image data. If particular components can be related to forest parameters, inversion procedures may be developed to estimate these parameters from the scattering components. Several methods have been used to decompose an averaged Stoke's matrix or covariance matrix into three components representing odd (surface), even (double-bounce) and diffuse (volume) scatterings. With these decomposition techniques, phenomena, such as canopy-ground interactions, randomness of orientation, and size of scatters can be examined from SAR data. In this study we applied the method recently reported by van Zyl (1992) to decompose averaged backscattering covariance matrices extracted from JPL SAR images over forest stands in Maine, USA. These stands are mostly mixed stands of coniferous and deciduous trees. Biomass data have been derived from field measurements of DBH and tree density using allometric equations. The interpretation of the decompositions and relationships with measured stand biomass are presented in this paper.

Sun, Guoqing↗

High-resolution high-throughput spatiotemporal strain imaging reveals loss mechanisms in a surface acoustic wave device

Surface acoustic wave devices are key components for processing radio frequency signals in wireless communication because these devices offer simultaneously high performance, compact size and low cost. The optimization of the device structure requires a quantitative understanding of energy conversion and loss mechanisms. Here we use stroboscopic full-field diffraction x-ray microscopy to reveal an unanticipated acoustic loss in a prototypical one-port resonator device. A non-uniform acoustic excitation in the active area was responsible for the substantial end and side leakages observed at the design frequency. Quantitative analysis of the strain amplitude using a wave decomposition method allowed the determination of several key device parameters. This high-resolution high-throughput spatiotemporal strain imaging technique is more generally applicable to the study of dynamic strain modulation in nanoscale acoustic, electronic, optical and quantum devices. The high sensitivity allows precise measurement of the strain modulation with picometer-scale amplitude.

Zhou, Tao [ESRF - The European Synchrotron, Grenob↗

Application of Spectral Analysis Techniques in the Intercomparison of Aerosol Data. Part II: Using Maximum Covariance Analysis to Effectively Compare Spatiotemporal Variability of Satellite and AERONET Measured Aerosol Optical Depth

Moderate Resolution Imaging SpectroRadiometer (MODIS) and Multi-angle Imaging Spectroradiomater (MISR) provide regular aerosol observations with global coverage. It is essential to examine the coherency between space- and ground-measured aerosol parameters in representing aerosol spatial and temporal variability, especially in the climate forcing and model validation context. In this paper, we introduce Maximum Covariance Analysis (MCA), also known as Singular Value Decomposition analysis as an effective way to compare correlated aerosol spatial and temporal patterns between satellite measurements and AERONET data. This technique not only successfully extracts the variability of major aerosol regimes but also allows the simultaneous examination of the aerosol variability both spatially and temporally. More importantly, it well accommodates the sparsely distributed AERONET data, for which other spectral decomposition methods, such as Principal Component Analysis, do not yield satisfactory results. The comparison shows overall good agreement between MODIS/MISR and AERONET AOD variability. The correlations between the first three modes of MCA results for both MODIS/AERONET and MISR/ AERONET are above 0.8 for the full data set and above 0.75 for the AOD anomaly data. The correlations between MODIS and MISR modes are also quite high (greater than 0.9). We also examine the extent of spatial agreement between satellite and AERONET AOD data at the selected stations. Some sites with disagreements in the MCA results, such as Kanpur, also have low spatial coherency. This should be associated partly with high AOD spatial variability and partly with uncertainties in satellite retrievals due to the seasonally varying aerosol types and surface properties.

aerosols↗

Introduction to Vector Field Visualization

Vector field visualization techniques are essential to help us understand the complex dynamics of flow fields. These can be found in a wide range of applications such as study of flows around an aircraft, the blood flow in our heart chambers, ocean circulation models, and severe weather predictions. The vector fields from these various applications can be visually depicted using a number of techniques such as particle traces and advecting textures. In this tutorial, we present several fundamental algorithms in flow visualization including particle integration, particle tracking in time-dependent flows, and seeding strategies. For flows near surfaces, a wide variety of synthetic texture-based algorithms have been developed to depict near-body flow features. The most common approach is based on the Line Integral Convolution (LIC) algorithm. There also exist extensions of LIC to support more flexible texture generations for 3D flow data. This tutorial reviews these algorithms. Tensor fields are found in several real-world applications and also require the aid of visualization to help users understand their data sets. Examples where one can find tensor fields include mechanics to see how material respond to external forces, civil engineering and geomechanics of roads and bridges, and the study of neural pathway via diffusion tensor imaging. This tutorial will provide an overview of the different tensor field visualization techniques, discuss basic tensor decompositions, and go into detail on glyph based methods, deformation based methods, and streamline based methods. Practical examples will be used when presenting the methods; and applications from some case studies will be used as part of the motivation.

Kao, David↗

Some practicable applications of quadtree data structures/representation in astronomy

Development of quadtree as hierarchical data structuring technique for representing spatial data (like points, regions, surfaces, lines, curves, volumes, etc.) has been motivated to a large extent by storage requirements of images, maps, and other multidimensional (spatially structured) data. For many spatial algorithms, time-efficiency of quadtrees in terms of execution may be as important as their space-efficiency concerning storage conditions. Briefly, the quadtree is a class of hierarchical data structures which is based on the recursive partition of a square region into quadrants and sub-quadrants until a predefined limit. Beyond the wide applicability of quadtrees in image processing, spatial information analysis, and building digital databases (processes becoming ordinary for the astronomical community), there may be numerous further applications in astronomy. Some of these practicable applications based on quadtree representation of astronomical data are presented and suggested for further considerations. Examples are shown for use of point as well as region quadtrees. Statistics of different leaf and non-leaf nodes (homogeneous and heterogeneous sub-quadrants respectively) at different levels may provide useful information on spatial structure of astronomical data in question. By altering the principle guiding the decomposition process, different types of spatial data may be focused on. Finally, a sampling method based on quadtree representation of an image is proposed which may prove to be efficient in the elaboration of sampling strategy in a region where observations were carried out previously either with different resolution or/and in different bands.

Pasztor, L.↗

Subband And Transform Compression Of Video Signals

Class of hierarchical subband coders developed primarily for compression of image data at video rates. Offers good performance with limited computational complexity and with flexibility inherent in subband decomposition. Particular subband decomposition chosen for coders appears to hide large quantitative errors effectively, largely because decomposition occurs along two-dimensional spatial-frequency-domain boundaries resembling spatial-frequency-domain curves of constant sensitivity of human visual system. Curves found approximately diamond-shaped: thus, low-pass filtering for reduction of data ideally involves nonrectangular passbands.

Sauer, Ken↗

Laboratory Study of Water Surface Roughness Generation by Wave-Current Interaction

Within the framework of this project, the blocking of waves by inhomogeneous currents was studied. A laboratory experiment was conducted in collaboration with Steven R. Long at the linear wave tank of the NASA Air-Sea Interaction Facility, Wallops Island, VA during May 1999. Mechanically-generated waves were blocked approximately 3m upstream from the wave paddle by an opposing current. A false bottom was used to obtain a spatially varying flow field in the measurement section of the wave tank. We used an imaging slope gauge, which was mounted directly underneath the sloping section of the false tank bottom to observe the wave field. For a given current speed, the amplitude and the frequency of the waves was adjusted so that the blocking occurred within the observed footprint. Image sequences of up to 600 images at up 100 Hz sampling rate were recorded for an area of approximately 25cm x 25cm. Unlike previous measurements with wave wire gauges, the captured image sequences show the generation of the capillary waves at the blocking point and give detailed insight into the spatial and temporal evolution of the blocking process. The image data were used to study the wave-current interaction for currents from 5 to 25 cm/s and waves with frequencies between 1 and 3 Hz. First the images were calibrated with regard to size and slope. Then standard Fourier techniques as well the empirical mode decomposition method developed by Dr. Norden Huang and Dr. Steven R. Long were employed to quantify the wave number downshift from the gravity to the capillary regime.

Klinke, Jochen↗