多级树集合分裂（SPIHT）算法的过程详解和Matlab实现（2）数学表述

上一篇文章我们讨论了SPIHT算法与EZW算法的关系，介绍了SPIHT算法的树结构、分集规则和有序表的构建。在此基础上，我们接下来讨论算法的编码原理。下文给出了比较详细的数学描述，吃透了这一过程，就比较容易写出程序代码了。

SPIHT算法的编码过程如下：
（1）初始化
      输出初始阈值T的指数 N = floor ( log2 ( max{| Cr,c |} ) ) （Matlab函数 floor( num ) 给出不大于数值 num 的最大整数）
      定义：  LSP 为空集
                  LIP = {(r,c) | (r,c)∈H }
                  LIS = {D(r,c) | (r,c)∈H 且(r,c)具有非零子孙}
      初始的LIS中各表项类型均为‘D’， LIS 和 LIP 中 (r,c) 的排列顺序与EZW算法零树结构的扫描顺序相同（即按从上到下、从左到右的“Z”型次序排列）。
（2）排序扫描
      1）扫描LIP队列