grad.m

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

function y =grad(x,y,v_x,v_y,nb,lambda)

  y=sum(2*pi*(x-v_x(1:nb)+y-v_y(1:nb))+lambda*(x-v_x(1:nb)-round(x-v_x(1:nb))));
  return