inv2.eps

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

)B/Po!Z`:&q&]__+<DXJrZD+\!$D.[!$(tW!?)jT(_d`$'GV>r'+tog&J#E`%LijU$OdFB$4-tE#RCY@
#6b;:r<3B."U+u0!Wi6%$NU;1!s/T,"9JZ-!r`90!<WB*!WrN+!s/Q+q?$m&!WrN-!rN-%!WrT(!s&N+
rW<0("Si-/!WrQ-"9\f0"U"o2"U+f."9o&6#6t,4"pkSC$OI.JrX8l=%M'*L&-NUY',2,Y'aGNk)&O50
)&iJn!$(tY0HM5T*ZlIF*?H7E+!2LC*$-(?*?6":)]BY5(`4&*()@W!'+tli&.oHb%h9$X$k3+C'F>*R
#mghB"pbD;"U529"pG,5!s8E'$3C;2!WrK,!WiH+rW45F!X&T,!<N?*!WrH*!WiH*!s/Q,!WrK*!<N?+
!s/N+"Si-&"9S`-!r`9+!WrN+!s8W("9AT,rWEH0"9\l2"pG20"pYA<"pYA=#mgM;!t,DD%/p\C$kEgZ
&J5Zg'GD2['a,9h(Ddqf)$_'!)]BY8*$,"u)BBe=*$$"?*#ot=)]Te:)]BV6)&X5,(`4#(rY5MK'+tlU
&H`OS%1N4F#7M"K$OR4H$2b&-#5SK(!sAN(!W`?)rW3Z4!s/N+!s/Q+!WrK+!WiK,!r`9.!WrK,!X&T.
!s8E'!WiH*rW<9+!s/Q-r<*$#rWE-$"p4u0!sA`1r<!'(#6"`<"U"o1"U"u3"pG5:#R:YB#RLk;#ltA@
rX95D$k<dY&.]<`&eblk'E/^V(DR`&)&XA3)B0S6)B0Uo)Z^otrYu(\*#fe8)u^]\)B'G1)&O/+(B"sR
'+toi&J,Nc&,m%X%1NXQ$O[:J$4."D#RCY>#6tD:"pG/5"pG26rW<'$!s/Q+!rE'"!rW3+!WrN+!sA]-
q?%!)!s/N,!sAZ,r<!3)!s/Q,!rW3,"9JZ."9SZ*"9&?'"9S`)!sSl2!sAc1"T/E1"9o&5"U5,6#6tP6
#7q.J#mgqF$OR=P$k*RT%M00M&.0$^'GD/q'G_H#(B,*R(`=2.)B&Mm!#tnUrYlCe)]KY4)B'J2(Dmr*
(Ddi&rY,PK&eYcf&.e[K!XoJHrX/]5rX&N/$3phA#6kA9"p>)4rWW9'!!<-%!W`?(qZ@')!s/Q,!<WE%
!t554!s/T,!sA]-!s8W'!s/K*"9&?%"8r9$"8i3&!s/Q+r<!0(!s8Q*rWE-&$3LD6!sJi2"9Sc0rWNl:
"U529#7(P>#RLbC$47.J$iUP@%1NdW%hB6`&J>!TrY5;HrYH:d(`=/.(`F80)AsG2)Aj;/)Aj>.)#b<P
()?fa%2BQk',(rj&J>Wc&.]6L%LNXR$OdCN#mgkE$2k,0#Q=oL"U5/7"p>&4"9el0!WiH+!WrQ+!s/Q+
!s8T+"9JZ,rW<r<!WrK*!s/N+!W`?)!WrN+!sAZ,!s/T(!rW3%!X/K($3L>3!WiK-!s&K*rW<0&!r`<3
!sA]/!sJi2"U,#5#6k)2!!W?,)@-WT#7(YC$OI.J$4I@P%1NdV%M'*^&J>ZgrY#DK'bhB!r>#JP(Ddkd
(]YBM(^:lp(`+#)(]>*V'GV;s'b_2n&J>]T&-*7?%fQn<$Om"ArWrW3#R^P9!<rf."p+r3#6"c0"9eo1
!s8E'$3C;4!s/N+!s&N,rW30)!s8]*!X/Z-!sA],quciu#6P)2!W`B*"TAK-!WiH*!s/W)!s&H+rW<<,
!WrK+"8r?)!s/Q,rWEN0"9el1"9eu6"oSW/#6tG<#71;8!=B25$iCG4%L3LT&.oHc&eY'U!Yc=arY>DJ
#8\0o(D[f()#P0N()Hf_#8@md'+trk')`CN&.T3]%1WgU%1<RO$iUJB#RUkC#mghA"pY>:":"f.$jHh<
"9el/!s/N)!sJT*r;m6,!s&K*"9S],rW<3)!sAK'!s8T+"TAH7!s/Q+!sAZ-!<WB*!WiH+"9\T(!W`E,
rW<0&!rN--"9Sc/!sA]."9Si*"9\l4"U>50"ptS?#m^bA#RLhD$N:D8$k*RTrXB&B%hK6_&J=sS"r%ad
'GD2qrY7.(()I`$()@W$(DR]$'bh>t'G_Ar'+trl&eGWc&.oE_%M''Z%1E^U$Od@L#mghC$2t26#6kD<
"pOu1#Qt;7"9eo3"9e]+!<N?%"9&B$"9&?+!s8T,!s&K&!rrB$!W`B*r;m-)!s&K+!r`9,!WrN+!s/Q,
quQm$!rW3+!WrN-!s/Q-pB16p%L!"?"9er5"UG89#6tS@r<No?$4.%I$471L%1EXS%h9'J%g3LV%hK?O
&e,Ke',))p'bqK!'bh>srY,AJ(&emO',2,o&ePcU&I&dY&.]<_&,m"D$k!IO$jm@K$N184#6tP6#8@@H
"U5,6"pG,6"9\l2"9el.!s8W'!s\l0!W`B)!s/B'!WiH+rW30(!WrQ'!WW?%!WW?%!XSr1!WrK+!WrK*
r<!$$rW<6(!sA]&!s&N+r<"8I!WrQ,"9eo1!sJf2!s8Z/"U,#4"U5/7#7(S?"pbM@#RLe9$N:G4%/p_8
%Kd4Q%hJUN#SI^^',(uk',:BYq@iuE',(3V#SI^]&eGTb&J+jNr=&l<%13INrX/W2rWi]5#6kD;#R:21
!X&Z0rWNK/!sJc/!X&E'"9\c.!s/B&r<!'%"9&?3!sA]-!s/Q*!s8Q+!X&Q'!W`B)rW<N0!s/Q,!WrN+
!s&K"!sJ`.!WrN-qu[9/"9Sc/"9J]/!sSZ,!X/]0r<30+rWWl>#6tG<#71V@#RLhF#mq"<$NpnI%1E^F
%KI%=&.&jX&.oKe&eP`ir=])C!##8E#o"!b&JGcj&J,QP&-i^U%LipX%1`@H#R^tI$4.%I#n$Y;!=&i-
#6>#4r<<f9!s&K,"9Sc0!s/N+!s/Q,r;n2G"9\f.!s8T*!s8T,!X&Q*"9JZ,!s/N+!s/Q,!WrK)!<WE$
"!.LF"9AT,!WiH*!WrQ-!s/N,!WrN+!sAZ-!WrQ-"9&E*!sA`.rWE?-"9o&5r<3K3"U>8<#R:V@#Q4u.
$3LYC$P!(D!t>_L%K6k@%1a$]&H37H&eP`g&J>ZfrY#/A!>6"C&H!%C%M''Z%LimH%/pV;#n%(H$N:;2
$2t5/#Q=o."U=r0"9\o2"pOu/"Tni.!X&W)!WrN+!s8B&"9JZ,!s8H("9JW,!sJT)#m:>4!X&T,!s/N&
!XSr3!s/N+!s/N,rW<9)"9S`,q?$a$"T/?5!WrN-"Tni0!s8]/!X/]1"o\W+"U"u-"UkS>#RCV?#R1M>
#RLh9$3LYE$k*%C$4RFT%LipZ&.]<_q@ElD&.oHc&.oNQ&.&s\&.]9`%hK0\rXJr<%K-\@$4@4K#mptG
$2Xu3#7(P?"pY&1!!`E*!WrT.rWE6*"oAB&!X&E'"Tni-!sA])"TAK("9SQ)!<N?#"9&B!!sJ`.!WiK*
rW<6(!s8W'!s/N*"9&B$!so#2"9S].!s8Z-rW<E-"9Sf1"9Sf0r<*Z8"U"u4"pG28"U#&8"pYD=rWrH/
!"/Z4%h&jS$k<[R%1NgW$k<aX&H34<&-3:Pr=B8G&.]<`%h9*Z&,ctJ%1WdT%1EUO$jm@N#mgt<#lk81
#6G/8#lY&-"UG;8"pG/5!sJc*":>56!sJc-!s/N+!rW3/"9JW,"9S`-!WrK*r<!-&!s8W(!t#)3!s/N+
!WrK+!rN-&!WrN+rW<9+"9S`-r<!*&!sAN($3L>4!WrN+!sA]/qu[*(!sA]/"o\W)"oSQ."U5,6"pP>3
#R^qD#RC_B#mgqH$4?b?q[=&G%LijW%M''[&.]6]&.]6^%fQt>%M/:DrXAi9rX/l=$O[7K$N(2H#71Y@
#7(P?"pY;9#6k>7"U"u3"U,#3rW<K1!s/Q,!W`?(!s8<$"9JZ-!Wi9&rW<'$!!<$"!<NE%!sJf0!WrN+
rW<E-!s8Z.!s/Q,r;mQ5!WiH*!s/T-"9\c."9Sc.!rW3&"9Sc)";:qA"U+u3"9er5"U,)6"U>8<#7(56
!XT/=r<`l<$4I:M$k*IO%1N.DrXJo=q[N`>%h9'\rXSo:!Y,\Lr=&Z5!=B56$3pkE#mgkC#7:A8!!`E-
%L*.A"pP58"U"u3"p4u3rWEN2!sA].!WiB&!W<*4!WiH*!s8T-"9Sc.!X&T,!s8E'!!N3$$3L>3!<WE*
!<WE,rW3f<!WrN-!<WH,!<WH*"9Sc/!X&T.!rE'+!X&Q,!sAZ.!X/K)"U"u3"9er+"T\f2r!*'(('k-O
#m^eD#mptF$O[:K$O[@N$k!IPr!`c>%1WgVrXSo:"UkkN%1E^E$lTQ`$OdCM$4@4J$47(G#RC_B#7(P>
"pYA;"TAW&"VLqA"9Sf0!sAc1!sA`/!s/T-!r`9*"9JZ.!X&E'!WrN*rW3<,!<N?)!s/Q&!sAZ,!s/Q!
!rrK'!s8Q*!X&B&!sA].!r`9(!WrN+q#^X#"9&H%"9&E1"U"u4"9\l2"9er*"TSf+"pG59#R1P>rWj;H
$47(H$4@4K$OR7K$OmIQ%1<RR$k3^G$kj'Z$k*RR%L`dU$O[CO$OR:M$4Hh?rWrW3$46\:!sf/;#Q+c.
"U,)6rWX2B"p>&2"9el1!sA`/"9JZ-!WrN-!sAW+r<!*&!s8H(!WiH*rW<B,!<WE*!sA])"8r9/!s/N+
!sA]-!s8T-rW3-'!s8H(quQm$"Sr3%!WrT(!sA`/!sAc)!t>A8"9eo1!sA`1"9o#4rWN<,"oSZ)#6Y;<
#7(V6#n%(H$47.J$OR7K$k3(B!t5VH%K6bK%1E[T$k*LP$OdCN$k!FM$O[:>$3UYB#mq";#lY&6"pG/7
#6b57q?6p'!sSZ,!<WK("9&E,!WiH*!s8E'quR!'!W`B)qZ?j#!rN0!!WW?$!rW6"!Xo/5!WiH*!sA]/
!WrN&!s/N*!<*'&"9AT+r<",E!sA`/"9\f/!X/Z,!X&Z/"9eo1"9nu4"U5/7"U>51"p5)9#71;6$jd4G
#mptG$4$qE$4?b?":>PF$4HhA!XfAEq[3H6$Ock>":GPE#n$V;!=0#2#QtA;#RCY5#6=u4rWNK1"U"u3
"pOu/r<*!$!!<-%#Qb)/!WrK+!X&E'"TSW+!s/W)!rrE%!s&N+r<!$&r<!$&rW3'%quR0,!WrN+!X&Q,
!rW39!WrK)!s8T+!s8Z.!WrQ-!sA]-!sJT)!X&T.r<!<.!sA`1"9Sf1rWNf8"U,)7"U5/6"U5,7"pYG4
#65)/#S7:J#mgqF#mq%I$4@.J$4@4>$O$kD#mq%Jq[*Z<$OI.H#RLbB#m^G8!so8=#5\Z)"pG,5"9eo3
rWET4!sA`1"9Sc1"9Sf+!rrK'!ZV=F!sA].!<WB(!s8T,"9J]-!WrK+!<WB(!<WK'!t#)2!s8T,!WrK*
"TAH-!s/N+!<N<)rW33)!WrN+rW<<,!WiH*!rW3*!WrN-!s8](!sA`/"9\o,"9J`0"T/E*#6b56r<E3*
!!`6)rWs)@#mgqF#mq"G#RUkF$4.+:$3:M@rX/T1"pYG?#m^eCrWiE.#6kA;"pY>:#Q4i3"U5)5"U"r3
r<*6-"9S`/r<!3+!s/N+!r`9*!<N?)"9\W)!WrN+rW<H-!sA]-"9AT,!rN-1!WiH+!WrK*!s8T-!sAc+
!rrK&"9&?2!s/T,!<WB*"9S`-!s/B&rW<-'rWE-$!X&W.rWE<*"9\i*"9Sf1"U+f.'F"^E"pG58"U52:
"pbG>#6tP@#7(26$OR1H$4."G$4$nE$Mt,2$4$q9#lk>.$2b#0#R:V5"qV"D#6k>9#6k>9"Tno2"9er-
"U>/5!sA]."9SK'#m(21"9S`-!s8])!t555"9Sc-!WrK+!s/Q%!u:q?!s/Q,"9JZ,!WrN-!WrK+!WrN+
!r`<$!WW9#!ser0!s/N,!s/Q!!s&N+rW<<*!X&Q-"TAK'"TAK/"Tnl0"9\l2"T/E/"pP58"pG59#6"i.
#R:85r<NT4#R:V?#RL>7&dSaL#m^hD#RLhE#RLbC#RCYArWiT2#R:S2#71P<#R:M:"pG,4r<3W5"U"r3
!sAc1"9Sc/"o\T."9S`/!s/Q'"9&?1!s/Q,!X&T-"9Sc/!r`9)"9J].!W)p%!s/N,!Wr<&rW3'%r;ls"
rW<<*!WrK+"T/<("9S`/"TAK'!r`9%!r`91!WiK,"9\f0!sAZ-"9&E*"Tno2rWN<*"TAQ)#Q=l,"T/K'
#64u-#64u-#Qb5:#Q=u/#Q>#,#poub#71YA$4-tE#RCYA#RC_B#6tM>#R:S="pY>9"pY;8"p>)5#6"c.
"9er5r<**(q?$g&!sA`*"TAK+!WrK+!r`9,!WrN+!s8Z-r<!*&!sAN(rW<B.!s/T,!s/Q&!tbS9!<N?*
"9S]+!WiH)!<N<)r;mK3!WiH,!X&T-!WrK,!WrQ'"TAK*!WrN-qu[-+!sA]/!sJT,!!N-%%0ct="U529
"9o&8#6b>2#8.:K#7(S?#RCY@#R:VA$4-tE$2t2<#71YB#7(V@#mLY=r<EN2"pG/7"pOr0!X/`3rWW9(
"9er5"9SN)!WrQ-rWMuu"9J].!Wr9%!WiH+quR*(!WrN-!<WE&!s&H)r<!*&"9SN(rW39+!WiH*"9\W)
"p>#/!s/Q,qZ7!*!s/N+!s/T(!Y,;8!X&T.!s&K,"9S`-"o\W1!sA`/!sJf1"o\W*"U4l."9er3"U=r0
rWW?+rWWo="pP;:#R:V@#71Y@#71V?rWiK1rWrK.!XK&9rW`Q3"pY>;r<<H0"pP58"U=o/!!`B*!X&W0
rWEf:!s8Z0"9Sc."9Sc.!WrQ-!r`9&!Wr?'r;mc;!<N<)!WrN-!<<0%!WrN+!<WB*!WrT)!rrB#!WiH+
!rW33!WiK,!WrK+"9J].!WiH*rW<B,!WrN-!s8W'!rrDt!rrK$!s&Q.qu[-+"U"o0"U+]+"U5,6"pG50
"Tnu6#7(54r<Ef:#6kD<"pP>="pYA;#7184":#/9"pk25'a=mI#6k>9"pG28"U,)6"p>&4#7(M4"98W'
"98W)"98N&!so)5!s/Q-!WiH+J,
~>
gr
4 w
DO
SO
6 w
0 sg
 899 4616 mt 6261 4616 L
 899  389 mt 6261  389 L
