10_01_10000.eps

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

%!PS-Adobe-3.0 EPSF-3.0
%%Creator: MATLAB, The Mathworks, Inc.
%%Title: C:\Documents and Settings\PB\Mes documents\Mes images\figures\metrpolis\nouveau\100110000.eps
%%CreationDate: 08/15/2007  12:57:22
%%DocumentNeededFonts: Helvetica
%%DocumentProcessColors: Cyan Magenta Yellow Black
%%LanguageLevel: 2
%%Pages: 1
%%BoundingBox:    73   212   543   583
%%EndComments

%%BeginProlog
% MathWorks dictionary
/MathWorks 160 dict begin
% definition operators
/bdef {bind def} bind def
/ldef {load def} bind def
/xdef {exch def} bdef
/xstore {exch store} bdef
% operator abbreviations
/c  /clip ldef
/cc /concat ldef
/cp /closepath ldef
/gr /grestore ldef
/gs /gsave ldef
/mt /moveto ldef
/np /newpath ldef
/cm /currentmatrix ldef
/sm /setmatrix ldef
/rm /rmoveto ldef
/rl /rlineto ldef
/s {show newpath} bdef
/sc {setcmykcolor} bdef
/sr /setrgbcolor ldef
/sg /setgray ldef
/w /setlinewidth ldef
/j /setlinejoin ldef
/cap /setlinecap ldef
/rc {rectclip} bdef
/rf {rectfill} bdef
% page state control
/pgsv () def
/bpage {/pgsv save def} bdef
/epage {pgsv restore} bdef
/bplot /gsave ldef
/eplot {stroke grestore} bdef
% orientation switch
/portraitMode 0 def /landscapeMode 1 def /rotateMode 2 def
% coordinate system mappings
/dpi2point 0 def
% font control
/FontSize 0 def
/FMS {/FontSize xstore findfont [FontSize 0 0 FontSize neg 0 0]
  makefont setfont} bdef
/ISOLatin1Encoding where {pop /WindowsLatin1Encoding 256 array bdef
ISOLatin1Encoding WindowsLatin1Encoding copy pop
/.notdef/.notdef/quotesinglbase/florin/quotedblbase/ellipsis/dagger
/daggerdbl/circumflex/perthousand/Scaron/guilsinglleft/OE/.notdef/.notdef
/.notdef/.notdef/quoteleft/quoteright/quotedblleft/quotedblright/bullet
/endash/emdash/tilde/trademark/scaron/guilsinglright/oe/.notdef/.notdef
/Ydieresis WindowsLatin1Encoding 128 32 getinterval astore pop}
{/WindowsLatin1Encoding StandardEncoding bdef} ifelse
/reencode {exch dup where {pop load} {pop StandardEncoding} ifelse
  exch dup 3 1 roll findfont dup length dict begin
  { 1 index /FID ne {def}{pop pop} ifelse } forall
  /Encoding exch def currentdict end definefont pop} bdef
/isroman {findfont /CharStrings get /Agrave known} bdef
/FMSR {3 1 roll 1 index dup isroman {reencode} {pop pop} ifelse
  exch FMS} bdef
/csm {1 dpi2point div -1 dpi2point div scale neg translate
 dup landscapeMode eq {pop -90 rotate}
  {rotateMode eq {90 rotate} if} ifelse} bdef
% line types: solid, dotted, dashed, dotdash
/SO { [] 0 setdash } bdef
/DO { [.5 dpi2point mul 4 dpi2point mul] 0 setdash } bdef
/DA { [6 dpi2point mul] 0 setdash } bdef
/DD { [.5 dpi2point mul 4 dpi2point mul 6 dpi2point mul 4
  dpi2point mul] 0 setdash } bdef
% macros for lines and objects
/L {lineto stroke} bdef
/MP {3 1 roll moveto 1 sub {rlineto} repeat} bdef
/AP {{rlineto} repeat} bdef
/PDlw -1 def
/W {/PDlw currentlinewidth def setlinewidth} def
/PP {closepath eofill} bdef
/DP {closepath stroke} bdef
/MR {4 -2 roll moveto dup  0 exch rlineto exch 0 rlineto
  neg 0 exch rlineto closepath} bdef
/FR {MR stroke} bdef
/PR {MR fill} bdef
/L1i {{currentfile picstr readhexstring pop} image} bdef
/tMatrix matrix def
/MakeOval {newpath tMatrix currentmatrix pop translate scale
0 0 1 0 360 arc tMatrix setmatrix} bdef
/FO {MakeOval stroke} bdef
/PO {MakeOval fill} bdef
/PD {currentlinewidth 2 div 0 360 arc fill
   PDlw -1 eq not {PDlw w /PDlw -1 def} if} def
