Search NASA⌕ Search

SEARCH · Search NASA

Results for “Deconvolution”

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 37 records · Page 2

Image Deconvolution and Point-spread Function Reconstruction with STARRED: A Wavelet-based Two-channel Method Optimized for Light-curve Extraction

We present starred, a point-spread function (PSF) reconstruction, two-channel deconvolution, and light-curve extraction method designed for high-precision photometric measurements in imaging time series. An improved resolution of the data is targeted rather than an infinite one, thereby minimizing deconvolution artifacts. In addition, starred performs a joint deconvolution of all available data, accounting for epoch-to-epoch variations of the PSF and decomposing the resulting deconvolved image into a point source and an extended source channel. The output is a high-signal-to-noise-ratio, high-resolution frame combining all data and the photometry of all point sources in the field of view as a function of time. Of note, starred also provides exquisite PSF models for each data frame. We showcase three applications of starred in the context of the imminent LSST survey and of JWST imaging: (i) the extraction of supernovae light curves and the scene representation of their host galaxy; (ii) the extraction of lensed quasar light curves for time-delay cosmography; and (iii) the measurement of the spectral energy distribution of globular clusters in the "Sparkler," a galaxy at redshift z = 1.378 strongly lensed by the galaxy cluster SMACS J0723.3-7327. starred is implemented in jax, leveraging automatic differentiation and graphics processing unit acceleration. This enables the rapid processing of large time-domain data sets, positioning the method as a powerful tool for extracting light curves from the multitude of lensed or unlensed variable and transient objects in the Rubin-LSST data, even when blended with intervening objects.

79 ASTRONOMY AND ASTROPHYSICS↗

Digital pulse deconvolution with adaptive shaping for real-time high-resolution high-throughput gamma spectroscopy

A high-resolution high-throughput real-time adaptive digital pulse processing system is developed for high count rate gamma-ray spectroscopy applications. The adaptive digital pulse processing algorithms are implemented on a reconfigurable FPGA and include a pulse deconvolver, adaptive shaping filter, timing filter, baseline restorer, and pile-up rejecter. Digital pulse deconvolution is implemented to reconstruct the original detector signal from the preamplifier signal, which reduces the resolution deterioration due to pulse pile-up. The deconvoluted signal is shaped with a trapezoid filter and the shaping parameter is selected adaptively based on the time separation between successive input pulses. Experimental measurements are performed with a 137Cs source under varying count-rate conditions and using germanium detectors equipped with resistive feedback and transistor reset preamplifiers (TRP). The results herein demonstrate that when using a TRP, adaptive digital signal processing allows handling 106 counts/s. In addition, the implementation of a deconvolution approach limits resolution deterioration for throughput rates that are 4 to 10 times better than achievable in typical digital and analog gamma-ray spectroscopy systems.

47 OTHER INSTRUMENTATION↗

Deconvolution of dynamic heterogeneity in protein structure

Heterogeneity is intrinsic to the dynamic process of a chemical reaction. As reactants are converted to products via intermediates, the nature and extent of heterogeneity vary temporally throughout the duration of the reaction and spatially across the molecular ensemble. The goal of many biophysical techniques, including crystallography and spectroscopy, is to establish a reaction trajectory that follows an experimentally provoked dynamic process. It is essential to properly analyze and resolve heterogeneity inevitably embedded in experimental datasets. We have developed a deconvolution technique based on singular value decomposition (SVD), which we have rigorously practiced in diverse research projects. In this review, we recapitulate the motivation and challenges in addressing the heterogeneity problem and lay out the mathematical foundation of our methodology that enables isolation of chemically sensible structural signals. We also present a few case studies to demonstrate the concept and outcome of the SVD-based deconvolution. Finally, we highlight a few recent studies with mechanistic insights made possible by heterogeneity deconvolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Qubit noise deconvolution

We present a noise deconvolution technique to remove a wide class of noises when performing arbitrary measurements on qubit systems. In particular, we derive the inverse map of the most common single qubit noisy channels, and exploit it at the data processing step to obtain noise-free estimates of observables evaluated on a qubit system subject to known noise. We illustrate a self-consistency check to ensure that the noise characterization is accurate providing simulation results for the deconvolution of a generic Pauli channel, as well as experimental evidence of the deconvolution of decoherence noise occurring on Rigetti quantum hardware.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Improved deconvolution of very weak confocal signals

