Search NASA⌕ Search

SEARCH · Search NASA

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 289 records · Page 16

Advancing the Frontiers of Deep Learning for Low-Dose 3D Cone-Beam CT Reconstruction

X-ray computed tomography (CT) is an important noninvasive medical imaging modality for studying the structural details of internal organs. Image reconstruction in CT is an inverse problem of recovering an object's internal structure from the absorption profile of X-ray beams (sinogram) measured using a detector. The classical variational approach for CT reconstruction minimizes an energy functional using an appropriate iterative algorithm. Motivated by the success of deep learning (DL), researchers have begun to leverage training data and enhanced computing capabilities in recent years to produce high-fidelity reconstructed images. Nonetheless, much of the academic research in DL algorithms for CT has focused primarily on the two-dimensional setting (with simplified forward operators and noise model) for proofs-of-concept, and a comprehensive benchmarking of various classical and data-driven CT reconstruction approaches has not beenundertaken. The key objective of our CT reconstruction grand challenge was to promote methodological advancements for both classical and DL-based approaches for clinical CT with a reasonably accurately simulated 3D CT forward operator and noise model. We have utilized the publicly available LIDC-IDRI dataset and simulated sinograms and FDK images corresponding to two dose levels (clinical- and low-dose, constituting two tracks of the challenge) starting from the normal-dose images as the ground truth. In this paper, we summarize the motivation, context, and results of our challenge, and highlight the future research directions in DL for clinical CT.

X-ray tomography↗

Total Variation Majorization Minimization (TV-MM) Approach to Radiometer Brightness Temperature Gridding and Reconstruction

This paper presents the implementation of an algorithm to enhance the image resolution of the Earth's surface brightness temperature (T B ) data measured by radiometers such as the one onboard of the Soil Moisture Active Passive (SMAP) mission. A key step in radiometer T B processing is the conversion of the swath-based calibrated antenna temperature (T A ) measurements to the Level 3 Earth-centered grid. The simplest algorithm to transform this data from swath to gridded format is called drop-in-the-bucket which simply averages surrounding noisy T A samples to form a T B value at the gridded location. This method reduces noise, however produces low resolution products. To obtain a higher resolution product, SMAP uses other techniques such the Backus-Gilbert (BG) algorithm, which is the conventional method used in microwave radiometry. Although this method performs the required interpolation, it is not effective in denoising and removing blurring effects due to antenna filtering of the radiometer image data. Our motivation for this development is to further improve the resolution through post-processing of the radiometer T B image, a highly cost-effective method of image enhancement. The approach adapted in this work is based on the minimization of the Total Variation (TV) regularized objective function that is used extensively in solving general ill-posed linear inverse problems in image processing. Since the TV-based objective function is convex but not everywhere differentiable, there exists many numerical algorithms that can estimate the solution and the one selected for this work is called Majorization- Minimization (MM). By applying this algorithm, simulation experiments were performed based on synthetic data from the Geophysical model as well as real SMAP data to demonstrate the effectiveness of the technique. Results were then compared against the BG method.

Wing Lee↗

Tomography by iterative convolution - Empirical study and application to interferometry

An algorithm for computer tomography has been developed that is applicable to reconstruction from data having incomplete projections because an opaque object blocks some of the probing radiation as it passes through the object field. The algorithm is based on iteration between the object domain and the projection (Radon transform) domain. Reconstructions are computed during each iteration by the well-known convolution method. Although it is demonstrated that this algorithm does not converge, an empirically justified criterion for terminating the iteration when the most accurate estimate has been computed is presented. The algorithm has been studied by using it to reconstruct several different object fields with several different opaque regions. It also has been used to reconstruct aerodynamic density fields from interferometric data recorded in wind tunnel tests.

Vest, C. M.↗

Real-Time Artificial Intelligence for Particle Reconstruction and Higgs Physics