/FA {newpath tMatrix currentmatrix pop translate scale
  0 0 1 5 -2 roll arc tMatrix setmatrix stroke} bdef
/PA {newpath tMatrix currentmatrix pop	translate 0 0 moveto scale
  0 0 1 5 -2 roll arc closepath tMatrix setmatrix fill} bdef
/FAn {newpath tMatrix currentmatrix pop translate scale
  0 0 1 5 -2 roll arcn tMatrix setmatrix stroke} bdef
/PAn {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
  0 0 1 5 -2 roll arcn closepath tMatrix setmatrix fill} bdef
/vradius 0 def /hradius 0 def /lry 0 def
/lrx 0 def /uly 0 def /ulx 0 def /rad 0 def
/MRR {/vradius xdef /hradius xdef /lry xdef /lrx xdef /uly xdef
  /ulx xdef newpath tMatrix currentmatrix pop ulx hradius add uly
  vradius add translate hradius vradius scale 0 0 1 180 270 arc 
  tMatrix setmatrix lrx hradius sub uly vradius add translate
  hradius vradius scale 0 0 1 270 360 arc tMatrix setmatrix
  lrx hradius sub lry vradius sub translate hradius vradius scale
  0 0 1 0 90 arc tMatrix setmatrix ulx hradius add lry vradius sub
  translate hradius vradius scale 0 0 1 90 180 arc tMatrix setmatrix
  closepath} bdef
/FRR {MRR stroke } bdef
/PRR {MRR fill } bdef
/MlrRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lry uly sub 2 div def
  newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
  rad rad scale 0 0 1 90 270 arc tMatrix setmatrix lrx rad sub lry rad
  sub translate rad rad scale 0 0 1 270 90 arc tMatrix setmatrix
  closepath} bdef
/FlrRR {MlrRR stroke } bdef
/PlrRR {MlrRR fill } bdef
/MtbRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lrx ulx sub 2 div def
  newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
  rad rad scale 0 0 1 180 360 arc tMatrix setmatrix lrx rad sub lry rad
  sub translate rad rad scale 0 0 1 0 180 arc tMatrix setmatrix
  closepath} bdef
/FtbRR {MtbRR stroke } bdef
/PtbRR {MtbRR fill } bdef
/stri 6 array def /dtri 6 array def
/smat 6 array def /dmat 6 array def
/tmat1 6 array def /tmat2 6 array def /dif 3 array def
/asub {/ind2 exch def /ind1 exch def dup dup
  ind1 get exch ind2 get sub exch } bdef
/tri_to_matrix {
  2 0 asub 3 1 asub 4 0 asub 5 1 asub
  dup 0 get exch 1 get 7 -1 roll astore } bdef
/compute_transform {
  dmat dtri tri_to_matrix tmat1 invertmatrix 
  smat stri tri_to_matrix tmat2 concatmatrix } bdef
/ds {stri astore pop} bdef
/dt {dtri astore pop} bdef
/db {2 copy /cols xdef /rows xdef mul dup 3 mul string
  currentfile 
  3 index 0 eq {/ASCIIHexDecode filter}
  {/ASCII85Decode filter 3 index 2 eq {/RunLengthDecode filter} if }
  ifelse exch readstring pop
  dup 0 3 index getinterval /rbmap xdef
  dup 2 index dup getinterval /gbmap xdef
  1 index dup 2 mul exch getinterval /bbmap xdef pop pop}bdef
/it {gs np dtri aload pop moveto lineto lineto cp c
  cols rows 8 compute_transform 
  rbmap gbmap bbmap true 3 colorimage gr}bdef
/il {newpath moveto lineto stroke}bdef
currentdict end def
%%EndProlog

%%BeginSetup
MathWorks begin

0 cap

end
%%EndSetup

%%Page: 1 1
%%BeginPageSetup
%%PageBoundingBox:    73   212   543   583
MathWorks begin
bpage
%%EndPageSetup

%%BeginObject: obj1
bplot

/dpi2point 12 def
portraitMode 0216 7344 csm

  664   341  5647  4448 rc
