discrete wavelet transform

Ivo Alxneit ivo.alxneit@psi.ch
Sun Mar 7 15:56:00 GMT 2004

Hash: SHA1


i have put together functions to calculate the discrete wavelet transform in 
1D (gsl_vector_double) and 2D (gsl_matrix_double). i have already implemented 
different wavelets (beylkin, coiflet, daubechies, haar).

i have verified the output against paul addisons illustrated wavelet transform 
book and matlabs wavekit. the numerical output is the same (minus 
rounding/numerics) but note that the position of the T submatrices occuring 
at each levels is different in the book compared with the output of wavekit. 
the position of Th and Tv (the off-diagonal blocks are exchanged). the output 
of my implemetation is the same as in the book.

is this something that could be included in gsl? if yes i would still need to 
write some tests (suggestions welcome) and documentation.

the attached archive contains the code and illustrations for 1D and 2D.
- -- 
Dr. Ivo Alxneit
Laboratory for Solar Technology   phone: +41 56 310 4092
CH-5232 Villigen                    fax: +41 56 310 2624
Paul Scherrer Institute          http://solar.web.psi.ch
Switzerland                        gnupg key: 0x515E30C7
Version: GnuPG v1.2.1 (GNU/Linux)

-------------- next part --------------
A non-text attachment was scrubbed...
Name: gsl_wavelet.tar.gz
Type: application/x-tgz
Size: 6310 bytes
Desc: not available
URL: <http://sourceware.org/pipermail/gsl-discuss/attachments/20040307/04bf777b/attachment.bin>

More information about the Gsl-discuss mailing list