Re: wavelet using wv_dwt [message #76509] |
Fri, 17 June 2011 22:07  |
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   |
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  |
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
|
|
|