SEARCH · Engineering Papers
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.
Image data compression schemes are discussed with regard to achieving maximum compression ratio. It is demonstrated that the compression ratio must be traded-off against the allowable residual distortion in the reconstructed image.
Active millimeter and microwave imaging techniques can be used to create a high-resolution 3D image volume of a target’s reflectivity. Millimeter-wave imaging is commonly used for personnel security screening and numerous other applications. Backprojection based image reconstruction techniques form a 3D complex-valued volume. The complex-valued volume is commonly reduced to a real-valued volume by taking the magnitude. For anomaly detection and cross image registration of an object it is beneficial to generate an accurate representation of an object’s surface. Extracting a smooth and accurate surface from a magnitude only image is challenging. One difficulty is the magnitude image has limited resolution in the depth direction which normally limits precision to a moderate fraction of the depth resolution. Another difficulty is that the magnitude of the image depends heavily on the geometry and orientation of the object being imaged. The phase information in a complex-valued image volume provides a means to decouple the magnitude of the image from the geometry of an object and provide precision much finer than the depth resolution would indicate. This enables the generation of a smooth and accurate point cloud representation of the surface of an imaged object. A method to extract a point cloud from the phase information in a 3D complex-valued millimeter-wave image volume is developed and results with simulated and experimental data are presented.
Atmospheric muons are typically high energy, highly penetrating charged particles. They interact with matter primarily through multiple Coulomb scatterings. Muon scattering intensities can be used to characterize the density and atomic number of the matter that they pass through. Previously, the Los Alamos National Laboratory (LANL) muon tomography team performed muon imaging of the partially filled MC-10 spent nuclear fuel (SNF) cask at Idaho National Laboratory (INL). This experiment demonstrated the feasibility of muon imaging for the verification of spent fuel container contents. That original effort used the mini muon tracker array, consisting of two arrays of drift tubes on either side of the SNF cask. The reconstructed image quality was limited by statistics, largely due to low muon flux at high zenith angles. A LANL led team will perform new measurements with a larger array, the Giant Muon Tracker (GMT), to improve data collection rates and statistics. In this work, simulations were performed with the GMT near the partially filled INL MC-10 cask. For more general fuel diversion detection, a full MC-10 cask and casks with a singular missing fuel bundle were also simulated. To understand minimum measurement times needed for missing bundle identification, 100 000 to millions of tracked muons (corresponding to 1.4 days to several weeks measurement time) were analyzed. Simulated images were then analyzed visually and numerically to explore techniques designed to minimize the collection time needed to identify the diversion of fuel in each scenario.
Calculations are given for constructing the three dimensional distribution of optically thin EUV emission sources associated with solar active regions, from two dimensional observations (projections) recorded by the spectroheliograph on the OSO 7 satellite. The relation of the method to other image reconstruction methods is briefly discussed as well as the special requirements imposed in the solar case such as a knowledge of the true solar rotation function. A useful correlation criterion for establishing the physical validity of solutions is given.
Abstract Small animal research is an essential tool in studying both pharmaceutical biodistribution and disease progression over time. Furthermore, through the rapid development of in vivo imaging technology over the last few decades, small animal imaging (also referred to as preclinical imaging) has become a mainstay for all fields of biologic research and a center point for most preclinical cancer research. Preclinical imaging modalities include optical, MRI and MRS, microCT, small animal PET, ultrasound, and photoacoustic, each with their individual strengths. The strong points of small animal PET are its translatability to the clinic; its quantitative imaging capabilities; its whole-body imaging ability to dynamically trace functional/biochemical processes; its ability to provide useful images with only nano- to pico- molar concentrations of administered compounds; and its ability to study animals serially over time. This review paper gives an overview of the development and evolution of small animal PET imaging. It provides an overview of detector designs; system configurations; multimodality PET imaging systems; image reconstruction and analysis tools; and an overview of research and commercially available small animal PET systems. It concludes with a look toward developing technologies/methodologies that will further enhance the impact of small animal PET imaging on medical research in the future.
In the era of big data, materials science workflows need to handle large-scale data distribution, storage, and computation. Any of these areas can become a performance bottleneck. We present a framework for analyzing internal material structures (e.g., cracks) to mitigate these bottlenecks. We demonstrate the effectiveness of our framework for a workflow performing synchrotron X-ray computed tomography reconstruction and segmentation of a silica-based structure. Our framework provides a cloud-based, cutting-edge solution to challenges such as growing intermediate and output data and heavy resource demands during image reconstruction and segmentation. Specifically, our framework efficiently manages data storage, scaling up compute resources on the cloud. The multi-layer software structure of our framework includes three layers. A top layer uses Jupyter notebooks and serves as the user interface. A middle layer uses Ansible for resource deployment and managing the execution environment. A low layer is dedicated to resource management and provides resource management and job scheduling on heterogeneous nodes (i.e., GPU and CPU). At the core of this layer, Kubernetes supports resource management, and Dask enables large-scale job scheduling for heterogeneous resources. The broader impact of our work is four-fold: through our framework, we hide the complexity of the cloud’s software stack to the user who otherwise is required to have expertise in cloud technologies; we manage job scheduling efficiently and in a scalable manner; we enable resource elasticity and workflow orchestration at a large scale; and we facilitate moving the study of nonporous structures, which has wide applications in engineering and scientific fields, to the cloud. While we demonstrate the capability of our framework for a specific materials science application, it can be adapted for other applications and domains because of its modular, multi-layer architecture.
The goal of this project was to develop and demonstrate an innovative, low cost, field deployable, stationary 3D x-ray computed tomography (CT) system that will image total root phenotypes with a micron size resolution at a throughput of hundreds of plants per cycle. This system is based on UHV’s unique low cost linear x-ray tube technology and sophisticated reconstruction & image segmentation algorithms developed at University of Massachusetts, Lowel and University of Nottingham; and was tested for several types of soils at University of Wisconsin and Texas A&M University. Currently, no technologies exist that have been designed to image roots in complex media such as agricultural field conditions. Due to its small size, high resolution & fast imaging of fine roots, low power consumption, large penetration depth (i.e. ability to see through several feet of soil) and ease of field deployability, this CT system will increase the speed and efficacy of discovery, field translation, and deployment of improved crops and systems that improve soil carbon accumulation and storage, decrease N2O emissions, and improve water efficiency leading towards advancements that could mitigate 10% of the total US Greenhouse gases. This degree of imaging in the field has never been available and would be invaluable to scientists in understanding how environmental conditions and phenotypic variations contribute to carbon deposition through root development.
A simple adaptive algorithm is applied to the line by line encoding of images. Among its desirable properties the algorithm includes a control of encoding errors independently of the image statistics. This block transform encoder, because of its adaptive nature gives rise to nonstationary errors and end of block effects. Two methods are proposed to reduce these effects and thus improve the subjective quality of reconstructed images. The performance of the algorithm is image dependent and several illustrative examples are given.
This presentation describes how a vision model incorporating contrast sensitivity, contrast masking, and light adaptation is used to design visually optimal quantization matrices for Discrete Cosine Transform image compression. The Discrete Cosine Transform (DCT) underlies several image compression standards (JPEG, MPEG, H.261). The DCT is applied to 8x8 pixel blocks, and the resulting coefficients are quantized by division and rounding. The 8x8 'quantization matrix' of divisors determines the visual quality of the reconstructed image; the design of this matrix is left to the user. Since each DCT coefficient corresponds to a particular spatial frequency in a particular image region, each quantization error consists of a local increment or decrement in a particular frequency. After adjustments for contrast sensitivity, local light adaptation, and local contrast masking, this coefficient error can be converted to a just-noticeable-difference (jnd). The jnd's for different frequencies and image blocks can be pooled to yield a global perceptual error metric. With this metric, we can compute for each image the quantization matrix that minimizes the bit-rate for a given perceptual error, or perceptual error for a given bit-rate. Implementation of this system demonstrates its advantages over existing techniques. A unique feature of this scheme is that the quantization matrix is optimized for each individual image. This is compatible with the JPEG standard, which requires transmission of the quantization matrix.
We present a novel machine learning CNN architecture that can learn from limited data combined appropriately with physics and statistical priors (e.g., forward models and noise models). To address the limited availability of training data we adopt a 3D patch-based approach for our models. Patch-based learning is central to several image reconstruction methods and demands fewer training data than DL approaches, as a single data volume can be broken into several millions of overlapping 3D sub-volumes or patches. This creates a very large number of training sub-volumes from a limited number of overall image volumes. A 3D Generative Adversarial Networks (GAN) is then trained to remove artifacts at the sub-volume level. The combination of a sub-volume-based approach with DL allows us to exploit the richness of the latter in extracting and representing image features, while avoiding risks associated with overfitting due to limited training data.
X-ray ptychography provides exceptional nanoscale resolution and is widely applied in materials science, biology, and nanotechnology. However, its full potential is constrained by the critical challenge of accurately reconstructing images when the illuminating probe is unknown. Conventional iterative methods and deep learning approaches are often suboptimal, particularly under the low-signal conditions inherent to low-dose and high-speed experiments. These limitations compromise reconstruction fidelity and restrict the broader adoption of the technique. In this work, we introduce the Ptychographic Implicit Neural Representation (PtyINR), a self-supervised framework that simultaneously addresses the object- and probe-recovery problem. By parameterizing both as continuous neural representations, PtyINR performs end-to-end reconstruction directly from raw diffraction patterns without requiring any pre-characterization of the probe. Extensive evaluations demonstrate that PtyINR achieves superior reconstruction quality on both simulated and experimental data, with remarkable robustness under challenging low-signal conditions. Furthermore, PtyINR offers a generalizable, physics-informed framework for addressing probe-dependent inverse problems, making it applicable to a wide range of computational microscopy problems.
Global observation of the dayside auroral region by the Ultraviolet Imager (UVI) during transient solar wind pressure pulse events on October 1, 1997 has revealed unusual features in the auroral precipitation. The auroral arc structure on the dayside, possibly connected with the LLBL, split into 2 arc structures; one moving poleward and fading over a 5 min period, and the other stationary or slightly shifted equatorward (by changes in the x component). The y component was large and positive, and the z component was small and negative. The splitting of the arc structure extended from 9 to 15 MLT and was concurrent with an enhancement of the convection in the cusp region identified by SuperDARN observations. The convection reversal on the morningside was adjacent to and poleward of the weak lower latitude band of precipitation. The sensitivity of the UVI instrument enabled observation of arc structures down to about 0.2 erg electron energy flux, as confirmed by comparison with particle measurements from the FAST satellite for other dayside events. Removal of the spacecraft wobble by PIXON image reconstruction restored the original resolution of the UVI of about 40 km from apogee. This event is being analyzed in connection with a larger study of global dynamics of dayside energy and momentum transfer related to changes in IMF conditions using UVI images in conjunction with observations from FAST and SuperDARN.
Four-dimensional scanning transmission electron microscopy (4D-STEM) is a versatile analytical tool for characterizing materials structural properties. However, extending such analysis to disordered materials is challenging, especially in technologically important samples with mixed ordered and disordered phases. Here, in this work, we present a new 4D-STEM method, called fluctuation cepstral STEM (FC-STEM), based on the fluctuation analysis of cepstral transform of diffraction patterns. The peaks in the associated transformation relate to inter-atomic distances in a thin sample. By varying the real-space range over which fluctuations are calculated, distinct ordered and disordered phases can be mapped in a diffractive image reconstruction. We demonstrate the principles of FC-STEM by characterizing a silicon anode, harvested from a cycled lithium-ion battery. A mixture of amorphous and nanocrystalline silicon, graphitic carbon, and electrolyte by-products is identified and mapped. Comparisons with conventional electron imaging and energy-dispersive X-ray spectroscopy show that FC-STEM is highly effective for the structure determination of mixed-phase amorphous materials.
Increases in X-ray brightness from synchrotron light sources lead to a requirement for higher frame rates from hybrid pixel array detectors (HPADs), while also favoring charge integration over photon counting. However, transfer of the full uncompressed data will begin to constrain detector design, as well as limit the achievable continuous frame rate. Here a data compression scheme that is easy to implement in a HPAD's application-specific integrated circuit (ASIC) is described, and how different degrees of compression affect image quality in ptychography, a commonly employed coherent imaging method, is examined. Using adaptive encoding quantization, it is shown in simulations that one can digitize signals up to 16383 photons per pixel (corresponding to 14 bits of information) using only 8 or 9 bits for data transfer, with negligible effect on the reconstructed image.
Measuring the radiation dose reaching a patient’s body is difficult. Here we report a technique for the tomographic reconstruction of the location of photon pairs originating from the annihilation of positron–electron pairs produced by high-energy X-rays travelling through tissue. We used Monte Carlo simulations on pre-recorded data from tissue-mimicking phantoms and from a patient with a brain tumour to show the feasibility of this imaging modality, which we named ‘pair-production tomography’, for the monitoring of radiotherapy dosing. We simulated three image-reconstruction methods, one applicable to a pencil X-ray beam scanning through a region of interest, and two applicable to the excitation of tissue volumes via broad beams (with temporal resolution sufficient to identify coincident photon pairs via filtered back projection, or with higher temporal resolution sufficient for the estimation of a photon’s time-of-flight). In addition to the monitoring of radiotherapy dosing, we show that image contrast resulting from pair-production tomography is highly proportional to the material’s atomic number. The technique may thus also allow for element mapping and for soft-tissue differentiation.
Wavelength-resolved neutron tomography (WRNT) is an emerging technique for characterizing samples relevant to the materials sciences in 3D. WRNT studies can be carried out at beam lines in spallation neutron or reactor-based user facilities. Because of the limited availability of experimental time, potential imperfections in the neutron source, or constraints placed on the acquisition time by the type of sample, the data can be extremely noisy resulting in tomographic reconstructions with significant artifacts when standard reconstruction algorithms are used. Furthermore, making a full tomographic measurement even with a low signal-to-noise ratio can take several days, resulting in a long wait time before the user can receive feedback from the experiment when traditional acquisition protocols are used. In this paper, we propose an interlaced scanning technique and combine it with a model-based image reconstruction algorithm to produce high-quality WRNT reconstructions concurrent with the measurements being made. The interlaced scan is designed to acquire data so that successive measurements are more diverse in contrast to typical sequential scanning protocols. The model-based reconstruction algorithm combines a data-fidelity term with a regularization term to formulate the wavelength-resolved reconstruction as minimizing a high-dimensional cost-function. Using an experimental dataset of a magnetite sample acquired over a span of about two days, we demonstrate that our technique can produce high-quality reconstructions even during the experiment compared to traditional acquisition and reconstruction techniques. In summary, the combination of the proposed acquisition strategy with an advanced reconstruction algorithm provides a novel guideline for designing WRNT systems at user facilities.
Sampling generally causes the response of a digital imaging system to be locally shift-variant and not directly amenable to Modulation Transfer Function (MTF) analysis. However, this paper demonstrates that a meaningful system response can be calculated by averaging over an ensemble of point-source system inputs to yield an MTF which accounts for the combined effects of image formation, sampling, and image reconstruction. As an illustration, the MTF of the Landsat MSS system is analyzed to reveal an average effective instantaneous field of view which is significantly larger than the commonly accepted value, particularly in the along-track direction where undersampling contributes markedly to an MTF reduction and resultant increase in image blur.
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.