Results for “reconstruction algorithms” — Page 23

DOE OSTI2025

A Quantum Leap for Dynamic Radiography

Dynamic radiography techniques have been instrumental in advancing the National Nuclear Security Administration's (NNSA) mission since WWII. This paper explores the transformative potential of quantum information science (QIS) to revolutionize dynamic radiography through enhanced image processing, tomographic reconstruction, statistics and uncertainty quantification (UQ), and artificial intelligence integration. By leveraging quantum algorithms to extract previously inaccessible information from existing datasets, this interdisciplinary approach promises unprecedented insights at the intersection of dynamic imaging, artificial intelligence, and quantum information technologies.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS1986

High resolution observations of the L1551 bipolar outflow

The nearby dark cloud Lynds 1551 contains one of the closest examples of a well-collimated bipolar molecular outflow. This source has the largest angular size of any known outflow and was the first bipolar outflow to be detected. The outflow originates from a low-luminosity young stellar object, IRS-5. Optical and radio continuum observations show the presence of a highly collimated, ionized stellar wind orginating from close to IRS-5 and aligned with the molecular outflow. However, we have little information on the actual mechanism that generates the stellar wind and collimates it into opposed jets. The Very Large Array (VLA) observations indicate that the winds originate within 10(15) cm of IRS-5, unfortunately at a size scale difficult to resolve. For these reasons, observations of the structure and dynamics of the hypersonic molecular gas may provide valuable information on the origin and evolution of these outflows. In addition, the study of the impact of the outflowing gas on the surrounding molecular material is essential to understand the consequence these outflows have on the evolution and star formation history of the entire cloud. Moriarty-Schieven et al. (1986) obtained a oversampled map of the CO emission of a portion of both the blueshifted and redshifted outflows in LI551 using Five College Radio Astronomy Observatory 14 m telescope. The oversampled maps have been reconstructed to an effective angular resolution of 20 arcsec using a maximum entropy algorithm. A continuation of the study of Moriarty-Schieven et al. is presented. The entire L1551 outflow has now been mapped at 12 arcsec sampling requiring roughly 4000 spectra. This data has been constructed to 20 arcsec resolution to provide the first high resolution picture of the entire L1551 outflow. This new data has shown that the blueshifted lobe is more extended than previously thought and has expanded downstream sufficiently to break out of the dense molecular cloud, but the redshifted outflow is still confined within the molecular cloud. Details of the structure and kinematics of the high velocity gas are used to test the various models of the origin and evolution of outflows.

Snell, R.↗

DOE OSTI2025

Characterization of Pinhole Collimators for High-Resolution Gamma Imaging of Irradiated Fuel

Post-irradiation examination (PIE) of nuclear fuels requires imaging tools capable of resolving isotopic and spatial features with high throughput. This project contributes to a proof-of-concept effort aimed at advancing gamma emission tomography (GET) by evaluating novel fine-aperture pinhole collimators. Two Rose’s metal collimators, 100 µm (20° acceptance angle) and 350 µm (30° acceptance angle), were prototyped and characterized for their effectiveness in transporting gamma rays through the pinhole aperture. To support data collection, a Python-based data acquisition system was developed to coordinate a rotation stage, linear stage, and CZT detector, reducing latency in high-rate gamma event logging to one second per acquisition. Queue-based file writing enabled seamless real-time data capture for count rates up to 35,000 counts per second (cps). List-mode parsing algorithms were implemented to differentiate single and simultaneous gamma interactions for future tomographic reconstruction. Detector response was evaluated in both spectroscopy and list mode acquisition methods across varying source distances to confirm absolute and collimator efficiencies. Preliminary efficiency figures suggest effective collimation of gamma-rays with energies below 700 keV, with ~4% residual intensity through the aperture for Cs-137. The impact of collimator geometry on image quality is currently being evaluated. This groundwork supports the ongoing development of a sub mm resolution cone-beam CT system for imaging fuel phantoms, an essential step toward improving GET efficiency and accelerating nuclear fuel qualification efforts.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

NASA NTRS1988

New insights in the numerical reliability properties of existing Kalman filter implementations

