function [x] = wavedecompress2(inputfile)
%   wavedecompress2  - Decompression for 2D data 
%
%   [x] = WAVEDECOMPRESS2(inputfile) takes the
%   compressed inputfile and returns a reconstructed
%   2-D signal x
%
%   To compress the data, see: wavecompress2
%
%   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,m,mode,s,wave] = loaddata(inputfile,level);

s = reshape(s,level+2,2);

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] = waverec2(wave,s,filter);

function [n,m,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');
    m = fread(fd,1,'double');
    mode = fread(fd,1,'double');
    s = fread(fd,(level+2)*2,'double');
    wave = fread(fd,m*n,'double');
    fclose(fd);

end

end