Quantitative phase reconstruction from input-output intensities using a shift-invariant linear imaging system

We formulate an iteration scheme for quantitative two-dimensional phase retrieval of complex scalar wave fields from input-output intensity profiles using a coherent rotationally symmetric shift-invariant linear optical imaging system. This method is a generalization of our previous method for quantitative single-image pure phase reconstruction using an imperfect shift-invariant linear imaging system. The method is somewhat analogous to the Gerchberg-Saxton iterative phase-retrieval algorithm, but typically converges much more rapidly. We have demonstrated the efficacy of the method using simulated data in the presence of noise.