main.asv

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

% premier esai pour developper le code de l'algorithme de l'article

clear all
close all
lambda=50;
h=0.05;
nbiter=5000000;
x=1:128;
y=x;
pixelx=length(x);
pixely=length(y);
[X,Y]=meshgrid(x,y);
PhaDerou=100*exp(-(X-64).^2/10000-(Y-64).^2/10000);
G=mod(PhaDerou,2*pi);
F=zeros(pixelx,pixely);
Grad=zeros(pixelx,pixely);
phi=zeros(pixelx,pixely);
F=2*pi*G;
tic
%=======================================================
  for k=1:nbiter
%     for i=1:pixelx
%     for j=1:pixely
i=ceil(pixelx*rand(1));% choix all閍tiore de pixels i et j
j=ceil(pixely*rand(1));
          f_s=F(i,j);
          g_s=G(i,j);
          voisin_f_s=voisinage(F,i,j);% acquisition des valeurs des site voisins selon un syst鑝e 4 connex par rapport