icmstep.m

Matlab code for encoding an unwrapping phase InSAR image based on Markov Random field

function  psim = icmstep(eta,mph,beta)
%   psim = icmstep(eta,mph,beta);
%   Compute the maximizer of
%   beta*cos(psi-eta) - (psi-mph)^2
%   for psi \in [-pi,pi];


N = 300; % course grid
delta = pi/N;
psi = linspace(-pi-delta,pi+delta,N);
[dummy i] = max(beta*cos(psi-eta) - (psi-mph).^2);
psim = psi(i);