Phase Retrieval with Application to Optical Imaging

This review article provides a contemporary overview of phase retrieval in optical imaging, linking the relevant optical physics to the information processing methods and algorithms. Its purpose is to describe the current state of the art in this area, identify challenges, and suggest vision and areas where signal processing methods can have a large impact on optical imaging and on the world of imaging at large, with applications in a variety of fields ranging from biology and chemistry to physics and engineering

Stability estimates for linearized near-field phase retrieval in X-ray phase contrast imaging

Propagation-based X-ray phase contrast enables nanoscale imaging of biological tissue by probing not only the attenuation, but also the real part of the refractive index of the sample. Since only intensities of diffracted waves can be measured, the main mathematical challenge consists in a phase-retrieval problem in the near-field regime. We treat an often used linearized version of this problem known as contract transfer function model. Surprisingly, this inverse problem turns out to be well-posed assuming only a compact support of the imaged object. Moreover, we establish bounds on the Lipschitz stability constant. In general this constant grows exponentially with the Fresnel number of the imaging setup. However, both for homogeneous objects, characterized by a fixed ratio of the induced refractive phase shifts and attenuation, and in the case of measurements at two distances, a much more favorable algebraic dependence on the Fresnel number can be shown. In some cases we establish order optimality of our estimates.Comment: Published version: Minor revision according to the referees' comments. + Edit 08-09-2017: Corrected typos due to deprecated notation in the proof of Lemma

Expansion-maximization-compression algorithm with spherical harmonics for single particle imaging with X-ray lasers

In 3D single particle imaging with X-ray free-electron lasers, particle orientation is not recorded during measurement but is instead recovered as a necessary step in the reconstruction of a 3D image from the diffraction data. Here we use harmonic analysis on the sphere to cleanly separate the angu- lar and radial degrees of freedom of this problem, providing new opportunities to efficiently use data and computational resources. We develop the Expansion-Maximization-Compression algorithm into a shell-by-shell approach and implement an angular bandwidth limit that can be gradually raised during the reconstruction. We study the minimum number of patterns and minimum rotation sampling required for a desired angular and radial resolution. These extensions provide new av- enues to improve computational efficiency and speed of convergence, which are critically important considering the very large datasets expected from experiment

MULTIBAT: Unified workflow for fast electrochemical 3D simulations of lithium-ion cells combining virtual stochastic microstructures, electrochemical degradation models and model order reduction

We present a simulation workflow for efficient investigations of the interplay between 3D lithium-ion electrode microstructures and electrochemical performance, with emphasis on lithium plating. Our approach addresses several challenges. First, the 3D microstructures of porous electrodes are generated by a parametric stochastic model, in order to significantly reduce the necessity of tomographic imaging. Secondly, we integrate a consistent microscopic, 3D spatially-resolved physical model for the electrochemical behavior of the lithium-ion cells taking lithium plating and stripping into account. This highly non-linear mathematical model is solved numerically on the complex 3D microstructures to compute the transient cell behavior. Due to the complexity of the model and the considerable size of realistic microstructures even a single charging cycle of the battery requires several hours computing time. This renders large scale parameter studies extremely time consuming. Hence, we develop a mathematical model order reduction scheme. We demonstrate how these aspects are integrated into one unified workflow, which is a step towards computer aided engineering for the development of more efficient lithium-ion cells

Wavelet transform and terahertz local tomography

Copyright © 2007 SPIE - The International Society for Optical Engineering. Copyright 2007 Society of Photo-Optical Instrumentation Engineers. This paper was published in Novel Optical Instrumentation for Biomedical Applications III, edited by Christian D. Depeursinge Proc. of SPIE-OSA Biomedical Optics, SPIE Vol. 6631, 663113 and is made available as an electronic reprint with permission of SPIE. One print or electronic copy may be made for personal use only. Systematic or multiple reproduction, distribution to multiple locations via electronic or other means, duplication of any material in this paper for a fee or for commercial purposes, or modification of the content of the paper are prohibited.We use the theory of two dimensional discrete wavelet transforms to derive inversion formulas for the Radon transform of terahertz datasets. These inversion formulas with good localised properties are implemented for the reconstruction of terahertz imaging in the area of interest, with a significant reduction in the required measurements. As a form of optical coherent tomography, terahertz CT complements the current imaging techniques and offers a promising approach for achieving non-invasive inspection of solid materials, with potentially numerous applications in industrial manufacturing and biomedical engineering. © 2007 SPIE-OSA.Xiaoxia Yin and Brian W.-H. Ng and Bradley Fergusona and Derek Abbot