Hi,
I am currently looking/ comparing different wavelet packages
have to say yours is a very nice package.
however there may be an issue in the correct calculation of the scaling.
In the standard Torrence/Compo wavelet package the default minimum scale (j0) is set to 2
Now according to your code (transform.py line 321) you should do the same (if I understood correctly)
However if I set dt =1 i get s0 = 1.93602661839, NOT 2.0!
Is that a mistake/ a rounding error coming from how you calculate s0 ?
I am using the python code hosted here for comparison (and the idl code also)
http://atoc.colorado.edu/research/wavelets/software.html
I mean, the result is overall the same, but the scales are slightly shifted, which is annoying.
Hi,
I am currently looking/ comparing different wavelet packages
have to say yours is a very nice package.
however there may be an issue in the correct calculation of the scaling.
In the standard Torrence/Compo wavelet package the default minimum scale (j0) is set to 2
Now according to your code (transform.py line 321) you should do the same (if I understood correctly)
However if I set dt =1 i get s0 = 1.93602661839, NOT 2.0!
Is that a mistake/ a rounding error coming from how you calculate s0 ?
I am using the python code hosted here for comparison (and the idl code also)
http://atoc.colorado.edu/research/wavelets/software.html
I mean, the result is overall the same, but the scales are slightly shifted, which is annoying.