Search NASA⌕ Search

SEARCH · Search NASA

Results for “iterative phase 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.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Penalized maximum likelihood reconstruction for x-ray differential phase-contrast tomography

Purpose: The purpose of this work is to propose a cost function with regularization to iteratively reconstruct attenuation, phase, and scatter images simultaneously from differential phase contrast (DPC) acquisitions, without the need of phase retrieval, and examine its properties. Furthermore this reconstruction method is applied to an acquisition pattern that is suitable for a DPC tomographic system with continuously rotating gantry (sliding window acquisition), overcoming the severe smearing in noniterative reconstruction. Methods: We derive a penalized maximum likelihood reconstruction algorithm to directly reconstruct attenuation, phase, and scatter image from the measured detector values of a DPC acquisition. The proposed penalty comprises, for each of the three images, an independent smoothing prior. Image quality of the proposed reconstruction is compared to images generated with FBP and iterative reconstruction after phase retrieval. Furthermore, the influence between the priors is analyzed. Finally, the proposed reconstruction algorithm is applied to experimental sliding window data acquired at a synchrotron and results are compared to reconstructions based on phase retrieval. Results: The results show that the proposed algorithm significantly increases image quality in comparison to reconstructions based on phase retrieval. No significant mutual influence between the proposed independent priors could be observed. Further it could be illustrated that the iterative reconstruction of a sliding window acquisition results in images with substantially reduced smearing artifacts. Conclusions: Although the proposed cost function is inherently nonconvex, it can be used to reconstruct images with less aliasing artifacts and less streak artifacts than reconstruction methods based on phase retrieval. Furthermore, the proposed method can be used to reconstruct images of sliding window acquisitions with negligible smearing artifacts.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2018

Dynamic diffraction artefacts in Bragg coherent diffractive imaging

This article reports a theoretical study on the reconstruction artefacts in Bragg coherent diffractive imaging caused by dynamical diffraction effects. It is shown that, unlike the absorption and refraction effects that can be corrected after reconstruction, dynamical diffraction effects have profound impacts on both the amplitude and the phase of the reconstructed complex object, causing strong artefacts. At the dynamical diffraction limit, the reconstructed shape is no longer correct, as a result of the strong extinction effect. Simulations for hemispherical particles of different sizes show the type, magnitude and extent of the dynamical diffraction artefacts, as well as the conditions under which they are negligible.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Atomically resolved structural determination of graphene and its point defects via extrapolation assisted phase retrieval

Previously reported crystalline structures obtained by an iterative phase retrieval reconstruction of their diffraction patterns seem to be free from displaying any irregularities or defects in the lattice, which appears to be unrealistic. We demonstrate here that the structure of a nanocrystal including its atomic defects can unambiguously be recovered from its diffraction pattern alone by applying a direct phase retrieval procedure not relying on prior information of the object shape. Individual point defects in the atomic lattice are clearly apparent. Conventional phase retrieval routines assume isotropic scattering. We show that when dealing with electrons, the quantitatively correct transmission function of the sample cannot be retrieved due to anisotropic, strong forward scattering specific to electrons. We summarize the conditions for this phase retrieval method and show that the diffraction pattern can be extrapolated beyond the original record to even reveal formerly not visible Bragg peaks. Such extrapolated wave field pattern leads to enhanced spatial resolution in the reconstruction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Projected Multi-Agent Consensus Equilibrium (PMACE) With Application to Ptychography

Multi-Agent Consensus Equilibrium (MACE) formulates an inverse imaging problem as a balance among multiple update agents such as data-fitting terms and denoisers. However, each such agent operates on a separate copy of the full image, leading to redundant memory use and slow convergence when each agent affects only a small subset of the full image. In this article, we extend MACE to Projected Multi-Agent Consensus Equilibrium (PMACE), in which each agent updates only a projected component of the full image, thus greatly reducing memory use for some applications. We describe PMACE in terms of an equilibrium problem and an equivalent fixed-point problem and show that in most cases the PMACE equilibrium is not the solution of an optimization problem. To demonstrate the value of PMACE, we apply it to the problem of ptychography, in which a sample is reconstructed from the diffraction patterns resulting from coherent X-ray illumination at multiple overlapping spots. In our PMACE formulation, each spot corresponds to a separate data-fitting agent, with the final solution found as an equilibrium among all the agents. In conclusion, our results demonstrate that the PMACE reconstruction algorithm generates more accurate reconstructions at a lower computational cost than existing ptychography algorithms when the spots are sparsely sampled.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Defect identification in simulated Bragg coherent diffraction imaging by automated AI