6261  389 mt 6261 4616 L
 899  389 mt  899 4616 L
 899 4616 mt 6261 4616 L
 899  389 mt  899 4616 L
1387 4616 mt 1387 4562 L
1387  389 mt 1387  442 L
%%IncludeResource: font Helvetica
/Helvetica /WindowsLatin1Encoding 120 FMSR

1321 4761 mt 
(20) s
1888 4616 mt 1888 4562 L
1888  389 mt 1888  442 L
1822 4761 mt 
(40) s
2389 4616 mt 2389 4562 L
2389  389 mt 2389  442 L
2323 4761 mt 
(60) s
2890 4616 mt 2890 4562 L
2890  389 mt 2890  442 L
2824 4761 mt 
(80) s
3392 4616 mt 3392 4562 L
3392  389 mt 3392  442 L
3292 4761 mt 
(100) s
3893 4616 mt 3893 4562 L
3893  389 mt 3893  442 L
3793 4761 mt 
(120) s
4394 4616 mt 4394 4562 L
4394  389 mt 4394  442 L
4294 4761 mt 
(140) s
4895 4616 mt 4895 4562 L
4895  389 mt 4895  442 L
4795 4761 mt 
(160) s
5396 4616 mt 5396 4562 L
5396  389 mt 5396  442 L
5296 4761 mt 
(180) s
5897 4616 mt 5897 4562 L
5897  389 mt 5897  442 L
5797 4761 mt 
(200) s
 899  774 mt  952  774 L
