2D Weighted Phase Unwrapping

In phase retrieval applications, such as interferometry, people are facing the problem of phase being wrapped (i.e. it has value only from -pi to pi). To make the results more useful, one needs to unwrap the phase. This scripts implements the algorithm from Ghiglia and Romero (1994) that uses least-square methods to solve the problem. The algorithm runs relatively quickly and the algorithm is also robust to noise.