fig32.eps

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

%1EUO#m^b>"UkbF$4$nA"U"o.!W305#RUtM&/,fs*$6@N-n6`!r[na6.k2qn+s.dB()%5h$OR.FrWcC.
!!33)"U>8=#mq"J"9o);#mq%K%1WUK#mq(M%M&gO$OdLU&J>K]%hK<`%MTQg&.T*V#nRUW$k!CJ#7h7P
$OR1G#6kVG$47(G#6k;6!W`Q4rW`u?#mq(O&eu6'+!N!Z.kND.r\4p;0.e\&,pFEM(_mYo$jm==#>tj/
!s8]2#RCbE$O[CI#7(VC$OdIS#7(VC$k3[W#RUqK%M03c$k3[X&.]3b'+k`b$k!@Q%h9!U$4$kI%1EUO
#m^b>%1<ON$4$k?"Tnf+#lOu>#RLkK&/5p"*[2mY.kWM22#]?D1bpa;.ju\f*Z>h.&.K!R#lP!+#QXu0
"9o);#mq"J%0d"B#mq%K%1WUK#n%.N%M&gP$OmRW&JGQ^%hK<`%MTQg&.T*V#nRUW$k!CK#7h7P$OR1G
"pPMF$47(F"pP24!<EK)#7q.J$k<mb(`XYD-n?o)1Gq,j3"S#+1GLL4-mToT(_mVn$jm:=#>tg-!WrT0
#7(YC$O[CI#7(VC$OdIS#7(VD$k3[W#RUqK%M03c$k<d[&.f<d'+kcc$k!@Q%h9!U$4-qJ%L`^P#m^b>
%1<OM$4$k?"Tnf4#lOu>#RV"P',DH,+sn`j0JbUH3rV2P3B/oQ0J+b$+rqR;'+bT[$2k*,#RLM6"9o&9
#mq"J$jHnA#mq%K%1`[L#n%.N%M/mQ$OmRW&eb]a&.oKb%M]Zi&.T*V&eGQ`$k!CK#7h7Q$OR1G"pPJE
$47(F"pP24!=9&1#7q1K%1j3i)]p:P.k`Y83''/#4V9b:3&NHE.juVa)Aa%u%1<IJrWc=,#m(/2"U>8>
#n%+M"U>;>$4@7P%gN@H$OdIS%hT-V%1Wp^'FYH^&J5Qa'b_/j%L`^X&.]6Z$OR1F%h9!U$4-tC"UthI
#mghA"U"l.#RLD7&I8[N%1s9k*$?LU/2/k=3BKD'4qKk=3]AlL/L_tg*#KA$%LWU@#>b^-!sAf4#RLhF
$k*=E#RLhG%1NgQ#RLkI%1Wp]$4@:Q&/#Wc%M03a%h9<g&J,E\$PF!]%L`^P#mUtM%1<LL#R1JA$jm=J
#R1J9!s&`+#7q1L%M9Bl*$?OV/ho4C4$5_+57p(B4#f)P/h/.j*#KA%%LWRKrWc=,#Qb&1"U>8=#mq%L
"U>;>$4@7P%gN@I$OdIS%hT-V%1Wp^'FbQa&J5Qa'b_/j%L`^X&.]6Z$OR1G%h9!U$4-tC"UthH#mghA
"U"l.#RLD7&I8[N%M9Bl*$HUX/ho4C4$>e,57g"A4#f)P/h/.j*#KA%%LWU@#>b^-!s8`3#RLhF$k*=E
#RLhG%1NgQ#mgqI%1Wp]$4@:Q&/#Wc%M03a%h9<g&J,E\$PF!]%L`^P#mUtM%1<LL#R1JA$jm=J#R1J9
!s&`+#7q1L%M9Bl*$?OV/ho4C4$5_+57p(B4#f)P/h/.j*#KA%%LWRKrWc=,#Qb&1"U>8=#mq%L"U>;>
$4@7P%gN@I$OdIS%hT-V%1Wp^'FbQa&J5Qa'b_/j%L`^X&.]6Z$OR1F%h9!U$4-tC"UthI#mgkB"U"l.
#mgM8&I8[N%1s9k*$?LU/2/k=3BKD'4qKk=3]AlL/L_qf*#KA$%LWU@#>kd.!sAf4#RLhF$k*=E#RLhG
%1NgQ#mgqI%1Wp]$4@:Q&/#Wc%M03a%h9<g&J,E\$PF!]%L`^P#mUtM%1<LL#R1JA$jm=J#R1J9!s&]6
r<FVT$4RO\(E4G@-S$f)2)mWY4?Yeg3]JuP0.SFq*uc"1&.K!R#lY',#RLP6"9f#8$4.(J%0ZtB#mq%K
%M&dM$4@7O%M'-Y$OmRW&ek`a%hT?a%M]Wh&.T*V#nRUW$k!CK#7h:Q$OR1G#6bPE$OI.G#6b85!XT/3
#8%4J$4RO[()\/:-7LN#1H%3Qr]:ZO3&WQH/Li(k*Z>e-&.JsQr<H4+#lt,2"pP>>$471M"pP>>$4@7O
%gN@H$OdIS%gWLM%1Wp^'+GE_&J5N`'GD&i%L`^O&.]6Z$OR1F%h/pT$4-tB"q1nI$4-tB"p>#0!=8W/
&-rUO%h]To*$?LT.kWP42)kLo&2uI&0eOt),Tn'D'bLrb$4-S8;@!7(!sJi5#RUkG$Om7E#RLhG$k3[O
#RLkI%1WmS$4@:Q&.oQa%M'-_%h06e&J#?[$OI@T%L`^P#R:kL$k!CK#R:MB$O[:K#R:M:!s/H/r<El>
#n%4S'G_T.+seWf/hf%:r\Y6C1,(=2-m^#X)Aa%u%1<IJr<H7,!<NB,"pYG@$471N"pP>>$4@7O%L37G
$OdIS%gWLM%1Wp]'+><\&.oE_'GD&i%L`^O&.]6Z$OR.E%h/pT$4-tC"q:tJ$4."D"p>&1!=/u1#9sK\
$OdOY'c7o4+seWf/MAe41,1I9/Lr4q+s.dA'bV#c$46V8;[EI,!sJl7#RUnH$Om7E#RLhG$k3[O#RLhH
%1WmS$4@:Q&.oQa%1Wp\%h03c&J#?[$OI@T%L`^P#R:hK$k!CK#R:M:$k!@K#RCV<!s8Q)#lOu=#mq(M
&/,fs*$6@N-n-Ytr[na5.k)km+s.dB()%5h$OR.FrWcC/!!36*"U>8=#mq"J"9o):#mq%K%1WUK#mq(M
%M&gO$OdLU&J>K]%hK9^%29Hf%h9!U#nRUW$k!CJ#7h7P$OR1G#6k;>$OR1H#6tD8!s/H/rWaYR#RLhH
%M9Bk)&s_B,UOlf.4Qbq-RL)_*ul.6'G1f`$OI+<#?D*9!<E9*"pP>=$4.(K"U52;#mq%K%1WUK#mq(M
%M&gP$OdIT&/#?[%M0-\$k*a\%h9!U#nRUW$k!CJ#6kVG$OR1G#R1G@$OI.G#mUY;"9JW1#lY&>#RCbF
%1a'c(E"26+<_sTrZr%",9\0N*#TM+&J#<Y#mgJ7<sT*5!s8]2#7(VB$4@:H"pYG@$4@7O"pYD?$4I@R
%gWIL$k3[X#n%+M%M0*Z$PF!]%1ERM#S.@Q$OR.E"q:tK$4."D#6bPE$4."E#R1J8!s&E.#lY&O#RCbF
%M'-c()Rr0*ZlLI+WqmL*ZQ(8().>k%1<LL#lP!/#QOo-"9f#8#mgqG$k*:D#RC_C$OdIK#7(VD$k3[O
#RLhH%1`sU$4I@R%h9!U&eGQ`$k!@I%h9!U$4-qA%1EUO#m^e@"p>>A$4."D#6b85!W`W,#7_"E#RCbF
%M'-b'c%W()uLZa)]9J/'b_/i%1ERM#lP!1#RCJ5"9\o5#RCbD$O[(@#6tPA$4@7O"pYD?$4I@R#7(VD
$OmRW#RUqJ%1`mV$5*m\%1ERM#S.CR$OR1F"q:tK$4."D#6b8=$jd7H#mU\<"Te`*$2t/=#R:VA$4@:Q
&/#Zl()Qlc$PsQn'+tid%1<LL#lFp1#QXr."9eu7#RLeE$4I%A"pbJ@#n%.N"pYD?$OdIS#7(VC$k*UV
#RCeG%1WgT$P<p[%1<LL#S.CR$OR1F"pGDC$4-tD"pP25$4@.I#RCV>"9ei.!"&W0$O?tC#RUqJ%1a!_
'D`CL'+tie%1EUO#mgG6>6k30!s8]1#6tP@$4.(B"U,,:#mgqH$j?h?#RUqJ%1WUK#mq%K%1`[L#n%.O
$k!@I%h9$W$4-tC%LidQ#m^_=$k!CK#m^e?"p5;A$4$qC#R1J9"9JW2$2t/;#R:V@$471M%M'*^r"05H
&.T-X$jm=I#lFp5#RLM5!sJi5#7(SB#n%+E":#/;#RLeF$j?h?#mq%K%0d%C#mq%K%0Zq@#mgtK$jm:H
%h9$V$4-tB%L`aQ#m^_=$k!CK#RCY?"U,8B#mpqE#7(M<"9\`,#mpS9#RCY@#RLhG$Om"DrXSo:"q1kH
$4-tDr!-R7#lt)0"9eu6#RC_C$474F"pP>=#mgqH$j?h@#mq%K%0d%D#mq"J$j?e>#71\F$OR.F%h9!U
$4-tC%LidQ#m^b>":GPC#m^e@#6Y2;$OI.G#mLY="p>&0!XT24#6P5:#RLD9rX/Q2rX&Z4$2t5.#@RlE
!<N?+"U>5;#71YC#n-n?"U>;>#mq"I"9eu8#mq%K%L*1E#mq%J%0Qk>#6tSE$OI(D%h/pT$4-tC"q:tJ
$4$k@"UkbG#RL_A"pP25$47%F#71S?"U5)4!<NQ6rW`N2#7(S6#l+f*#l4d9#RLM5!sAc4"pYA=#RLeF
!X/`4#7(SA#n%.E"pYD?$O[CI#6tPA$OR7M"9\o5#RCeF#mUY<%1EUO$4$k@%1EUO#mUY<":GPC#RCY@
"U5&;#mpqE#7(M>"U5)4!<NQ6p]h!/#7(P?#71#.B*\eE!WrQ/"pP;;#RC_C$OR%?#6kG>#mgqH"9eu7
#RUnI%0d"B#mgqH$NpS9"U52<$4-qA"pGGE$OR1H#R1bI$OR1G#6b8=$OI+F#R:S;"Tf)=$4$qC#R1P<
"p>&1!WWN3#j_e3#R:V8!<N?+"U5/9#7(SA#mq"B!sSo7#RLeE$4I%A"pbMB$OdIK#7(VB$OR:E"9\l4
#RCbC#6kA9%1EUO$4$nB"UtkJ#m^b?"U#/>#mgkC"pY87$Od:J#RL_A"pbA:"9\c-!!ro:mK[Y9#RLJ3
!WrQ/"pP;;#RC_C$471D"pG8;#RC_C$NgP:#RC_D$k3@E#RLeE$4@4E"9eo4#71\B#6kA9":PYG$4-tD
"qD%L$4-tB"U#5@$4$qC#R1M:"UYVC$3pkB#R1P<"p>&1!W`<-#mUb.#7CeB#m^P4!WiK-"U=r21'e-l
#mq"J!sJi5#7(SA#n%+E"U>;>$4@7O"pYD>$47.K!s8Z/"U55<$4$k@"U,#;rX1=b#RCY=$k!CK#R:P;
!t,GC#RL_A"pY55$4@+H#71V@"pbA;"9ei/!<E0##RUG8!XAu9qZmQ0!!30'!X&W0"pYD3#^Qd`$4I"?
"U>8<#RLeE$3UJ:#7(VC$Od1C#RC_C$471D"9\l3#6tPA#R:P;"U"r:$OR4I#R:PC$k!CK#6tD9$k!@J
#RL_A"pP/4$4@+H#71V@"pbD;"9el0!WiB'!!*B4#RLbC#RLe<!!**%!WrN-"9o&8#Q4pS#mgqI!sJi5
#7(SA#RUnI"9o&9#RLkI%0d%C#mgqG$j-V8"U,)8#lk;<#6kA9"U#5A$OR1H#6k>?$OR1G#6tD8$O[7I
#RL_A"U5&3$4@+H#RL_B"pbD<"U+u2!WrK*!<N9'!<N9'!<N9'!WrN-"9er5"pYD3#]L(V$4?q="U>8<
#RLbD$4I%A"pbJ@$4@7G"pYD?$47.L"U##5#6tM@!X]5="pP57"9],@$OR1G#6kVF$4-tD#6tD8$O[7I
#RCY@"U5&3$4@+H#RL_B#7(M="U5)5!sA].!WrK+!WrK+!WrQ'"B>X0"pP;;#RC_B#mgtI!sAc4#7(SA
#RUnI"9eu8#71\D$Om:F#RLeE$Od.A"pP;;#RLk?!XT,;#6kA9"UtkI$4-tC"pGDB$4-tD#6k>7$OR1H
r<O2C"p4u8$O@(F#m^hB#R:S="pG/5"U##-"<.FG"9eo3"9er5"pYA=#7(SA#RLeF$NgM8"pb)4,mand
"9eu7#7(VC$Od1C#7(VC$OdIJ#6kG=#RLh>"9Jr;#R:V6"r@FP$OR1H#6tDA$OR1H#RCY>"UkbG#mpP9
2$F0d":>PD$4$tE#m^hB#R:S="pG27"pG27"pG27"pG28#6tM>#m^hC#mgqF$O[(?"pP>2#Ta9Z$Od1C
#6tM?$47.M"U52;#mq%K"U52:#RLeE$iU82!Xf;?rWa;F"q:tJ$4-tC"pGDC$4."E#R:P;$k!CK$2k,c
"U5&4$Od=K#mpqF#RLbC#7(P>"pY>;"pY>;"pY>;#7(P@#RLbC#RUkF#n%+D"9o&8rWkgs#mq%K"U52:
#RC_C$Od1C#7(VB$OdIJ#6tM@#mq%C"9o#4$OR1H#R:S="UkbG#mgkA"p>>A$4."E#6tD9$k*%@'aP'M
#R(G9"Tf)>$OI.H$4."E#m^k9#QtA<#7(S5#WE%q#m^hD$4.%G$OI4C"U##6#6tM?#mgqH$j?h?#7(SA
#n%+E"U>8<#mq%K"9o)9#RLhG$j?P6*sDcT$4."E#6tDA$OR1H#RCV=$k*IM#m^e@#6b8>$iUGK#RCY@
"pY87!sAu?$4@.I#mptG#mptGq?\5N#mptG#mptG#mq"J$3UG7"9o&8#7(SA#mq%K"U552#[I`C$O[CI
#6tPA$47.L"U52;#RUqJ"9o&9#6k>?$OR1G#R:P:$jm=J#mU_?"q1nI$4-tD#6tA8$k!CK#m^hB#R1M:
"p>);%1<RO$4.%G$4.%G$4.%G$4.%GrX(jr$4.%H$O[@N"Tnr3#6kG=#RC_C$4@7G"pYA=#RLeF$4R+C
#71YB#n%+D"U>8<#mq%K"U52;#R:P:$k!FA#Uor]$k*IM#mgkC"pGDC$4."D#R:S;"UkbG#mgkD#7(M=
"U5)5!sAu?r<ic9$OR7K$iUJ:$OR7K$OZh?4:DT#!sJf3"U>8<#7(SA#mq%K"U52:#RC_C$OR=G#6tP@
#mgqH%0ZtA#RC_D$NpV;#71YA"pPMF$46\;.0]qa$OR4I#m^b?"UkbG#RL_A#6tD9$k!CK#m^hB#R:V>
#6b;7"U#;5$k*LP$OdCO$OdCO$OdFA$igb9"F:7U#6kG=#R:YA#mgqH$j?h?#7(SA#mq"K"9o&:#RLeE
$Om4C#7(SA$471E"U55=#m^b?"UkbF#m^eA"U,;B$OR4I#R:P:$jm=I#mU_@#6b8>$O[7I#RLbB#7(M=
"pP58"U,#4%1N^T$k3US$k3US$k3XL"T8K9"pP;:#R:V@#m^hD$O[@H"pb)4"pkVD$j6b=rWrc9$4@:G
"pb,5:^II:!sJi5#RLeC#6k>7$OR1H#RCY=":PYG$47%E"pPJE#mpqE#7(M="pG,;$OI.G#m^eA#R:V>
#6kD;"pG27"pG27"pG27"pG27"pOr0!s]#8#Q4o7#m^hD$O[(@"pb)4&.&XM$k3CF#RLeE$4@7F"pb,5
2[Kg!!sA`2#7(VC#R:P<":PYF#mgkC#6bPG$Od@L#RCS<$Od=K#RLbB#6tG:":GSE$4."E$2k,D#7(P?
#6tG<"pY>;"pY>;"pY>;#5ncR#R:V@#m^hC#mgqF$O[(@"pYD>#m^hD$4I@I"pbMA$OR:N%L*1;#V$,e
$Od.@"U52;#mptD#6b86$k!CK#m^b@"q;%L$jd7G#R(D?$jm=J$2t2D#6k>8$k*IM#mgkD#RLbC#RLbC
mK`t&#RLbC#lP&e$471L"9eu7#RC_B#mgtI%0ZqA#RLkG$k*UL"pYG?#mgtI"9\l3#6tMA$4-tC#6b86
$jut?+:&2V"q;%L$jd7G#R(D@$k!CL#mgkC#6tD9!t5ME$4>Vr=UGHC$4I"?"U52;#RLbC#n%+M"U52;
$4@4N$k<FG#71YB#n%+D"9er5"pbMB$4-tD"pY;8$k*IM#mgkC"pYSI$4I4J#7(J;":PYF$4."E#R:S<
"U"u-$iUJ9#mptG$1&!##mq"<$=!rH$j6\:"U>8<#RLbC#n.4G"U>8=$4@7O%1`UI#7(SA#mq%L"9o#6
"pYG@$O[:J#R:P="U,;B#mptF#RCV>%1W^R#mpnC"pP5>$jm=J#m^eA#6kA9"U#7j$RH2`"U,)7#6tM?
#m^kE$k*:C#6tPB$OdIR%gE:E#RC_C$474Fr<<f=#RLk?!t,GB#6tG<"U#5Ar<X)B#R1MC%LNUN$4$nB
"pGDD$4?b=":,5;"pY&1!!Mj$!Y#PIndGL$rWP7c"pbG?#RLeE$4I=H"U>8=$4@7O%M&aL#7(SA#mq"J
!sJl6#7(SA$N^G6$k!CK#RCV="U#5@r<X)B#R1MC%LW[O$4$nB"pGDD$OZk>rWiT2#6kD2"o&N4%LimW
%LimW%epM2"U##7#7(V7#nRFN$OdFJ"UGA?$471M%M&gC#Ta9Y$47.C"U52:#RC_E$NpV:":PYE#RCY?
"U,8B$N18C#RCV>%1W^R$47(F#6kA8$k*"@!so8>#lY)*#5o#4%hAFH!=fY=%f$;/#Q>#.#o!^R$O[@G
"pP>=#mq"J$k<aP#7(V7#p'E\$Od.B"pYD>#mq%C"U5/7":PYF#m^b@"9f/@$N18C#RCV>%1W^R$47(F
#R1J:%1N4D!=B25$2t4h#6>)9rX&Q2'+50V%0Qk>#R:YB$4@4O%0m+D#lY/0$6BTX"U52:#RLhG"9\o6
#6k>7$jm=I#mLV;"Ub\F$N:>F#RCV>%1NXP$4."E#R:P;"UtkKrX/]5r<`'$!XT,=o*G:$(CCKW$OdFQ
"U,)7#R:YB$4@4N%0m(C#lY/0$6T`Z"U,,9#RLhG!sAf5#RCY>"UkeH$3pk@#6P,:$OR7?#n7.G"pPMG
$OR4I#m^eArWWB,%K6b:$Ock>p'Cj/#mptGpB^j-$2k23$OcqB!=]53"q:eC#RLeE$4@:P"pP>=rX&Q2
+UeGV"U>5<#RUnJ"9o&9#RL_@"U,8B#n$tE"U>,5$4@4>#n.(F"pPMF$OR4I#m^h9#6tG;"UtnL$Oln?
q[33.!"/Z4r<`Q5$Mt29$k!IQ%0Zq?rW`i;#RLeE$4@7O"pP>=rWsnY$4@7E"U,,8#m^nG%0Hb;#6tPA
#R:P:"Ub_E$O-n@#6P,:$OZh=#RCY>"U#5A$4?b=rWrK."U5,?%1E^($jmFP%1NIF"pYA=#lP&8$47.L
"U,)8#lP&W$4@4E"U52:#m^nG%0Hb;#6tP@#m^e@#6Y2;%1!=G#mCP9"Ub_Gr<WZ6#6tD9$k3(Ar<W`8
#R:S="U5DGiXH)%%0ct?"pb,5rX&N1"q1nD"U>8<r<Xn[$O[(@"pYA>#R^tK!sSo6"pbJ@$4-qB"U>/7
$4R7K"pkG<!sK&@$N18:#RCY="UkeI$iUM4$2t22#7(23j:;+m!=&l/#lY/0$6T]a$k*:D#7(SA#mptG
$O[(?"pP>=#n%(L!sSo6"pYD3#os9V#6b>8"pteF$O-n@#6Y2<$k!@K#mgkC#6k>7rX8f8$iUJ:#mgkD
#RL83q@<cA%hK6_%hSRJq?Qs)!=9)3$3:MBrX:^g"U52;#RLeE#n%+L"9\o5#RC_D$Om1B"U5/9#7(SB
$4-tD"pY>;"9o2B#R^hD"U>/6":P\;#m:M>#6tJ4"U5DF$k!FM$N183#RU,/q@W?*!=9&1$2t86$OdFQ
rWPIi#6tM@#mq"I$NpS9"pG8;$4.(K"9o#6"pY><#RUS5$OI(E"pY>;"9o2C#R^kD"pY;8":PYF#mgkC
#7(23"9f5C$Ocn?p'L3qp'L[)":GSG%0c\7'aG!M#mgqG$k*7B"U>5;#7:bG$O$_5"sO9W#RC_;!X&i;
#m^eA#6kD9"q(kG$O6qA#6k>8$k*IL$2k,<#6tG;"U#8B$OR7K$i122$OZ;0!=B52$3pqG$OR7L%0ZY7
!=&l.#S@@M$4I=H"U52:#R:\C$k!4ArWXVQ#7(SA$3C;3":>JB#RCY?"pY88$Od:J#RCY>"pG/5$k!CK
r!3]9#6kA8":GVG$OR7Kr<rK2oF:[+r<j#@$OR7L$k!1@"pP;;r<NZ8#mq%K!sJi4rW`W5#mh"J"U+f/
r<FJP$3C;4"p4u8$4-tD#6tG<"9f/A#mgkC#6kA9"U#5A$N:A2#S%(E"U,#3!t5PG$OdCO%/gY1%KI%7
%fHhH$k*LP$OdCO!sA`1"U>;2#S[RO$4@7E"U,)7#7(SA#n.1F"9o).#<E+t!<N?+"U5&3$47(F#R1M<
#6Y2<$jm:H#R:P<"U"r1$O[:K#mgkC#6tD9!sA]/!t>YJ%1WgW%fHq1&H*+N%LrpX%1N^S!sA]/!sJi5
#Q4rE$47.K!sA`1"pP;<#RLhG$O$Y;"pYA3#9X9Z!<N<)":,/8!sAr<#m^e@#6b;7"UkbG#m^b@"pG20
"98l0$O-nD#RCV="oAE-"U"u=%hB3_r=K&C&c*+@&JG$S#7qFW%h8aG"oAE,"U5/:#RLk:$Q'!L"U,)7
#6tPA$4@4E"9o#7"pYA=#RLP2!Z_@G"pY88!sAr<#m^eA#6kD:"To,>#m^b@"pP57rWF)I$k!CK#R:P<
"9eo3"9er5"pbG>#RM+W')`IH'GD/o'`A[_&ebla#RCY?"pP57"9eo3"9o&9#RLhG$k3+<!<iZ+"q_(G
#mq(C"9\o4#6kD<#RC_;rW4#@!sAl5#6Y55":GPC#RCY?"pY;8!t,GC#R:V5#6"c2":PYF#m^b?rWN]5
"U,)8#RC_C$47.Kr<r]9%/p\:$k<.D"q(eG$4."ErW`l9"U"u3"U52;#mq%K"TAT'"qh.G#RLhH!sJi4
"pY><#7(VC!WE-=!s8Z/"pbA;"U5)4$OR1G#R:V?#6k>7$k!FA#lOu7"pP57!t,GB#RC83$O6kA#RUkG
$OdFQ%JpV8&,QhE%1N^S$O[7J#RCV>r!!H4#RLhG!sJi4#5nf+#S%1K$j-Y;#6kD<#R:VA$N:,*!uh@G
#mLY>#6b;7"9Jr<#mgkD#6tG:"9T#=#mpS:r<F)A"Tni6$4$nB#R:V?#RCbE$k*RS%M&FJ!>#e@&-3:Q
r=9YR%LijU$jm=I#R:V?#R:V@$46h:"U4i0rWj)B#RUn@!sJl6#7(P@#RLe<!rN0""<%FK#mLY>#R1J9
"9Jo;#mptF#R:P;"9T#=$N1;1#T3jN"Tni-!=/u;#m^hC$4@7N%Ls![&.fEbp^nAT&J,Ka%h9$W$k!FK
#m^hC#m^P4!s8]1#Q>#,#nRFN!sAc3#7(VA#mgqF!<NB""TJXU
~>
gr
gr

0 sg
3348  189 mt 
(reconstruction) s
4 w
DO
SO
6 w
2877 3366 mt 4566 3366 L
2877  284 mt 4566  284 L
4566  284 mt 4566 3366 L
2877  284 mt 2877 3366 L
2877 3366 mt 4566 3366 L
2877  284 mt 2877 3366 L
3134 3366 mt 3134 3335 L
3134  284 mt 3134  314 L
3068 3511 mt 
(20) s
3398 3366 mt 3398 3335 L
3398  284 mt 3398  314 L
3332 3511 mt 
(40) s
3662 3366 mt 3662 3335 L
3662  284 mt 3662  314 L
3596 3511 mt 
(60) s
3926 3366 mt 3926 3335 L
3926  284 mt 3926  314 L
3860 3511 mt 
(80) s
4189 3366 mt 4189 3335 L
4189  284 mt 4189  314 L
4089 3511 mt 
(100) s
4453 3366 mt 4453 3335 L
4453  284 mt 4453  314 L
4353 3511 mt 
(120) s
2877  753 mt 2907  753 L
4566  753 mt 4535  753 L
2709  797 mt 
(20) s
2877 1235 mt 2907 1235 L
4566 1235 mt 4535 1235 L
2709 1279 mt 
(40) s
2877 1716 mt 2907 1716 L
4566 1716 mt 4535 1716 L
2709 1760 mt 
(60) s
2877 2198 mt 2907 2198 L
4566 2198 mt 4535 2198 L
2709 2242 mt 
(80) s
2877 2679 mt 2907 2679 L
4566 2679 mt 4535 2679 L
2642 2723 mt 
(100) s
2877 3161 mt 2907 3161 L
4566 3161 mt 4535 3161 L
2642 3205 mt 
(120) s
2877 3366 mt 4566 3366 L
2877  284 mt 4566  284 L
4566  284 mt 4566 3366 L
2877  284 mt 2877 3366 L

end %%Color Dict

eplot
%%EndObject

epage
end

showpage

%%Trailer
%%EOF