85 dict begin %Colortable dictionary
/c0 { 0 0 0 sr} bdef
/c1 { 1 1 1 sr} bdef
/c2 { 1 0 0 sr} bdef
/c3 { 0 1 0 sr} bdef
/c4 { 0 0 1 sr} bdef
/c5 { 1 1 0 sr} bdef
/c6 { 1 0 1 sr} bdef
/c7 { 0 1 1 sr} bdef
c0
1 j
1 sg
   0    0 6919 5187 rf
6 w
0 -4227 5362 0 0 4227 899 389 4 MP
PP
-5362 0 0 -4227 5362 0 0 4227 899 389 5 MP stroke
gs 899 389 5363 4228 rc
gr

/mwscm { [/Indexed /DeviceRGB 63
  <
    000000 040404 080808 0c0c0c 101010 141414 181818 1c1c1c 202020 242424 
    282828 2c2c2c 303030 343434 383838 3c3c3c 404040 444444 484848 4c4c4c 
    505050 555555 595959 5d5d5d 616161 656565 696969 6d6d6d 717171 757575 
    797979 7d7d7d 818181 858585 898989 8d8d8d 919191 959595 999999 9d9d9d 
    a1a1a1 a5a5a5 aaaaaa aeaeae b2b2b2 b6b6b6 bababa bebebe c2c2c2 c6c6c6 
    cacaca cecece d2d2d2 d6d6d6 dadada dedede e2e2e2 e6e6e6 eaeaea eeeeee 
    f2f2f2 f6f6f6 fafafa ffffff 
  >
] setcolorspace } bdef mwscm
gs np 899 389 mt 0 4227 rl 5362 0 rl 0 -4227 rl cp c np
[5362 0 0 4227 899 389] cc
<< % Image dictionary
  /ImageType 1
  /Width 128 /Height 128
  /BitsPerComponent 8
  /Decode [0 255]
  /ImageMatrix [128.000000 0 0 128.000000 0 0]
  /DataSource currentfile
  /ASCII85Decode filter
  /RunLengthDecode filter
