lis_decode.m

spiht编码

function [LSP,LIP,LIS,LisFlag,DecodeMat,N]=lis_decode(DecodeMat,N,Sn,LSP,LIP,LIS,LisFlag)
% 函数 LIS_DECODE() 根据排序位流 Sn，更新列表LIP、LSP和重构系数矩阵 DecodeMat
% 输入参数：DecodeMat —— 本级 LIP 扫描后更新的重构系数矩阵
%          N —— 本级解码阈值的指数
%          Sn —— 经本级 LIP 扫描截取后的排序位流
%          LSP —— 经本级 LIP 扫描后更新的重要系数列表
%          LIP —— 经本级 LIP 扫描后更新的不重要系数列表
%          LIS —— 上一级编码生成的不重要子集列表
%          LisFlag —— 上一级编码生成的不重要子集表项类型列表
% 输出参数：LSP —— 本级 LIS 扫描后更新的重要系数列表
%          LIP —— 经本级 LIS 扫描处理后更新的不重要系数列表
%          LIS —— 本级 LIS 扫描后更新的不重要子集列表
%          LisFlag —— 本级 LIS 扫描后更新的不重要子集表项类型列表
%          DecodeMat —— 本级 LIS 扫描后更新的重构系数矩阵
%          N —— 下一级解码阈值的指数

global rMat cMat
% rMat、cMat是Mat的行、列数，作为全局变量，在编码、解码的相关程序中使用

% 读入当前 LIS 的表长
rlis=size(LIS,1);
% ls 是指向 LIS 当前表项位置的指针，初始位置为1
ls=1;
while ls<=rlis
    % 读入当前 LIS 表项的类型
	switch LisFlag(ls)
        % ‘D’类表项，包含孩子和非直系子孙
        case 'D'
            % 读入该表项的坐标值
            rP=LIS(ls,1);
            cP=LIS(ls,2);
            % 根据 Sn 判断该表项‘D’型子孙树是否重要
            if Sn(1)==1
                % 每次判断都是读入 Sn 的首位数，判断后立即删除这一位数
                Sn(1)=[];
                % 生成该表项的孩子树
                chO=coef_DOL(rP,cP,'O');
                % 分别判断每个孩子的重要性
                for r=1:4
                    % 读入孩子的坐标值
                    rO=chO(r,1);
                    cO=chO(r,2);
                    % 判断该孩子的重要性
                    if Sn(1)==1
                        Sn(1)=[];
                        % 判断该孩子的正负符号
                        if Sn(1)==1
                            Sn(1)=[];
                            % 生成该孩子的系数值
                            DecodeMat(rO,cO)=1.5*2^N;
                        else
                            Sn(1)=[];
                            DecodeMat(rO,cO)=-1.5*2^N;
                        end
                        % 将该孩子添加到重要系数列表 LSP
                        LSP=[LSP;chO(r,:)];
                    else
                        % 如果不重要，则这个孩子的系数值为 0
                        DecodeMat(rO,cO)=0;
                        Sn(1)=[];
                        % 将该孩子添加到不重要系数列表 LIP
                        % 本级阈值下不重要的系数在下一级解码中可能是重要的
                        LIP=[LIP;chO(r,:)];
                    end
                end
                % 生成该表项的非直系子孙树
                chL=coef_DOL(rP,cP,'L');
                if ~isempty(chL)
                    % 如果‘L’型树非空，则将该表项添加到列表 LIS 的尾端等待扫描
                    LIS=[LIS;LIS(ls,:)];
                    % 表项类型更改为‘L’型
                    LisFlag=[LisFlag,'L'];
                    % 至此，该表项的‘D’型LIS扫描结束，在LIS中删除该项及其类型符
                    LIS(ls,:)=[];
                    LisFlag(ls)=[];
                else
                    % 如果‘L’型树为空集
                    % 则该表项的‘D’型LIS扫描结束，在LIS中删除该项及其类型符
                    LIS(ls,:)=[];
                    LisFlag(ls)=[];
                end
            else
                % 如果该表项的‘D’型子孙树不重要，则先删除读入的 Sn 值
                Sn(1)=[];
                % 然后将指针指向下一个 LIS 表项
                ls=ls+1;
            end
            % 更新当前 LIS 的表长，转入下一表项的扫描
            rlis=size(LIS,1);
        case 'L'
            % 对‘L’类表项，不需判断孩子的重要性
            % 读入该表项的坐标值
            rP=LIS(ls,1);
            cP=LIS(ls,2);
            % 判断该表项的‘L’型子孙树是否重要
            if Sn(1)==1
                % 如果该子孙树重要
                Sn(1)=[];
                % 则生成该表项的孩子树
                chO=coef_DOL(rP,cP,'O');
                % 将该‘L’类表项从 LIS 中删除
                LIS(ls,:)=[];
                LisFlag(ls)=[];
                % 将表项的四个孩子添加到 LIS 的结尾，标记为‘D’类表项
                LIS=[LIS;chO(1:4,:)];
                LisFlag(end+1:end+4)='D';
            else
                % 如果该表项的‘L’型子孙树不重要，则先删除读入的 Sn 值
                Sn(1)=[];
                % 然后将指针指向下一个 LIS 表项
                ls=ls+1;
            end
            % 更新当前 LIS 的表长，转入下一表项的扫描
            rlis=size(LIS,1);
	end
end
% 对 LIS 的扫描结束，将本级阈值的指数减1，准备进入下一级编码
N=N-1;