% RES = blurDn(IM, LEVELS, FILT)
%
% Blur and downsample an image.  The blurring is done with filter
% kernel specified by FILT (default = 'binom5'), which can be a string
% (to be passed to namedFilter), a vector (applied separably as a 1D
% convolution kernel in X and Y), or a matrix (applied as a 2D
% convolution kernel).  The downsampling is always by 2 in each
% direction.
%
% The procedure is applied recursively LEVELS times (default=1).

% Eero Simoncelli, 3/97.

function res = blurDn(im, nlevs, filt)

%------------------------------------------------------------
%% OPTIONAL ARGS:

if (exist('nlevs') ~= 1) 
  nlevs = 1;
end

if (exist('filt') ~= 1) 
  filt = 'binom5';
end

%------------------------------------------------------------

if isstr(filt)
  filt = namedFilter(filt);
end  

filt = filt/sum(filt(:));

if nlevs > 1
  im = blurDn(im,nlevs-1,filt);
end

if (nlevs >= 1)
  if (any(size(im)==1))
    if (~any(size(filt)==1))
      error('Cant  apply 2D filter to 1D signal');
    end
    if (size(im,2)==1)
      filt = filt(:);
    else
      filt = filt(:)';
    end
    res = corrDn(im,filt,'reflect1',(size(im)~=1)+1);
  elseif (any(size(filt)==1))
    filt = filt(:);
    res = corrDn(im,filt,'reflect1',[2 1]);
    res = corrDn(res,filt','reflect1',[1 2]);
  else
    res = corrDn(im,filt,'reflect1',[2 2]);
  end
else
  res = im;
end