Sparse-view tomographic reconstruction of MagLIF emission volumes using learned basis functions.
Abstract not provided.
Fein, Jeffrey↗
SEARCH · Engineering Papers
Results for “tomographic reconstruction”
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.
Single-View Tomographic Reconstruction Using Learned Primal Dual
Explore the source record for details and available documents.
97 MATHEMATICS AND COMPUTING↗
Seeing the whole picture: Methods for getting the most from micro X-ray computed tomography of TRISO nuclear fuel particles
Tristructural isotropic (TRISO) coated fuel particles are a nuclear fuel form under extensive study for use in advanced nuclear reactor concepts. TRISO fuels are subjected to high temperature neutron irradiations and then examined to assess their performance by determining fission product retention and studying morphological changes. Micro X-ray computed tomography is one method of nondestructively studying the effects of TRISO performance. This work addresses the need for image processing to remove X-ray tomographic reconstruction artifacts that prevent the study of TRISO features, as the TRISO particles’ high Z kernel can introduce metal artifacts that degrade the image quality in the surrounding low Z coating layers. These metal artifacts were reduced by imaging the TRISO particles with both high- and low-energy X-rays and applying a mask to the radiographs obtained with low-energy X-rays to digitally remove the dense fuel kernel region. These masked radiographs were then used to produce a tomographic reconstruction which was combined with the tomographic reconstruction of the high-energy data. This enabled the relatively-low-density TRISO buffer layer to be examined in more detail, providing information on irradiation induced dimensional changes of the coatings. This methodology, which helps see the full picture of a TRISO particle, is not limited to nuclear fuels but can be applied to systems that contain highly attenuating material surrounded by less dense materials.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Accurate real space iterative reconstruction (RESIRE) algorithm for tomography
Tomography has made a revolutionary impact on the physical, biological and medical sciences. The mathematical foundation of tomography is to reconstruct a three-dimensional (3D) object from a set of two-dimensional (2D) projections. As the number of projections that can be measured from a sample is usually limited by the tolerable radiation dose and/or the geometric constraint on the tilt range, a main challenge in tomography is to achieve the best possible 3D reconstruction from a limited number of projections with noise. Over the years, a number of tomographic reconstruction methods have been developed including direct inversion, real-space, and Fourier-based iterative algorithms. Here, we report the development of a real-space iterative reconstruction (RESIRE) algorithm for accurate tomographic reconstruction. RESIRE iterates between the update of a reconstructed 3D object and the measured projections using a forward and back projection step. The forward projection step is implemented by the Fourier slice theorem or the Radon transform, and the back projection step by a linear transformation. Our numerical and experimental results demonstrate that RESIRE performs more accurate 3D reconstructions than other existing tomographic algorithms, when there are a limited number of projections with noise. Furthermore, RESIRE can be used to reconstruct the 3D structure of extended objects as demonstrated by the determination of the 3D atomic structure of an amorphous Ta thin film. We expect that RESIRE can be widely employed in the tomography applications in different fields. Finally, to make the method accessible to the general user community, the MATLAB source code of RESIRE and all the simulated and experimental data are available at https://zenodo.org/record/7273314.
97 MATHEMATICS AND COMPUTING↗
Longitudinal Differences of Ionospheric Vertical Density Distribution and Equatorial Electrodynamics
Accurate estimation of global vertical distribution of ionospheric and plasmaspheric density as a function of local time, season, and magnetic activity is required to improve the operation of space-based navigation and communication systems. The vertical density distribution, especially at low and equatorial latitudes, is governed by the equatorial electrodynamics that produces a vertical driving force. The vertical structure of the equatorial density distribution can be observed by using tomographic reconstruction techniques on ground-based global positioning system (GPS) total electron content (TEC). Similarly, the vertical drift, which is one of the driving mechanisms that govern equatorial electrodynamics and strongly affect the structure and dynamics of the ionosphere in the low/midlatitude region, can be estimated using ground magnetometer observations. We present tomographically reconstructed density distribution and the corresponding vertical drifts at two different longitudes: the East African and west South American sectors. Chains of GPS stations in the east African and west South American longitudinal sectors, covering the equatorial anomaly region of meridian approx. 37 deg and 290 deg E, respectively, are used to reconstruct the vertical density distribution. Similarly, magnetometer sites of African Meridian B-field Education and Research (AMBER) and INTERMAGNET for the east African sector and South American Meridional B-field Array (SAMBA) and Low Latitude Ionospheric Sensor Network (LISN) are used to estimate the vertical drift velocity at two distinct longitudes. The comparison between the reconstructed and Jicamarca Incoherent Scatter Radar (ISR) measured density profiles shows excellent agreement, demonstrating the usefulness of tomographic reconstruction technique in providing the vertical density distribution at different longitudes. Similarly, the comparison between magnetometer estimated vertical drift and other independent drift observation, such as from VEFI onboard Communication/Navigation Outage Forecasting System (C/NOFS) satellite and JULIA radar, is equally promising. The observations at different longitudes suggest that the vertical drift velocities and the vertical density distribution have significant longitudinal differences; especially the equatorial anomaly peaks expand to higher latitudes more in American sector than the African sector, indicating that the vertical drift in the American sector is stronger than the African sector.
total electron content↗
Quantitative phenotyping of crop roots with spectral electrical impedance tomography: a rhizotron study with optimized measurement design
Background: Root systems are key contributors to plant health, resilience, and, ultimately, yield of agricultural crops. To optimize plant performance, phenotyping trials are conducted to breed plants with diverse root traits. However, traditional analysis methods are often labour-intensive and invasive to the root system, therefore limiting high-throughput phenotyping. Spectral electrical impedance tomography (sEIT) could help as a non-invasive and cost-efficient alternative to optical root analysis, potentially providing 2D or 3D spatio-temporal information on root development and activity. Although impedance measurements have been shown to be sensitive to root biomass, nutrient status, and diurnal activity, only few attempts have been made to employ tomographic algorithms to recover spatially resolved information on root systems. In this study, we aim to establish relationships between tomographic electrical polarization signatures and root traits of different fine root systems (maize, pinto bean, black bean, and soy bean) under hydroponic conditions. Results: Our results show that, with the use of an optimized data acquisition scheme, sEIT is capable of providing spatially resolved information on root biomass and root surface area for all investigated root systems. We found strong correlations between the total polarization strength and the root biomass (R 2 = 0.82) and root surface area (R 2 = 0.8). Our findings suggest that the captured polarization signature is dominated by cell-scale polarization processes. Additionally, we demonstrate that the resolution characteristics of the measurement scheme can have a significant impact on the tomographic reconstruction of root traits. Conclusion: Our findings showcase that sEIT is a promising tool for the tomographic reconstruction of root traits in high-throughput root phenotyping trials and should be evaluated as a substitute for traditional, often time-consuming, root characterization methods.
59 BASIC BIOLOGICAL SCIENCES↗
Physics-assisted generative adversarial network for X-ray tomography
X-ray tomography is capable of imaging the interior of objects in three dimensions non-invasively, with applications in biomedical imaging, materials science, electronic inspection, and other fields. The reconstruction process can be an ill-conditioned inverse problem, requiring regularization to obtain satisfactory results. Recently, deep learning has been adopted for tomographic reconstruction. Unlike iterative algorithms which require a distribution that is known a priori , deep reconstruction networks can learn a prior distribution through sampling the training distributions. In this work, we develop a Physics-assisted Generative Adversarial Network (PGAN), a two-step algorithm for tomographic reconstruction. In contrast to previous efforts, our PGAN utilizes maximum-likelihood estimates derived from the measurements to regularize the reconstruction with both known physics and the learned prior. Compared with methods with less physics assisting in training, PGAN can reduce the photon requirement with limited projection angles to achieve a given error rate. The advantages of using a physics-assisted learned prior in X-ray tomography may further enable low-photon nanoscale imaging.
47 OTHER INSTRUMENTATION↗
X-ray scattering based scanning tomography for imaging and structural characterization of cellulose in plants
X-ray and neutron scattering have long been used for structural characterization of cellulose in plants. Due to averaging over the illuminated sample volume, these measurements traditionally overlooked the compositional and morphological heterogeneity within the sample. Here, a scanning tomographic imaging method is described, using contrast derived from the X-ray scattering intensity, for virtually sectioning the sample to reveal its internal structure at a resolution of a few micrometres. This method provides a means for retrieving the local scattering signal that corresponds to any voxel within the virtual section, enabling characterization of the local structure using traditional data-analysis methods. This is accomplished through tomographic reconstruction of the spatial distribution of a handful of mathematical components identified by non-negative matrix factorization from the large dataset of X-ray scattering intensity. Joint analysis of multiple datasets, to find similarity between voxels by clustering of the decomposed data, could help elucidate systematic differences between samples, such as those expected from genetic modifications, chemical treatments or fungal decay. The spatial distribution of the microfibril angle can also be analyzed, based on the tomographically reconstructed scattering intensity as a function of the azimuthal angle.
36 MATERIALS SCIENCE↗
High Fidelity Tomographic Density Reconstructions
Explore the source record for details and available documents.
Klasky, Marc Louis↗
X-ray nano-holotomography reconstruction with simultaneous probe retrieval
In conventional tomographic reconstruction, the pre-processing step includes flat-field correction, where each sample projection on the detector is divided by a reference image taken without the sample. When using coherent X-rays as a probe, this approach overlooks the phase component of the illumination field (probe), leading to artifacts in phase-retrieved projection images, which are then propagated to the reconstructed 3D sample representation. The problem intensifies in nano-holotomography with focusing optics, which, due to various imperfections creates high-frequency components in the probe function. Here, we present a new iterative reconstruction scheme for holotomography, simultaneously retrieving the complex-valued probe function. Implemented on GPUs, this algorithm results in 3D reconstruction resolving twice thinner layers in a 3D ALD standard sample measured using nano-holotomography.
Nikitin, Viktor↗
3D imaging of magnetic domains in Nd 2 Fe 14 B using scanning hard X-ray nanotomography
Nanoscale structural and electronic heterogeneities are prevalent in condensed matter physics. Investigating these heterogeneities in 3D has become an important task for understanding material properties. To provide a tool to unravel the connection between nanoscale heterogeneity and macroscopic emergent properties in magnetic materials, scanning transmission X-ray microscopy (STXM) is combined with X-ray magnetic circular dichroism. A vector tomography algorithm has been developed to reconstruct the full 3D magnetic vector field without any prior noise assumptions or knowledge about the sample. Two tomographic scans around the vertical axis are acquired on single-crystalline Nd 2 Fe 14 B pillars tilted at two different angles, with 2D STXM projections recorded using a focused 120 nm X-ray beam with left and right circular polarization. Image alignment and iterative registration have been implemented based on the 2D STXM projections for the two tilts. Dichroic projections obtained from difference images are used for the tomographic reconstruction to obtain the 3D magnetization distribution at the nanoscale.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Calibration method for a spectral computerized tomography system
A calibration method for an x-ray computerized tomography system and a method of tomographic reconstruction are provided. The calibration method includes steps of measuring at least one point spread function (PSF) at each of a plurality of points, compressing each PSF, and in one or more storing operations, storing the compressed PSFs in a computer-accessible storage medium. The PSF measurements are made in a grid of calibration points in a field of view (FOV) of the system. In the measuring step, an absorber is positioned at each of the calibration points, and an x-ray projection is taken at least once at each of those absorber positions. In the method of tomographic image reconstruction, projection data from an x-ray tomographic projection system are input to an iterative image reconstruction algorithm. The algorithm retrieves and utilizes a priori system information (APSI) The APSI comprises comprising point spread functions (PSFs) of all voxels in a voxelization of the field of view that are compressed in the form of vectors of parameters. For utilization, each retrieved vector of parameters is decompressed so as to generate a discretized PSF.
Jimenez, Jr., Edward Steven↗
Experimental Investigation of Pool Boiling Heat Transfer Enhancement in Microgravity in the Presence of Electric Fields
In boiling high heat fluxes are possible driven by relatively small temperature differences, which make its use increasingly attractive in aerospace applications. The objective of the research is to develop ways to overcome specific problems associated with boiling in the low gravity environment by substituting the buoyancy force with the electric force to enhance bubble removal from the heated surface. Previous studies indicate that in terrestrial applications nucleate boiling heat transfer can be increased by a factor of 50, as compared to values obtained for the same system without electric fields. The goal of our research is to experimentally explore the mechanisms responsible for EHD heat transfer enhancement in boiling in low gravity conditions, by visualizing the temperature distributions in the vicinity of the heated surface and around the bubble during boiling using real-time holographic interferometry (HI) combined with high-speed cinematography. In the first phase of the project the influence of the electric field on a single bubble is investigated. Pool boiling is simulated by injecting a single bubble through a nozzle into the subcooled liquid or into the thermal boundary layer developed along the flat heater surface. Since the exact location of bubble formation is known, the optical equipment can be aligned and focused accurately, which is an essential requirement for precision measurements of bubble shape, size and deformation, as well as the visualization of temperature fields by HI. The size of the bubble and the frequency of bubble departure can be controlled by suitable selection of nozzle diameter and mass flow rate of vapor. In this approach effects due to the presence of the electric field can be separated from effects caused by the temperature gradients in the thermal boundary layer. The influence of the thermal boundary layer can be investigated after activating the heater at a later stage of the research. For the visualization experiments a test cell was developed. All four vertical walls of the test cell are transparent, and they allow transillumination with laser light for visualization experiments by HI. The bottom electrode is a copper cylinder, which is electrically grounded. The copper block is heated with a resistive heater and it is equipped with 6 thermocouples that provide reference temperatures for the measurements with HI. The top electrode is a mesh electrode. Bubbles are injected with a syringe into the test cell through the bottom electrode. The working fluids presently used in the interferometric visualization experiments, water and PF 5052, satisfy requirements regarding thermophysical, optical and electrical properties. A 30kV power supply equipped with a voltmeter allows to apply the electric field to the electrodes during the experiments. The magnitude of the applied voltage can be adjusted either manually or through the LabVIEW data acquisition and control system connected to a PC. Temperatures of the heated block are recorded using type-T thermocouples, whose output is read by a data acquisition system. Images of the bubbles are recorded with 35mm photographic and 16mm high-speed cameras, scanned and analyzed using various software packages. Visualized temperature fields HI allows the visualization of temperature fields in the vicinity of bubbles during boiling in the form of fringes. Typical visualized temperature distributions around the air bubbles injected into the thermal boundary layer in PF5052 are shown. The temperature of the heated surface is 35 C. The temperature difference for a pair of fringes is approximately 0.05 C. The heat flux applied to the bottom surface is moderate, and the fringe patterns are regular. In the image a bubble penetrating the thermal boundary layer is visible. Because of the axial symmetry of the problem, simplified reconstruction techniques can be applied to recover the temperature field. The thermal plume developing above the heated surface for more intensive heating is shown. The temperature distribution in the liquid is clearly 3D, and tomographic techniques have to be applied to recover the temperature distribution in such a physical situation. A sequence of interferometric images showing the temperature distribution around the rising bubble, recorded with a high-speed camera is shown. Again, the temperature distribution is 3D, and a more complex approach to the evaluation, the tomographic reconstruction has to be taken. Measurement of the temperature distribution from the fringe pattern temperature distributions that yield important information regarding heat transfer are determined. Two algorithms that allow the quantitative evaluation of interferometric fringe patterns and the reconstruction of temperature fields during boiling have been developed at the Heat Transfer Laboratory of the Johns Hopkins University. In the first algorithm the bubble is assumed to be axially symmetrical, which significantly reduces the computational effort for quantifying the temperature distribution around the bubble. For this purpose the thermal boundary layer around the bubble is divided into equidistant concentric shells, and the refractive index is assumed to be constant in each of the shells. Since large temperature gradients are expected in the vicinity of the bubble during boiling, the deflection of the light beam cannot be neglected in boiling experiments. Since the exit angle of the light beam is known, this allows to account for the deflections and phase shifts outside the boundary layer (in the bulk fluid and in the windows of the test cell). Three dimensional temperature distributions in the vicinity of the bubble are reconstructed using tomographic techniques. In tomography, the measurement volume is sliced into 2D planes. In the present study these planes are parallel to the heated surface. The objective is to determine the values of the field parameter of interest in form of the field function in these 2D planes. The field parameter is the change of the refractive index of the liquid in the measurement volume caused by temperature changes. By superimposing data for many 2D planes recorded at the same time instant, the 3D temperature distribution in the measurement volume is recovered.
Herman, Cila↗
Joint iterative reconstruction and 3D rigid alignment for X-ray tomography
X-ray tomography is widely used for three-dimensional structure determination in many areas of science, from the millimeter to the nanometer scale. The resolution and quality of the 3D reconstruction is limited by the availability of alignment parameters that correct for the mechanical shifts of the sample or sample stage for the images that constitute a scan. In this paper we describe an algorithm for marker-free, fully automated and accurately aligned and reconstructed X-ray tomography data. Our approach solves the tomographic reconstruction jointly with projection data alignment based on a rigid-body deformation model. We demonstrate the robustness of our method on both synthetic phantom and experimental data and show that our method is highly efficient in recovering relatively large alignment errors without prior knowledge of a low resolution approximation of the 3D structure or a reasonable estimate of alignment parameters.
36 MATERIALS SCIENCE↗
FIRM: federated image reconstruction using multimodal tomographic data
Here, we propose a federated algorithm for reconstructing images using multimodal tomographic data sourced from dispersed locations, addressing the challenges of traditional unimodal approaches that are prone to noise and reduced image quality, as well as the limitations of centralized multimodal approaches that require extensive data transfer, leading to significant communication overhead, storage demands, and potential data privacy concerns. Our approach formulates a joint inverse optimization problem incorporating multimodality constraints and solves it in a federated framework through local gradient computations complemented by lightweight central operations, thereby ensuring data decentralization. Leveraging the connection between our federated algorithm and the quadratic penalty method, we introduce an adaptive step-size rule with guaranteed sublinear convergence. Numerical results demonstrate superior computational efficiency and improved image reconstruction quality compared to existing approaches.
federated algorithm↗
A Feasibility Study of 4-D Tomography of Soft X-ray Magnetosheath Emissivities Using Multi-Spacecraft Measurements
Upcoming heliophysics missions utilize state-of-the-art wide field-of-view (FOV) imaging technology to measure and investigate the space plasma environment on a global scale. At Earth, remote sensing of soft X-ray emissions, which are generated via the charge exchange interaction between heavy solar wind ions and exospheric neutral atoms, is a promising means to investigate the global magnetosheath structure, its response to varying solar wind conditions, and the spatiotemporal properties of the dayside magnetic reconnection. Data analysis techniques such as optical tomography can provide additional structural and time-varying information from the observed target and thus enhance the mission’s scientific return. In this work, we simulate multiple and simultaneous observations of the dayside magnetosphere using soft X-ray imagers located at long-distance vantage points to reconstruct the time-dependent, three-dimensional (3-D) structure of the magnetosheath using a dynamic tomographic approach. The OpenGCCM MHD model is used to simulate the time-varying response of the magnetosheath to solar wind conditions and, subsequently, generate synthetic soft X-ray images from multiple spacecraft vantage points separated along a common orbit. A detailed analysis is then performed to identify the nominal set of spacecraft that produces the highest fidelity tomographic reconstruction of the magnetopause. This work aims to (i) demonstrate, for the first time, the use of dynamic tomography to retrieve the time-varying magnetosheath structure and (ii) identify a nominal mission design for multi-spacecraft configurations aiming for optical tomography.
3-D/4_d tomography↗
3D reconstruction of an inertial-confinement fusion implosion with neural networks using multiple heterogeneous data sources
3D asymmetries are major degradation mechanisms in inertial-confinement fusion implosions at the National Ignition Facility (NIF). These asymmetries can be diagnosed and reconstructed with the neutron imaging system (NIS) on three lines of sight around the NIF target chamber. Conventional tomographic reconstructions are used to reconstruct the 3D morphology of the implosion using NIS [Volegov et al., J. Appl. Phys. 127, 083301 (2020)], but the problem is ill-posed with only three imaging lines of sight. Asymmetries can also be diagnosed with the real-time neutron activation diagnostics (RTNAD) and the neutron time-of-flight (nToF) suite. Since the NIS, RTNAD, and nToF each sample a different part of the implosion using different physical principles, we propose that it is possible to overcome the limitations of too few imaging lines of sight by performing 3D reconstructions that combine information from all three heterogeneous data sources. This work presents a new machine learning-based reconstruction technique to do just this. By using a simple physics model and group of neural networks to map 3D morphologies to data, this technique can easily account for data of multiple different types. A simple proof-of-principle is presented, demonstrating that this technique can accurately reconstruct a hot-spot shape using synthetic primary neutron images and a hot-spot velocity vector. In particular, the hot-spot’s asymmetry, quantified as spherical harmonic coefficients, is reconstructed to within ±4% of the radius in 90% of test cases. In the future, this technique will be applied to actual NIS, RTNAD, and nToF data to better understand 3D asymmetries at the NIF.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Evaluating multi-slice ptychography tomography for X-ray imaging
X-ray ptychography typically relies on single-slice reconstruction approaches that assume minimal beam propagation through the specimen. However, this approximation breaks down for thick samples or when using high numerical aperture optics, leading to reconstruction artifacts due to multiple scattering. Multi-slice ptychographic tomography (MSPT) addresses this limitation by explicitly modeling beam propagation through the specimen while reducing angular sampling requirements for tomography, but its performance depends critically on reconstruction parameters and experimental conditions. Here, we systematically evaluate MSPT using numerical simulations and experimental X-ray datasets to understand how factors such as object thickness, optical parameters, and initialization strategies affect reconstruction accuracy. We demonstrate that informed object initialization significantly improves depth resolution and enables accurate projection extension, reducing the angular sampling requirements for tomographic reconstruction. Experimental demonstration confirms that MSPT produces reconstructions with reduced artifacts, establishing it as a robust approach for nanoscale imaging at modern synchrotron facilities.
Luktuke, Amey [Argonne National Laboratory (ANL), ↗