Gerchberg-Saxton Phase Retrieval

Iterative algorithm: recover phase from amplitude constraints in image & Fourier domains

Iter: 0
Target Amplitude |T(x,y)|
Current Estimate
Fourier Amplitude |F(u,v)|
Retrieved Phase φ(x,y)
Error vs Iteration
GS Algorithm: Alternate projections — enforce target amplitude in image domain, free phase in Fourier domain