The convergence properties of Kalman filter algorithms are investigated analytically. A theoretical error analysis is performed on four types of algorithms, as defined by Anderson and Moore (1979): (1) a conventional Kalman filter, (2) a square-root (SR) covariance filter, (3) the Chandrasekhar SR filter, and (4) an SR information filter. The derivations are given in detail, and numerical results for the flight-path reconstruction problem studied by Verhaegen (1987) are presented in tables and graphs. It is shown that error propagation in algorithms (1) and (2) is sensitive to the condition number of the innovation-signal covariance matrix and the spectral norm of the filter state-transition matrix, whereas other parameters are dominant in (3) and (4). Filter (2) is found to be the most reliable for the class of problems studied.

Verhaegen, M. H.↗

NASA NTRS1996

Machine-Vision Aids for Improved Flight Operations

The development of machine vision based pilot aids to help reduce night approach and landing accidents is explored. The techniques developed are motivated by the desire to use the available information sources for navigation such as the airport lighting layout, attitude sensors and Global Positioning System to derive more precise aircraft position and orientation information. The fact that airport lighting geometry is known and that images of airport lighting can be acquired by the camera, has lead to the synthesis of machine vision based algorithms for runway relative aircraft position and orientation estimation. The main contribution of this research is the synthesis of seven navigation algorithms based on two broad families of solutions. The first family of solution methods consists of techniques that reconstruct the airport lighting layout from the camera image and then estimate the aircraft position components by comparing the reconstructed lighting layout geometry with the known model of the airport lighting layout geometry. The second family of methods comprises techniques that synthesize the image of the airport lighting layout using a camera model and estimate the aircraft position and orientation by comparing this image with the actual image of the airport lighting acquired by the camera. Algorithms 1 through 4 belong to the first family of solutions while Algorithms 5 through 7 belong to the second family of solutions. Algorithms 1 and 2 are parameter optimization methods, Algorithms 3 and 4 are feature correspondence methods and Algorithms 5 through 7 are Kalman filter centered algorithms. Results of computer simulation are presented to demonstrate the performance of all the seven algorithms developed.

Menon, P. K.↗

NASA NTRS1997

An Adaptive Flow Solver for Air-Borne Vehicles Undergoing Time-Dependent Motions/Deformations

This report describes a concurrent Euler flow solver for flows around complex 3-D bodies. The solver is based on a cell-centered finite volume methodology on 3-D unstructured tetrahedral grids. In this algorithm, spatial discretization for the inviscid convective term is accomplished using an upwind scheme. A localized reconstruction is done for flow variables which is second order accurate. Evolution in time is accomplished using an explicit three-stage Runge-Kutta method which has second order temporal accuracy. This is adapted for concurrent execution using another proven methodology based on concurrent graph abstraction. This solver operates on heterogeneous network architectures. These architectures may include a broad variety of UNIX workstations and PCs running Windows NT, symmetric multiprocessors and distributed-memory multi-computers. The unstructured grid is generated using commercial grid generation tools. The grid is automatically partitioned using a concurrent algorithm based on heat diffusion. This results in memory requirements that are inversely proportional to the number of processors. The solver uses automatic granularity control and resource management techniques both to balance load and communication requirements, and deal with differing memory constraints. These ideas are again based on heat diffusion. Results are subsequently combined for visualization and analysis using commercial CFD tools. Flow simulation results are demonstrated for a constant section wing at subsonic, transonic, and a supersonic case. These results are compared with experimental data and numerical results of other researchers. Performance results are under way for a variety of network topologies.

Singh, Jatinder↗

DOE OSTI2025

Streaming Compression of Scientific Data via Weak-SINDy

