comp.lang.idl-pvwave archive: archive » Re: unwrap modulo 2pi

Home » Public Forums » archive » Re: unwrap modulo 2pi

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: unwrap modulo 2pi [message #23634 is a reply to message #23629]

Wed, 07 February 2001 13:07

Craig Markwardt
Messages: 1869
Registered: November 1996

Senior Member

"Pavel A. Romashkin" <pavel.romashkin@noaa.gov> writes:

> I just don't see that you can do this at all. Modulo operator discards
> information about the number of multiples of B from A, leaving only the
> remainder of the division operation. For example, if:
>
... deleted ...
>
> Unless I am missing something, I see no way how a unique solution can be
> obtained from the *remainder of division* and *divisor*. But I studied
> arithmetics a long time ago :-)

Pavel, and others--

He is looking for an algorithm that can *reconstruct* the full number.
Of course this will use local information from the surrounding points.

His example is perhaps a little too simple. Try this one:

ph = findgen(100)
ph1 = atan(sin(ph), cos(ph))

Now, PH and PH1 represent the same phase angle on the circle, but PH1
has the disadvantage of being a discontinuous sawtooth function.
Sometimes you want to reconstruct PH based only on knowledge of PH1.
Assuming that the function is nearly monotonic and there is never a
phase jump more than !dpi, I believe this can be done.

The naive solution (which I've never gotten beyond) is to look for
discontinuities in PH1 of (say) more than !dpi. When such a
discontinuity is found, assume we have wound around once, so add
another 2*!dpi to the number.

Craig

--
------------------------------------------------------------ --------------
Craig B. Markwardt, Ph.D. EMAIL: craigmnet@cow.physics.wisc.edu
Astrophysics, IDL, Finance, Derivatives | Remove "net" for better response
------------------------------------------------------------ --------------

Report message to a moderator

[Message index]

		Re: unwrap modulo 2pi By: Pavel A. Romashkin on Wed, 07 February 2001 15:25
		Re: unwrap modulo 2pi By: Craig Markwardt on Wed, 07 February 2001 13:07
		Re: unwrap modulo 2pi By: Pavel A. Romashkin on Wed, 07 February 2001 12:39
		Re: unwrap modulo 2pi By: Randall Skelton on Wed, 07 February 2001 11:48
		Re: unwrap modulo 2pi By: tam on Wed, 07 February 2001 11:18
		Re: unwrap modulo 2pi By: Dave Greenwood on Wed, 07 February 2001 10:48
		Re: unwrap modulo 2pi By: Dave Greenwood on Wed, 07 February 2001 10:29
		Re: unwrap modulo 2pi By: graham_wilson on Wed, 07 February 2001 09:52
		Re: unwrap modulo 2pi By: Pavel A. Romashkin on Wed, 07 February 2001 09:02
		Re: unwrap modulo 2pi By: graham_wilson on Wed, 07 February 2001 07:30
		Re: unwrap modulo 2pi By: Wim Bouwman on Fri, 09 February 2001 00:42
		Re: unwrap modulo 2pi By: Craig Markwardt on Thu, 08 February 2001 07:43
		Re: unwrap modulo 2pi By: Wim Bouwman on Thu, 08 February 2001 01:13

Previous Topic:	IDL programmers wanted
Next Topic:	Re: True Color, but !d.table_size eq 64?

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sat Oct 11 00:57:29 PDT 2025

Total time taken to generate the page: 0.08143 seconds