Deconvolution is typically used to sharpen fluorescence images, but when the signal-to-noise ratio is low, the primary benefit is reduced noise and a smoother appearance of the fluorescent structures. 3D time-lapse (4D) confocal image sets can be improved by deconvolution. However, when the confocal signals are very weak, the popular Huygens deconvolution software erases fluorescent structures that are clearly visible in the raw data. We find that this problem can be avoided by prefiltering the optical sections with a Gaussian blur. Analysis of real and simulated data indicates that the Gaussian blur prefilter preserves meaningful signals while enabling removal of background noise. This approach is very simple, and it allows Huygens to be used with 4D imaging conditions that minimize photodamage .

Day, Kasey J.↗

Improved deconvolution of very weak confocal signals

Deconvolution is typically used to sharpen fluorescence images, but when the signal-to-noise ratio is low, the primary benefit is reduced noise and a smoother appearance of the fluorescent structures. 3D time-lapse (4D) confocal image sets can be improved by deconvolution. However, when the confocal signals are very weak, the popular Huygens deconvolution software erases fluorescent structures that are clearly visible in the raw data. We find that this problem can be avoided by prefiltering the optical sections with a Gaussian blur. Analysis of real and simulated data indicates that the Gaussian blur prefilter preserves meaningful signals while enabling removal of background noise. This approach is very simple, and it allows Huygens to be used with 4D imaging conditions that minimize photodamage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The application of a time-domain deconvolution technique for identification of experimental acoustic-emission signals

A method is presented for the signature analysis of pulses by reconstructing in the time domain the shape of the pulse prior to its passing through the measurement system. This deconvolution technique is first evaluated using an idealized system and analytical pulse models and is shown to provide improved results. An experimental situation is then treated; system-component models are developed for the digitizer, tape recorder, filter, transducer and mechanical structure. To accommodate both calibration results and manufacturer's data, and to provide stable mathematical models entails considerable effort: some 30 parameters must be identified to model this system - which is still a substantial approximation - albeit of very high order. Experimental pulses generated by a ball drop, spark discharge and a tearing crack are then deconvoluted 'back through' the system as modeled, using this technique. These results are compared and indicate (a) that consistent shapes may be expected from a given type of source and (b) that some sources can be identified with greater clarity using the deconvolution approach.

Houghton, J. R.↗

The direct deconvolution of X-ray spectra

The method of deconvolution of proportional counter X-ray spectra as outlined by Blissett and Cruise is reviewed with particular emphasis on low-energy spectra. This method involves the expansion of the incident spectrum in terms of a set of orthonormal singular functions which propagate independently through the detector matrix. When applied to low-energy detectors which typically exhibit absorption features in their efficiency curves, the initial Blissett and Cruise prescription is shown to be unstable. Alternative methods for handling the efficiency are described and evaluated. In the deconvolution of steep spectra, additional distortions are shown to arise as a result of sidelobe oscillations in the effective response functions. A selective weighting procedure is thus introduced to suppress spurious features of this type. Finally, the role of filtering in the deconvolution procedure is discussed and several suitable forms for the filter are suggested.

Kahn, S. M.↗

Data enhancement and analysis through mathematical deconvolution of signals from scientific measuring instruments

Mathematical deconvolution of digitized analog signals from scientific measuring instruments is shown to be a means of extracting important information which is otherwise hidden due to time-constant and other broadening or distortion effects caused by the experiment. Three different approaches to deconvolution and their subsequent application to recorded data from three analytical instruments are considered. To demonstrate the efficacy of deconvolution, the use of these approaches to solve the convolution integral for the gas chromatograph, magnetic mass spectrometer, and the time-of-flight mass spectrometer are described. Other possible applications of these types of numerical treatment of data to yield superior results from analog signals of the physical parameters normally measured in aerospace simulation facilities are suggested and briefly discussed.

Wood, G. M.↗

Optimal application of Morrison's iterative noise removal for deconvolution

Morrison's iterative method of noise removal can be applied for both noise removal alone and noise removal prior to deconvolution. This method is applied to noise of various noise levels added to determine the optimum use of the method. The phase shift method of migration and modeling is evaluated and the results are compared to Stolt's approach. A method is introduced by which the optimum iterative number for deconvolution can be found. Statistical computer simulation is used to describe the optimum use of two convergent iterative techniques for seismic data. The Always-Convergent deconvolution technique was applied to data recorded during the quantitative analysis of materials through NonDestructive Evaluation (NDE) in which ultrasonic signals were used to detect flaws in substances such as composites.