With the discovery of the Higgs boson at the CERN LHC, the world's highest-energy particle accelerator complex, scientists have acquired an important tool to study the fundamental building blocks of the universe. Precision measurements of Higgs bosons produced with large momentum allow for unique insights into the structure of the interactions of the Higgs boson with other particles that may shed light on physics beyond the standard model. While experimentally challenging, exploring such interactions with novel artificial intelligence (AI) methods can advance our understanding of the Higgs sector, including the Higgs boson's self-interaction. Moreover, the LHC is undergoing a major upgrade to further increase its particle collision rate and thereby operate for an additional decade. The experimental detectors at the upgraded facility must process at least a factor of ten more data at rates of hundreds of terabytes per second all under challenging conditions. New AI techniques are required to reconstruct and select, or trigger on, the most physics-sensitive events in real-time to handle the resulting avalanche of data. The proposed research will achieve the goals of the LHC program at the CMS experiment by developing a sub-microsecond event reconstruction system using real-time AI algorithms that employ field-programmable gate array technologies. By harnessing sophisticated AI techniques, this research focuses on measuring the production of Higgs bosons at large momentum while enhancing particle reconstruction methods in the trigger and beyond. Overall, the proposed research has broader implications for the use of AI in resource-constrained, low-latency embedded applications across all fields of science.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Direct Numerical Simulation of Acoustic Waves Interacting with a Shock Wave in a Quasi-1D Convergent-Divergent Nozzle Using an Unstructured Finite Volume Algorithm

Numerical simulation of a very small amplitude acoustic wave interacting with a shock wave in a quasi-1D convergent-divergent nozzle is performed using an unstructured finite volume algorithm with a piece-wise linear, least square reconstruction, Roe flux difference splitting, and second-order MacCormack time marching. First, the spatial accuracy of the algorithm is evaluated for steady flows with and without the normal shock by running the simulation with a sequence of successively finer meshes. Then the accuracy of the Roe flux difference splitting near the sonic transition point is examined for different reconstruction schemes. Finally, the unsteady numerical solutions with the acoustic perturbation are presented and compared with linear theory results.

Bui, Trong T.↗

Energy functions for regularization algorithms

Regularization techniques are widely used for inverse problem solving in computer vision such as surface reconstruction, edge detection, or optical flow estimation. Energy functions used for regularization algorithms measure how smooth a curve or surface is, and to render acceptable solutions these energies must verify certain properties such as invariance with Euclidean transformations or invariance with parameterization. The notion of smoothness energy is extended here to the notion of a differential stabilizer, and it is shown that to void the systematic underestimation of undercurvature for planar curve fitting, it is necessary that circles be the curves of maximum smoothness. A set of stabilizers is proposed that meet this condition as well as invariance with rotation and parameterization.

Delingette, H.↗

A Parallel, Finite-Volume Algorithm for Large-Eddy Simulation of Turbulent Flows

A parallel, finite-volume algorithm has been developed for large-eddy simulation (LES) of compressible turbulent flows. This algorithm includes piecewise linear least-square reconstruction, trilinear finite-element interpolation, Roe flux-difference splitting, and second-order MacCormack time marching. Parallel implementation is done using the message-passing programming model. In this paper, the numerical algorithm is described. To validate the numerical method for turbulence simulation, LES of fully developed turbulent flow in a square duct is performed for a Reynolds number of 320 based on the average friction velocity and the hydraulic diameter of the duct. Direct numerical simulation (DNS) results are available for this test case, and the accuracy of this algorithm for turbulence simulations can be ascertained by comparing the LES solutions with the DNS results. The effects of grid resolution, upwind numerical dissipation, and subgrid-scale dissipation on the accuracy of the LES are examined. Comparison with DNS results shows that the standard Roe flux-difference splitting dissipation adversely affects the accuracy of the turbulence simulation. For accurate turbulence simulations, only 3-5 percent of the standard Roe flux-difference splitting dissipation is needed.

Bui, Trong T.↗

Dual X-ray computed tomography-aided classification of melt pool boundaries and flaws in crept additively manufactured parts

