Adaptive-additive algorithm

Last updated

In the studies of Fourier optics, sound synthesis, stellar interferometry, optical tweezers, and diffractive optical elements (DOEs) it is often important to know the spatial frequency phase of an observed wave source. In order to reconstruct this phase the Adaptive-Additive Algorithm (or AA algorithm), which derives from a group of adaptive (input-output) algorithms, can be used. The AA algorithm is an iterative algorithm that utilizes the Fourier Transform to calculate an unknown part of a propagating wave, normally the spatial frequency phase (k space). This can be done when given the phase’s known counterparts, usually an observed amplitude (position space) and an assumed starting amplitude (k space). To find the correct phase the algorithm uses error conversion, or the error between the desired and the theoretical intensities.

Contents

The algorithm

History

The adaptive-additive algorithm was originally created to reconstruct the spatial frequency phase of light intensity in the study of stellar interferometry. Since then, the AA algorithm has been adapted to work in the fields of Fourier Optics by Soifer and Dr. Hill, soft matter and optical tweezers by Dr. Grier, and sound synthesis by Röbel.

Algorithm

  1. Define input amplitude and random phase
  2. Forward Fourier Transform
  3. Separate transformed amplitude and phase
  4. Compare transformed amplitude/intensity to desired output amplitude/intensity
  5. Check convergence conditions
  6. Mix transformed amplitude with desired output amplitude and combine with transformed phase
  7. Inverse Fourier Transform
  8. Separate new amplitude and new phase
  9. Combine new phase with original input amplitude
  10. Loop back to Forward Fourier Transform

Example

For the problem of reconstructing the spatial frequency phase (k-space) for a desired intensity in the image plane (x-space). Assume the amplitude and the starting phase of the wave in k-space is and respectively. Fourier transform the wave in k-space to x space.

Then compare the transformed intensity with the desired intensity , where

Check against the convergence requirements. If the requirements are not met then mix the transformed amplitude with desired amplitude .

where a is mixing ratio and

.

Note that a is a percentage, defined on the interval 0 ≤ a ≤ 1.

Combine mixed amplitude with the x-space phase and inverse Fourier transform.

Separate and and combine with . Increase loop by one and repeat.

Limits

  • If then the AA algorithm becomes the Gerchberg–Saxton algorithm.
  • If then .

See also

Related Research Articles

Additive synthesis is a sound synthesis technique that creates timbre by adding sine waves together.

<span class="mw-page-title-main">Wave interference</span> Phenomenon resulting from the superposition of two waves

In physics, interference is a phenomenon in which two waves combine by adding their displacement together at every single point in space and time, to form a resultant wave of greater, lower, or the same amplitude. Constructive and destructive interference result from the interaction of waves that are correlated or coherent with each other, either because they come from the same source or because they have the same or nearly the same frequency. Interference effects can be observed with all types of waves, for example, light, radio, acoustic, surface water waves, gravity waves, or matter waves.

<span class="mw-page-title-main">Astronomical seeing</span> Atmospheric distortions of light

In astronomy, seeing is the degradation of the image of an astronomical object due to turbulence in the atmosphere of Earth that may become visible as blurring, twinkling or variable distortion. The origin of this effect are rapidly changing variations of the optical refractive index along the light path of the object. Seeing is a major limitation to the angular resolution in astronomical observations with telescopes that would otherwise be limited through diffraction by the size of the telescope aperture. Today, many large scientific ground-based optical telescopes include adaptive optics to overcome seeing.

Fourier optics is the study of classical optics using Fourier transforms (FTs), in which the waveform being considered is regarded as made up of a combination, or superposition, of plane waves. It has some parallels to the Huygens–Fresnel principle, in which the wavefront is regarded as being made up of a combination of spherical wavefronts whose sum is the wavefront being studied. A key difference is that Fourier optics considers the plane waves to be natural modes of the propagation medium, as opposed to Huygens–Fresnel, where the spherical waves originate in the physical medium.

<span class="mw-page-title-main">Mach–Zehnder interferometer</span>

The Mach–Zehnder interferometer is a device used to determine the relative phase shift variations between two collimated beams derived by splitting light from a single source. The interferometer has been used, among other things, to measure phase shifts between the two beams caused by a sample or a change in length of one of the paths. The apparatus is named after the physicists Ludwig Mach and Ludwig Zehnder; Zehnder's proposal in an 1891 article was refined by Mach in an 1892 article. Demonstrations of Mach–Zehnder interferometry with particles other than photons had been demonstrated as well in multiple experiments.

<span class="mw-page-title-main">Superposition principle</span> Fundamental physics principle stating that physical solutions of linear systems are linear

The superposition principle, also known as superposition property, states that, for all linear systems, the net response caused by two or more stimuli is the sum of the responses that would have been caused by each stimulus individually. So that if input A produces response X and input B produces response Y then input produces response.

