Improved phase imaging from intensity measurements in multiple planes.

Problems stemming from quantitative phase imaging from intensity measurements play a key role in many fields of physics. Techniques based on the transport of intensity equation require an estimate of the axial derivative of the intensity to invert the problem. Derivation formulas in two adjacent planes are commonly used to experimentally compute the derivative of the irradiance. Here we propose a formula that improves the estimate of the derivative by using a higher number of planes and taking the noisy nature of the measurements into account. We also establish an upper and lower limit for the estimate error and provide the distance between planes that optimizes the estimate of the derivative.