perform_atrou_transform.m

小波提升算法的实现

function y = perform_atrou_transform(x,Jmin,options)

% perform_atrou_transform - compute the "a trou" wavelet transform, 
%   i.e. without subsampling.
%
%   w_list = perform_atrou_transform(M,Jmin,options);
%
%   'w_list' is a cell array, w_list{ 3*(j-Jmin)+q }
%   is an imagette of same size as M containing the full transform
%   at scale j and quadrant q.
%   The lowest resolution image is in w_list{3*(Jmax-Jmin)+4} =
%   w_list{end}.
%
%   The ordering is :
%       { H_j,V_j,D_j, H_{j-1},V_{j-1},D_{j-1}, ..., H_1,V_1,D_1, C }
%
%   'options' is an (optional) structure that may contains:
%       options.wavelet_type: kind of wavelet used (see perform_wavelet_transform)
%       options.wavelet_vm: the number of vanishing moments of the wavelet transform.
%
%   Copyright (c) 2004 Gabriel Peyr