Dual-Tree Complex Wavelet Transform Pack - version 4.3
Nick Kingsbury, Cambridge University, June 2003.
This pack (version 4.2) includes the following functions for performing
the Dual Tree Complex wavelet Transform on 1-D and 2-D signals:
The main DT CWT functions:
dtwavexfm => 1D DTCWT decomposition
dtwaveifm => 1D DTCWT reconstruction
dtwavexfm2 => 2D DTCWT decomposition
dtwaveifm2 => 2D DTCWT reconstruction
Lower level functions:
colfilter => Column filtering of a matrix with symmetric extension
coldfilt => Column filtering with decimation by 2, using 2 filters.
colifilt => Column filtering with interpolation (upsampling) by 2, using 2 filters.
reflect => Reflect a vector about max and min limits (used for sym extension).
There are also equivalent functions for the usual Discrete Wavelet Transform
(included mainly for comparison purposes):
wavexfm => 1D DWT decomposition
waveifm => 1D DWT reconstruction
wavexfm2 => 2D DWT decomposition
waveifm2 => 2D DWT reconstruction
Lower level DWT functions:
coldwtfilt => Column filtering with decimation by 2.
coliwtfilt => Column filtering with interpolation (upsampling) by 2.
There are 2 test M-files which show how the main functions
can be used, and also demonstrate the shift invariance properties
of the DT CWT, compared with the DWT, in both 1-D and 2-D:
shift_test_1D => Demonstrate shift invariance in 1-D
shift_test_2D => Demonstrate shift invariance in 1-D
These M-files use the following ancilliary functions:
draw => Draw an image in a correctly sized figure window.
cimage5 => Draw a complex subimage using a colour palette for the complex numbers.
drawcirc => Generate a circular disc image.
setfig => Set up a predefined figure window.
settitle => Set the title of a figure window.
Various .MAT files contain the wavelet filter coefficients.
To test the 1-D DTCWT package use the following Matlab commands:
X = rand(512,1);
figure; plot(X); drawnow
[Yl,Yh] = dtwavexfm(X,5,'near_sym_b','qshift_b');
Z = dtwaveifm(Yl,Yh,'near_sym_b','qshift_b');
figure; plot(Z)
dtcwt_error = max(abs(Z-X)) % Error should be < 1e-12
To test the 2-D DTCWT package use the following Matlab commands:
load lenna
figure; draw(X); drawnow
[Yl,Yh] = dtwavexfm2(X,4,'near_sym_b','qshift_b');
Z = dtwaveifm2(Yl,Yh,'near_sym_b','qshift_b');
figure; draw(Z)
dtcwt_error = max(abs(Z(:)-X(:))) % Error should be < 1e-12
For more interesting tests, run the M-files shift_test_1D and
shift_test_2D.
********************************
Further information on the DT CWT can be obtained from papers
downloadable from my website (given below). The best tutorial is in
the 1999 Royal Society Paper. In particular this explains the conversion
between 'real' quad-number subimages and pairs of complex subimages.
The Q-shift filters are explained in the ICIP 2000 paper and in more detail
in the May 2001 paper for the Journal on Applied and Computational
Harmonic Analysis.
This code is copyright and is supplied free of charge for research
purposes only. In return for supplying the code, all I ask is that, if
you use the algorithms, you give due reference to this work in any
papers that you write and that you let me know if you find any good
applications for the DT CWT. If the applications are good, I would be
very interested in collaboration. I accept no liability arising from use
of these algorithms.
Nick Kingsbury,
Cambridge University, June 2003.
***********************************************************
Dr N G Kingsbury,
Dept. of Engineering, University of Cambridge,
Trumpington St., Cambridge CB2 1PZ, UK.
or
Trinity College, Cambridge CB2 1TQ, UK.
Phone: (0 or +44) 1223 338514 / 332647; Home: 1954 211152;
Fax: 1223 338564 / 332662; E-mail: ngk@eng.cam.ac.uk
Web home page: http://www.eng.cam.ac.uk/~ngk/
***********************************************************
Version 4.3: Bug corrected at lines 45 and 48 of waveifm2.m
(LoLo instead of Lo). June 2003.
没有合适的资源?快使用搜索试试~ 我知道了~
资源推荐
资源详情
资源评论
收起资源包目录
dtcwt_toolbox4_3.rar (34个子文件)
dtcwt_toolbox4_3
draw.m 2KB
drawcirc.m 505B
coldwtfilt.m 2KB
qshift_c.mat 2KB
qshift_a.mat 1KB
antonini.mat 608B
lenna.mat 76KB
dtwavexfm.m 3KB
qshift_32.mat 3KB
wavexfm.m 2KB
shift_test_1D.m 3KB
near_sym_b.mat 864B
colifilt.m 3KB
colfilter.m 919B
setfig.m 221B
qshift_d.mat 2KB
reflect.m 813B
coldfilt.m 2KB
qshift_b.mat 1KB
coliwtfilt.m 2KB
wavexfm2.m 3KB
qshift_06.mat 1KB
cimage5.m 2KB
waveifm.m 2KB
dtwaveifm2.m 5KB
near_sym_a.mat 544B
shift_test_2D.m 3KB
waveifm2.m 3KB
SETTITLE.M 133B
dtwavexfm2.m 6KB
shiftmovie.m 5KB
dtwaveifm.m 3KB
legall.mat 480B
README.TXT 4KB
共 34 条
- 1
资源评论
- 吉龙2012-07-26二维的就好了!
- LeonDu_MUST2012-07-08写得很详细,很透彻。 学习到了。
sh0thl
- 粉丝: 0
- 资源: 4
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功