ibafrit.m

为stanford大学donoho教授所编写的redgit变换源代码。是用c编写的

function x = ibafrit(r, m, ind, l, wpar)
% IBAFRIT  Inverse block adaptive finite ridgelet transform
%	x = ibafrit(r, m, ind, l, wpar)
%
% Input:
%	r:	block ridgelet coefficients in a matrix
%	m:	block mean value
%	ind:	block list of projections in which DWT is replaced by DCT
%	l:	number of wavelet decomposition levels
%	wpar:	parameters for spline biorthogonal filters, e.g. [1,3]
%
% Output:
%	x:	reconstructed image.
%
% Note:
%	This uses folded FRAT with symmetric extention and biorthogonal
%	symmetric wavelet transform (WaveLab function)
%
% See also:	IAFRITFD

[a, b] = size(m);
n = size(r, 1) / a;

x = zeros(a * n, b * n);

% Inverse transform bock-by-block
for i = 1:a
    for j = 1:b
	x((1+(i-1)*n):(i*n), (1+(j-1)*n):(j*n)) = ...
	    iafritfd(r((1+(i-1)*n):(i*n), (1+(j-1)*(n+1)):(j*(n+1))), ...
		     m(i, j), ind{i,j}, l, wpar);
    end
end