Half-quadratic cost functions for phase unwrapping

We present a generic regularized formulation, based on robust half-quadratic regularization, for unwrapping noisy and discontinuous wrapped phase maps. Two cases are presented: the convex case and the nonconvex case. The unwrapped phase with the convex formulation is unique and robust to noise; however, the convex function solution deteriorates as a result of real discontinuities in phase maps. Therefore we also present a nonconvex formulation that, with a parameter continuation strategy, shows superior performance. (C) 2004 Optical Society of America.