function [x] = wavedecompress(inputfile)
%   wavedecompress  - Decompression for 1D data 
%
%   [x] = WAVEDECOMPRESS(inputfile) takes the
%   compressed inputfile and returns a reconstructed
%   1-D signal x
%
%   To compress the data, see: wavecompress
%
%   C. Fleizach 16-May-2006
%   cfleizac@cs.ucsd.edu


filter = 'bior5.5';
limit = .0000001;
level = 5;

unzip(inputfile);
delete(inputfile);

inputfile = strrep(inputfile,'.wcm','');

[n,mode,s,wave] = loaddata(inputfile,level);

switch(mode)
    case 1
        filter = 'bior3.1';
        limit = .0001;
        level = 5;
    case 2
        filter = 'bior3.1';
        limit = .00001;
        level = 5;
    case 3
        filter = 'bior5.5';
        limit = .0000001;
        level = 5;
end

[x] = waverec(wave,s,filter);

function [n,mode,s,wave] = loaddata(filename,level)
    fd = fopen(filename , 'r');
    if fd < 0
        disp(['could not open the frame ' filename]);
    end;
    
    n = fread(fd,1,'double');
    mode = fread(fd,1,'double');
    s = fread(fd,(level+2),'double');
    wave = fread(fd,n,'double');
    fclose(fd);

end

end