Iterative algorithms for phase retrieval from intensity data are compared to gradient search methods. Both the problem of phase retrieval from two intensity measurements (in electron microscopy or wave front sensing) and the problem of phase retrieval from a single intensity measurement plus a non-negativity constraint (in astronomy) are considered, with emphasis on the latter. It is shown that both the error-reduction algorithm for the problem of a single intensity measurement and the Gerchberg-Saxton algorithm for the problem of two intensity measurements converge. The error-reduction algorithm is also shown to be closely related to the steepest-descent method. Other algorithms, including the input-output algorithm and the conjugate-gradient method, are shown to converge in practice much faster than the error-reduction algorithm. Examples are shown.

Phase retrieval algorithms: a comparison.

1. Basic Concepts. 2. Nonparametric Methods. 3. Parametric Methods for Rational Spectra. 4. Parametric Methods for Line Spectra. 5. Filter Bank Methods. 6. Spatial Methods. Appendix A: Linear Algebra and Matrix Analysis Tools. Appendix B: Cramer-Rao Bound Tools. Appendix C: Model Order Selection Tools. Appendix D: Answers to Selected Exercises. Bibliography. References Grouped by Subject. Subject Index.

/pdf/spectral-analysis-of-signals-4dbupf9m4n.pdf

Spectral analysis of signals

Extending the methodology of X-ray crystallography to allow imaging of micrometre-sized non-crystalline specimens

Foundering of mafic lower continental crust into underlying convecting mantle has been proposed as one means to explain the unusually evolved chemical composition of Earth's continental crust, yet direct evidence of this process has been scarce. Here we report that Late Jurassic high-magnesium andesites, dacites and adakites (siliceous lavas with high strontium and low heavy-rare-earth element and yttrium contents) from the North China craton have chemical and petrographic features consistent with their origin as partial melts of eclogite that subsequently interacted with mantle peridotite. Similar features observed in adakites and some Archaean sodium-rich granitoids of the tonalite-trondhjemite-granodiorite series have been interpreted to result from interaction of slab melts with the mantle wedge. Unlike their arc-related counterparts, however, the Chinese magmas carry inherited Archaean zircons and have neodymium and strontium isotopic compositions overlapping those of eclogite xenoliths derived from the lower crust of the North China craton. Such features cannot be produced by crustal assimilation of slab melts, given the high Mg#, nickel and chromium contents of the lavas. We infer that the Chinese lavas derive from ancient mafic lower crust that foundered into the convecting mantle and subsequently melted and interacted with peridotite. We suggest that lower crustal foundering occurred within the North China craton during the Late Jurassic, and thus provides constraints on the timing of lithosphere removal beneath the North China craton.

Recycling lower continental crust in the North China craton

We summarize the problem of measuring an ultrashort laser pulse and describe in detail a technique that completely characterizes a pulse in time: frequency-resolved optical gating. Emphasis is placed on the choice of experimental beam geometry and the implementation of the iterative phase-retrieval algorithm that together yield an accurate measurement of the pulse time-dependent intensity and phase over a wide range of circumstances. We compare several commonly used beam geometries, displaying sample traces for each and showing where each is appropriate, and we give a detailed description of the pulse-retrieval algorithm for each of these cases.

Measuring ultrashort laser pulses in the time-frequency domain using frequency-resolved optical gating

We present a digital method for solving the phase-retrieval problem of optical-coherence theory: the reconstruction of a general object from the modulus of its Fourier transform. This technique should be useful for obtaining high-resolution imagery from interferometer data.

Reconstruction of an object from the modulus of its Fourier transform.

Three new algorithms for 2D translation image registration to within a small fraction of a pixel that use nonlinear optimization and matrix-multiply discrete Fourier transforms are compared. These algorithms can achieve registration with an accuracy equivalent to that of the conventional fast Fourier transform upsampling approach in a small fraction of the computation time and with greatly reduced memory requirements. Their accuracy and computation time are compared for the purpose of evaluating a translation-invariant error metric.

Efficient subpixel image registration algorithms

Previously it was shown that one can reconstruct an object from the modulus of its Fourier transform (solve the phase-retrieval problem) by using the iterative Fourier-transform algorithm if one has a nonnegativity constraint and a loose support constraint on the object. In this paper it is shown that it is possible to reconstruct a complex-valued object from the modulus of its Fourier transform if one has a sufficiently strong support constraint. Sufficiently strong support constraints include certain special shapes and separated supports. Reconstruction results are shown, including the effect of tapered edges on the object’s support.

Reconstruction of a complex-valued object from the modulus of its Fourier transform using a support constraint

The iterative Fourier-transform algorithm has been demonstrated to be a practical method for reconstructing an object from the modulus of its Fourier transform (i.e., solving the problem of recovering phase from a single intensity measurement). In some circumstances the algorithm may stagnate. New methods are described that allow the algorithm to overcome three different modes of stagnation: those characterized by (1) twin images, (2) stripes, and (3) truncation of the image by the support constraint. Curious properties of Fourier transforms of images are also described: the zero reversal for the striped images and the relationship between the zero lines of the real and imaginary parts of the Fourier transform. A detailed description of the reconstruction method is given to aid those employing the iterative transform algorithm.

/pdf/phase-retrieval-stagnation-problems-and-solutions-aeeh4hay04.pdf

James R. Fienup

Papers

Phase retrieval algorithms: a comparison.

Reconstruction of an object from the modulus of its Fourier transform.

Efficient subpixel image registration algorithms

Reconstruction of a complex-valued object from the modulus of its Fourier transform using a support constraint

Phase-retrieval stagnation problems and solutions