PCG_unwrap_2D(ph_ni​i, ph_file, mask_nii, max_iter, epsi_con, N)

2D weighted Laplace unwrapping with preconditioned conjugate gradient algorithm

This phase unwrapping algorithm is based on the theory and c++ code from: "Two-Dimensional Phase Unwrapping, Theory, Algorithms, and Software" written by Dennis C. Ghiglia and Mark D. Pritt
This program uses weighted least-squares phase unwrapping technique (i.e. weighted Laplace unwrapping) with a solution using iterative algorithm, namely preconditioned conjugate gradient (PCG).
In order to get rid of the residual phase ambiguities after unwrapping process (which occur for all least-squares approaches) a congruence operation is applied.

In order to run this script the following files from MATLAB exchange server are required:

1) http://de.mathworks.com/matlabcentral/fileexchange/8797-tools-for-nifti-and-analyze-image
2) https://de.mathworks.com/matlabcentral/fileexchange/30853-field-mapping-toolbox/content/vector.m
3) http://www.mathworks.com/matlabcentral/fileexchange/25634-fast--n-easy-smoothing/content/smoothn.m
4) http://de.mathworks.com/matlabcentral/fileexchange/26040-dct-and-dst--+-inverse--in-arbitrary-dimension/content/dctn.m