基于累计残差图进行的相位解包裹算法

% Code 1 for testing the Unwrapping method ARM corresponding to
%
% Mariano Rivera, Francisco Hernandez-Lopez and Adonai Gonzalez, "Phase 
% unwrapping by accumulation of residual maps," to appear in Optics 
% and Lasers in Engineering, 2014.
%
% Mariano Rivera
% July 2014

clear all;
close all;

[GT_UnWrap, Wrap, ~, ~] = syntheticphase(0.0, pi);
Mask = ones(size(Wrap));

%%

%parameters for syntetic experiment 
lambda = 1e-1;      % smoothness
mu      = pi/10;    % edge detection sensitivity
numIter = 4000;     % number of iterations per level
numlevels = 5;      % num of levels in the multigrid pyramid

% UnwrapPhase
UnWrap = nestedMultigridUnwrapARM(double(Wrap),double(Mask),...
                                           lambda, mu, numIter, numlevels);
%%
%Mean Square Error
mnGTuwp = min(GT_UnWrap(:));
mnuwp = min(UnWrap(:));

GT_UnWrap = 2*pi./(max(GT_UnWrap(:))-mnGTuwp) .*(GT_UnWrap-mnGTuwp)-pi;
E_UnWrap  = 2*pi./(max