X-ray Bragg coherent diffraction imaging is a powerful technique for operando and in situ materials characterization and provides a unique means of quantifying the influence of one-dimensional (1D) and two-dimensional (2D) material defects on material response. However, obtaining full images from raw x-ray diffraction data is nontrivial and computationally intensive, precluding real-time experimental feedback. Here, we present a machine learning approach to identify the presence of crystalline line defects (edge and screw) in samples from the raw, 2D, coherent diffraction data without the need for image reconstruction through iterative phase retrieval. Further, we compare different approaches to designing neural networks for this application and demonstrate the potential of automated ML (autoML) approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Reducing radiation dose for multi-phase contrast-enhanced dual energy renal CT: pilot study evaluating prior iterative reconstruction

Purpose: Prior iterative reconstruction (PIR) uses spatial information from one phase of enhancement to reduce image noise in other phases. We sought to determine if PIR could reduce radiation dose while preserving observer performance and CT number at multi-phase dual energy (DE) renal CT. Methods: CT projection data from multi-phase DE renal CT examinations were collected. Images corresponding to 40% radiation dose were reconstructed using validated noise insertion and PIR. Three genitourinary radiologists examined routine and 40% dose PIR images. Probability of malignancy was assessed [from 0 to 100] with malignancy assumed at probability ≥ 75. Observer performance was compared on a per patient and per lesion level. CT number accuracy was measured. Results: Twenty-three patients had 49 renal lesions (11 solid renal neoplasms). CT number was nearly identical between techniques (mean CT number difference: unenhanced 2 ± 2 HU; enhanced 4 ± 4 HU). AUC for malignancy was similar between multi-phase routine dose DE and lower dose PIR images [per patient: 0.950 vs. 0.916 (p = 0.356); per lesion: 0.931 vs. 0.884 (p = 0.304)]. Per patient sensitivity was also similar (78% routine dose vs. 82% lower dose [p ≥ 0.99]), as was specificity (91% routine dose vs. 93% lower dose PIR [p > 0.99]), with similar findings on a per lesion level. Subjective image quality was also similar (p = 0.34). Conclusions: Prior iterative reconstruction is a new reconstruction method for multi-phase CT examinations that promises to facilitate radiation dose reduction by over 50% for multi-phase DE renal CT exams without compromising CT number or observer performance.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Direct single-shot phase retrieval from the diffraction pattern of separated objects

The non-crystallographic phase problem arises in numerous scientific and technological fields. An important application is coherent diffractive imaging. Recent advances in X-ray free-electron lasers allow capturing of the diffraction pattern from a single nanoparticle before it disintegrates, in so-called ‘diffraction before destruction’ experiments. Presently, the phase is reconstructed by iterative algorithms, imposing a non-convex computational challenge, or by Fourier holography, requiring a well-characterized reference field. Here we present a convex scheme for single-shot phase retrieval for two (or more) sufficiently separated objects, demonstrated in two dimensions. In our approach, the objects serve as unknown references to one another, reducing the phase problem to a solvable set of linear equations. We establish our method numerically and experimentally in the optical domain and demonstrate a proof-of-principle single-shot coherent diffractive imaging using X-ray free-electron lasers pulses. Lastly, our scheme alleviates several limitations of current methods, offering a new pathway towards direct reconstruction of complex objects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Denoised Wigner distribution deconvolution via low-rank matrix completion

Wigner distribution deconvolution (WDD) is a decades-old method for recovering phase from intensity measurements. Although the technique offers an elegant linear solution to the quadratic phase retrieval problem, it has seen limited adoption due to its high computational/memory requirements and the fact that the technique often exhibits high noise sensitivity. Here, we propose a method for noise suppression in WDD via low-rank noisy matrix completion. Our technique exploits the redundancy of an object’s phase space to denoise its WDD reconstruction. We show in model calculations that our technique outperforms other WDD algorithms as well as modern iterative methods for phase retrieval such as ptychography. Here, our results suggest that a class of phase retrieval techniques relying on regularized direct inversion of ptychographic datasets (instead of iterative reconstruction techniques) can provide accurate quantitative phase information in the presence of high levels of noise.

