Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reconstruction algorithms”

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 703 records · Page 39

AdjointBackMapV2: Precise reconstruction of arbitrary CNN unit’s activation via adjoint operators

Adjoint operators have been found to be effective in the exploration of CNN’s inner workings (Wan and Choe, 2022). However, the previous no-bias assumption restricted its generalization. We overcome the restriction via embedding input images into an extended normed space that includes bias in all CNN layers as part of the extended space and propose an adjoint-operator-based algorithm that maps high-level weights back to the extended input space for reconstructing an effective hypersurface. Such hypersurface can be computed for an arbitrary unit in the CNN, and here we prove that this reconstructed hypersurface, when multiplied by the original input (through an inner product), will precisely replicate the output value of each unit. We show experimental results based on the CIFAR-10 and CIFAR-100 data sets where the proposed approach achieves near 0 activation value reconstruction error.

97 MATHEMATICS AND COMPUTING↗

Inverse Method for Noninvasive Material Discrimination with Multienergy X-Ray Radiography

Active x-ray inspections have great potential for noninvasive detection of illicit materials. Dual-energy x-ray radiography, in particular, has been used to determine material composition by utilizing the energy-dependence of the x-ray attenuation coefficients. However, current implementations of this method are limited in their ability to determine the material composition of composite objects. Here, in this paper, we present an inverse algorithm that uses multienergy x-ray radiography data to noninvasively reconstruct the contents of a container. A critical feature of the current contribution is that material identification can be performed using conventional detectors and x-ray spectra with different endpoint energies produced by varying nothing more than the tube voltage. Adaptive regularization is used to increase the accuracy of material estimations from multienergy data sets. The utility of these methods is demonstrated with experimentally acquired radiographs obtained using a tunable x-ray source that produces spectra with endpoint energies of 100–450 keV. The object inspected is a scale model of a nuclear materials storage container composed of three-dimensional printed plastic and stainless-steel spheres inside a thin-walled steel container. Reconstructions of the steel sphere thicknesses are within a root-mean-square error of 0.37 cm.

47 OTHER INSTRUMENTATION↗

Measurement of lepton mass squared moments in B → X c ℓ ν ¯ ℓ decays with the Belle II experiment

We present measurements of the first to fourth moments of the lepton mass squared $q$ 2 of $B$ → $X_cℓ\overline{ν}_ℓ$ decays for $ℓ = e, μ$ and with $X_c$ a hadronic system containing a charm quark. These results use a sample of electron-positron collisions at the $Υ$($4S$) resonance corresponding to 62.8 fb -1 of integrated luminosity and collected by the Belle II 2 experiment in 2019 and 2020. To identify the $X_c$ system and reconstruct $q$ 2 , one of the $B$ mesons from an $Υ$($4S$) → $B\overline{B}$ decay is fully reconstructed in a hadronic decay mode using a multivariate $B$ tagging algorithm. We report raw and central moments for $q$ 2 > 1.5 $GeV$ 2 =$c$ 4 up to $q$ 2 > 8.5 $GeV$ 2 =$c$ 4 , probing up to 77% of the accessible $B$ → $X_cℓ\overline{ν}_ℓ$ phase space. This is the first measurement of moments in the experimentally challenging range of [1.5; 2.5] $GeV$ 2 =$c$ 4 . The results can be used for a new determination of |$V_{cb}$| using inclusive $B$ → $X_cℓ\overline{ν}_ℓ$ decays.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

2.5D Super-Resolution Approaches for X-Ray Computed Tomography-Based Inspection of Additively Manufactured Parts

X-ray computed tomography (XCT) is a key tool in non-destructive evaluation of additively manufactured (AM) parts, allowing for internal inspection and defect detection. Despite its widespread use, obtaining high-resolution CT scans can be extremely time consuming. This issue can be mitigated by performing scans at lower resolutions; however, reducing the resolution compromises spatial detail, limiting the accuracy of defect detection. Super-resolution algorithms offer a promising solution for overcoming resolution limitations in XCT reconstructions of AM parts, enabling more accurate detection of defects. While 2D super-resolution methods have demonstrated state-of-the-art performance on natural images, they tend to under-perform when directly applied to XCT slices. On the other hand, 3D super-resolution methods are computationally expensive, making them infeasible for large-scale applications. To address these challenges, we propose a 2.5D super-resolution approach tailored for XCT of AM parts. Our method enhances the resolution of individual slices by leveraging multi-slice information from neighboring 2D slices without the significant computational overhead of full 3D methods. Specifically, we use neighboring low-resolution slices to super-resolve the center slice, exploiting inter-slice spatial context while maintaining computational efficiency. This approach bridges the gap between 2D and 3D methods, offering a practical solution for high-throughput defect detection in AM parts.

