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

Home » Public Forums » archive » Re: wavelet using wv_dwt
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: wavelet using wv_dwt [message #76509] Fri, 17 June 2011 22:07 Go to next message
Lavanya is currently offline  Lavanya
Messages: 7
Registered: August 2010
Junior Member
On Jun 16, 10:17 pm, Konstantinos <moonlightsha...@hotmail.gr> wrote:
> If i understand what you want I will propose yoy this piece of code
>
> pro wavelet_for_forum
>
>     x1 = read_tiff(filepath('image.tif',subdir =
> ['examples','data']))
>     z=128
>     x = congrid(x1, z,z ,/cubic)
>
>     info = wv_fn_daubechies(2, wavelet, scaling,ioff, joff)
>     wv_dwtpartial = wv_dwt(x,wavelet,scaling,ioff,joff, N_LEVELS = 2)
>
>     LL2=wv_dwtpartial[0:z/4-1,0:z/4-1]
>
>     tv, LL2
>
> end
>
> LL2 stands for the LL component of the second level
>
> IF you want wavelet code more than what ITTVIS provides (eg a trous,
> or Daubechies 7/9 wavelet etc) and without the limnitation of the
> image been of the type: image[2^n, 2^n] (eg image[240, 458] without
> the need of congrid) ask me
>
> Kostas

Thanks Kostas, that works perfectly for me.

Please tell me on how to make the wavelet code work when the image is
not ^2. Is there a better way to use congrid?


-
Thanks
Lavanya
Re: wavelet using wv_dwt [message #76540 is a reply to message #76509] Thu, 16 June 2011 10:17 Go to previous messageGo to next message
Konstantinos is currently offline  Konstantinos
Messages: 29
Registered: April 2011
Junior Member
If i understand what you want I will propose yoy this piece of code

pro wavelet_for_forum

x1 = read_tiff(filepath('image.tif',subdir =
['examples','data']))
z=128
x = congrid(x1, z,z ,/cubic)

info = wv_fn_daubechies(2, wavelet, scaling,ioff, joff)
wv_dwtpartial = wv_dwt(x,wavelet,scaling,ioff,joff, N_LEVELS = 2)

LL2=wv_dwtpartial[0:z/4-1,0:z/4-1]

tv, LL2


end

LL2 stands for the LL component of the second level

IF you want wavelet code more than what ITTVIS provides (eg a trous,
or Daubechies 7/9 wavelet etc) and without the limnitation of the
image been of the type: image[2^n, 2^n] (eg image[240, 458] without
the need of congrid) ask me

Kostas
Re: wavelet using wv_dwt [message #76657 is a reply to message #76509] Sat, 18 June 2011 00:58 Go to previous message
Konstantinos is currently offline  Konstantinos
Messages: 29
Registered: April 2011
Junior Member
Lavanya,

Wavelets are filters.
All wavelet filters do not work the same way. (they do not decompose
the image the same way)

IDL wavelet transform use the wavelet (=filter) that described in the
book "section 13.10 of Numerical Recipes in C: The Art of Scientific
Computing, 2nd ed."

Do you have in depth knowledge of wavelets or you just want to use
them?

Give more details of what you want to do and I will help you.

Kostas
  Switch to threaded view of this topic Create a new topic Submit Reply
Previous Topic: Looking for more elegant plot scaling...
Next Topic: Shapefile Opening Error

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

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

Total time taken to generate the page: 0.00669 seconds