In metal additive manufacturing (AM), understanding the process-structure-performance relationships requires a combination of multi-scale characterization techniques that allows for the measurement of the melt pool shape and boundary and classifying various defects and flaws in the AM parts. Such approaches can be destructive, only 2D in nature, or have a small field of view and can be complex to co-register and analyze. Here, in this work, we present a non-destructive 3D inspection technique that employs dual-energy X-ray computed tomography (XCT) along with a model-based iterative reconstruction (MBIR) and a new segmentation algorithm. The proposed approach and algorithm are not only capable of classifying and quantifying flaws such as pores, cracks, and inclusions, but they also allow for the extraction of microstructural features such as melt pool boundaries (MPB) and melt pool regions (MPR), that can help understand process-structure-performance relationships for alloys under study. As an exemplar application, we employed the method for characterization of an additively manufactured aluminum alloy crept under tensile stress at 300 °C for 1064 h. Our results demonstrate high quality segmentation and classification of various flaws and MPB and MPR, for the first time, using 3D X-ray CT inspection. The delineated MPB and MPR in the crept samples reveal the preferential growth paths of cracks that formed during creep deformation. The technique was used for successfully quantifying the characteristics (number of defects, their density, volume fraction, etc.) of the manufacturing-induced pores and creep-induced cracks, which is necessary to better understand the creep failure mechanisms of the material.

36 MATERIALS SCIENCE↗

Improving missing transverse momentum estimation with a deep neural network

At hadron colliders, the net transverse momentum of particles that do not interact with the detector (missing transverse momentum, $^→_𝑝$$^{miss}_{T}$) is a crucial observable in many analyses. In the standard model, $^→_𝑝$$^{miss}_{T}$ originates from neutrinos. Many beyond-the-standard-model particles, such as dark matter candidates, are also expected to leave the experimental apparatus undetected. This paper presents a novel deep neural network based $^→_𝑝$$^{miss}_{T}$ estimator, DeepMET, developed by the CMS Collaboration at the LHC. The DeepMET algorithm produces a weight for each reconstructed particle based on its properties. The estimator is based on the negative vector sum of the weighted transverse momenta of all reconstructed particles in an event. Compared with other estimators currently employed by CMS, DeepMET improves the $^→_𝑝$$^{miss}_{T}$ resolution by 10%–30%, shows improvement for a wide range of final states, is easier to train, and is more resilient against the effects of additional proton-proton interactions accompanying the collision of interest.

artificial neural networks↗

Observation of Enhanced Long-Range Elliptic Anisotropies Inside High-Multiplicity Jets in p p Collisions at s = 13 TeV

A search for collective effects inside jets produced in proton-proton collisions is performed via correlation measurements of charged particles using the CMS detector at the CERN LHC. The analysis uses data collected at a center-of-mass energy of s = 13 TeV , corresponding to an integrated luminosity of 138 fb − 1 . Jets are reconstructed with the anti- k T algorithm with a distance parameter of 0.8 and are required to have transverse momentum greater than 550 GeV and pseudorapidity | η jet | < 1.6 . Two-particle correlations among the charged particles within the jets are studied as functions of the particles’ azimuthal angle and pseudorapidity separations ( Δ ϕ * and Δ η * ) in a jet coordinate basis, where particles’ η * , ϕ * are defined relative to the direction of the jet. The correlation functions are studied in classes of in-jet charged-particle multiplicity up to N ch j ≈ 100 . Fourier harmonics are extracted from long-range azimuthal correlation functions to characterize azimuthal anisotropy for | Δ η * | > 2 . For low- N ch j jets, the long-range elliptic anisotropic harmonic, v 2 * , is observed to decrease with N ch j . This trend is well described by Monte Carlo event generators. However, a rising trend for v 2 * emerges at N ch j ≳ 80 , hinting at a possible onset of collective behavior, which is not reproduced by the models tested. This observation yields new insights into the dynamics of jet evolution in the vacuum. © 2024 CERN, for the CMS Collaboration 2024 CERN

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Simulated Feasibility of 3-D Lightning Mapping From Space

