main.asv

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

clear all 
close all
%------ construction d'une elevation synthtetique = une Gaussienne------- 
% x=[-5:0.0469:5]       ;
% y=x                 ;
x=1:256;
y=x;
pixelx=length(x)    ;
pixely=length(y)    ;
[X,Y]=meshgrid(x,y) ;
PhaDerou=100*exp(-(X-128).^2/10000-(Y-128).^2/10000);

%PhaDerou=60*exp(-0.3*(X+2).^2-0.3*(Y+2).^2)-100*exp(-0.3*(X-2).^2-0.3*(Y-2).^2);
G0  =   mod(PhaDerou,2*pi)       ;%phase confin閑 entre 0 et 2*pi
 
%G1= pi*cos(PhaDerou);

% %----ajout de bruit additif------------
moyenne   =0.6      ;
variance  =0.9     ;%pour une variance du bruit d閜asant 3 la recontruction n'est plus robuste
n = moyenne + (sqrt(variance) * randn(pixelx,pixely));%bruit de distribution normale
                                                %moyenne=0 variance=1 
G1  =   G0  +   n   ;% ajout d'un bruit additif

%v_nul=[0;0;0];

      G1(60:80,200:220,1)=zeos

     % G1(200:220,60:80,2)=0;
%     
     % G1(200:220,60:80,3)=0;
% %---- ajout de bruit multiplicatif ------------
% power=2;%modulation en amplitude du bruit
% n=sqrt(power/2)*random('norm',0,1,pixelx,pixely);%
% G2=G1.*n;% ajout d'un bruit multiplicatif

%%%%%%%%%%%%%%%  ZpiM %%%%%%%%%%%%%%%%%
%============== zstep ================


%phder = zstep(G1)   ;


%---- plot --------------
figure(1);
mesh(x,y,PhaDerou);%title('Gaussienne')
figure(2);
imagesc(G0);%title('interferogram')
figure(3);
imagesc(G1);%title('ajout de bruit multiplicatif')
% figure(4);
% imagesc(phder);%title('phase d閞oul閑')
% figure(5);
% mesh(x,y,PhaDerou);%title('Relief synth閠ique')
% figure(6);
% mesh(x,y,phder);%title('Relief reconstruit')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


% %============== pistep =================
% 
% clear phder
% phder=zpim(G2);
% 
% %---- plot --------------
% figure(3);
% subplot(121);imagesc(phder);title('phase d閞oul閑 et liss閑')
% subplot(122);mesh(x,y,phder);title('Relief reconstruit apr鑣 lissage')