inter2.eps

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

1FaOc&IB-P3\r<6*uP_%#m<IS2DQp9-m^&ZbQ.M@r;qHN!sA`1"U5/9#RLhG$k<dZ&el)u)B9hC,UY)p
1,_0W!sf>K(`jqP0fVB%%Mg*40/rt>'-&>P3<BU)+"8s-"V)=j-nmSE"q2.^*[<$_/ho.=2)[?M1bp^8
.4$)U((^cR4#AK6)\NGW2_H?p&I1KV.NoW="["Rk+rV(0$3N@H/ghe]'bC]U!B`=u0J"Y",9S&N!;cfr
!r`:!"9\l3"pYD?#n%.N%M03c'c.`,*?ZRR.4d//2`h!B#n@Xe+!iHn3<09p(EY"X2uj6t)^I!n57%rA
,V:o<#7hUm-7q&8!!N]?',MQ,*[(D)4Wb^a(_mVl#Qm4K0Ie7c'+4n'0dn.[%0SaJ.3B<7"["Rk,9%:(
"?eXr.3]`J&I\mG4ukJT/h/4o+Wh]I!;cfq!`K1'"9\l3"pYA=#mq%K%1a!_'GVE&*$-7K-Rp]&2*!f#
"q(tW)^$IZ1HGOB%i6<80/l*$&/cZA1HYgM(*4hV2a7KV'cS>F/Mf=J!!<E3$OltC4UVT%"TgkK2)$O,
*uGRs"$8:h,8q1$!&uVY*>Jtc4#8B3)\WV`4Z>&E,9@[7%1!&#4#\rK/1;eg+<;ED!;cfr!r`:O"9\l4
"pYD?#mq(M%M'-b'G_N)*?QIO-n?o*2EF#'#7M1[*$H[^1ck^D&/QB7/iGlt%Mp370K9%<&/ZN<0/brr
$4msl+=/Hg0JbRF3&s#\3B/oR1+t.*,9@^9%0lnp1Fj^j'FY+,1b0di&de\$1+FIe'+>%.2_cd-&g%u*
%LE:=4?,/O/L`"k+W_U@aoM;>r;qHN!sA`1"U52:#RLhG$k3^Y&ebur)&jV>,:4ij0/GOK56D0#&f2N3
-n[><!=02O*[NBp49Q',)^?jh3s,g%)'C:Y1cdi)$4ddd)BBqE,UFcb-RL,a+Wh[A'bLiZ!^/Ir.O,oJ
%L,-V0..kZ&-r7q0I\.`'+>%.3&)ot,p+$>%giLA4ZPAS/h/4o+s.gD(s)dU!W<*!"&B%&"U5/9#RC_D
$OmRV&J>cn(`FA9+XAEb/MK"@49,Kj%MK[",:G0$3roTt'cePM1-,FB&/ZK:/i>`n$5"'q,UkB'3B[BH
#R_+R&ek6Y4qeP@%h/jN"$\k*0e=[q)\`be!')ea,9%=*"[+as-Qj9>$O&^R0e4Om)ANek"p-tL2DQp9
-m^&[*?,k4aT22=r;ogu!sA`1"U5/9#RLhG$k3[X&JGlq)&jS<+seTe/ho4C495Tl%MK[",:G0$3WTHq
'H85F0K20!$kjO&-ndJ@!XK8M)BU4Q/2/n>3]oS+!WE-a!<>>I4#f,R/h/1k*#B1q#!Y+'/0u>R&.&D!
1Fsgn(_I)X4Z>)G-6O6C&I\mH5<:_Z0.\P!,9\*J)Aj72!;llr!r`:O"9\l3"pYD>#mq%K%M'*`'GVE%
)]g+H,q(8r1,V'T!<if?',V`7-n[>;!!WiF)'C7X1cbXC%M]p,-n[>;!!N]?'Gqf4,:4ii/M8_20eb:7
/hAJ#,pFBK((q&^!^/Is.jZ8S&dnh+2_Z^+*#9%k"$JRt.jZ;W((gr\&-G-Y2`*3?.OQJb*ul19()4gi
qZ6WsIffcU"9er5"pYD?$4@7O%M'-a'G_N(*$-7J-7LK!1,_0V!X8uB',V`6-nR584p23)(EF_M0K)&r
#nIai+X\fr2`h!B#n.@Y()Rr0*ZlIH+<MXF)]9G+&Io-Q!^8V$0.J4g(_R5]!'2nf-Qj<A%0lnq2D?U,
+<)(.$k`gK!'E4t0eFk&,pFEO)]9G,'?:+N!W<*!"&/n$"U,)8#7(VC$OdIT&.oQj()Rr0+!DmV.4d//
2EF#'"ptkS(`aeJ/2B.I!sfAN)B^@X1,qH##nI^g+!`9e1,_-U!!<E2$4I@Rr=;%$%1<II"TgnM3&E9=
-6X?F&e"sF4#JZ?,97R4$O$Mk1bU=(*uYn,$jZt94?,2Q/h8=r+s7pG)&F#%'#t"L!W<'#!sAN*Ig$#\
#7(SA$4@7P%hK?e'c.`,*?QIN-Rp]%1H.?Y!X8uA',MW3-7^c.3rfHm&JlE2-nR574p)'"&f2K/,:=uo
0f1gL4$5\k56*o45!;"h3&WQH/L_ti*>oP%$O$Ml2)$L+*uYh'#m3@N1G1+$*uYk+$jZt94?,2R/ds6Y
-6jWS*#]Y0'GD(u!;cfr!r`9p"9\l3"pP;<#mq%K%1Wm\&el&s)&jS=+seTd/MJt>3^*TL$4[[b)^$FV
0/P[Q!sfAM(`jnM/M]7I!<if?&f)B++<r3].P*/)r%L8e/hJS&-R9lX)]08#$jZr"3]/TA-6O6C&.8[C
4#SfE-R'TL'G(ZY"9:PD2)6g9.4-8_*ZQ(8()7Jp&B4_I!W<*!"%rb""U,)8#7(VC$OdIS&.oQi()Rr0
*ZuXQ-Rp]%1H.?Y!<ic=&Jc9,,:G,u2ELmB$4dgg*[<$`0f:sT!X/l=&/,cq)B9e?+<g_.4s;'l*ZZ1;
(DIGk$O6e74?#&K.O?2V(D7/_"$S_%0.J:l*#KA#$O6b64?,2R0.\P!,U"6M)]9G,'+tie`;oc9rW<'$
9*>A$"U5/9#RC_D$OdIT&.oQi()Rr0*Zu[R-n6f&1H%9W!<i`;&/6!&+XSZl1cRW""q(tV)''kI.P<J7
3B[?E"pbSG%M03b'E/\0',(uj%h9!S#6Y&k4#\uM/1;bc)\s%q#m<IS2DHd3,TdsA&e56Q"$\n-2)6g9
.4-;`*uu:<(DRVs&J#Aj!;cfq!r`9n"9\l3"pP;<#mq"I%1Wm[&ebuq(`FA8+<r0\.P3>12E<o#":#;F
',MT0,Ub2u2*!i%"q(qT(E4G@-Rp`'1cIEW4obTc"9n]-4Ttul!WbMK3]JuP0.\Lt+WVF8&If$L!'E1r
/h/1l*>oS(%1*4?5<Ch^1+t1--6s`V*?,k4().An&.T/f!;cfq!W`B*rWG^p"U52:#RLhG$k3[W&/#Wk
(Dn&1*Zu[R-n6f&1H%6U56D)s%MBKo*?c^X/M]4F4otoq%M9Em*$6@N-n?l&1,CdG2`Nhs3^,Pa3&``P
1,1C4.OQJb*Z>h.%h&^H!'E4t0InLq*uc"0&.AgJ!BrP'1bgU6.4-;a*uu:<(D[\u&J,H^%)`/C!W<*!
"%rb""U,)7#7(SA$4@7O%M'-a'GVE$)B9e@,:+]e/MAk:3'0;*"UPSK',MT0+sncl1,_-U!<`Z9%MBKn
)BBqE,UOog.P!&%/M@#S4=i!4-mp;c+W_U@()%5g$3gSs4#\uM/1;be*>oS)%LNFC!'NA%1bgU5.4-;`
*uu:<(`!f!&J,H_$k'rQqZ6WsrW>[o"9er5"pYD>#mq%K%1Wp]&el&s)&aJ9+<i*Z.P*5.1cIHY!<`W7
%M9Em*$6CQ.P<G53'0>+"UGJG&/,cp)&jP:+!;^Mr?=BJ+<V^G*#]Y0'G1ia#mLJ34ZYJW0.\Lt+W_O<
'+bQY"TgnN3ArZI/Lr4r,U"6M)]9J.'GD&i%h9$W$c2o@!W<*!"%NIs"U,)8#7(VB$4@7O%M03b'GVE$
)B9e@+s\Nb/1rY52E<la!X/i<%h]Wq*$?IQ.4m502E<la!<`T5$OmUX&ebro()Hib4VeeJ'G:ug%L`[M
"Tn`f4#f)P/h8=r+W_R='G1c]#6Y&k4#f)Q0.eY%-6jZU*?,k5().Dp&.]6[$k!EZ!;cfq!r`<$"@E=r
"pYA=#mq"I$k3^Y&J>cm(Dn&1*ZuXP-7LGs0JYLF4$EWH"pk_M',DH*+!N!Z.kWP42E3c]561ii"pbMB
$Och?4:;K&#R:P:!WYDH3]AoO0.\S#,U"3J(DIGk$O?n;5<Lqb1bgX7.OQMe+Wh^D)&O,''+tie%LigS
$44NKqZ6WsrW>Ok"9er5"pYA=#mq%J$k3^Y&J>cm(Dn&1*ZlOM-7C>p0/5:B3BTM,":#8C&/,fr)]g.I
-7LGt0JYIC3'',a5!K#JrW3'%rW51_!<>>I4?GP^2)?s?/1N"n+s.gD(DIGk$OI">!'ND'2)?p=.k)hl
+s7pI)]9G-'GD)j%h9$W$OR1G_#X<4r;us#6j*Vr"U5/9#7(VC$O[@P%hK<c'GVB#)B0\=+X89\.4d)*
1Gq-Q4[/uO#7:nN&eu3$*$6=K-7C>p/h\q71c70M3&pdq4?,8Y2Dd3F0J>"..46Dc+<DI>()%8i$OI">
!BrS*2Dd-A/Li.q,U+?P*#fb3().Dp&.]6[$k!CK#me<Hq>pNrrW>Ii"9er5"pYD>#mq%K$k3^Y&J>`l
()Ro.*?H=I,UOoh/MAh82E3c]!!3</$4ICV',DE'*$-4H,UOlf.P!&%/hd&T4>/??/M&D$-mg5b+Wh^D
)&Eu"&.K!Q"p4li4?>DY1,(=3.46Dd+WhaF)As>,'GD&i%h9$W$OR1H#RA*EqZ6WsrWE-&5Qq;q"pYD>
#mq%K$k3[X&J>`l()Ii-*$-4G,:+]d.kWM21c@?U4[/uN"pkYH%hTHi(E"/4*ZuUM,:"Q^-R\g?4XVU(
-7'o^+s7sJ*#]Y0'G:re$jm7D!s(SJ3]AoP0eY(--mg2`+Wh^D)Aj8*'GD&i%h9$W$OR1H#R:UK!;cfq
!r`<$"@!%n"pP;<#RLhG$OdIS%hTEf'G_K&)B9b>+X/0Y-n6c#0JYIC3BKAg!<WK1#n%1Q&J>cn(Dn&0
*$$(A+!:A&4<GXd*?6":)&O,('+tid$jm7E"9JQe4?>GZ1GLO7.k)hl,9e6O*#fb4()7Mr&J,H_%1EUO
#mgkB#6qmBq>pNrrWE-&49Ylm"pYD>#mq%J$k3[W&.oQi'c.]))]]tB+X89\.4Zr%0JYIC3''/c!!39,
#71_G%M'*`',2/s()Qrer>@a8(`*u('bh;o&J,E]$OR.D"9SZh4ZbY_2)?s?/Lr;!-6scY*uu=>)&F&'
'+tlg%h9$W$OR1H#R:S=#/($5!W<*!"9&Ed"U,)7#7(SA$47.L%1Wm[&J>cm()Ii,*$$+D+s\H^.4Zr%
0JYIC2`Wr_4obWd"U52<$4@7O%M'*^&J=pRrXh1&%hB-Z$k!@I#6b53!C&_03]K#T1GLO7/1N%p,pOQU
*ZZ1<(`*o$'+tlg%h9$W$OR1H#R:S="pMX>q>pNrrWE-&3s>cl"pYD>#mq"I$k3[W&.oNg'G_K%)B0Y;
+!DgR-7:5m/1rV21Gq'M3BK>e561ff"9eu7#RLeE$Mb&k$47(F#R:P;"9JT(5!;"h3B&fP1,(@5.k2qo
,pOTV*uu=>)&O,('GD&i%h9$W$k!CK#m^b@"pP57])_^/quZj"rWG7c"U5/9#RC_C$4@7O%M'']&ebro
()Ii,)]]tB+X86Z-RgMr/MAe51Gq'L3BB5a5!K#J!<N?*"9eN(3s5Zg!s/K'5<_4m4#o8Z2)I*D0J4q-
.4?Mh,9e6P*ZQ+:(`*r%'+tlg%h9$W$O[:J#RCY?"pP57"Fp
~>
gr
4 w
DO
SO
6 w
0 sg
 899 4616 mt 5483 4616 L
 899  389 mt 5483  389 L
5483  389 mt 5483 4616 L
 899  389 mt  899 4616 L
 899 4616 mt 5483 4616 L
 899  389 mt  899 4616 L
1316 4616 mt 1316 4570 L
1316  389 mt 1316  434 L
%%IncludeResource: font Helvetica
/Helvetica /WindowsLatin1Encoding 120 FMSR

1250 4761 mt 
(20) s
1745 4616 mt 1745 4570 L
1745  389 mt 1745  434 L
1679 4761 mt 
(40) s
2173 4616 mt 2173 4570 L
2173  389 mt 2173  434 L
2107 4761 mt 
(60) s
2601 4616 mt 2601 4570 L
2601  389 mt 2601  434 L
2535 4761 mt 
(80) s
3030 4616 mt 3030 4570 L
3030  389 mt 3030  434 L
2930 4761 mt 
(100) s
3458 4616 mt 3458 4570 L
3458  389 mt 3458  434 L
3358 4761 mt 
(120) s
3887 4616 mt 3887 4570 L
3887  389 mt 3887  434 L
3787 4761 mt 
(140) s
4315 4616 mt 4315 4570 L
4315  389 mt 4315  434 L
4215 4761 mt 
(160) s
4743 4616 mt 4743 4570 L
4743  389 mt 4743  434 L
4643 4761 mt 
(180) s
5172 4616 mt 5172 4570 L
5172  389 mt 5172  434 L
5072 4761 mt 
(200) s
 899  774 mt  944  774 L
5483  774 mt 5437  774 L
 731  818 mt 
(20) s
 899 1169 mt  944 1169 L
5483 1169 mt 5437 1169 L
 731 1213 mt 
(40) s
 899 1564 mt  944 1564 L
5483 1564 mt 5437 1564 L
 731 1608 mt 
(60) s
 899 1959 mt  944 1959 L
5483 1959 mt 5437 1959 L
 731 2003 mt 
(80) s
 899 2354 mt  944 2354 L
5483 2354 mt 5437 2354 L
 664 2398 mt 
(100) s
 899 2749 mt  944 2749 L
5483 2749 mt 5437 2749 L
 664 2793 mt 
(120) s
 899 3144 mt  944 3144 L
5483 3144 mt 5437 3144 L
 664 3188 mt 
(140) s
 899 3539 mt  944 3539 L
5483 3539 mt 5437 3539 L
 664 3583 mt 
(160) s
 899 3934 mt  944 3934 L
5483 3934 mt 5437 3934 L
 664 3978 mt 
(180) s
 899 4329 mt  944 4329 L
5483 4329 mt 5437 4329 L
 664 4373 mt 
(200) s
 899 4616 mt 5483 4616 L
 899  389 mt 5483  389 L
5483  389 mt 5483 4616 L
 899  389 mt  899 4616 L
1 sg
0 4227 403 0 0 -4227 5751 4616 4 MP
PP
-403 0 0 4227 403 0 0 -4227 5751 4616 5 MP stroke
gs 5751 389 404 4228 rc
gr

mwscm
gs np 5751 389 mt 0 4227 rl 403 0 rl 0 -4227 rl cp c np
[403 0 0 4227 5751 389] cc
<< % Image dictionary
  /ImageType 1
  /Width 1 /Height 64
  /BitsPerComponent 8
  /Decode [0 255]
  /ImageMatrix [1.000000 0 0 64.000000 0 0]
  /DataSource currentfile
  /ASCII85Decode filter
  /RunLengthDecode filter
>>
image
5<h=p4?GSa3&``Q1c$mA0J>%1/1W2!-mp>f,U4KV+<MXF*#fe6(`*r&'GD)k&.]6[$k!CK#R:P;"9S]+
!.Y
~>
gr
4 w
DO
SO
6 w
0 sg
5751  389 mt 6154  389 L
5751 4616 mt 6154 4616 L
6154 4616 mt 6154  389 L
5751 4616 mt 5751  389 L
5751 4616 mt 6154 4616 L
6154 4616 mt 6154  389 L
6154 3943 mt 6111 3943 L
5751 3943 mt 5793 3943 L
6188 3987 mt 
(1) s
6154 3270 mt 6111 3270 L
5751 3270 mt 5793 3270 L
6188 3314 mt 
(2) s
6154 2597 mt 6111 2597 L
5751 2597 mt 5793 2597 L
6188 2641 mt 
(3) s
6154 1924 mt 6111 1924 L
5751 1924 mt 5793 1924 L
6188 1968 mt 
(4) s
6154 1252 mt 6111 1252 L
5751 1252 mt 5793 1252 L
6188 1296 mt 
(5) s
6154  579 mt 6111  579 L
5751  579 mt 5793  579 L
6188  623 mt 
(6) s
5751  389 mt 6154  389 L
5751 4616 mt 6154 4616 L
6154 4616 mt 6154  389 L
5751 4616 mt 5751  389 L

end

eplot
%%EndObject

epage
end

showpage

%%Trailer
%%EOF