基于状态位图的一种改进算法SPIHT（2）——编码与解密程序

根据前文所述的改进算法流程，编写了以下核心代码，其余相关函数如childMat()、coef_DOL()、SnOut()、fracnum2bin()等的实现可参考我此前发布的关于SPIHT算法过程详解与Matlab实现的文章，主要代码如下：

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

global rMat cMat

% ----- 起始时间 ----- %
strtime=cputime;

[rMat,cMat]=size(Im);
% ----- 小波分解 ----- %
DecIm=mywavedec2(Im,imDim);
% ----- M_SPIHT 编码过程 ----- %
[T,SnList,RnList]=spihtcoding(DecIm,codeDim);
% ----- M_SPIHT 解码过程 ----- %
DecodeMat=spihtdecoding(T,SnList,RnList,decodeDim);
% ----- 小波重构 -----%
RecIm=mywaverec2(DecodeMat,imDim,decodeDim);

% ----- 运行时间统计 -----%
runtime=cputime-