Here, in this paper, a streaming weak-SINDy algorithm is developed specifically for compressing streaming scientific data. The production of scientific data, either via simulation or experiments, is undergoing a stage of exponential growth, which makes data compression important and often necessary for storing and utilizing large scientific data sets. As opposed to classical “offline” compression algorithms that perform compression on a readily available data set, streaming compression algorithms compress data “online” while the data generated from simulation or experiments is still flowing through the system. This feature makes streaming compression algorithms well suited for scientific data compression, where storing the full data set offline is often infeasible. This work proposes a new streaming compression algorithm, streaming weak-SINDy, which takes advantage of the underlying data characteristics during compression. The streaming weak-SINDy algorithm constructs feature matrices and target vectors in the online stage via a streaming integration method in a memory efficient manner. The feature matrices and target vectors are then used in the offline stage to build a model through a regression process that aims to recover equations that govern the evolution of the data. For compressing high-dimensional streaming data, we adopt a streaming proper orthogonal decomposition (POD) process to reduce the data dimension and then use the streaming weak-SINDy algorithm to compress the temporal data of the POD expansion. We propose modifications to the streaming weak-SINDy algorithm to accommodate the dynamically updated POD basis. By combining the built model from the streaming weak-SINDy algorithm and a small amount of data samples, the full data flow could be reconstructed accurately at a low memory cost, as shown in the numerical tests.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1973

ERTS image data compression technique evaluation

The author has identified the following significant results. Tapes of compressed ERTS data were obtained to permit later reconstruction and to prove that in general four ERTS MSS tapes can be put onto a single compressed tape. A compressed tape was reconstructed and imagery made. The data were compressed using the essentially information preserving SSDIAM algorithm, with mappings of from 1 to 3 levels and imagery was made of the result. This imagery shows that no visual degradation results from the one level mapping while compression is significantly increased. Mappings of up to three levels shows negligible deterioration in areas of moderate to high data activity, but contouring is noticeable in areas of uniform data such as the plains region.

Spencer, D. J.↗

NASA NTRS1983

Sensory redundancy management: The development of a design methodology for determining threshold values through a statistical analysis of sensor output data

Sensor redundancy management (SRM) requires a system which will detect failures and reconstruct avionics accordingly. A probability density function to determine false alarm rates, using an algorithmic approach was generated. Microcomputer software was developed which will print out tables of values for the cummulative probability of being in the domain of failure; system reliability; and false alarm probability, given a signal is in the domain of failure. The microcomputer software was applied to the sensor output data for various AFT1 F-16 flights and sensor parameters. Practical recommendations for further research were made.

Scalzo, F.↗

NASA NTRS1989

Effect Of Noise In The Ideal State Reconstructor

More measurements yield better estimate. Report discusses effects of measurement noise on system including linear, time-invariant plant of measurable inputs and known parameters and deterministic digital control subsystem governed by algorithm called "ideal state reconstructor." So named because in absence of noise, it exactly reconstructs vector represting state of plant, even without knowledge of initial state. Ideal state reconstructor adds no new states or eigenvalues to system; affects measurement equation only.

Polites, Michael E.↗

NASA NTRS1989

Kalman estimation for SEDS measurements

The first on-orbit experiment of the Small Expendable Deployer System (SEDS) for tethered satellites will collect telemetry data for tethered length, rate of deployment, and tether tension. The post-flight analysis will use this data to reconstruct the deployment history and determine payload position and tether shape. Two Kalman estimator algorithms were written, and output using simulated measurement data was compared. Both estimators exhibited the same estimated state histories, indicating that numerical instability in the traditional algorithm was not the cause of filter divergence. Estimation of acceleration biases was added, which reduced the error but did not correct the divergence. An add-a-bead estimator that adds lumped masses as the tether is deployed was written, which provides a state model that matches the BEADSIM simulation providing the true measurements and states. This twenty-one bead estimator produced state histories similar to those of the two-bead estimator, indicating that the filter divergence was not caused by a reduced-order model. The noise models used to date are relatively simple and may be the source of estimator divergence. The investigation of colored noise models, cross-correlated measurement and process covariances, and noise-adaptive filter techniques is recommended.

Carrington, Connie K.↗

NASA NTRS2013

Semi-Supervised Eigenbasis Novelty Detection

