Search NASA⌕ Search

SEARCH · Search NASA

Results for “image 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.

160 records · Page 9

An assessment of image reconstruction from balloon-borne and the IRAS data

Angular resolution and structural information from the far-infrared mapping of astronomical sources (Galactic star forming regions, spiral galaxies, etc.) made using the TIFR 1 m balloon-borne telescope and the IRAS have been compared. The effective wavelengths of the TIFR two-band photometer are 58 and 150 microns. From IRAS, the survey COADD data, additional observations (AO's) made with the survey detectors with different Macros (DPS, DSD, DPM), as well as the chopped photometric channel (CPC) data have been considered here. The observed signals have been processed using different deconvolution strategies, either based on a maximum entropy method (MEM) developed at TIFR or the HiRes package developed at IPAC. Relative merits of each of these, under different conditions of signal to noise ratio, are highlighted. The following sources have been selected for illustration: Carina complex, W31 region, IRAS 10361-5830 (all Galactic), M101 and M81 (extragalactic). The main conclusions are: far-infrared maps from MEM deconvolution of balloon-borne data have the best angular resolution; MEM deconvolution of IRAS AO's gives resolution comparable to HiRes but with less amount of computation, though the dynamic range in MEM maps is less than in HiRes maps.

Ghosh, S. K.↗

Strongly intermittent far scrape-off layer fluctuations in Alcator C-Mod plasmas close to the empirical discharge density limit

Abstract Intermittent plasma fluctuations in the boundary region of the Alcator C-Mod device were comprehensively investigated using data time-series from gas puff imaging and mirror Langmuir probe diagnostics. Fluctuations were sampled during stationary plasma conditions in ohmically heated, lower single null diverted configurations with scans in both line-averaged density and plasma current, with Greenwald density fractions up to 0.85. Utilizing a stochastic model, we describe the plasma fluctuations as a super-position of uncorrelated pulses, with large-amplitude events corresponding to blob-like filaments moving through the scrape-off layer (SOL). A deconvolution method is used to estimate the pulse arrival times and amplitudes. The analysis reveals a significant increase of pulse amplitudes and waiting times as the line-averaged density approaches the empirical discharge density limit. Broadened and flattened average radial profiles are thus accompanied by strongly intermittent and large-amplitude fluctuations. Although these filaments are arriving less frequently at high line-averaged densities, we show that there are significant increases in radial far-SOL particle and heat fluxes which will further enhance plasma–wall interactions. The stochastic model has been used as a framework for study of the scalings in the intermittency parameter, flux and mean amplitude and waiting times, and is being used to inform predictive capability for the effects of filamentary transport as a function of Greenwald fraction.

Physics↗

Evidence on Sizes and Fragmentation of the Nuclei of Comet Shoemaker-Levy 9 from Hubble Space Telescope Images

Central regions on the digital maps of 13 nuclear condensations of Comet Shoemaker-Levy 9, obtained with the Planetary Camera of the Hubble Space Telescope on January 24-25, March 28-30, and July 4, 1994, have been analyzed with the aim to identify the presence of distinct, major fragments in each condensation, to deconvolve their contributions to the signal that also includes the contribution from a surrounding cloud of dust (modeled as an extended source, using two different laws), to estimate the dimensions of the fragments and to study their temporal variations, and to determine the spatial distributions of the fragments as projected onto the plane of the sky. The deconvolution method applied is described and the results of the analysis are summarized, including the finding that sizable fragments die survive until the time of atmospheric entry. This result does not contradict evidence of the comet's continuing, apparently spontaneous fragmentation, which still went on long after the extremely close approach to Jupiter in July 1992 and which, because of the jovian tidal effects, may even have intensified in the final days before the crash on Jupiter. On plausible assumptions, the largest fragments are found to have had effect diameters of 4 km as late as March and event early July 1994. In most condensations, several sizable companions ( 1 km or more across) have been detected within 1000 km of the projected location of the brightest fragment and the surrounding dust cloud has been found to be centered on a point that is shifted in the general direction of the tail, probably due to effects of solar radiation pressure. Since the developed approach is based on certain premises and involves approximations, the results should be viewed as preliminary and the problem should be a subject of further investigation.

Comet Shoemaker-Levy 9↗

A hard X-ray imaging instrument for solar and cosmic sources

A hard X-ray imaging instrument is described which is capable of high-resolution imaging of solar and cosmic hard X-ray sources between 2 and 80 keV during Shuttle sortie flights. The properties of solar burst sources and the resulting instrument requirements are discussed. The instrument envelope of 1.2 x 1.2 x 3.0 meters includes a tungsten multigrid collimator which has 4-arcsec resolution, a 40-arcmin response envelope and a point-source effective area of 26 sq cm. A combination of periodic fan beams and nonperiodic pencil beams enables a unique deconvolution to be achieved within a 128 x 128 arcsec field without mechanical scanning. The detector system is a set of direct-readout 40 atm-cm xenon-filled proportional counters designed to minimize background. The instrument is capable of refurbishment to optimize the collimator configuration for specific solar or cosmic scientific objectives, to upgrade the angular resolution, or to extend the high-energy response.

Hurford, G. J.↗

From minerals to rocks: Toward modeling lithologies with remote sensing

High spectral resolution imaging spectroscopy will play an important role in future planetary missions. Sophisticated approaches will be needed to unravel subtle, super-imposed spectral features typically of natural systems, and to maximize the science return of these instruments. Carefully controlled laboratory investigations using homogeneous mineral separates have demonstrated that variations due to solid solution, changes in modal abundances, and the effects of particle size are well understood from a physical basis. In many cases, these variations can be modeled quantitatively using photometric models, mixing approaches, and deconvolution procedures. However, relative to the spectra of individual mineral components, reflectance spectra of rocks and natural surfaces exhibit a reduced spectral contrast. In addition, soils or regolith, which are likely to dominate any natural planetary surface, exhibit spectral properties that have some similarities to the parent materials, but due to weathering and alteration, differences remain that cannot yet be fully recreated in the laboratory or through mixture modeling. A significant challenge is therefore to integrate modeling approaches to derive both lithologic determinations and include the effects of alteration. We are currently conducting laboratory investigations in lithologic modeling to expand upon the basic results of previous analyses with our initial goal to more closely match physical state of natural systems. The effects of alteration are to be considered separately.

Mustard, John F.↗

Lens Model Accuracy in the Expected LSST Lensed AGN Sample

Strong gravitational lensing of active galactic nuclei (AGN) enables measurements of cosmological parameters through time-delay cosmography (TDC). With data from the upcoming LSST survey, we anticipate using a sample of O(1000) lensed AGN for TDC. To prepare for this dataset and enable this measurement, we construct and analyze a realistic mock sample of 1300 systems drawn from the OM10 (Oguri & Marshall 2010) catalog of simulated lenses with AGN sources at $z<3.1$ in order to test a key aspect of the analysis pipeline, that of the lens modeling. We realize the lenses as power law elliptical mass distributions and simulate 5-year LSST i-band coadd images. From every image, we infer the lens mass model parameters using neural posterior estimation (NPE). Focusing on the key model parameters, $θ_E$ (the Einstein Radius) and $γ_{lens}$ (the projected mass density profile slope), with consistent mass-light ellipticity correlations in test and training data, we recover $θ_E$ with less than 1% bias per lens, 6.5% precision per lens and $γ_{lens}$ with less than 3% bias per lens, 8% precision per lens. We find that lens light subtraction prior to modeling is only useful when applied to data sampled from the training prior. If emulated deconvolution is applied to the data prior to modeling, precision improves across all parameters by a factor of 2. Finally, we combine the inferred lens mass models using Bayesian Hierarchical Inference to recover the global properties of the lens sample with less than 1% bias.

Venkatraman, Padmavathi [Illinois U., Urbana; KIPA↗

RIFTA: A Robust Iterative Fourier Transform-based dwell time Algorithm for ultra-precision ion beam figuring of synchrotron mirrors

With the rapid evolution of synchrotron X-ray sources, the demand for high-precision X-ray mirrors has greatly increased. Single nanometer profile error is required to keep imaging capability at the diffraction limit. Ion Beam Figuring (IBF), as a highly deterministic surfacing technique, has been used for ultra-precision finishing of mirrors. One crucial step that guides the IBF process is dwell time calculation. A valid dwell time solution should be non-negative and duplicate the shape of the desired removal map. Another important aspect is to minimize the total dwell time. In this study, we propose a Robust Iterative Fourier Transform-based dwell time Algorithm (RIFTA) that automatically fulfills these requirements. First, the thresholded inverse filtering in Fourier transform-based deconvolution is stabilized and automated by optimizing the threshold value using the Nelder-Mead simplex algorithm. Second, a novel two-level iterative scheme is proposed to guarantee the minimized total dwell time with its non-negativity at each dwell point. Third, a bicubic resampling is employed to flexibly adapt the calculated dwell time map to any IBF process intervals. The performance of RIFTA is first studied with simulation, followed by a comparison with the other state-of-the-art dwell time algorithms. We then demonstrate with an experiment that, using the dwell time calculated by the RIFTA, the total dwell time is shortened by a factor of two and the RMS in a 5 × 50 mm clear aperture was reduced from 3.4 nm to 1.1 nm after one IBF run, which proves the effectiveness and the efficiency of the proposed algorithm.

42 ENGINEERING↗

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↗

The GPM Combined Algorithm

In this paper, the operational Global Precipitation Measurement (GPM) mission combined radar-radiometer algorithm is thoroughly described. The operational combined algorithm is designed to reduce uncertainties in GPM Core Observatory precipitation estimates by effectively integrating complementary information from the GPM Dual-Frequency Precipitation Radar (DPR) and the GPM Microwave Imager (GMI) into an optimal, physically consistent precipitation product. Although similar in many respects to previously developed combined algorithms, the GPM combined algorithm has several unique features that are specifically designed to meet the GPM objectives of deriving, based on GPM Core Observatory information, accurate and physically consistent precipitation estimates from multiple spaceborne instruments, and ancillary environmental data from reanalyses. The algorithm features an optimal estimation framework based on a statistical formulation of the Gauss-Newton method, a parameterization for the nonuniform distribution of precipitation within the radar fields of view, a methodology to detect and account for multiple scattering in Ka-band DPR observations, and a statistical deconvolution technique that allows for an efficient sequential incorporation of radiometer information into DPR precipitation retrievals.

Grecu, Mircea↗

Machine Learning Benchmarks for the Classification of Equivalent Circuit Models from Electrochemical Impedance Spectra

Analysis of Electrochemical Impedance Spectroscopy (EIS) data for electrochemical systems often consists of defining an Equivalent Circuit Model (ECM) using expert knowledge and then optimizing the model parameters to deconvolute various resistance, capacitive, inductive, or diffusion responses. For small data sets, this procedure can be conducted manually; however, it is not feasible to manually define a proper ECM for extensive data sets with a wide range of EIS responses. Automatic identification of an ECM would substantially accelerate the analysis of large sets of EIS data. We showcase machine learning methods to classify the ECMs of 9,300 impedance spectra provided by QuantumScape for the BatteryDEV hackathon. The best-performing approach is a gradient-boosted tree model utilizing a library to automatically generate features, followed by a random forest model using the raw spectral data. A convolutional neural network using boolean images of Nyquist representations is presented as an alternative, although it achieves a lower accuracy. We publish the data and open source the associated code. The approaches described in this article can serve as benchmarks for further studies. A key remaining challenge is the identifiability of the labels, underlined by the model performances and the comparison of misclassified spectra.

25 ENERGY STORAGE↗

Comparison between mirror Langmuir probe and gas-puff imaging measurements of intermittent fluctuations in the Alcator C-Mod scrape-off layer

Statistical properties of the scrape-off layer plasma fluctuations are studied in ohmically heated plasmas in the Alcator C-Mod tokamak. For the first time, plasma fluctuations as well as parameters that describe the fluctuations are compared across measurements from a mirror Langmuir probe (MLP) and from gas-puff imaging (GPI) that sample the same plasma discharge. This comparison is complemented by an analysis of line emission time-series data, synthesized from the MLP electron density and temperature measurements. The fluctuations observed by the MLP and GPI typically display relative fluctuation amplitudes of order unity together with positively skewed and flattened probability density functions. Such data time series are well described by an established stochastic framework that models the data as a superposition of uncorrelated, two-sided exponential pulses. The most important parameter of the process is the intermittency parameter, $γ=τ_d/τ_w$, where τd denotes the duration time of a single pulse and $τ_w$ gives the average waiting time between consecutive pulses. In this study we show, using a new deconvolution method, that these parameters can be consistently estimated from different statistics of the data. We also show that the statistical properties of the data sampled by the MLP and GPI diagnostic are very similar. Finally, a synthetic GPI signal using only plasma parameters sampled by the MLP shows qualitatively different fluctuation statistics from the measured GPI signal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fourier-based three-dimensional multistage transformer for aberration correction in multicellular specimens

High-resolution tissue imaging is often compromised by sample-induced optical aberrations that degrade resolution and contrast. Although wavefront sensor-based adaptive optics (AO) can measure these aberrations, such hardware solutions are typically complex, expensive to implement and slow when serially mapping spatially varying aberrations across large fields of view. Here we introduce AOViFT (adaptive optical vision Fourier transformer)—a machine learning-based aberration sensing framework built around a three-dimensional multistage vision transformer that operates on Fourier domain embeddings. AOViFT infers aberrations and restores diffraction-limited performance in puncta-labeled specimens with substantially reduced computational cost, training time and memory footprint compared to conventional architectures or real-space networks. We validated AOViFT on live gene-edited zebrafish embryos, demonstrating its ability to correct spatially varying aberrations using either a deformable mirror or postacquisition deconvolution. By eliminating the need for the guide star and wavefront sensing hardware and simplifying the experimental workflow, AOViFT lowers technical barriers for high-resolution volumetric microscopy across diverse biological samples.

Alshaabi, Thayer [Howard Hughes Medical Institute,↗

Radiation Mapping for an Unmanned Aerial Vehicle: Development and Simulated Testing of Algorithms for Source Mapping and Navigation Path Generation

Image reconstruction algorithms were developed for radiation source mapping and used for generating the search path of a moving radiation detector, such as one onboard an unmanned aerial vehicle. Simulations consisted of first assuming radioactive sources of varying complexity and estimating the radiation fields that would then be produced by that source distribution. Next, the "measurements" that would result from a pair of adjacent spatial locations were computed. A crude estimate of the source distribution likely to have produced such "measurements" was reconstructed based upon the limited measurements. Location of the next "measurement" was then determined as halfway between the location of the estimated source and the current "measurement." With each additional sample, improved source distribution reconstructions were made and used to inform the immediate direction of detector motion. Source reconstruction or mapping was formulated as an inverse problem solved with either maximum a posteriori or least squares (LS) regression deconvolution methods. Different amounts of noise were added to the simulated "measurements," allowing evaluation of the methods' performances as functions of signal-to-noise ratio of the measured map. As expected, methods that promote sparsity were better suited in reconstructing point sources. Reliable prior information of the source distribution also improved the reconstruction results, especially with distributed sources. With a non-negative least square algorithm and the suggested paths it generated, location of sources was successfully estimated to an accuracy of 0.014 m within nine iterations in a single-source scenario and 12 iterations in a two-source scenario, given a 10% error on the integrated counts and a Poisson distribution of the noise associated with the measured counts.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Three-Dimensional Morphological Analysis of ALH84001 Magnetite Using Electron Tomography

We report here the crystal morphologies of MV-1 and ALH84001 magnetites as calculated by back-projection using electron tomography. In the present study, we used a 300 keV TEM with a field emission gun (Tecnai F-30 from FEI Inc.), equipped with a 2048 x 2048 pixel CCD camera from Gatan Inc. to image magnetite crystals over tilt ranges of approx. +/- 72 deg in 2 deg tilt intervals. The images were aligned for back-projection, either manually, or through the use of fiducial 5 nm Au spheres affixed to the specimen prior to microscopy. Three-dimensional (3-D) reconstructions were computed using weighted back-projection of the tilted views. The tomograms were viewed and analyzed as a series of slices 1.0 nm thick, taken parallel to the specimen-supporting grid, using the IMOD software package. The shape of each magnetite crystal was determined by defining the external contour of a given magnetite in each slice and assembling a stack of these contours in 3-D. To aid in visualization, the stacked contour array was reduced to an optimal mesh by Delaunay triangulation. The surface normal to each of the triangles in the mesh was calculated and the triangle faces colored according to the orientation of that surface normal relative to the principal crystallographic axis of magnetite. Green surfaces correspond to {111} orientations, blue surfaces to {100} orientations, and red surfaces to {110} orientations. Triangles whose surface normal did not correspond to one of the principal axes were colored gray. Within the experimental and numerical uncertainties of the deconvolution, the tomographic reconstruction of both MV-1 and ALH84001 magnetites are equivalent and correspond to a truncated hexa-octahedral morphology.

Thomas-Keprta, Kathie L.↗

Fe and O EELS Studies of Ion Irradiated Murchison CM2 Carbonaceous Chondrite Matrix

Introduction: The physical and chemical response of hydrated carbonaceous chondrite materials to space weathering processes is poorly understood. Improving this understanding is a key part of establishing how regoliths on primitive carbonaceous asteroids respond to space weathering processes, knowledge that supports future sample return missions (Hayabusa 2 and OSIRISREx) that are targeting objects of this type. We previously reported on He+ irradiation of Murchison matrix and showed that the irradiation resulted in amorphization of the matrix phyllosilicates, loss of OH, and surface vesiculation. Here, we report electron energy-loss spectroscopy (EELS) measurements of the irradiated material with emphasis on the Fe and O speciation. Sample and Methods: A polished thin section of the Murchison CM2 carbonaceous chondrite was irradiated with 4 kilovolts He(+) (normal incidence) to a total dose of 1 x 10(exp 18) He(+) per square centimeter. We extracted thin sections from both irradiated and unirradiated regions in matrix using focused ion beam (FIB) techniques with electron beam deposition for the protective carbon strap to minimize surface damage artifacts from the FIB milling. The FIB sections were analyzed using a JEOL 2500SE scanning and transmission electron microscope (STEM) equipped with a Gatan Tridiem imaging filter. EELS spectra were collected from 50 nanometer diameter regions with an energy resolution of 0.7 electronvolts FWHM at the zero loss. EELS spectra were collected at low electron doses to minimize possible artifacts from electron-beam irradiation damage. Results and Discussion: Fe L (sub 2,3) EELS spectra from matrix phyllosilicates in CM chondrites show mixed Fe(2+)/Fe(3+) oxidation states with Fe(3+)/Sigma Fe approximately 0.5. Fe L(sub 2,3) spectra from the irradiated/ amorphized matrix phyllosilicates show higher Fe(2+)/Fe(3+) ratios compared to spectra obtained from pristine material at depths beyond the implantation/amorphization layer. We also obtained O Ka spectra from phyllosilicates in both regions of the sample. The O Ka spectra show a pre-edge feature at approximately 530.5 electronvolts that is related to O 2p states hybridized with Fe 3d states. The intensity ratio of the O Ka pre-edge peak relative to the main part of the O Ka edge (that results from transitions of O 1s to 2p states) is lower in the irradiated layer compared to the pristine material and may reflect the loss of O (as OH) as was observed by IR spectroscopy. Conclusions: In addition to amorphization and OH loss, EELS spectra of He(+) irradiated matrix phyllosilicates in Murchison show that some of the Fe(3+) is reduced to Fe(2+). Spectral deconvolution is underway to extract quantitative ratios from the EELS spectra.

Keller, L. P.↗

The ALEXIS data processing package: An IDL based system

The Array of Low Energy X-ray Imaging Sensors (ALEXIS) experiment consists of a mini-satellite containing six wide angle EUV/ultrasoft x-ray telescopes. Its purpose is to map out the sky in three narrow (approximately 5 percent) bandpasses around 66, 71, and 93 eV. The 66 and 71 eV bandpasses are centered on intense Fe emission lines which are characteristic of million degree plasmas such as the one thought to produce the soft x-ray background. The 93 eV bandpass is not near any strong emission lines and is more sensitive to continuum sources. The mission will be launched on the Pegasus Air Launched Vehicle in the second half of 1992 into a 400-nautical-mile, high inclination orbit and will be controlled entirely from a small ground station located at Los Alamos. The project is a collaborative effort between Los Alamos National Laboratory, Sandia National Laboratory, and the University of California-Berkeley Space Sciences Laboratory. The six telescopes are arranged in three pairs. As the satellite spins twice a minute they scan the entire anti-solar hemisphere. Each f/1 telescope consists of a spherical, multilayer-coated mirror with a curved, microchannel plate detector located at the prime focus. The multilayer coatings determine the bandpasses of the telescopes. The field of view of each telescope is 30 degrees with a spatial resolution of 0.5 degree, limited by spherical aberration. The data processing requirements for ALEXIS are large. Each event is one of the six telescopes is telemetered to the ground with its time of arrival and position on the detector. This information must be folded with the aspect solution for the satellite to reconstruct the direction on the sky from which the photon came. Because of the way the six telescopes scan the sky, the effective exposure calculation is also very computationally intensive. ALEXIS may generate up to 100 megabytes of raw data per day, which are converted into a gigabyte per day of processed data. While the processing job for ALEXIS is sizable, the programming staff is small. To maximize programming efficiency, and to make the best use of tools available in the public domain, we chose IDL as our software development platform. IDL was used from the start of instrument development through flight. We use IDL as a top-level executive for the processing tasks (replacing Unix shell scripts), as a device independent graphics engine, as a database manager, and as a final data manipulator. IDL routines spawn special purpose C programs to perform detailed telemetry deconvolution and other specialized functions. We discuss the use of IDL and C within the processing and archiving strategy for the ALEXIS data anlaysis system as implemented on a SPARCstation platform. We also show results from our End-to-End software simulation capability as processed by our analysis codes.

Bloch, J. J.↗