voisinage.m

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

function syst_voisi= voisinage(I,i,j)
[n,m]=size(I);
    pos=(i-1)*n+j;
    neighboors=pos+[-1 1 -n m];
    neighboors(find([j==1 j==m i==1 i==m]))=[];
    syst_voisi=I(neighboors);
    return