function [QT,Theta] = compute_quadtree(M,T,j_min,j_max,s)

% bandelet_transform_fwd - compute the quadtree that optimize the Lagrangian.
%
%   [QT,Theta] = compute_quadtree(M,T,j_min,j_max,s);
%
%   M is the original image
%   T is the selected threshold (the higher, the most compressed the data)
%   j_min is the depth minimum of the QT 
%       (ie 2^j_min is the size minimum of the square).
%   j_max is the depth maximum of the QT        [default : min(5,log2(n))]
%       (ie 2^j_max is the size maximum of the square).
%   s is the super-resolution for the geometry [default 2]
%
%   QT is an image representing the levels of the quadtree.
%   Theta is an image representing the optimal angle choosed on each
%       quadtree (Inf token for no geometry).
%
%   Copyright (c) 2005 Gabriel Peyr