Recent discoveries in high-time-resolution radio astronomy data have focused attention on a new class of events. Fast transients are rare pulses of radio frequency energy lasting from microseconds to seconds that might be produced by a variety of exotic astrophysical phenomena. For example, X-ray bursts, neutron stars, and active galactic nuclei are all possible sources of short-duration, transient radio signals. It is difficult to anticipate where such signals might appear, and they are most commonly discovered through analysis of high-time- resolution data that had been collected for other purposes. Transients are often faint and difficult to detect, so improved detection algorithms can directly benefit the science yield of all such commensal monitoring. A new detection algorithm learns a low-dimensional linear manifold for describing the normal data. High reconstruction error indicates a novel signal that does not match the patterns of normal data. One unsupervised portion of the manifold model adapts its representation in response to recent data. A second supervised portion of the model is made of a basis trained in advance using labeled examples of RFI; this prevents false positives due to these events. For a linear model, an orthonormalization operation is used to combine these bases prior to the anomaly detection decision. Another novel aspect of the approach lies in combining basis vectors learned in an unsupervised, online fashion from the data stream with supervised basis vectors learned in advance from known examples of false alarms. Adaptive, data-driven detection is achieved that is also informed by existing domain knowledge about signals that may be statistically anomalous, but are not interesting and should therefore be ignored. The method was evaluated using data from the Parkes Multibeam Survey. This data set was originally collected to search for pulsars, which are astronomical sources that emit radio pulses at regular periods. However, several non-pulsar anomalies have recently been discovered in this dataset, making it a compelling test case. By explicitly filtering known false alarm patterns, the approach yields significantly better performance than current transient detection methods.

Wagstaff, Kiri L.↗

NASA NTRSDate not supplied

Investigating Photogrammetric Accuracy of a Lunar-lander-induced Crater Measurement System

Laboratory measurements have been made to validate the performance of the Stereo CAmeras for Lunar Plume-Surface Studies (SCALPSS) stereo photogrammetry systems which will be flying to the moon on two of NASA’s upcoming Commercial Lunar Payload Services (CLPS) missions. Until recently, the system’s accuracy had only been studied using idealized geometric shapes as measurement targets. A realistic crater model of representative scale and an idealized ‘staircase’ target have been used to compare measurement accuracy of ideal versus lunar-like objects, with the commercial V-STARS® system being used to provide the known reference values for comparison. In a parametric study, altitude, lens focal length, and camera separation are varied to assess each parameter’s impact on photogrammetric accuracy in relation to the scaling law prediction developed previously. The SCALPSS 1.0 and 1.1 configurations have been validated on the crater model within acceptable accuracy for the missions, performing significantly better than the scaling law prediction in some cases. A semi-automated post-processing routine was developed in MATLAB® and proved successful for the cross-correlation of features between two stereo images. For some cases of extreme convergence angles between a camera pair, manual feature detection and matching was required. By using this manual process, the crater depth map was reconstructed but with worse accuracy than the idealized staircase measurements; refinements to the processing algorithm are expected to improve future results. Also examined in this work is the impact of illumination environments, both natural (e.g., Sun angles) and artificial (diffuse or structured illumination sources), on the camera system’s ability to measure the erosion of the lunar terrain.

Plume-surface interaction↗

NASA NTRSDate not supplied

Investigating Photogrammetric Accuracy of a Lunar-lander-induced Crater Measurement System

Laboratory measurements have been made to validate the performance of the Stereo CAmeras for Lunar Plume-Surface Studies (SCALPSS) stereo photogrammetry systems which will be flying to the moon on two of NASA’s upcoming Commercial Lunar Payload Services (CLPS) missions. Until recently, the system’s accuracy had only been studied using idealized geometric shapes as measurement targets. A realistic crater model of representative scale and an idealized ‘staircase’ target have been used to compare measurement accuracy of ideal versus lunar-like objects, with the commercial V-STARS® system being used to provide the known reference values for comparison. In a parametric study, altitude, lens focal length, and camera separation are varied to assess each parameter’s impact on photogrammetric accuracy in relation to the scaling law prediction developed previously. The SCALPSS 1.0 and 1.1 configurations have been validated on the crater model within acceptable accuracy for the missions, performing significantly better than the scaling law prediction in some cases. A semi-automated post-processing routine was developed in MATLAB® and proved successful for the cross-correlation of features between two stereo images. For some cases of extreme convergence angles between a camera pair, manual feature detection and matching was required. By using this manual process, the crater depth map was reconstructed but with worse accuracy than the idealized staircase measurements; refinements to the processing algorithm are expected to improve future results. Also examined in this work is the impact of illumination environments, both natural (e.g., Sun angles) and artificial (diffuse or structured illumination sources), on the camera system’s ability to measure the erosion of the lunar terrain.