6261  774 mt 6207  774 L
 731  818 mt 
(20) s
 899 1169 mt  952 1169 L
6261 1169 mt 6207 1169 L
 731 1213 mt 
(40) s
 899 1564 mt  952 1564 L
6261 1564 mt 6207 1564 L
 731 1608 mt 
(60) s
 899 1959 mt  952 1959 L
6261 1959 mt 6207 1959 L
 731 2003 mt 
(80) s
 899 2354 mt  952 2354 L
6261 2354 mt 6207 2354 L
 664 2398 mt 
(100) s
 899 2749 mt  952 2749 L
6261 2749 mt 6207 2749 L
 664 2793 mt 
(120) s
 899 3144 mt  952 3144 L
6261 3144 mt 6207 3144 L
 664 3188 mt 
(140) s
 899 3539 mt  952 3539 L
6261 3539 mt 6207 3539 L
 664 3583 mt 
(160) s
 899 3934 mt  952 3934 L
6261 3934 mt 6207 3934 L
 664 3978 mt 
(180) s
 899 4329 mt  952 4329 L
6261 4329 mt 6207 4329 L
 664 4373 mt 
(200) s
 899 4616 mt 6261 4616 L
 899  389 mt 6261  389 L
6261  389 mt 6261 4616 L
 899  389 mt  899 4616 L

end

eplot
%%EndObject

epage
end

showpage

%%Trailer
%%EOF