% PUMA (Phase Unwrapping via MAxflow) 
%   demo4 - High phase rate Gaussian Hill with a non vertical and non horizontal 
%   alined sector set to zero. 

%   The clique potential is given by  V(x)= |Q2p(x)|^2 where Q2p is the 
%   quantization function; the employed cliques are
%   given by vectors [1 0], [1 1], and [-1 1].

%   Notice: the error image rendered at the end of the script represents
%   the difference between the original (not wrapped) image and the PUMA
%   unwrapped solution.

%   For further details about using PUMA please see puma_ho.m file. See 
%   also: J. Bioucas-Dias and G. Valad鉶, "Phase Unwrapping via Graph Cuts"
%   IEEE Transactions Image Processing, 2007 (to appear).

% Jos