An algorithm of unordered wavefront propagation in terahertz phase retrieval with dense multiplane data acquisition

Iterative phase retrieval algorithms from multiple diffraction patterns in the terahertz (THz) frequency range are a promising tool of computational imaging capable of providing high spatial resolution of reconstructed phase images. One of the commonly used algorithms is SBMIR, which employs multiple intensity distributions of the diffraction object wavefield as input data. Compared with single-frame methods, the multi-plane approach allows for a faster convergence, but requires time-consuming data acquisition from a receiver positioned at a variety of distances from the object. Previously, we proposed a method for THz data acquisition in a single scan mode, which allows one to quickly obtain an exhaustive set of diffraction distributions. In this paper we evaluate an up-to-date phase retrieval algorithm based on the SBMIR/R-SBMIR method (which utilizes stochastic wavefront propagation) on the experimental data captured by a single-scan technique. Unlike a number of conventional phase retrieval algorithms, which may require a series of numerical experiments for determining optimal intensity distributions from a large dataset, disordered propagation of the estimation wavefront guarantees the high-contrast and high-resolution image reconstruction without pre-setting the parameters. It is shown that the package use of the single-scan technique with the subsequent data processing using the R-SBMIR algorithm has application potential for automation of the multi-plane phase retrieval in the THz range.