function RecIm=spiht(Im,imDim,codeDim,decodeDim)

global rMat cMat
strtime=cputime;
[rMat,cMat]=size(Im);
% ----- Image Wavedec ----- %
DecIm=mywavedec2(Im,imDim);

% ----- SPIHT Coding ----- %
[T,SnList,RnList,ini_LSP,ini_LIP,ini_LIS,ini_LisFlag]=spihtcoding(DecIm,imDim,codeDim);

% ----- SPIHT Decoding ----- %
DecodeMat=spihtdecoding(T,SnList,RnList,ini_LSP,ini_LIP,ini_LIS,ini_LisFlag,decodeDim);
% ----- Waverec -----%
m=2^imDim;
DecImLL=DecIm(1:rMat/m,1:cMat/m);
DecodeLL=DecodeMat(1:rMat/m,1:cMat/m);
RecIm=mywaverec2(DecodeMat,imDim,decodeDim);
Runtime=cputime-strtime


% function DecodeMat=spiht(Im,imDim,codeDim,decodeDim)
% global rMat cMat
% [rMat,cMat]=size(Im);
% [T,SnList,RnList,ini_LSP,ini_LIP,ini_LIS,ini_LisFlag]=spihtcoding(Im,imDim,codeDim);
% DecodeMat=spihtdecoding(T,SnList,RnList,ini_LSP,ini_LIP,ini_LIS,ini_LisFlag,decodeDim);
% Im;

% DecIm=importdata('D:\MATLAB\R2007b\work\Wavelet\Imcomprs\Images and Datas\testMat.txt');
% Im=imread('D:\MATLAB\R2007b\work\Wavelet\Imcomprs\Images and Datas\lena.bmp');