47 OTHER INSTRUMENTATION↗
DOE OSTI2014

SU-D-17A-03: 5D Respiratory Motion Model Based Iterative Reconstruction Method for 4D Cone-Beam CT

Purpose: The purpose of this work is to develop a new iterative reconstruction method for 4D cone-beam CT (CBCT) based on a published time-independent 5D respiratory motion model. The proposed method will offer a single high-resolution image at a user-selected breathing phase and the 5D motion model parameters, which could be used to generate the breathing pattern during the CT acquisition. Methods: 5D respiratory motion model was proposed for accurately modeling the motion of lung and lung tumor tissues. 4D images are then parameterized by a reference image, measured breathing amplitude, breathing rate, two time-independent vector fields that describe the 5D model parameters, and a scalar field that describes the change in HU as a function of breathing amplitude. In contrast with the traditional method of reconstructing multiple temporal image phases to reduce respiratory artifact, 5D model based method simplify the problem into the reconstruction of a single reference image and the 5D motion model parameters. The reconstruction formulation of the reference image and scalar and vector fields is a nonlinear least-square optimization problem that consists of solving the reference image and fields alternately, in which the reference image is regularized with the total variation sparsity transform and the vector fields are solved through linearizations regularized by the H1 norm. 2D lung simulations were performed in this proof-of-concept study. Results: The breathing amplitude, its rate, and the corresponding scalar and vector fields were generated from a patient case. Compared with filtered backprojection method and sparsity regularized iterative method for the phase-by-phase reconstruction, the proposed 5D motion model based method yielded improved image quality. Conclusion: Based on 5D respiratory motion model, we have developed a new iterative reconstruction method for 4D CBCT that has the potential for improving image quality while providing needed on-table motion information.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

3D reconstruction of the magnetic vector potential using model based iterative reconstruction

Lorentz transmission electron microscopy (TEM) observations of magnetic nanoparticles contain information on the magnetic and electrostatic potentials. Vector field electron tomography (VFET) can be used to reconstruct electromagnetic potentials of the nanoparticles from their corresponding LTEM images. The VFET approach is based on the conventional filtered back projection approach to tomographic reconstructions and the availability of an incomplete set of measurements due to experimental limitations means that the reconstructed vector fields exhibit significant artifacts. In this paper, we outline a model-based iterative reconstruction (MBIR) algorithm to reconstruct the magnetic vector potential of magnetic nanoparticles. We combine a forward model for image formation in TEM experiments with a prior model to formulate the tomographic problem as a maximum a-posteriori probability estimation problem (MAP). The MAP cost function is minimized iteratively to determine the vector potential. Here, a comparative reconstruction study of simulated as well as experimental data sets show that the MBIR approach yields quantifiably better reconstructions than the VFET approach.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Communication: X-ray coherent diffractive imaging by immersion in nanodroplets

Lensless x-ray microscopy requires the recovery of the phase of the radiation scattered from a specimen. Here, we demonstrate a de novo phase retrieval technique by encapsulating an object in a superfluid helium nanodroplet, which provides both a physical support and an approximate scattering phase for the iterative image reconstruction. The technique is robust, fast-converging, and yields the complex density of the immersed object. As a result, images of xenon clusters embedded in superfluid helium droplets reveal transient configurations of quantum vortices in this fragile system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

SU-D-12A-05: Iterative Reconstruction Techniques to Enable Intrinsic Respiratory Gated CT in Mice

Purpose: Longitudinal studies of lung function in mice need the ability to image different phases of ventilation in free-breathing mice using retrospective gating. However, retrospective gating often produces under-sampled and uneven angular samples, resulting in severe reconstruction artifacts when using traditional FDK based reconstruction algorithms. We wanted to demonstrate the utility of iterative reconstruction method to enable intrinsic respiratory gating in small-animal CT. Methods: Free-breathing mice were imaged using a Siemens Inveon PET/micro-CT system. Evenly distributed projection images were acquired at 360 angles. Retrospective respiratory gating was performed using an intrinsic marker based on the average intensity in a region covering the diaphragm. Projections were classified into 4 and 6 phases (finer temporal resolution) resulting in 138 and 67 projections respectively. Reconstruction was carried out using 3 Methods: conventional FDK, iterative penalized least-square (PWLS) with total variation (TV), and PWLS with edge-preserving penalty. The performance of the methods was compared using contrast-to-noise (CNR) in a region of interest (ROI). Line profile through a specific region was plotted to evaluate the preserving of edges. Results: In both the cases with 4 and 6 phases, inadequate and non-uniform angular sampling results in artifacts using conventional FDK. However, such artifacts are minimized using both the iterative methods. Using both 4 and 6 phases, the iterative techniques outperformed FDK in terms of CNR and maintaining sharp edges. This is further evidenced especially with increased artifacts using FDK for 6 phases. Conclusion: This work indicates fewer artifacts and better image details can be achieved with iterative reconstruction methods in non-uniform under-sampled reconstruction. Using iterative methods can enable free-breathing intrinsic respiratory gating in small-animal CT. Further studies are needed to compare the computational complexity for large datasets.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