Plume-surface interaction↗

DOE OSTI2024

PAIReD jet: A multi-pronged resonance tagging strategy across all Lorentz boosts

We propose a new approach of jet-based event reconstruction that aims to optimally exploit correlations between the products of a hadronic multi-pronged decay across all Lorentz boost regimes. The new approach utilizes clustered small-radius jets as seeds to define unconventional jets, referred to as PAIReD jets. The constituents of these jets are subsequently used as inputs to machine learning-based algorithms to identify the flavor content of the jet. We demonstrate that this approach achieves higher efficiencies in the reconstruction of signal events containing heavy-flavor jets compared to other event reconstruction strategies at all Lorentz boost regimes. Classifiers trained on PAIReD jets also have significantly better background rejections compared to those based on traditional event reconstruction approaches using small-radius jets at low Lorentz boost regimes. The combined effect of a higher signal reconstruction efficiency and better classification performance results in a two to four times stronger rejection of light-flavor jets compared to conventional strategies at low Lorentz-boosts, and rejection rates similar to classifiers based on large-radius multi-pronged jets at high Lorentz-boost regimes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS2022

Quantifying Errors in 3D CME Parameters Derived from Synthetic Data Using White-Light Reconstruction Techniques

Current efforts in space weather forecasting of CMEs have been focused on predicting their arrival time and magnetic structure. To make these predictions, methods have been developed to derive the true CME speed, size, position, and mass, among others. Difficulties in determining the input parameters for CME forecasting models arise from the lack of direct measurements of the coronal magnetic fields and uncertainties in estimating the CME 3D geometric and kinematic parameters after eruption. White-light coronagraph images are usually employed by a variety of CME reconstruction techniques that assume more or less complex geometries. This is the first study from our International Space Science Institute (ISSI) team “Understanding Our Capabilities in Observing and Modeling Coronal Mass Ejections”, in which we explore how subjectivity affects the 3D CME parameters that are obtained from the Graduated Cylindrical Shell (GCS) reconstruction technique, which is widely used in CME research. To be able to quantify such uncertainties, the “true” values that are being fitted should be known, which are impossible to derive from observational data. We have designed two different synthetic scenarios where the “true” geometric parameters are known in order to quantify such uncertainties for the first time. We explore this by using two sets of synthetic data: 1) Using the ray-tracing option from the GCS model software itself, and 2) Using 3D magnetohydrodynamic (MHD) simulation data from the Magnetohydrodynamic Algorithm outside a Sphere code. Our experiment includes different viewing configurations using single and multiple viewpoints. CME reconstructions using a single viewpoint had the largest errors and error ranges overall for both synthetic GCS and simulated MHD white-light data. As the number of viewpoints increased from one to two, the errors decreased by approximately 4° in latitude, 22° in longitude, 14° in tilt, and 10° in half-angle. Our results quantitatively show the critical need for at least two viewpoints to be able to reduce the uncertainty in deriving CME parameters. We did not find a significant decrease in errors when going from two to three viewpoints for our specific hypothetical three spacecraft scenario using synthetic GCS white-light data. As we expected, considering all configurations and numbers of viewpoints, the mean absolute errors in the measured CME parameters are generally significantly higher in the case of the simulated MHD white-light data compared to those from the synthetic white-light images generated by the GCS model. We found the following CME parameter error bars as a starting point for quantifying the minimum error in CME parameters from white-light reconstructions: Δθ (latitude)=6° +2° -3° , Δϕ (longitude)=11° +18° -6° , Δγ (tilt)=25° +8° -7° , Δx (half-angle)=10° +12° -6° , Δh (height)=0.6 +1.2 -0.4 R ⨀ , and Δκ (ratio)=0.1 +0.03 -0.02 .

Coronal mass ejections↗

NASA NTRS2022

Quantifying Errors in 3D CME Parameters Derived From Synthetic Data Using White-Light Reconstruction Techniques

