i2.eps

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

&eb0UrXf&>rXSo:rXA`5rX/T1r<W?,r<<6*r!**)r<WB/r<iN3rXA`7rXSo<rX])B&cNFB'E/[H'c$]`
rYPPN!>udS)ZC]T*<$uX*r[8\+T<P`,5ibb,lAtc-M]%U./>7f-N#1g,lJte,5rbb+T<J^*r[2Z*<$oV
)?C]nrYYVNrY>ML(&epH'E/UF&J>!RrX\u<r=/`7rX8Z3rX&N/r<N6)!<iZ*#5nf*#lP)/$N1A3%/p_8
%fR"<&H3:@')iRD'EJp`rYGJLrYPYR)?(QR)u^iV*W@,Z+9!D]+oW\a,Q&hb-2T%K-i57g-2]%e,Q8nc
+oWV`+9!>\*W@&X)u^`V)&`DjrYPPL!>Q@K'`JdF')iLB&H34>%fQq:%/gS5$N1;0#lP#,#6"f(#5nf*
#lP)/$N1A3%/p_8%fR"<&H3:@')iRD'EJp`rYGJLrYPYR)?(QR)u^iV*W@,Z+8m>]+oNV_,Q/nb-28hT
-ho%c-2f+e,Q/hc+oNP_+9!>\*W@&X)u^`V)&`DjrYPPL!>Q@K'`JdF')iLB&H34>%fQq:%/gS5$N1;0
#lP#+#5n`(#5e`*#lP).$N:G3%/p_8%fR"<&H3:@')iOF'GUN]rYGJLrYYVP!?)mU)u^iU*W@,Z+8m>]
+oNV_,Q&h`-0-EL-2]%d,Q/hc+oNP_+9!>[*W@&X)Z^iqrYb\PrYPPLrY5GJ'`JdF')iLB&H34>%fQq9
%/pY5$N1;1#lFr+#5eZ&#5nf*#lP)/$N1A3%/p_7%fR"<&H3:@&ciXZrY5>HrYGJLrYYVPrYkbTrZ(nX
rZ;"[rZM._r$)%`q]tn^ng=)Wq]u(ar$)(_rZM.]rZ;%ZrZ(nVrYkbRrYYVNrYGJJrY,AH'E/XD&cN@@
&,d";%K6e7$iUM3$2k/.#Q+f'"oSZ'#Q4r-$2k50$iUS6%K6k:&,d(=&HNLWrY,8FrY>DJrYPPNrYb\R
rYthVrZ1qYr?)"]r#u"_q]kn^kU$$Mq]l%`r#u%^r?(tZrZ1tXrYthTrYb\PrYPPLrY>DHrY#;F')iLA
&H34>%fQq:%/gS4$N:A1#lFr+#5JH$#5nf*#lP).$N:G3%/p_8%fR"<&H3:@')iRD'`JjH(B,-L)#bEP
)ZC]S*<$uW*r[8[+T*D\,5NPW,kr\Y,l/b`,5`V_+T<J]*r[2Y*<$oV)ZCWR)#b?N(B,'J'`JdF')iLB
&H34>%fQq9%/pY5$N1;0#lP#,#5AB"#5nf*#lP)/$N1A3%/gY7%fR"<&H3:@')iRD'`JjH(B,-L)#bEO
)ZC]T*;poW*rR2Y+T*D[,5EJK,l&\^,5`V^+T3D]*rR,Y*<$oU)ZCWR)#b?N(B,'J'`JdF')iLB&H34>
%fHk9%/gS5$N1;0#lP#+#58<"#5e`)#lP)/$N1A3%/gY7%fR"<&H3:@')iRD'`JjH(B,-K)#bEP)ZC]S
*<$uW*rI,Y+T!>Y,2sjD,5WP^+T*>[*r[2Y*<$oV)ZCWQ)#b?N(B,'J'`JdF')iLB&H34>%fHk9%/gS5
$N1;0#lFr+#5eW)"T/H$#5nf*#lP).$N1A3%/p_7%fR"<&H3:@')iRD'`JjG(B,-L)#bEP)Z:WR*<$uW
*rI,X+T!>W,4$QL,5WP]+T*>[*r[2Y*;piU)ZCWR)#b?M(B,'J'`JdF')iLB&H34=%fQq:%/gS4$N1;0
#lP#+#5eW*"9eW+r!**)r<W?.r<iN3r=&Z7rXSl;rXf&@rY#2Dr=o5GrYGJLrYYSOrYk_SrZ(kWr#YeY
q]P\Xj<FCCq]PhZr#YeWrZ(kUrYk_QrYYVNrYGJJr=o5ErY#2BrXf#=rXSo:r=&Z5r<iK0r<W?,r!*'&
rWN*%r<E0)r<W?.rX/T3r=&Z7rXSl;rXf&@rY#2Dr=o5GrYGJLr>>MOr>PYSr>bbVr#YbXqB5MUn07TM
qB5\Xr#YeWr>beUr>PYQr>>MMrYGJJr=o5ErY#2BrXf#=rXSo:r=&W4rX/T1r<W<+r<E-&r<3$%r!**)
r<W?.r<iK2rXA`7rXSl;rXf&@rY#/CrY5>HrYGGKrYYSOrYk_Sr>b_Ur#Y_Wq&nfCq&oPVr#YbVr>bbT
rYk_QrYYSMrYGJJrY5;ErY#2BrXf#=rXSl9rXA`5r<iK0r<W?,r!*'&qZQg#r!**)r<W?.r<iK2rXA`7
rXSl;rXf&@rY#/CrY5>Hr>,AKr>>MOr>PVRr>b_Uq]>SUocWrOocX)Qq]>YUr>bbTr>PYQr>>MMr>,AI
rY5;ErY#2BrXf#=rXSl9rXA`5r<iK0r<W?,r!*'&q?6[!r<E0)r<W?.r<iK2rXA`7rXSl;rXf&@r=])C
rY5;GrYGGKrYYSOr>PVRr#GVTq]>MSkTKXBq]>YUr#GYSr>PVPrYYSMrYGGIrY5>Fr=])ArXf#=rXSl9
rXA`5r<iK0r<W<+r<E-&q#pU!r!**)r!<6-r<iN3r=&W6rXSl;rXf&@r=])Cr=o5Gr>,AKr>>JNr>PVR
r#GSSqB"K8qB#MSr#GYSr>PVPr>>MMr>,AIr=o5Er=])ArXf#=rXSl9r=&Z5r<iK0r!<6+r!*'&pB:Bt
r!**)r<W<-r<iN3r=&W6rXSl;rXf#?rY#2Dr=o5Gr>,>Jr>>JNr>PSQr#GPRp`A]Bp`B8Pr#GVRr>PVP
r>>JLr>,AIr=o5ErY#/ArXf#=rXSl9r=&Z5r<iH/r<W?,r!*'&p&t6rr<E0)r<W?.r<iK2r=&W6rXSl;
rXf#?rY#/CrY5;Gr>,AKr>>GMr>PSQq],DPi#_\7q],MQr>PSOr>>MMr>,>HrY5;ErY#/ArXf#=rXSl9
r=&W4r<iK0r<W<+r<E-&qud!%qZQg#r!*'(r<W?.r<iK2r=&W6rXSl;rXf#?rY#/Cr=o5Gr>,>Jr>>GM
r#5JPqAf5MllPmAqAfDPr#5JNr>>JLr>,AIr=o2DrY#/ArXf#=rXSl9r=&W4r<iK0r<W<+r!*'&qZHp%
"8`6!"oJT&#Q4r,$2k50$iLM4%K6k9&,d(=&cE@A'E&XD(&\pH(]>3K)>kEL)u1K5*Vg]O)uLWP)>tEN
(]>-J(&\jG'E&RC&cE:>&,m(;%K-_6$iLG2$2k/-#Q+f("oAH$"8i<!"oJT'#Q+l+$2k50$iLM4%K6k9
&,d(=&cE@@'E&XD(&\pH(]5-J)>kEK)tk9@*VLKK)uLWP)>k?M(]>-J(&\jF'E&RC&cE:>&,m(;%K-_6
$iLG2$2b)-#Q+f'"oJN$"8`5u"oJT'#Q+l+$2k50$iLM4%K6k9&,d(<&cE@A'E&XD(&SjG(]5-J)>Y9H
)rr"6)u:KN)>k?M(]5'I(&\jG'E&RB&cE:>&,m(;%K-_6$iLG2$2b)-#Q+f'"oAH""8`6!"oJT&#Q+l+
$2k50$iLM4%K-e9&,d(<&cE@@'E&XD(&SjG(],'H)>Y9D)tY-B)u:KM)>b9L(]5'I(&\jF'E&RB&cE:?
&,d":%K-_6$iLG2$2b),#Q+f("oAGu"8`6!"oJT&#Q4r+$2k50$iLM4%K-e8&,d(<&cNF@'E&XD(&SjF
(]5-H)>G-2)u(?J)>k?L(]5'I(&\jE'E/XC&cE:>&,d":%K-_6$iLG1$2k/-#Q+f("oAGt"8`5u"oJT&
#Q4r+$2k50$iLM4%K-e8&,d(<&cE@@'DrRC(&SjF(]#!F)=nd?)tP!E)>Y3J(]5'I(&SdE'E&RB&cE:>
&,d":%K-_6$iLG1$2k/-#Q+f'"oAGs"8`5u"oJT&#Q4r+$2k50$iLM4%K-e8&,d(<&c<:?'E&XC(&SjE
(]#!C)<Mk3)>Y3I(]5'H(&\jF'DrLA&cE:>&,d":%K-_6$iLG1$2k/-#Q+f'"oAGq"8`5u"oJT&#Q4r+
$2k50$iLM3%K-e8&,d(<&cE@?'DrRB(&SjE(\np;)>kEB)>P-H(]5'H(&SdD'E&RB&cE:>&,d"9%K-_6
$iLG1$2k/-#Q+f'"oAGo"8`6!"oAN%#Q4r+$2k50$iCG3%K-e8&,d(;&cE@?'DrRB(&JdC(\\d/)>>!E
(],!G(&SdD'E&RA&cE:>&,d":%K$Y5$iLG1$2k/-#Q"`'"oAH!!s&H)qZHZtr!!!&r!3-*r!E</r!WH3
r=/]8r"&`;r=Sr?r"JuBq\AuDq%qa%q%roDq\B#Cr"Ju@r=Su>r"&`9r=/]6r!WH1r!E9,r!3-(quup"
qZHg!q?-TtqZZm%r!3-*r!E</r!WH3r=/Z7r=Af;r=Sr?r"JuBqA&iBpD;s/pD<ZAqA&oBr"Ju@r=Sr=
r=Af9r=/]6r!WH1r!E9,r!3-(qZZj"q?-Ztq?-TtqZZm%r!3-*r!E9.r<rQ4r!iT7r"&`;r"8i>r"JrA
qA&fAg)'Q#qA&lAr"Ju@r"8l=r"&`9r!iT5r<rN1r!E9,r!3-(qZZj"q?-Trq?-TtqZZm%r!3-*r!E9.
r<rN3r=/Z7r=Af;r"8i>q\/i@q%`W>jqmb-q%`c@q\/l?r"8i<r=Af9r=/Z5r<rN1r!E9,r!3-(qZZj"
q?-Npq?-TtqZZm%r!3-*r!E9.r!WH3r!iT7r"&]:r"8i>q\/f?p_DKup_EW>q\/l?r"8i<r"&`9r!iT5
r!WE0r!E9,r!3-(qZZj"q?-Hnq?-TtqZZm%qZm$)r<`B/r!WE2r=/Z7r"&]:r"8f=q\/c>p(c^*p(dB;
q\/i>r"8i<r"&]8r=/Z5r!WE0r<`B-qZm$'qZZj"q?-Emq?-QsqZZm%qZm$)r!E</r!WE2r!iQ6r"&]:
q[r]<q@iW<g(jAsq@i`=q[r`;r"&]8r!iQ4r!WH1r!E9,qZm$'qZZg!q?-Blq?-QsqZZm%qZm$)r!E9.
r!WH3r!iQ6q[`T9q[rZ;q@iQ:j;%A&q@i]<q[r`;q[`T7r!iT5r!WE0r!E9,qZm$'qZZg!q?-<jq?-Qs
qZZm%qZm$)r!E9.r!WE2r!iQ6q[`T9q@WQ:p_29op_3K:q@WW:q[`T7r!iQ4r!WE0r!E9,qZm$'qZZg!
q?-6hq?-QsqZZm%qZm$)r!E9.q[<<1r!iQ6q[`Q8q@WQ:ob6C#ob707q@WT9q[`T7r!iQ4q[<</r!E9,
qZm$'qZZg!q?-Np!<E5s!r<#q"T&B"#5e`(#l=r+$N(;0%/UM3%f6e6&GZq5'(?S0'^ue2')<.:&Gm"9
%f6_5%/^M2$Mt/-#lFr)#5\T#"Sr5r!r`5q!r<#q"T&B"#5e`'#lG#,$Mt5/%/UM3%f6e5&GZq3'&XGs
')<.9&Gm"9%f6_5%/UG1$N(5-#lFr)#5\T#"Sr5q!rW/p!r<#q"T&B"#5\Z'#l=r+$N(;/%/UM2%f6e5
&GHe-'(Q_,')*"7&Gm"8%f6_4%/^M2$Mt/-#l=l(#5\T#"Sr5q!rE#o!r2rp"T&B"#5\Z'#l=r+$Mt5.
%/UM2%f-_3&G?^m'(uq5&Gcq7%f6_4%/UG1$Mt/-#l=l(#5\T#"Si/q!r2lm!r2rp"T&B"#5\Z&#lG#+
$Mt5.%/UM2%f$Y1&FpG$'(QY0&GZk6%f6_4%/UG0$N(5-#l=l(#5\T#"Si/q!qu`k!r2rp"T&B!#5e`'
#l=r*$Mt5.%/UM1%f$Y/&D@`i&GZk5%f6_4%/UG0$Mt/,#lFr(#5\T#"Si/q!qcTi!r2rp"Sr<!#5\Z&
#l=r*$Mt5.%/LG0%epS+&Ej_t&GQe4%f-Y3%/UG0$Mt/,#l=l(#5SN""Si/q!qQHg!r2rp"Sr<!#5\Z&
#l=r*$Mk/-%/LG.%egLb&GH_1%f-Y3%/LA/$Mt/,#l=l(#5SN""Si/q!q?<e!r2rp"Sr<!#5ST%#l=r*
$Mk/,%/LG.%eL:k&G-M.%f-Y2%/LA/$Mt/,#l4f'#5SN""Si/q!q-0c!r2rp"Sr;u#5\Z&#l4l)$Mk/,
%/CA+%djl&&FL)&%f$S1%/LA/$Mk)+#l=l'#5SN""Si/q!pp$a!r2ro"Sr<!#5ST%#l4l)$Mk/+%/:;)
%bqTf%epM/%/LA/$Mk)+#l4f'#5SN!"Si/q!p]m_!r2ro"Sr;u#5\Z%#l=r)$Mb)*%/15$%dX_r%egG.
%/C;-$Mt/+#l=l'#5SN!"Si/q!qlWq!!)]nq#^<nq?6Wuq?Hd$q?Zs)q$Qs+p^Hg)edh?^p^I!,q$R!*
q?Zp&q?Hd"q?6Qqq#^6jr;cNlp]C6nq?6Ttq?Hd$q[!!(q?m$+pC-X&iXYPhpC-j*q?m$)q[!$'q?Ha!
q?6Trp]C-iqZ-<jp]C6nq#pNtq?Hd$q?Zm'q$Qp*p'g3qr!r!$p'ga)q$Qp(q?Zp&q?Hd"q#pKqp]C-i
q#L*hp]C6nq#pKsq?Hd$q?Zm'q$Qj(oF0F_oF1I%q$Qp(q?Zp&q?Ha!q#pKqp]C-ipAjmfp]C6nq#pKs
q?Ha#q?Zm'p^6^&nI4UjnI5+!p^6g'q?Zm%q?Ha!q#pKqp]C-io`4[dp]C3mq#pKsq?Ha#q$?d&p^6[%
cj]RTp^6g'q$?d$q?Ha!q#pHpp]C-io)SIbp]C3mq#pKsq$-X"q$?a%pBpL"g^Nc^pBp[%q$?d$q$-Wu
q#pHpp]C-inGr7`pB(-mp]UBrq$-X"q$?^$p'U3qmL8Lkp'UO#q$?d$q$-Wup]UBppB($hmf<%^pB(*l
q#pHrq$-X"p^$U#oEs"SoEt=!p^$[#q$-Ttq#pHppB($hm/Zh\pB(*lp]UBrp]gO!pB^I!nd=4]nd>'s
pB^R"p]gNtp]U?opB($hlN$VZpB('kq#pHrp]gKupB^EuljE:kljEClpB^O!p]gKsq#pEopB($hklCAW
pB(*lp]U?qp]gKup'C9scjKCOp'CEup]gKsp]U?opB(!gk5b/UpB('kp]U?qp]gHtoa('oh?rcZoa(9s
p]gKsp]U<npB(!go)O:
~>
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
1785 4616 mt 1785 4570 L
1785  389 mt 1785  434 L
%%IncludeResource: font Helvetica
/Helvetica /WindowsLatin1Encoding 120 FMSR

1719 4761 mt 
(50) s
2680 4616 mt 2680 4570 L
2680  389 mt 2680  434 L
2580 4761 mt 
(100) s
3575 4616 mt 3575 4570 L
3575  389 mt 3575  434 L
3475 4761 mt 
(150) s
4471 4616 mt 4471 4570 L
4471  389 mt 4471  434 L
4371 4761 mt 
(200) s
5366 4616 mt 5366 4570 L
5366  389 mt 5366  434 L
5266 4761 mt 
(250) s
 899 1206 mt  944 1206 L
5483 1206 mt 5437 1206 L
 731 1250 mt 
(50) s
 899 2031 mt  944 2031 L
5483 2031 mt 5437 2031 L
 664 2075 mt 
(100) s
 899 2857 mt  944 2857 L
5483 2857 mt 5437 2857 L
 664 2901 mt 
(150) s
 899 3683 mt  944 3683 L
5483 3683 mt 5437 3683 L
 664 3727 mt 
(200) s
 899 4508 mt  944 4508 L
5483 4508 mt 5437 4508 L
 664 4552 mt 
(250) 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 4342 mt 6111 4342 L
5751 4342 mt 5793 4342 L
6188 4386 mt 
(10) s
6154 3903 mt 6111 3903 L
5751 3903 mt 5793 3903 L
6188 3947 mt 
(20) s
6154 3463 mt 6111 3463 L
5751 3463 mt 5793 3463 L
6188 3507 mt 
(30) s
6154 3024 mt 6111 3024 L
5751 3024 mt 5793 3024 L
6188 3068 mt 
(40) s
6154 2585 mt 6111 2585 L
5751 2585 mt 5793 2585 L
6188 2629 mt 
(50) s
6154 2145 mt 6111 2145 L
5751 2145 mt 5793 2145 L
6188 2189 mt 
(60) s
6154 1706 mt 6111 1706 L
5751 1706 mt 5793 1706 L
6188 1750 mt 
(70) s
6154 1267 mt 6111 1267 L
5751 1267 mt 5793 1267 L
6188 1311 mt 
(80) s
6154  828 mt 6111  828 L
5751  828 mt 5793  828 L
6188  872 mt 
(90) 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