Re: Unwrapping routine

["R.G.S." <dejastockwell(at)my-dejanews.com>]

Here is a simple phase unwrapping code. (attached)
No attempt is made to fix noisy phase.
A 1D and 2D version is included.

Cheers,
bob




<mvogel@hppcunfmudet.foo.bar> wrote in message
news:8c22jg$juk$1@popeye.eur.nl...
> Dear reader,
>
> I wonder whether there is any 'unwrapping' routine available in IDL 5.3
> The unwrapping routine is thought to correct transitions (e.g Pi) with an
> extendend dynamic range and unwrapping the transition to a value in this
new
> range.
>
> Thank you for your interest,
>
> regards,
>
> Vogel


begin 666 phase_unroll.pro
M.R!U;G)O;&QS('!H87-E(&9U;F-T:6]N(&%M9"!R971U<FYS(&ET#0H[(&-H
M96-K<R!F;W(@8V]M<&QE>"!A;F0@871A;G,@:70-"@T*.R!C;W!Y<FEG:'0@
M4BY'+B!3=&]C:W=E;&P@2&%C:W=A<F4-"@T*#0IF=6YC=&EO;B!P:&%S95]U
M;G)O;&PL=&5M< T*#0IP(#T@=&5M< T*#0II;F9O(#T@<VEZ92AP*0T*:68@
M:6YF;RAI;F9O*# I*S$I(&5Q(#8@=&AE;B!B96=I;@T*"0EP<FEN="PG8V]M
M<&QE>"!D871A+"!C86QC=6QA=&EN9R!A<F=U;65N="<-"@D)<" ](&%T86XH
M:6UA9VEN87)Y*' I+&9L;V%T*' I*0T*96YD:68-"@T*:68@:6YF;R@P*2!E
M<2 R('1H96X@8F5G:6X-"@D[<')I;G0L("=T=V\@9&EM96YS:6]N86P@87)R
M87DG#0H):68@:6YF;R@Q*2!E<2 Q('1H96X@8F5G:6X-"@D).W!R:6YT+'-T
M<F-O;7!R97-S*"=!("<K<W1R:6YG*&EN9F\H,2DI*R<@>" G*W-T<FEN9RAI
M;F9O*#(I*2LG(&%R<F%Y+B<I#0H)"3MP<FEN="PG4F5F;W)M871T:6YG('1O
M(&$@;VYE(&1I;65N<VEO;F%L(&%R<F%Y+B<-"@D)<" ](')E9F]R;2AP*0T*
M"65N9&EF#0IE;F1I9@T*:6YF;R ]('-I>F4H<"D-"FEF(&EN9F\H,"D@97$@
M,2!T:&5N(&)E9VEN#0H);&5N9W1H/6Y?96QE;65N=',H<"D-"@ET<" ](#(J
M(5!I#0H)<&H]," [(&EF(&IU;7 @/" H/BD@4$D@=&AE;B!A9&0@*'-U8G)A
M8W0I('1H92!I;G9A<FEA;G0@,E!)#0H)9F]R(&,],4PL;&5N9W1H+3$@9&\@
M8F5G:6X-"@D@"6EF("@@86)S*' H8RDK<&HM<"AC+3$I*2!G=" A<&DI('1H
M96X@) T*(" @( D):68@*" H<"AC*2MP:BUP*&,M,2DI(&QT(# I('1H96X@
M<&H@/2!P:B K('1P("0-"B @(" @("!E;'-E(&EF("@H<"AC*2MP:BUP*&,M
M,2DI(&=T(# I('1H96X@<&H@/2!P:B M('1P#0H@"7 H8RD@/2!P*&,I("L@
M<&H-"@EE;F1F;W(-"@ER971U<FXL< T*96YD:68@#0H-"F5N9 T*#0H-"@T*
M#0H[(#(@9&EM96YS:6]N86P@=F5R<VEO;@T*9G5N8W1I;VX@<&AA<V5?=6YR
M;VQL7S)D+'1E;7 -"@T*.R!R;W<@,"!D;W=N(&%L;"!C;VQU;6YS(&%N9"!R
M;W<@;VYE(&-O;" Q('1H96X@86QL(&-O; T*#0IT96UP*"HL,"D@/2!P:&%S
M95]U;G)O;&PH=&5M<"@J+# I*0T*#0IC;VP@/2!P:&%S95]U;G)O;&PH=&5M
M<"@P+"HI*0T*9&EF9B ](&-O;" M('1E;7 H,"PJ*0T*:6YF;R ]('-I>F4H
M=&5M<"D-"F9O<B!I(#T@,"QN7V5L96UE;G1S*'1E;7 H*BPP*2DM,2 @9&\@
M8F5G:6X-"@ET96UP*&DL*BD@/2!T96UP*&DL*BD@*R!D:69F#0IE;F1F;W(-
M"@T*9F]R(&D@/2 P+&Y?96QE;65N=',H=&5M<"@P+"HI*2TQ(&1O(&)E9VEN
M#0H)=&5M<"@J+&DI(#T@<&AA<V5?=6YR;VQL*'1E;7 H*BQI*2D-"F5N9&9O
;<@T*#0IR971U<FXL=&5M< T*96YD#0H-"@T*
`
end