Ioup, George E.↗

Nimbus 7 earth radiation budget wide field of view climate data set improvement. I - The earth albedo from deconvolution of shortwave measurements

A deconvolution method for extracting the top of the atmosphere (TOA) mean, daily albedo field from a set of wide-FOV (WFOV) shortwave radiometer measurements is proposed. The method is based on constructing a synthetic measurement for each satellite observation. The albedo field is represented as a truncated series of spherical harmonic functions, and these linear equations are presented. Simulation studies were conducted to determine the sensitivity of the method. It is observed that a maximum of about 289 pieces of data can be extracted from a set of Nimbus 7 WFOV satellite measurements. The albedos derived using the deconvolution method are compared with albedos derived using the WFOV archival method; the developed albedo field achieved a 20 percent reduction in the global rms regional reflected flux density errors. The deconvolution method is applied to estimate the mean, daily average TOA albedo field for January 1983. A strong and extensive albedo maximum (0.42), which corresponds to the El Nino/Southern Oscillation event of 1982-1983, is detected over the south central Pacific Ocean.

Hucek, Richard R.↗

Optimal application of Morrison's iterative noise removal for deconvolution. Appendices

Morrison's iterative method of noise removal, or Morrison's smoothing, is applied in a simulation to noise-added data sets of various noise levels to determine its optimum use. Morrison's smoothing is applied for noise removal alone, and for noise removal prior to deconvolution. For the latter, an accurate method is analyzed to provide confidence in the optimization. The method consists of convolving the data with an inverse filter calculated by taking the inverse discrete Fourier transform of the reciprocal of the transform of the response of the system. Various length filters are calculated for the narrow and wide Gaussian response functions used. Deconvolution of non-noisy data is performed, and the error in each deconvolution calculated. Plots are produced of error versus filter length; and from these plots the most accurate length filters determined. The statistical methodologies employed in the optimizations of Morrison's method are similar. A typical peak-type input is selected and convolved with the two response functions to produce the data sets to be analyzed. Both constant and ordinate-dependent Gaussian distributed noise is added to the data, where the noise levels of the data are characterized by their signal-to-noise ratios. The error measures employed in the optimizations are the L1 and L2 norms. Results of the optimizations for both Gaussians, both noise types, and both norms include figures of optimum iteration number and error improvement versus signal-to-noise ratio, and tables of results. The statistical variation of all quantities considered is also given.

Ioup, George E.↗

Extracting ocean surface information from altimeter returns - The deconvolution method

An evaluation of the deconvolution method for estimating ocean surface parameters from ocean altimeter waveforms is presented. It is shown that this method presents a fast, accurate way of determining the ocean surface parameters from noisy altimeter data. Three parameters may be estimated by using this method, including the altimeter-height error, the ocean-surface standard deviation, and the ocean-surface skewness. By means of a Monte Carlo experiment, an 'optimum' deconvolution algorithm and the accuracies with which the above parameters may be estimated using this algorithm are determined. Then the influence of instrument effects, such as errors in calibration and pointing-angle estimation, on the estimated parameters is examined. Finally, the deconvolution algorithm is used to estimate height and ocean-surface parameters from Seasat data.

Rodriguez, Ernesto↗

Determination of design and operation parameters for upper atmospheric research instrumentation to yield optimum resolution with deconvolution

The final report for work on the determination of design and operation parameters for upper atmospheric research instrumentation to yield optimum resolution with deconvolution is presented. Papers and theses prepared during the research report period are included. Among all the research results reported, note should be made of the specific investigation of the determination of design and operation parameters for upper atmospheric research instrumentation to yield optimum resolution with deconvolution. A methodology was developed to determine design and operation parameters for error minimization when deconvolution is included in data analysis. An error surface is plotted versus the signal-to-noise ratio (SNR) and all parameters of interest. Instrumental characteristics will determine a curve in this space. The SNR and parameter values which give the projection from the curve to the surface, corresponding to the smallest value for the error, are the optimum values. These values are constrained by the curve and so will not necessarily correspond to an absolute minimum in the error surface.

Ioup, George E.↗

Wavespace-Based Coherent Deconvolution