>>
image
IloqJ&260e+XJ0a3AD6T+rgRk-TE^]1`\\'#:MAE(-W<O)(%!^5:8BB(b[KA'GVSo0+pGP1)CBQ1+Xb(
(DRo')BobP)&XD,,rHZ%,:a?Q!X9#^+Uo8(&hY.X1b:@(.1A$K1*mqs(Bkd2!$W0m/.W4D3!h]"#r`7q
%Y$&m5;3W^-6`j0'IO&'-8PiV!":1Y'Eedu-8l`#-6Np?.2=H_1(G9Y,RuRU.2iF40.%c#,rH/W""tZ.
3u0@Z*\&EW.2sor)&4ki-Or0n0EX*X#o47&'.=n>(aCae!?r?i.LIs^)AtLT*BGi+*>T.s,Us!H'cIrE
/0fWk&e-WA-R:3($8ami3B[rh3YrSH,6L(L*uH4I(-4&h3#W#J+V?+.'+5Hn0bQ&G$kj+8/M&+^'*A:n
+=IX.!?!a33?SS@1cb[m2%)2m-k8Ht-Q3R7-n$9&(FftX.4lMG+:nl//f#<J+>+ls4rtLj/Jg&_'JC^U
!"oNd-5SBg.NK0`%4aIi'e(m`,o.pQ+VZ4T*Yo25$90gC%kJe=*Z?LA*$?Ub)F#GV&e$uE/g`S3";MLl
(c=Ac#U'Hm%h^9I"Xaa#(+(.V)C@'c-R/m73Ys,#)DN]P+ql=8-NFnt2CAef/.3^=2(T^N'Gr#.&2jbG
))*-EIhr\?#pD26$4.=o+tONd2@;VU*ASfF&LT(c$P3po0aeX44pO4U.hs?^2$Xp)#UU2o(,%N^(c;jM
,p4'L3>3Sq0,k3N+qY\G.kWFS%h:KM&3VcL/N#+')%JMF3Y2fQ)&+P42%VW-*!\80!=1Fl'+cH756s7Y
!?r0\+b(IL2C:7D,onfO#UUr[3uoXp'0R3].m=IE3$/Yk.00Jj)Aa&4%Q4kI*%;Ob#m2M1"#^N4.KMCG
*=``;*YpC0"Y;534qn/G%150&(,Z=I+X.7E%kAkT2$P'V%NHB'.6TEp.j-o$2CJq\)`&$U$R?oE'GiSe
1bB1N.Oop*#sTRa'IbXV+W(VL+Z2A?+>ONG$TBgS*=4kf$nEbN)]L+I,n1GH"<S@.!@TKE&j%3H-7KZ8
+sS`N1b^$t*"s4o-jV($.L6D,'-mf$)':OA+Uf&I*YKn+*#'VK&2Pph2$,d%#otrQ1DqAV*AeK2.OtNR
4r>^d*]59o+YP/_#!X73+uM854V.f4"V)mb2^CU5)%n/8)Cc.1"r8gH""b&k(c27&&0_bl#6cCY)B:k#
"W]<@!sLIp(aT@j.0(J72$,WP4W#j]$R-Mr+<3Ks.3p2>*>^Lc%Pp!W$T]I^(bm`S0JO7M3$SnW/3+R[
'fm$i+Z^o$IfU;h*?d3[&LeMJ&K`5A+uq(j2@LQW%5o1T4Ym`Z+<hsb(_7Q<!]:2u!t5GS&/,NW.O$i.
.OkNS,8D4Z*Z?"G4tmfm%43J.)[cm0+!(dn/i=OW+sn?U0d.hI)a,9)!C%GF1+Xae.O,NR)Boh!!uhY5
'HSDS58tRo(Om5.58=8?)\"M+,qU,H-S,Q>*^(bM)EJEJ+#"ZR.4HJk1F4.Z0G4Bi#mhY+)EA`p"p-+m
&dBL-.kVY\&-sI--5-k""t2#4+XdmZ#SK$4!?WRL":m@#+=.^1(c<Ei+W_LI0/G-b$6o`g3>!5A1+tC%
,<-AS$nhT3-8'd.,:XZ<+#OER2CpmP,S:.d,<$2W)&j53*B4iB(`F"u)_N<:3%Z?h!ZjQ^-94+:(apgI
*#K_A0cglY/KkcL4sVO(%3[DQ)'^XX2'!VS+rMC="XklF/J/(51*[tV&KWtS-8H8[)B^L<%QPOV(I&KA
"Y_nL+r1Up57dgE,8hRE/.NOC+UTk;)$U?^,o[^..Ld+?'JCpY-ijSi+X7LX$R7e^,Sh:=*$R!g.2"<j
3u0j\2[9@'(+]e:-7Lu0)B];@)]'_H'I,+=+#=<G'0e>s+;m(")E8j$+Y#!#$ofgp'F#`d('c&c3$&h_
,9mm[#o4iu)&XY@Ijbs`(GH^N#UV5I&h5UV-79*S-4U_%(H+2l&eGF"1*\Cu-5@UM.3BNH0,Pl6/M&U^
+q-4:"=c2l+"%+92?l5E(^i5F2D-`u+VHLG&iC=10J"U\"$R)E(^r#8*>9,!+VH[Q*>Kt&%iGu`0Jts@
,U+$I-ncMa/NZ!2%=]$>""+E_0e+aq""kQV3slQ-//L,Q+X&Wp+<^gr*&9;q/M@hQ"Z?i>)](@])*0&h
$9'+>/3O=`.PN&#/-R+H"W%Om1a!Dj.5)ko,Rk^r(E4Fu1'K]V3%bXS!?;de2ErQF"YLi8!%9!A.N9rG
":,Gj+V6Im)_!cV*BATn4s;XD(,m`H-k$;1#o+U@+XK&j+s[g7-R']e//BEk'*B%&(^`#J.j-bm(*GOD
">U8k*>^UQ.4biX3]/2p-TNh+(G%Eo.i&O!$l9s!2Zb3P1D_Au)(>V*#qZK1$4d^^%Q5ai*=O/%2^LR"
*"3$(+=.O0&crpl+>O%K3YX%b(-XAg(]ag^1)qW(#T4s-,UEs6)B^Ua'e1jX+;,520F\p_#:_DG,V0<D
'dah;&cj+!)`])d4;9>-(BHBN)DWZk'e^"W#nnL?&/t^$![\j:2Eqm&-jg>02]F+]#pLW34TH!?+>X!E
(,@Ek"=Fs!2&HuN/e8:,Ij#[[)*Bf$*[1q$/f>0%)$W,J#;ZQC#Tt37#s'Xu!Asm0%N5lc+s\?O.Le$N
-l=Z@0+\j@-T*=q2@(i=-7Koh0Jb%-3\hTW&gBD$3X$09(a&km.k;Ad0e,C54#\K="u.S9+TNTD)EJ'+
,X+@+-8?Yp3BJc;3s6EL.=W!>&g[ZE)$h5h"="Eq%M045*Z,"_-nuGd4<>I[2)>sf*%N9]0.nj\&M#FF
1C5?T'eTqA+#Xo^0a]BA1Goas/cdUb1-H$u&.T+,&0)lK.Q.nc$7-WU&N9t8(E"_N!%ed6*\\$:2_#s_
+ql^D.kN/(0g6pG-p7\n57E#C'-\V?().c?0*ELH2*)$k56iqg/Lhto(G6FW4$Pe13uK5'!uqIe()&\N
&f3Yk*],)C#q$]D2&ngl,n^V3$l1cN$3LbY-mT0T'/1[G*@hU)(,>jm$jSHp(bdQV$m,Ku(Fh$g)@-BP
!A"a:4tnN9)'1UT'g<lh1clV*/MRkD0GbfZ(DeqX#q7)O,m#VH)_W3T,mQCR-6OT\&J>'W56rYc.kias
+p&-6*Y(@Y$9/+i-P.j9![\X/.P<7^/IEjG-O1eK+X\?7)%JS3(d&6B'Fb]d$8aOA%l4M=0cV)I(b7'a
!Z_gi56`nE.MY#I4VoFO/f?\jIk_Wl-8QW#)^HdV$p>mV&d:B=.l9ID+?MS!-7'?S2`1Ud+"&6R(F^.C
!B)bN-QF'L&i:[Y+uptc+XScf.g@g[-O9f(/K#Tb*'5f/3t48+/MS2),9]8f/gClJ#Sf?(/0lVu(C:6R
-T`:J'bq#n3#Mc3,SVL61DMJQ#Ccgm'J]kL(cEKR"q;al0d?c%*@r]r"W%j^5;OiH)E9kh*BG5W#:;JP
#U9il+<;"6-jr-W-NtP8.L[1;4WG@b&hsnZ2Ali(,n;U]!ZM1`(`ro;)A3cA$4nEs4<u'b(-a<)2$5B4
-nR.f+?0o\+Wr^,/fQ_R(cs,d&hAG[!A>3N#TY-2+uUbq%08Hl.i]cj1E6uO(Ejb90GuJd*#o#-.g$VK
,;L,f(`O2L"WJsC0/=k1#W)]I1EeJ(.if!?'.=PH4pW&i(,m9X"?$o6*WS#T#RDRh3#<>\-m9f[2D-:#
/KQYe$Rdnn0.81h+9`\a%6G^n"=kaI&1\.l0dmk[+ZorB(aBe'4TllH+smj<)&XA6)_;gK&i(pb+=7j=
'K[HV-7ofF%iI)]+<<X$2`3TX'JhKo',E);0I&M,+TVQQ&e6uE'fQI21+4"Q%MU!""t1;g&f!DN3uTh&
-kAKe3"HNY)_W?9-m/L1#:pZI3<p6HIhN\$,UalV3%?d$/2970![f`J,q9fB(+gLJ"uRP='+Plh,:=`E
$6prI4;8JH'/^"4+<3BK$4A7*&fqnl$Q1EA*]#DZ-mh#.'b^uj.K]bm!$siI4Y&H10/=Fs)_)jZ"=GKW
1B93U"q;nC/e1*$)C#_=3"A.s.PN"^%t?&Q-6O9H*>(Fp"V2_?#!"(<1c7&n+Y"ut+Y4ud3#OCm.lo14
&.Smb!sf>;.2iR4%l,+Z(E4;-,pkAZ)'Ke;,WZkY!uMpp$lK*m+X\fP'FQ9*.O>f;+?U5g.3p_u"Uth?
1D_5F!=Jl])]8tj.4$Gt'/_NI,rZ#8#RZ&"4#nc<&Ha$o&326W+q-%J2'W5W,r$Sd+r*<[*!e2X"uIqa
%Po7D/.E76(_.WC"r8d6(`*0$1b_'D'GNSH*?,Um&f)Q7'e9Il/0d,53'KGU//B9D-o`.u2(&kQ*?6dl
.1@R0!#Zq4+$2@i&f2Q5"=+F'%kee9(Bt%m3?'FK&fMc;1+aXS1E.]!4t7L"/grOh0a8a;)[6Tg/1)Pg
)]L"F&HW^l(+U:R'*/Cn3uAes,o.jO/g:TN"#i"M2@(ZF2@0[L1^sX%1Cb9C3^465&-rLW-ib)%*=rQ=
.kED(!>GeI*"2j1)?i)J%k/tW3[5Rg#QY)LIj?NR(]u?.)'L"B'fn5`+s8Ec1F<2F1a"aY*siW#2Cfk3