spiht.m

一种SPIHT的改进算法!!!mallat算法实现,还有待提高!多联系[e-mail]

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]=spihtcoding(DecIm,codeDim);

% ----- SPIHT Decoding ----- %
DecodeMat=spihtdecoding(T,SnList,RnList,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);
Mod_Runtime=cputime-strtime


% function DecodeMat=spiht(Im,codeDim,decodeDim)
% global rMat cMat
% [rMat,cMat]=size(Im);
% [T,SnList,RnList]=spihtcoding(Im,codeDim);
% DecodeMat=spihtdecoding(T,SnList,RnList,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');