In addition to the awe it inspires, lightning can illuminate the microphysical processes hidden away within deep convection. The current generation of space-based lightning mapping uses mostly 2-D optical imaging to connect overall flash characteristics to their parent storm dynamics, but are missing a dimension’s worth of information. With lightning now classified as an essential climate variable, future spaceborne mappers will need improved capabilities to take advantage of the 3-D structure of lightning flashes to support meteorological and climate modeling. We report here on a study of the feasibility of high-resolution 3-D lightning mapping using a radio frequency (RF)-based network of satellites from low-Earth orbit (LEO). Lightning sources are simulated using existing lightning mapping array (LMA) tools, modified for orbital detection, and spatially reconstructed using a Levenberg–Marquardt geolocation algorithm to assess sources of uncertainty in these solutions. We analyze the benefits and limitations of this approach compared to existing orbital and ground-based methods. Results of this study show that lightning can be mapped in 3-D with a vertical location accuracy better than 2 km using as few as five satellites in LEO capable of measuring the time-of-arrival of impulsive RF signals in the very high-frequency (VHF) band. The consequence of this study is that high-resolution, spaceborne 3-D mapping of lightning is achievable across most of the globe, having crucial implications for our understanding of not only lightning, but also severe weather development, climate science, and more.

47 OTHER INSTRUMENTATION↗

Measurement of multidifferential cross sections for dijet production in proton–proton collisions at $\sqrt{s}$ = 13 TeV

A measurement of the dijet production cross section is reported based on proton–proton collision data collected in 2016 at $\sqrt{s}$ = 13 TeV by the CMS experiment at the CERN LHC, corresponding to an integrated luminosity of up to 36.3 fb –1 . Jets are reconstructed with the anti-$k$ T algorithm for distance parameters of R = 0.4 and 0.8. Cross sections are measured double-differentially (2D) as a function of the largest absolute rapidity |y| max of the two jets with the highest transverse momenta $p$ T and their invariant mass $m$ 1,2 , and triple-differentially (3D) as a function of the rapidity separation $y$*, the total boost $y$ b , and either $m$ 1,2 or the average $p$ T of the two jets. The cross sections are unfolded to correct for detector effects and are compared with fixed-order calculations derived at next-to-next-to-leading order in perturbative quantum chromodynamics. The impact of the measurements on the parton distribution functions and the strong coupling constant at the mass of the Z boson is investigated, yielding a value of $α$ S ($m$ Z ) = 0.1179 ± 0.0019.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Speckle imaging for planetary research

The present study of speckle imaging technique effectiveness encompasses image reconstruction by means of a division algorithm for Fourier amplitudes, and the Knox-Thompson (1974) algorithm for Fourier phases. Results which have been obtained for Io, Titan, Pallas, Jupiter and Uranus indicate that spatial resolutions lower than the seeing limit by a factor of four are obtainable for objects brighter than Uranus. The resolutions obtained are well above the diffraction limit, due to inadequacies of the video camera employed. A photon-counting camera has been developed to overcome these difficulties, making possible the diffraction-limited resolution of objects as faint as Charon.

Nisenson, P.↗

A laboratory demonstration of high-resolution hard X-ray and gamma-ray imaging using Fourier-transform techniques

A laboratory imaging system has been developed to study the use of Fourier-transform techniques in high-resolution hard X-ray and gamma-ray imaging, with particular emphasis on possible applications to high-energy astronomy. Considerations for the design of a Fourier-transform imager and the instrumentation used in the laboratory studies is described. Several analysis methods for image reconstruction are discussed including the CLEAN algorithm and maximum entropy methods. Images obtained using these methods are presented.

Palmer, David↗

Non-Spherical Aerosol Phase Functions Derived from MODIS and AERONET Observations

