sbind2.m

基于小波变换的特征检索算法

function ind = sbind2(s, n)
% SBIND2 Return indices for accessing a subband from wavelet decomposition
%	IND = SBIND2(S, N)
%
% Input:
%	S:	wavelet decomposition structure (see WAVEDEC2)
%	N:	subband number, 0 for approximation, 
%		1 for H(K), 2 for V(K), 3 for D(K), 4 for H(K-1), ...
%		where K is number of wavelet decomposition levels
%
% Output:
%	IND:	subband indices

if n == 0
    ind = 1:s(1,1)*s(1,2);
else
   q = floor((n-1) / 3);
   r = n - 1 - 3*q;

   length = s(q+2,1)*s(q+2,2);
   first  = s(1,1)*s(1,2) + 3*sum(s(2:q+1,1).*s(2:q+1,2)) + r*length + 1;
   last   = first + length - 1;

   ind = first:last;
end