“…Many effective numerical methods have been proposed to solve the TIE, such as the Green's function method [3,7], the Zernike polynomial expansion method [8], the fast Fourier transform (FFT) method [9,10], and the multi-grid method [10,11]. TIE has further been applied in a variety of applications including adaptive optics [11,12], microscopy [13,14], x-ray optics [15,16], etc, but the fundamental problems inherent in these methods are the non-linear effect of the higher-order intensity derivatives, the noise effect [17], and the finite optical coherence [16,18,19] which are essential considerations in industrial applications. Multiple intensity measurements [20][21][22] can mitigate the non-linear effect by estimating the higher-order intensity derivatives and thus can improve the accuracy of the retrieved phase.…”