Array deconvolution is commonly used in aeroacoustic analysis to remove the influence of a microphone array's point spread function from a conventional beamforming map. Unfortunately, the majority of deconvolution algorithms assume that the acoustic sources in a measurement are incoherent, which can be problematic for some aeroacoustic phenomena with coherent, spatially-distributed characteristics. While several algorithms have been proposed to handle coherent sources, some are computationally intractable for many problems while others require restrictive assumptions about the source field. Newer generalized inverse techniques hold promise, but are still under investigation for general use. An alternate coherent deconvolution method is proposed based on a wavespace transformation of the array data. Wavespace analysis offers advantages over curved-wave array processing, such as providing an explicit shift-invariance in the convolution of the array sampling function with the acoustic wave field. However, usage of the wavespace transformation assumes the acoustic wave field is accurately approximated as a superposition of plane wave fields, regardless of true wavefront curvature. The wavespace technique leverages Fourier transforms to quickly evaluate a shift-invariant convolution. The method is derived for and applied to ideal incoherent and coherent plane wave fields to demonstrate its ability to determine magnitude and relative phase of multiple coherent sources. Multi-scale processing is explored as a means of accelerating solution convergence. A case with a spherical wave front is evaluated. Finally, a trailing edge noise experiment case is considered. Results show the method successfully deconvolves incoherent, partially-coherent, and coherent plane wave fields to a degree necessary for quantitative evaluation. Curved wave front cases warrant further investigation. A potential extension to nearfield beamforming is proposed.

Bahr, Christopher J.↗

Mass Spectra Deconvolution of Gaseous Mixtures Containing Volatile Organic Compounds

The Spacecraft Atmosphere Monitor (S.A.M.) analyzes all gaseous pollutants in the 1-150Th mass-to-charge range based on their positive-ion mass spectra. Of particular interest is the Trace Gas Analysis (TGA) mode of operation in which S.A.M detects minute amounts of volatile organic compounds (VOC). Pollutant component of the ambient air is sampled on demand in twenty full mass spectra per second. All spectra is accumulated in two-second intervals by an on-board Field Programmable Gate Arrays (FPGA) unit. Mass spectra represents the number of detected fragment ions in a given mass-to-charge ratio channel. Electron impact ionization of neutral molecules results in creation of numerous fragment ions and most of organic compounds will contribute several identical fragments. Despite the presence of these molecular isobars, we developed the deconvolution algorithm capable of identifying target species based on their characteristic fragmentation patterns. We investigate the efficiency of deconvolution algorithm as a function of mass resolution with which mass spectrum is acquired. Higher the mass resolution, better the chances are to distinguish between stereoisomers using their fragmentation propensities, but also results in higher data volumes to be processed by a set of small footprint software stacks hosted by an on-board computer. Finding the balance between deconvolution accuracy and generated data volume under time constrains and limited computing resources is the main topic of this study.

Darrach, M.↗

Fast Bayesian blind deconvolution with Huber Super Gaussian priors

We report expectation Maximization (EM) based inference has already proven to be a very powerful tool to solve blind image deconvolution (BID) problems. Unfortunately, three important problems still impede the application of EM in BID: the undesirable saddle points and local minima caused by highly nonconvex priors, the instability around zero of some of the most interesting sparsity promoting priors, and the intrinsic high computational cost of the corresponding BID algorithm. In this paper we first show how Super Gaussian priors can be made numerically tractable around zero by introducing the family of Huber Super Gaussian priors and then present a fast EM based blind deconvolution method formulated in the image space. In the proposed computational approach, image and kernel estimation are performed by using the Alternating Direction Method of Multipliers (ADMM), which allows to exploit the advantages of FFT computation. For highly nonconvex priors, we propose a Smooth ADMM (SADMM) approach to avoid poor BID estimates. In conclusion, extensive experiments demonstrate that the proposed method significantly outperforms state-of-the-art BID methods in terms of quality of the reconstructions and speed.

97 MATHEMATICS AND COMPUTING↗

Multiqubit noise deconvolution and characterization

We present a noise deconvolution technique for obtaining noiseless expectation values of noisy observables at the output of multiqubit quantum channels. For any number of qubits or in the presence of correlations, our protocol applies to any mathematically invertible noise model, even when its inverse map is not physically implementable, i.e., when it is neither completely positive nor trace preserving. For a generic observable affected by Pauli noise it provides a quadratic speedup, always producing a rescaling of its Pauli basis components. We show that it is still possible to achieve the deconvolution while experimentally estimating the noise parameters, whenever these are unknown (bypassing resource-heavy techniques such as quantum process tomography). We provide a simulation, with examples for both Pauli and non-Pauli channels.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