In optics, an ultrashort pulse, also known as an ultrafast event, is an electromagnetic pulse whose time duration is of the order of a picosecond or less. Such pulses have a broadband optical spectrum, and can be created by mode-locked oscillators. Amplification of ultrashort pulses almost always requires the technique of chirped pulse amplification, in order to avoid damage to the gain medium of the amplifier.

<span class="mw-page-title-main">Spectral phase interferometry for direct electric-field reconstruction</span>

In ultrafast optics, spectral phase interferometry for direct electric-field reconstruction (SPIDER) is an ultrashort pulse measurement technique originally developed by Chris Iaconis and Ian Walmsley.

The Strehl ratio is a measure of the quality of optical image formation, originally proposed by Karl Strehl, after whom the term is named. Used variously in situations where optical resolution is compromised due to lens aberrations or due to imaging through the turbulent atmosphere, the Strehl ratio has a value between 0 and 1, with a hypothetical, perfectly unaberrated optical system having a Strehl ratio of 1.

Multiphoton intrapulse interference phase scan (MIIPS) is a method used in ultrashort laser technology that simultaneously measures, and compensates femtosecond laser pulses using an adaptive pulse shaper. When an ultrashort laser pulse reaches a duration of less than a few hundred femtosecond, it becomes critical to characterize its duration, its temporal intensity curve, or its electric field as a function of time. Classical photodetectors measuring the intensity of light are still too slow to allow for a direct measurement, even with the fastest photodiodes or streak cameras.

Molecular replacement is a method of solving the phase problem in X-ray crystallography. MR relies upon the existence of a previously solved protein structure which is similar to our unknown structure from which the diffraction data is derived. This could come from a homologous protein, or from the lower-resolution protein NMR structure of the same protein.

Phase retrieval is the process of algorithmically finding solutions to the phase problem. Given a complex signal , of amplitude , and phase :

In statistical signal processing, the goal of spectral density estimation (SDE) or simply spectral estimation is to estimate the spectral density of a signal from a sequence of time samples of the signal. Intuitively speaking, the spectral density characterizes the frequency content of the signal. One purpose of estimating the spectral density is to detect any periodicities in the data, by observing peaks at the frequencies corresponding to these periodicities.

Quantum mechanics was first applied to optics, and interference in particular, by Paul Dirac. Richard Feynman, in his Lectures on Physics, uses Dirac's notation to describe thought experiments on double-slit interference of electrons. Feynman's approach was extended to N-slit interferometers for either single-photon illumination, or narrow-linewidth laser illumination, that is, illumination by indistinguishable photons, by Frank Duarte. The N-slit interferometer was first applied in the generation and measurement of complex interference patterns.

The closure phase is an observable quantity in imaging astronomical interferometry, which allowed the use of interferometry with very long baselines. It forms the basis of the self-calibration approach to interferometric imaging. The observable which is usually used in most "closure phase" observations is actually the complex quantity called the triple product. The closure phase is the phase (argument) of this complex quantity.

An anamorphic stretch transform (AST) also referred to as warped stretch transform is a physics-inspired signal transform that emerged from time stretch dispersive Fourier transform. The transform can be applied to analog temporal signals such as communication signals, or to digital spatial data such as images. The transform reshapes the data in such a way that its output has properties conducive for data compression and analytics. The reshaping consists of warped stretching in the Fourier domain. The name "Anamorphic" is used because of the metaphoric analogy between the warped stretch operation and warping of images in anamorphosis and surrealist artworks.

<span class="mw-page-title-main">Phase stretch transform</span>

Phase stretch transform (PST) is a computational approach to signal and image processing. One of its utilities is for feature detection and classification. PST is related to time stretch dispersive Fourier transform. It transforms the image by emulating propagation through a diffractive medium with engineered 3D dispersive property. The operation relies on symmetry of the dispersion profile and can be understood in terms of dispersive eigenfunctions or stretch modes. PST performs similar functionality as phase-contrast microscopy, but on digital images. PST can be applied to digital images and temporal data. It is a physics-based feature engineering algorithm.

In crystallography, direct methods is a set of techniques used for structure determination using diffraction data and a priori information. It is a solution to the crystallographic phase problem, where phase information is lost during a diffraction measurement. Direct methods provides a method of estimating the phase information by establishing statistical relationships between the recorded amplitude information and phases of strong reflections.

Kernel-phases are observable quantities used in high resolution astronomical imaging used for superresolution image creation. It can be seen as a generalization of closure phases for redundant arrays. For this reason, when the wavefront quality requirement are met, it is an alternative to aperture masking interferometry that can be executed without a mask while retaining phase error rejection properties. The observables are computed through linear algebra from the Fourier transform of direct images. They can then be used for statistical testing, model fitting, or image reconstruction.

Spectral interferometry (SI) or frequency-domain interferometry is a linear technique used to measure optical pulses, with the condition that a reference pulse that was previously characterized is available. This technique provides information about the intensity and phase of the pulses. SI was first proposed by Claude Froehly and coworkers in the 1970s.

References