PtychoNet: Fast and High Quality Phase Retrieval for Ptychography

Ptychography is a coherent diffractive imaging method that captures multiple diffraction patterns of a sample with a set of shifted localized illuminations (“probes”). The reconstruction problem, known as “phase retrieval”, is typically solved by iterative algorithms. In this paper, we propose PtychoNet, a deep learning based method to perform phase retrieval for ptychography in a non-iterative manner. We devise a generative network to encode a full ptychography scan, reverse the diffractions at each scanning point and compute the amplitude and phase of the object. We demonstrate successful reconstructions using PtychoNet as well as recovering fine features in the case of extreme sparse scanning where conventional iterative methods fail to give recognizable features.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Randomized probe imaging through deep k-learning

Randomized probe imaging (RPI) is a single-frame diffractive imaging method that uses highly randomized light to reconstruct the spatial features of a scattering object. The reconstruction process, known as phase retrieval, aims to recover a unique solution for the object without measuring the far-field phase information. Typically, reconstruction is done via time-consuming iterative algorithms. In this work, we propose a fast and efficient deep learning based method to reconstruct phase objects from RPI data. The method, which we call deep k-learning, applies the physical propagation operator to generate an approximation of the object as an input to the neural network. This way, the network no longer needs to parametrize the far-field diffraction physics, dramatically improving the results. Deep k-learning is shown to be computationally efficient and robust to Poisson noise. The advantages provided by our method may enable the analysis of far larger datasets in photon starved conditions, with important applications to the study of dynamic phenomena in physical science and biological engineering.

Guo, Zhen (ORCID:0000000213473451)↗
DOE OSTI2016

MO-DE-207A-08: Four-Dimensional Cone-Beam CT Iterative Reconstruction with Time-Ordered Chain Graph Model for Non-Periodic Organ Motion and Deformation

Purpose: The purpose of this study is to propose a new concept of four-dimensional (4D) cone-beam CT (CBCT) reconstruction for non-periodic organ motion using the Time-ordered Chain Graph Model (TCGM), and to compare the reconstructed results with the previously proposed methods, the total variation-based compressed sensing (TVCS) and prior-image constrained compressed sensing (PICCS). Methods: CBCT reconstruction method introduced in this study consisted of maximum a posteriori (MAP) iterative reconstruction combined with a regularization term derived from a concept of TCGM, which includes a constraint coming from the images of neighbouring time-phases. The time-ordered image series were concurrently reconstructed in the MAP iterative reconstruction framework. Angular range of projections for each time-phase was 90 degrees for TCGM and PICCS, and 200 degrees for TVCS. Two kinds of projection data, an elliptic-cylindrical digital phantom data and two clinical patients’ data, were used for reconstruction. The digital phantom contained an air sphere moving 3 cm along longitudinal axis, and temporal resolution of each method was evaluated by measuring the penumbral width of reconstructed moving air sphere. The clinical feasibility of non-periodic time-ordered 4D CBCT reconstruction was also examined using projection data of prostate cancer patients. Results: The results of reconstructed digital phantom shows that the penumbral widths of TCGM yielded the narrowest result; PICCS and TCGM were 10.6% and 17.4% narrower than that of TVCS, respectively. This suggests that the TCGM has the better temporal resolution than the others. Patients’ CBCT projection data were also reconstructed and all three reconstructed results showed motion of rectal gas and stool. The result of TCGM provided visually clearer and less blurring images. Conclusion: The present study demonstrates that the new concept for 4D CBCT reconstruction, TCGM, combined with MAP iterative reconstruction framework enables time-ordered image reconstruction with narrower time-window.

60 APPLIED LIFE SCIENCES↗