Gerchberg–Saxton algorithm: Difference between revisions

Revision as of 13:47, 24 May 2006

Latest revision as of 16:19, 21 May 2025

(76 intermediate revisions by 48 users not shown)

Line 1:

{{Short description|Algorithm for phase retrieval}}

'''Gerchberg Saxton Algorithm'''

[[File:Gerchberg-Saxton algorithm.jpg|thumb|400px|The Gerchberg-Saxton algorithm. FT is Fourier transform.]]

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.<ref>{{Cite journal|last=Gerchberg|first=R. W.|last2=Saxton|first2=W. O.|date=1972|title=A practical algorithm for the determination of the phase from image and diffraction plane pictures|url=http://www.u.arizona.edu/~ppoon/GerchbergandSaxton1972.pdf|archive-url=https://web.archive.org/web/20160328053000/http://www.u.arizona.edu/~ppoon/GerchbergandSaxton1972.pdf|url-status=dead|archive-date=March 28, 2016|journal=Optik|language=EN|volume=35|pages=237–246}}</ref> 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 Gerchberg Saxton (GS) algorithm is an algorithm used for retrieving the phase of a pair of light distributions (or any other mathematically valid distribution) related via Fourier transform if their intensities at their respective optical planes are known.

ItThe is[[pseudocode]] oftenbelow necessary to know onlyperforms the phaseGS distributionalgorithm fromto oneobtain of the planes, since thea phase distribution onfor the other plane can"Source", besuch obtainedthat by performing aits Fourier transform onwould have the planeamplitude whosedistribution phase is known. Although often used for two dimensional signalsof the GS algorithm is also valid for 1-dplane signals"Target".

The Gerchberg-Saxton algorithm is one of the most prevalent methods used to create [[computer-generated hologram]]s.<ref>{{Cite journal |last=Memmolo |first=Pasquale |last2=Miccio |first2=Lisa |last3=Merola |first3=Francesco |last4=Paciello |first4=Antonio |last5=Embrione |first5=Valerio |last6=Fusco |first6=Sabato |last7=Ferraro |first7=Pietro |last8=Antonio Netti |first8=Paolo |date=2014-01-01 |title=Investigation on specific solutions of Gerchberg–Saxton algorithm |url=https://www.sciencedirect.com/science/article/pii/S0143816613001942 |journal=Optics and Lasers in Engineering |volume=52 |pages=206–211 |doi=10.1016/j.optlaseng.2013.06.008 |issn=0143-8166|url-access=subscription }}</ref>

The Pseudo-code below performs the GS algorithm to obtain a phase distribution for the plane, Source, such that its Fourier transform would have the amplitude distrribution of the plane, Target.

===Pseudo-codePseudocode algorithm===

'''Let:'''

===Pseudo-code algorithm===

FT -– Forwardforward Fourier Transformtransform

IFT -– Inverseinverse Fourier Transformtransform

Let:

''i'' – the imaginary unit, √&minus;1 (square root of &minus;1)

exp – exponential function (exp(x) = ''e''^''x'')

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- Amplitde extracting function:

ege.g. for complex ''z'' = ''x'' + ''iy'', Amplitudeamplitude(''z'') = sqrt(''x*''·''x'' + ''y*''·''y'')

for real ''x'', amplitude(''x'') = |''x''|

Phase - Phase extracting function:

Phase eg.– Phase(z) =extracting arctan(y/x)function:

e.g. Phase(z) = arctan(y / x)

FT - Forward Fourier Transform

'''end of Let:'''

IFT - Inverse Fourier Transform

end of Let:

'''algorithm''' Gerchberg–Saxton(Source, Target, Retrieved_Phase) '''is'''

A := IFT(Target)

Gercberg Saxton Algorithm(Source, Tareget)

'''while''' error criterion is not satifiedsatisfied

A = IFT(Target)

B := Amplitude(Source) × exp(i × Phase(A))

while error criterion is not satified

B = Amplitude(Source) * exp C := FT(Phase(A)B)

D := Amplitude(Target) × exp(i × Phase(C))

C = FT(B)

D = Amplitude(Target) * exp A := IFT(Phase(C)D)

A ='''end IFT(D)while'''

Retrieved_Phase = Phase(A)

end while

Retrieved Phase = Phase(A)

end of Gerchberg Saxton Algorithm

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.

==RelatedSee Articlesalso==

* [[FourierPhase opticsretrieval]]

* [[Fourier optics]]

* [[Holography]]

* [[Adaptive-additive algorithm]]

==External LinkReferences==

{{reflist}}

* [http://www.ysbl.york.ac.uk/~cowtan/fourier/coeff.html Graphical explanatory material by Kevin Cowtan]

==External links==

[[Category:Algorithms]]

* Dr W. Owen Saxton's pages [http://www-hrem.msm.cam.ac.uk/people/saxton/] {{Webarchive|url=https://web.archive.org/web/20080613024950/http://www-hrem.msm.cam.ac.uk/people/saxton/ |date=2008-06-13 }}, [https://www.murrayedwards.cam.ac.uk/fellows/dr-w-owen-saxton]

[[Category:Physics]]

* [http://www.optics.rochester.edu/workgroups/fienup/index.html Applications and publications on phase retrieval from the University of Rochester, Institute of Optics]

[[Category:Optics]]

* [https://www.creatgraphy.com/05/2020/allgemein/lightmodulation-gerchberg-saxton-algorithmus-gsa-660/ A Python-Script of the GS by Dominik Doellerer]

* MATLAB GS algorithms [https://ch.mathworks.com/matlabcentral/fileexchange/68647-gerchberg-saxton-phase-retrieval-algorithm/], [https://ch.mathworks.com/matlabcentral/fileexchange/65979-gerchberg-saxton-algorithm]

{{DEFAULTSORT:Gerchberg-Saxton algorithm}}

[[Category:Digital signal processing]]

[[Category:AlgorithmsPhysical optics]]

[[Category:Articles with example pseudocode]]