function feats = calcWaveFeats(path,filespec,maxn,L)
% function feats = calcWaveFeats(path,filespec,maxn)
% path is a directory
% filespec is a file spec within that directory (e.g., *.tif)
% maxn is the maximum number of images to process
% L is the number of levels of decomposition desired
%
h=daubcqf(4,'min');
feats=[];
dirls=dir(strcat(path,'/',filespec));
for i=1:min(length(dirls),maxn)
    img=imread(strcat(path,'/',dirls(i).name));
    a=double(img);
    b=(a-min(min(a)))./(max(max(a))-min(min(a)));
    %imshow(b);
    bsiz=min(fix(size(b)/(2^L))*(2^L));
    [y,L]=mdwt(b(1:bsiz,1:bsiz),h,L);
    imshow(y,[min(min(y)) max(max(y))]);
    colormap(prism);
    pause(1);
end