Sullivan, Haley↗

Provable Convergence of Plug-and-Play Priors With MMSE Denoisers

Plug-and-play priors (PnP) is a methodology for regularized image reconstruction that specifies the prior through an image denoiser. While PnP algorithms are well understood for denoisers performing maximum a posteriori probability (MAP) estimation, they have not been analyzed for the minimum mean squared error (MMSE) denoisers. Here we address this gap by establishing the first theoretical convergence result for the iterative shrinkage/thresholding algorithm (ISTA) variant of PnP for MMSE denoisers. We show that the iterates produced by PnP-ISTA with an MMSE denoiser converge to a stationary point of some global cost function. We validate our analysis on sparse signal recovery in compressive sensing by comparing two types of denoisers, namely the exact MMSE denoiser and the approximate MMSE denoiser obtained by training a deep neural net.

97 MATHEMATICS AND COMPUTING↗

The adaptation of a strap down formulation for processing inertial platform data

An estimator propagation formulation which ultilizes dynamic data (attitude and sensed acceleration information) from a gimballed inertial platform has been developed to aid in the Shuttle post-flight trajectory reconstruction process and aerodynamic coefficient determination studies. Unlike the classical inertial algorithms, this formulation yields a six degree-of-freedom fully coupled state and attitude estimate. Furthermore, this inertial version is shown to be independent of initial unknown platform misalignments. Results obtained using actual Inertial Measurement Unit (IMU) data and Aerodynamic Coefficient Identification Package (ACIP) strap down data from Shuttle flights are presented.

Heck, M. L.↗

Image data compression with vector quantization in the transform domain

In this paper, an algorithm is presented for image data compression based upon vector quantization of the two-dimensional discrete cosine transformed coefficients. The ac energies of the transformed blocks are used to classify them into eight different ac classes. The ac coefficients of the transformed blocks of class one are set to zero, while those of classes two through eight are transmitted by seven different code books. The dc coefficients of all eight classes are scalar quantized by an adaptive uniform quantizer. As a result, only 4.5 bits instead of eight bits are required to transmit the dc coefficient with negligible additional degradation. Overall, this algorithm requires approximately 0.75 bits per pixel and gives an average reconstruction error of 7.1.

Abdelwahab, A. A.↗

Evaluation of phase-diversity techniques for solar-image restoration

Phase-diversity techniques provide a novel observational method for overcomming the effects of turbulence and instrument-induced aberrations in ground-based astronomy. Two implementations of phase-diversity techniques that differ with regard to noise model, estimator, optimization algorithm, method of regularization, and treatment of edge effects are described. Reconstructions of solar granulation derived by applying these two implementations to common data sets are shown to yield nearly identical images. For both implementations, reconstructions from phase-diverse speckle data (involving multiple realizations of turbulence) are shown to be superior to those derived from conventional phase-diversity data (involving a single realization). Phase-diverse speckle reconstructions are shown to achieve near diffraction-limited resolution and are validated by internal and external consistency tests, including a comparison with a reconstruction using a well-accepted speckle-imaging method.

Paxman, Richard G.↗

Generating Exploration Mission-3 Trajectories to a 9:2 NRHO using Machine Learning

