3D Wave field phase retrieval from multi-plane observations

We reconstruct a spatially distributed three-dimensional (3D) wave field from a number of intensity observations, obtained in different sensor planes, parallel to the object plane. The proposed algorithm can be treated as a multiple plane iterative Gerchberg-Saxton algorithm [1]. It is obtained from the best linear estimate of the complex-valued object distribution derived for the complex-valued observations. This estimator is modified for the intensity measurements in the sensor planes. The algorithm is studied by numerical experiments performed for amplitude and phase object distributions. It is shown that the proposed method allows reconstructing the whole 3D wave fields for different setup parameters. This technique can be applied for 3D imaging. The comparison versus the successive iterative method [2] shows an accuracy advantage of the proposed algorithm provided that the type of modulation in the object plane is known.