comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: Help with WV_CWT
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Switch to threaded view of this topic Create a new topic Submit Reply
Re: Help with WV_CWT [message #81654] Thu, 11 October 2012 19:38
Balthasar Indermuehle is currently offline  Balthasar Indermuehle
Messages: 22
Registered: August 2012
Junior Member
On Monday, October 8, 2012 10:25:52 PM UTC-3, Phillip Bitzer wrote:


Awesome hints Phillip, many thanks!

Cheers

- Balt
Re: Help with WV_CWT [message #81680 is a reply to message #81654] Mon, 08 October 2012 18:25 Go to previous message
Phillip M. Bitzer is currently offline  Phillip M. Bitzer
Messages: 7
Registered: September 2012
Junior Member
Balt-

AFAIK, there is no "simple" connection between frequency and scale. It's
not quite as simple as an FFT (although you might argue to ordering of
frequencies returned by IDL isn't exactly simple :-) )

There are a couple things to help you along. First, the IDL Wavelet
toolkit helped me a lot when learning about wavelets:

vis.lbl.gov/NERSC/Software/idl/help/docs6.0/wavelet.pdf

I also suggest you take a peek at the site:

http://paos.colorado.edu/research/wavelets/

Chis Torrence wrote most (all?) of the wavelet routines for IDL. The
article "A Practical Guide to Wavelet Analysis" he is lead author on is
quite useful.

For the immediate problem of how scale relates to frequency, check out
this from the aforementioned web site:

The period (or inverse frequency) is the approximate Fourier period that
corresponds to the oscillations within the wavelet. For all wavelets,
there is a one-to-one relationship between the scale and period. The
relationship can be derived by finding the wavelet transform of a pure
cosine wave with a known Fourier period, and then computing the scale at
which the wavelet power spectrum reaches its maximum.

http://paos.colorado.edu/research/wavelets/faq.html#scale

I recall that I mucked around a lot "under the hood" with the wavelet
routines to really understand them. I also use wv_applet quite a bit to
check my own code.

Hope this helps!
Re: Help with WV_CWT [message #81681 is a reply to message #81680] Mon, 08 October 2012 15:53 Go to previous message
Michael Galloy is currently offline  Michael Galloy
Messages: 1114
Registered: April 2006
Senior Member
On 10/8/12 4:38 PM, Balthasar Indermuehle wrote:
> Further examination of this seems to suggest that the Y scale is in
> fact contained in the returned variable "scales", however, this is an
> assumption as the frequency components it finds in the test signal
> are now closer to where they should be but still not at the exact
> spot. I find no documentation that would evidence that this is what
> it is... wavelets with IDL, the great unknown...?

Maybe you could use the "iTools help system":

$IDL_ROOT/lib/wavelet/source/wv_cwt.pro

Mike
--
Michael Galloy
www.michaelgalloy.com
Modern IDL: A Guide to IDL Programming (http://modernidl.idldev.com)
Research Mathematician
Tech-X Corporation
Re: Help with WV_CWT [message #81682 is a reply to message #81681] Mon, 08 October 2012 15:38 Go to previous message
Balthasar Indermuehle is currently offline  Balthasar Indermuehle
Messages: 22
Registered: August 2012
Junior Member
Further examination of this seems to suggest that the Y scale is in fact contained in the returned variable "scales", however, this is an assumption as the frequency components it finds in the test signal are now closer to where they should be but still not at the exact spot. I find no documentation that would evidence that this is what it is... wavelets with IDL, the great unknown...?
  Switch to threaded view of this topic Create a new topic Submit Reply
Previous Topic: Re: Scale many regions in one image without using xroi
Next Topic: Re: reading netCDF 4.1.1 using IDL version 6.1

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Wed Oct 08 15:39:50 PDT 2025

Total time taken to generate the page: 0.00578 seconds