Bayesian Inference for Phase Unwrapping Using Conjugate Gradient Method in One and Two Dimensions

We investigated statistical performance of Bayesian inference using maximum entropy and MAP estimation for several models which approximated wave-fronts in remote sensing using SAR interferometry. Using Monte Carlo simulation for a set of wave-fronts generated by assumed true prior, we found that the method of maximum entropy realized the optimal performance around the Bayes-optimal conditions by using model of the true prior and the likelihood representing optical measurement due to the interferometer. Also, we found that the MAP estimation regarded as a deterministic limit of maximum entropy almost achieved the same performance as the Bayes-optimal solution for the set of wave-fronts. Then, we clarified that the MAP estimation perfectly carried out phase unwrapping without using prior information, and also that the MAP estimation realized accurate phase unwrapping using conjugate gradient (CG) method, if we assumed the model of the true prior appropriately.

• Yohei Saika
• Hiroki Sakaematsu
• Shota Akiyama

• Bayesian inference using maximum entropy
• MAP estimation using conjugate gradient method
• SAR interferometry.

[1] J. R. Jensen, Remote "Sensing of the Environment: An Earth Resource
Perspective (2nd Edition)," (Prentice Hall; 2 edition) 2006, ch. 2.
[2] D. C. Ghiglia and M. D. Pritt, "Two-Dimensional Phase Unwrapping
Theory, Algorithm and Software," New York: Wiley, 1998, ch. 3.
[3] R. M. Goldstein and H. A. Zebker, "Interferometric radar mapping of
ocean currents," Nature, vol. 328, pp. 707-709, Aug. 1987.
[4] D. L. Fried, "Least-square fitting of a wave-front distortion estimate to an
array of phase differences measurements," J. Opt. Soc. Am., vol. 67, pp.
370-375, Mar. 1977.
[5] R. H. Hudgin, "Wave-front reconstruction for compensated imaging", J.
Opt. Soc. Am., vol. 67, pp. 375-378, Mar. 1977.
[6] H. Takajyo and T. Takahashi, "Least squares phase estimation from phase
difference", J. Opt. Soc. Am. A, vol. 5, pp. 416-425, Nov. 1988.
[7] D. C. Ghiglia and L. A. Romero, "Robust two-dimensional weighted and
unweighted phase unwrapping that uses fast transforms and iterative
method," J. Opt. Soc. Am. A, vol. 11, pp. 107-117, Jan. 1994.
[8] L. Guerriero, G. Nico, G. Pasquariello and Stramaglia, "A new
regularization scheme for phase unwrapping," Appl. Opt. vol. 37, pp.
3058-3058, May. 1998.
[9] G. Nico, G. Palubinskas and M. Datcu, "Bayesian Approaches to Phase
Unwrapping: Theoretical Study," IEEE Trans. Signal Processing, vol.
48(4), pp. 2545-2556, Sep. 2000.
[10] Y. Saika and H. Nishimori, "Statistical-mechanical Approach to the
Problem of Phase Retrieval by the Q-Ising Model," Progress Theoretical
Physics Supplement, vol. 157, pp. 292-295, 2005.
[11] Y. Saika and H. Nishimori, "Statistical-mechanical approaches to the
problem of phase retrieval in adaptive optics in astronomy," J. Phys : Conf.
Ser., vol. 31, Mar. pp. 169-170, 2006.
[12] Y. Saika and T. Uezu, "Statistical Mechanics of Phase Unwrapping using
the Q-Ising Model," IEICE Technical Report, NLP2012-12(2012-4), pp.
61-65, Apr. 2012.
[13] J. L. Marroquin and M. Rivera, "Quadratic reguralization functionals for
phase unwrapping," J. Opt. Soc. Am. A, vol. 12, pp. 2393-2400, Sep.
1995.
[14] H. Sakaematsu and Y. Saika, "Statistical Performance of Conjugate
Gradient Method for Phase Unwrapping in Adaptive Optics," in Proc. of
2012 12th International Conference on Control, Automation and Systems,
Korea, 2012, pp. 1279-1284.
[15] Yohei Saika, Hiroki Sakaematsu and Shota Akiyama, "Maximum Entropy
and MAP Estimation using Conjugate Gradient Method for Phase
Unwrapping," the International Journal of Science and Engineering
Investigations, to be published.