We compare MODIS (Moderate Resolution Imaging Spectroradiometer) satellite aerosol retrievals of spectral optical thickness and size parameters over ocean with the same quantities derived from AERONET (Aerosol Robotic Network) observations made at island and coastal sites. Over much of the globe, the satellite-derived quantities agree well with the AERONET quantities. However, in regimes dominated by desert dust aerosol, the agreement is less robust. In the dust regimes, the MODIS retrievals show greater spectral dependence and report smaller particle sizes than do the AERONET derivations. We suggest that the reason for this discrepancy is the nonspherical nature of desert dust particles, which the initial MODIS algorithm is not able to handle. Using the discrepancy between MODIS and AERONET derived spectral optical thickness as an asset, instead of a detriment, we reconstruct the aerosol phase functions that the MODIS algorithm would have needed in order to match the AERONET retrievals. No assumptions of particle shape are used in the derivation of these functions and the results are empirical total column, ambient phase functions. We compare the empirically derived phase functions with phase functions calculated from spheres and spheroids, both situations in which assumptions about particle shape must be made. The resulting empirical nonspherical phase functions will be included in future updates of the MODIS algorithm.

Remer, L. A.↗

Low-Complexity Lossless and Near-Lossless Data Compression Technique for Multispectral Imagery

This work extends the lossless data compression technique described in Fast Lossless Compression of Multispectral- Image Data, (NPO-42517) NASA Tech Briefs, Vol. 30, No. 8 (August 2006), page 26. The original technique was extended to include a near-lossless compression option, allowing substantially smaller compressed file sizes when a small amount of distortion can be tolerated. Near-lossless compression is obtained by including a quantization step prior to encoding of prediction residuals. The original technique uses lossless predictive compression and is designed for use on multispectral imagery. A lossless predictive data compression algorithm compresses a digitized signal one sample at a time as follows: First, a sample value is predicted from previously encoded samples. The difference between the actual sample value and the prediction is called the prediction residual. The prediction residual is encoded into the compressed file. The decompressor can form the same predicted sample and can decode the prediction residual from the compressed file, and so can reconstruct the original sample. A lossless predictive compression algorithm can generally be converted to a near-lossless compression algorithm by quantizing the prediction residuals prior to encoding them. In this case, since the reconstructed sample values will not be identical to the original sample values, the encoder must determine the values that will be reconstructed and use these values for predicting later sample values. The technique described here uses this method, starting with the original technique, to allow near-lossless compression. The extension to allow near-lossless compression adds the ability to achieve much more compression when small amounts of distortion are tolerable, while retaining the low complexity and good overall compression effectiveness of the original algorithm.

Xie, Hua↗

Reconstruction of internal density distributions in porous bodies from laser ultrasonic data

It is presently shown that, for density-reconstruction problems in which information about the inhomogeneity is known a priori, the nonlinear least-squares algorithm yields satisfactory results on the basis of limited projection data. The back-projection algorithm, which obviates assumptions about the objective function to be reconstructed, does not recover the boundary of the inhomogeneity when the number of projections is limited and ray-bending is ignored.

Lu, Yichi↗

Calibration of a synthetic aperture radiometer

Calibration algorithms for a synthetic aperture microwave radiometer are presented. The calibration is geared to Earth remote sensing applications, and is demonstrated on an airborne prototype of thinned array imager. Two approaches to the system calibration are presented; the first utilizes commonly available reference brightness temperature scenes, such as open water, and the second utilizes data collected on the antenna range. Both algorithms yield spatial response information which is cast in matrix form and inverted to obtain the image reconstruction formula. Experimental results are examined, and errors in some reconstructed images are linked to the present prototype antenna design. Algorithms for improving the synthesized antenna pattern side lobe performance are also presented. In one solution, the pattern efficiency is optimized by minimizing the pattern outside a defined beam. In another solution, the patterns are matched to a desired model pattern by the method of least squared errors. Both techniques offer an attractive alternative to aperture weighting.

Tanner, Alan B.↗