多层次树算法的详细解析及其Matlab实现——编码流程主程序

上一篇内容已对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