Gerchberg–Saxton algorithm
Nowadays, Gerchberg–Saxton algorithm is a topic that has captured the attention of many people around the world. With the advancement of technology and globalization, Gerchberg–Saxton algorithm has become an issue of universal relevance that impacts different sectors of society. Whether on a personal, professional or social level, Gerchberg–Saxton algorithm has generated widespread debate and has sparked the interest of experts and fans alike. In this article, we will thoroughly explore the impact of Gerchberg–Saxton algorithm and discuss its implications on our daily lives. From its origins to its current evolution, Gerchberg–Saxton algorithm represents a relevant phenomenon that deserves to be understood in its entirety.
The Gerchberg–Saxton (GS) algorithm is an iterative phase retrieval algorithm for retrieving the phase of a complex-valued wavefront from two intensity measurements acquired in two different planes.[1] Typically, the two planes are the image plane and the far field (diffraction) plane, and the wavefront propagation between these two planes is given by the Fourier transform. The original paper by Gerchberg and Saxton considered image and diffraction pattern of a sample acquired in an electron microscope.
It is often necessary to know only the phase distribution from one of the planes, since the phase distribution on the other plane can be obtained by performing a Fourier transform on the plane whose phase is known. Although often used for two-dimensional signals, the GS algorithm is also valid for one-dimensional signals.
The pseudocode below performs the GS algorithm to obtain a phase distribution for the plane "Source", such that its Fourier transform would have the amplitude distribution of the plane "Target".
The Gerchberg-Saxton algorithm is one of the most prevalent methods used to create computer-generated holograms.[2]
Pseudocode algorithm
Let:
FT – forward Fourier transform
IFT – inverse Fourier transform
i – the imaginary unit, √−1 (square root of −1)
exp – exponential function (exp(x) = ex)
Target and Source be the Target and Source Amplitude planes respectively
A, B, C & D be complex planes with the same dimension as Target and Source
Amplitude – Amplitude-extracting function:
e.g. for complex z = x + iy, amplitude(z) = sqrt(x·x + y·y)
for real x, amplitude(x) = |x|
Phase – Phase extracting function:
e.g. Phase(z) = arctan(y / x)
end Let
algorithm Gerchberg–Saxton(Source, Target, Retrieved_Phase) is
A := IFT(Target)
while error criterion is not satisfied
B := Amplitude(Source) × exp(i × Phase(A))
C := FT(B)
D := Amplitude(Target) × exp(i × Phase(C))
A := IFT(D)
end while
Retrieved_Phase = Phase(A)
This is just one of the many ways to implement the GS algorithm. Aside from optimizations, others may start by performing a forward Fourier transform to the source distribution.
See also
References
- ^ Gerchberg, R. W.; Saxton, W. O. (1972). "A practical algorithm for the determination of the phase from image and diffraction plane pictures" (PDF). Optik. 35: 237–246. Archived from the original (PDF) on March 28, 2016.
- ^ Memmolo, Pasquale; Miccio, Lisa; Merola, Francesco; Paciello, Antonio; Embrione, Valerio; Fusco, Sabato; Ferraro, Pietro; Antonio Netti, Paolo (2014-01-01). "Investigation on specific solutions of Gerchberg–Saxton algorithm". Optics and Lasers in Engineering. 52: 206–211. doi:10.1016/j.optlaseng.2013.06.008. ISSN 0143-8166.
External links
- Dr W. Owen Saxton's pages Archived 2008-06-13 at the Wayback Machine, Archived 2018-03-16 at the Wayback Machine
- Applications and publications on phase retrieval from the University of Rochester, Institute of Optics Archived 2006-05-09 at the Wayback Machine
- A Python-Script of the GS by Dominik Doellerer
- MATLAB GS algorithms ,