lpr.m

This toolbox contains Matlab files that implement the Laplacian pyramid (LP) decomposition and the

function x = lpr(y, pfilt)
% LPR   Multi-level Laplacian pyramid reconstruction
%
%	x = lpr(y, pfilt)
%
% Input:
%   y:      output of a Laplacian pyramid in a cell vector
%   pfilt:  pyramid filter name (see PFILTERS)
%
% Output:
%   x:      reconstructed signal
%
% See also: LPD
%
% Note:     1-D input signals have to be in column vectors

nd = ndims(y{end});
% Consider column vectors as 1-D signals 
if nd == 2 & size(y{end}, 2) == 1
    nd = 1;
end

% Get the pyramidal filters from the filter name
[h, g] = pfilters(pfilt);

% Decide extension mode
switch pfilt
    case {'9-7', '9/7', '5-3', '5/3', 'Burt'}
        extmod = 'sym';
        
    otherwise
        extmod = 'per';
        
end

x = y{1};

for n = 2:length(y)
    x = lprec1(x, y{n}, nd, h, g, extmod);
end