The purpose of this thesis is to design a machine learning algorithm platform that provides expanded knowledge of mission availability through a launch season by improving trajectory resolution and introducing launch mission forecasting. The specific scenario addressed in this paper is one in which data is provided for four deterministic translational maneuvers through a mission to a Near Rectilinear Halo Orbit (NRHO) with a 9:2 synodic frequency. Current launch availability knowledge under NASA's Orion Orbit Performance Team is established by altering optimization variables associated to given reference launch epochs. This current method can bean abstract task and relies on an orbit analyst to structure a mission based on an established mission design methodology associated to the performance of Orion and NASA's Space Launch System. Introducing a machine learning algorithm trained to construct mission scenarios within the feasible range of known trajectories reduces the required interaction of the orbit analyst by removing the needed step of optimizing the orbit to fit an expected translational response required of the spacecraft. In this study, k-Nearest Neighbor and Bayesian Linear Regression successfully predicted classical orbital elements for the launch windows observed. However both algorithms had limitations due to their approaches to model fitting. Training machine learning algorithms off of classical orbital elements introduced a repetitive approach to reconstructing mission segments for different arrival opportunities through the launch window and can prove to be a viable method of launch window scan generation for future missions.

Guzman, Esteban↗

Depth-First Search Image Chasing

This presentation contains results from work on an active NNSS Site-directed research and development (SDRD) project: NLV-019-20. It describes a new algorithm for "few-angle tomography", which can be used to reconstruct an approximation to an object from limited quantities or radiographs under specific conditions. This presentation is an initial documentation of research results thus-far, and we would like to share these methods with colleagues at Lawrence Livermore National Lab who have an expertise in the area of few-angle tomography. Sharing this presentation will be done via email and or teleconference amongst DoE personnel. In the future, we would also like to share these results with collaborators from Academia as well.

97 MATHEMATICS AND COMPUTING↗

Heterogeneous Reconstruction of Tracks and Primary Vertices With the CMS Pixel Tracker

The High-Luminosity upgrade of the Large Hadron Collider (LHC) will see the accelerator reach an instantaneous luminosity of 7 × 10 34 cm −2 s −1 with an average pileup of 200 proton-proton collisions. These conditions will pose an unprecedented challenge to the online and offline reconstruction software developed by the experiments. The computational complexity will exceed by far the expected increase in processing power for conventional CPUs, demanding an alternative approach. Industry and High-Performance Computing (HPC) centers are successfully using heterogeneous computing platforms to achieve higher throughput and better energy efficiency by matching each job to the most appropriate architecture. In this paper we will describe the results of a heterogeneous implementation of pixel tracks and vertices reconstruction chain on Graphics Processing Units (GPUs). The framework has been designed and developed to be integrated in the CMS reconstruction software, CMSSW. The speed up achieved by leveraging GPUs allows for more complex algorithms to be executed, obtaining better physics output and a higher throughput.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Bragg coherent diffraction imaging by simultaneous reconstruction of multiple diffraction peaks

Bragg coherent diffractive imaging (BCDI) is a non-invasive microscopy technique that can visualize the morphology and internal lattice deviations of crystals with nanoscale spatial resolution and picometer deformation sensitivity. While BCDI has been successfully applied in various studies of materials, it is less successful for highly-strained crystals. Specifically, it is difficult to correctly reconstruct the electron density of a highly-strained object using conventional phase retrieval algorithms. Although various algorithms have been developed to overcome this challenge, most of them require a priori knowledge that is not always available in practice. Here we report on a new phase retrieval workflow that can invert diffraction patterns from multiple Bragg peaks simultaneously. The new workflow is explored via simulated diffraction from crystals with various strain conditions. Furthermore, reconstructions from the new workflow consistently demonstrate more accurate electron density maps, in comparison with the conventional method. For highly-strained crystals, the new workflow improves the reliability and consistency of BCDI phase retrieval significantly.

36 MATERIALS SCIENCE↗

Cognition at the Point of Sensing

Over the last 15 years, compressive sensing techniques have been developed which have the potential to greatly reduce the amount of data collected by systems while preserving the amount of information obtained. A cost of this efficiency is that a computationally-intensive optimization routine must be used to put the sensed data into a form that a person can interpret. At the same time, machine learning techniques have experienced tremendous growth as well. Machines have demonstrated the ability learn how to effectively perform tasks such as detection and classification at speeds much faster than humanly possible. Our goal in this project was to study the feasibility of using compressive sensing systems "at the edge." That is, how can compressive sensing sensors be deployed such that information is created at the remote sensor rather than sending raw data to a central processing location? Studies were performed to analyze whether machine learning could be done on the compressively sensed data in its raw form. If a machine is performing the task, is it possible to do so without putting the data into a human interpretable form? We show that this is possible for some systems, in particular a compressive sensing snapshot imaging spectrometer. Machine learning tasks were demonstrated to be more effective and more robust to noise when the machine learning algorithm worked on data in its raw form. This system is shown to outperform a traditional spectrometer. Techniques for reducing the complexity of the reconstruction routine were also analyzed. Techniques for such as data regularization, deep neural networks, and matrix completion were studied and shown to have benefits over traditional reconstruction techniques. In this project we showed that compressive sensing sensors are indeed feasible at the edge. As always, sensors and algorithms must be carefully tuned to work in the constrained environment. In this project we developed tools and techniques to enable those analyses.

