多级树集合分裂（SPIHT）算法的过程详解和Matlab实现（3）编码过程——主程序

上一篇文章已经详细介绍了SPIHT算法的编码过程，接下来有关编码和解码的部分就直接把代码写出来啦，我的代码里有详细的中文注释，基本上把程序的每个步骤都作了说明，呵呵，利人也利己！

1、首先给出编码的主程序

function [T,SnList,RnList,ini_LSP,ini_LIP,ini_LIS,ini_LisFlag]=spihtcoding(DecIm,imDim,codeDim)
% 函数 SPIHTCODING() 是SPIHT算法的编码主程序
% 输入参数：DecIm   ——小波分解系数矩阵；
%                     imDim   ——小波分解层数；
%                     codeDim ——编码级数。
% 输出参数：T —— 初始阈值，T=2^N，N=floor(log2(max{|c(i,j)|}))，c(i,j)为小波系数矩阵的元素
%                     SnList —— 排序扫描输出位流
%                     RnList —— 精细扫描输出位流
%                     ini_L* —— 初始系数（集合）表
%                     LSP：重要系数表
%                     LIP：不重要系数表
%                     LIS：不重要子集表，其中的表项是D型或L型表项的树根点
%                     LisFlag：LIS中各表项的类型，包括D型和L型两种

global Mat rMat cMat
% Mat是输入的小波分解系数矩阵，作为全局变量，在编码的相关程序中使用
% rMat、cMat是Mat的行、列数，作为全局变量，在编码、解码的相关程序中使用

%---------------------------%
% ----- Threshold ----- %
%---------------------------%
Mat=DecIm;
MaxMat=max(ma