Laplacian Pyramid

lappyramid.jl

# generates a laplacian pyramid decomposition of an image A, with nlevels + 1 residual level
function lappyramid(A, nlevels)
  residual = A
  pyramid = ()

  for n = 1:nlevels
    (N, Ns) = impyramid(residual, sizes="both")
    L = residual - N
    residual = Ns
    pyramid = tuple(pyramid..., L)
  end

  return tuple(pyramid..., residual)
end