47 OTHER INSTRUMENTATION↗

Ptychographic phase retrieval by proximal algorithms

We derive a set of ptychography phase-retrieval iterative engines based on proximal algorithms originally developed in convex optimization theory, and discuss their connections with existing ones. The use of proximal operator creates a simple frame work that allows us to incorporate the effect of noise from a maximum-likelihood (ML) principle. We focus on three particular algorithms, namely proximal minimization, alternating direction method of multiplier and accelerated proximal gradient (APG). We benchmark their performance with numerical simulations, and discuss their optimal conditions for convergence and accuracy. An experimental dataset is used to demonstrate their effectiveness as well, in which case an array of cubic Au nanoparticles with a size of 50 nm is imaged. We show that with the presence of Poisson noise, a dataset with photon counts up to 10 4 at one detector pixel already requires ML-based methods to achieve a stable convergence. Among the three algorithms derived in this work, APG method is reported first time for its application in ptychographic reconstruction and shows superior performance in terms of both accuracy and convergence rate with a noisy dataset.

36 MATERIALS SCIENCE↗

Shock capturing schemes for multidimensional flow

Progress made in the development of a genuinely multidimensional finite volume algorithm for problems in inviscid gas dynamics is presented. The approach entails: (1) the reconstruction of flowfield data using a planar wave pattern in which the strengths and orientations of the component waves are derived independently of the mesh geometry; and (2) the development of a flux formula which provides a numerical approximation to the flux at a finite volume cell face during the passage of waves which are in general oblique to the face. Several algorithms are outlined, and the most recent developments are included. The results of several numerical cases are also included.

Parpia, Ijaz H.↗

Hubble Space Telescope characterized by using phase-retrieval algorithms

We describe several results characterizing the Hubble Space Telescope from measured point spread functions by using phase-retrieval algorithms. The Cramer-Rao lower bounds show that point spread functions taken well out of focus result in smaller errors when aberrations are estimated and that, for those images, photon noise is not a limiting factor. Reconstruction experiments with both simulated and real data show that the calculation of wave-front propagation by the retrieval algorithms must be performed with a multiple-plane propagation rather than a simple fast Fourier transform to ensure the high accuracy required. Pupil reconstruction was performed and indicates a misalignment of the optical axis of a camera relay telescope relative to the main telescope. After we accounted for measured spherical aberration in the relay telescope, our estimate of the conic constant of the primary mirror of the HST was - 1.0144.

Fienup, J. R.↗

ForSE: A GAN-based Algorithm for Extending CMB Foreground Models to Subdegree Angular Scales

We present ForSE (Foreground Scale Extender), a novel Python package that aims to overcome the current limitations in the simulation of diffuse Galactic radiation, in the context of cosmic microwave background (CMB) experiments. ForSE exploits the ability of generative adversarial neural networks (GANs) to learn and reproduce complex features present in a set of images, with the goal of simulating realistic and non-Gaussian foreground radiation at subdegree angular scales. This is of great importance in order to estimate the foreground contamination to lensing reconstruction, delensing, and primordial B-modes for future CMB experiments. We applied this algorithm to Galactic thermal dust emission in both total intensity and polarization. Our results show how ForSE is able to generate small-scale features (at 12') having as input the large-scale ones (80'). The injected structures have statistical properties, evaluated by means of the Minkowski functionals, in good agreement with those of the real sky and which show the correct amplitude scaling as a function of the angular dimension. Furthermore, the obtained thermal dust Stokes Q and U full-sky maps as well as the ForSE package are publicly available for download.

79 ASTRONOMY AND ASTROPHYSICS↗