Current efforts in space weather forecasting of CMEs have been focused on predicting their arrival time and magnetic structure. To make these predictions, methods have been developed to derive the true CME speed, size, position, and mass, among others. Difficulties in determining the input parameters for CME forecasting models arise from the lack of direct measurements of the coronal magnetic fields and uncertainties in estimating the CME 3D geometric and kinematic parameters after eruption. White-light coronagraph images are usually employed by a variety of CME reconstruction techniques that assume more or less complex geometries. This is the first study from our International Space Science Institute (ISSI) team “Understanding Our Capabilities in Observing and Modeling Coronal Mass Ejections”, in which we explore how subjectivity affects the 3D CME parameters that are obtained from the Graduated Cylindrical Shell (GCS) reconstruction technique, which is widely used in CME research. To be able to quantify such uncertainties, the “true” values that are being fitted should be known, which are impossible to derive from observational data. We have designed two different synthetic scenarios where the “true” geometric parameters are known in order to quantify such uncertainties for the first time. We explore this by using two sets of synthetic data: 1) Using the ray-tracing option from the GCS model software itself, and 2) Using 3D magnetohydrodynamic (MHD) simulation data from the Magnetohydrodynamic Algorithm outside a Sphere code. Our experiment includes different viewing configurations using single and multiple viewpoints. CME reconstructions using a single viewpoint had the largest errors and error ranges overall for both synthetic GCS and simulated MHD white-light data. As the number of viewpoints increased from one to two, the errors decreased by approximately 4° in latitude, 22° in longitude, 14° in tilt, and 10° in half-angle. Our results quantitatively show the critical need for at least two viewpoints to be able to reduce the uncertainty in deriving CME parameters. We did not find a significant decrease in errors when going from two to three viewpoints for our specific hypothetical three spacecraft scenario using synthetic GCS white-light data. As we expected, considering all configurations and numbers of viewpoints, the mean absolute errors in the measured CME parameters are generally significantly higher in the case of the simulated MHD white-light data compared to those from the synthetic white-light images generated by the GCS model. We found the following CME parameter error bars as a starting point for quantifying the minimum error in CME parameters from white-light reconstructions: Δθ (latitude)=6° +2° -3° , Δϕ (longitude)=11° +18° -6° , Δγ (tilt)=25° +8° -7° , Δx (half-angle)=10° +12° -6° , Δh (height)=0.6 +1.2 -0.4 R ⨀ , and Δκ (ratio)=0.1 +0.03 -0.02 .

Coronal mass ejections↗

DOE OSTI2025

Self-supervised physics-informed generative networks for phase retrieval from a single X-ray hologram

X-ray phase contrast imaging significantly improves the visualization of structures with weak or uniform absorption, broadening its applications across a wide range of scientific disciplines. Propagation-based phase contrast is particularly suitable for time- or dose-critical in vivo/in situ/operando (tomography) experiments because it requires only a single intensity measurement. However, the phase information of the wave field is lost during the measurement and must be recovered. Conventional algebraic and iterative methods often rely on specific approximations or boundary conditions that may not be met by many samples or experimental setups. In addition, they require manual tuning of reconstruction parameters by experts, making them less adaptable for complex or variable conditions. Here we present a self-learning approach for solving the inverse problem of phase retrieval in the near-field regime of Fresnel theory using a single intensity measurement (hologram). A physics-informed generative adversarial network is employed to reconstruct both the phase and absorbance of the unpropagated wave field in the sample plane from a single hologram. Unlike most state-of-the-art deep learning approaches for phase retrieval, our approach does not require paired, unpaired, or simulated training data. This significantly broadens the applicability of our approach, as acquiring or generating suitable training data remains a major challenge due to the wide variability in sample types and experimental configurations. The algorithm demonstrates robust and consistent performance across diverse imaging conditions and sample types, delivering quantitative, high-quality reconstructions for both simulated data and experimental datasets acquired at beamline P05 at PETRA III (DESY, Hamburg), operated by Helmholtz-Zentrum Hereon. Furthermore, it enables the simultaneous retrieval of both phase and absorption information.

36 MATERIALS SCIENCE↗