main1.m

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=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);
G1  =   mod(PhaDerou,2*pi)       ;%phase sofin閑 entre 0 et 2*pi
%----ajout de bruit additif------------
moyenne   =0.6      ;
variance  =2        ;%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 
G2  =   G1  +   n   ;% ajout d'un bruit additif
%%%%%%%%%%%%%%%  ZpiM %%%%%%%%%%%%%%%%%
%============== zstep ================


phder = zstep(G2)   ;


%---- plot --------------
figure(1);
subplot(141);imagesc(PhaDerou);title('Gaussienne')
colorbar
subplot(142);imagesc(G1);title('interferogram')
colorbar
subplot(143);imagesc(G2);title('ajout de bruit additif')
subplot(144);imagesc(phder);title('phase d閞oul閑')
colorbar
figure(2);
subplot(121);mesh(x,y,PhaDerou);title('Relief synth閠ique')
subplot(122);mesh(x,y,phder);title('Relief reconstruit')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%