log2linchap4.m

斯坦福大学的小波的工具箱 比matlab自带的工具箱更完善

function linim = log2lin(logim,nvoice);
%
% function linim = log2lin(logim,nvoice);
%
% Description : takes a log-scale image l , and transforms it into
%               a linear scale image n
%
% Inputs :	logim : logscale image
%              	nvoice : number of voices (default = 1)
%
% Outputs :	linim : linear scale image
%
% See also : CWT, ImageCWT
%
   if nargin < 2
	nvoice =1;
   end
%
 
   [mn,n] = size(logim);
   m= mn/nvoice;
   s0 = n/(2^m);
   j0 = log2(s0);
    s1 = n;
   max = floor(n/2^(1/nvoice));
% normally, linim has size (max-s0+1,n);
% but this is too big (out of memory for n = 1024)
% so we replace it by (floor((max-s0+1)/10),n);

   sz = floor((max-s0)/4)+1;
%sz = floor((max-s0)/1)+1;
   linim = zeros(sz,n);
    for j = 1:sz,
        i = (j-1)*4 + s0*2^(3/nvoice);
%	i = (j-1) + s0;
%   for i = s0:max,
 	li = (log2(i)-j0) * nvoice;
	fli = floor(li);
	rli = li - fli;
%	if rli ~= 0,
%	    linim(j,:) = rli.*logim(fli+2,:);
%        end; % if
%	if rli ~= 1,
%    	    linim(j,:) =linim(j,:)+(1-rli).*logim(fli+1,:);
%	end; % if
 	    linim(j,:) = quadinterpol(logim(fli-1,:),logim(fli,:),logim(fli+1,:),rli);
   end % i loop 


% Written by Maureen Clerc and Jerome Kalifa, 1997
% [e-mail], [e-mail]
    
    
%   
% Part of WaveLab Version 802
% Built Sunday, October 3, 1999 8:52:27 AM
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%