lis_scan.asv

spiht编码

function [LSP,LIP,LIS,LisFlag,Sn,N]=lis_scan(N,Sn,LSP,LIP,LIS,LisFlag)
% 函数 LIS_SCAN() 检查LIS表的各个表项是否重要，更新列表LIP、LSP、LIS、LisFlag和排序位流 Sn
% 输入参数：N —— 本级编码阈值的指数
%          Sn —— 本级编码排序位流，为空表
%          LSP —— 上一级编码生成的重要系数列表
%          LIP —— 上一级编码生成的不重要系数列表
%          LIS —— 上一级编码生成的不重要子集列表
%          LisFlag —— 上一级编码生成的不重要子集表项类型列表
% 输出参数：LSP —— 本级LIS列表扫描后更新的重要系数列表
%          LIP —— 经本级LIS扫描处理后更新的不重要系数列表
%          LIS —— 本级LIS列表扫描后更新的不重要子集列表
%          LisFlag —— 本级LIS列表扫描后更新的不重要子集表项类型列表
%          Sn —— 本级LIS列表扫描后更新的排序位流
%          N —— 下一级编码阈值的指数

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

% 读入当前 LIS 的表长
rlis=size(LIS,1);
% ls 是指向 LIS 当前表项位置的指针，初始位置为1
ls=1;
% 由于循环过程中列表 LIS 的表长会变化，不适合用 for 循环，故采用 while 循环
while ls<=rlis
    % 读入当前 LIS 表项的类型
	switch LisFlag(ls)
        % ‘D’类表项，包含孩子和非直系子孙
        case 'D'
            % 读入该表项的坐标值
            rP=LIS(ls,1);
            cP=LIS(ls,2);
            % 生成‘D’型子孙树
            chD=coef_DOL(rP,cP,'D');
            % 判断该子孙树是否重要
            isImt=SnOut(chD,N);
            if isImt
                % 如果该子孙树重要，就输入‘1’到 Sn
                Sn=[Sn,1];
                % 生成该表项的孩子树
                chO=coef_DOL(rP,cP,'O');
                % 分别判断每个孩子的重要性
                for r=1:4
                    % 判断该孩子的重要性和正负符号
                    [isImt,Sign]=SnOut(chO(r,:),N);
                    if isImt
                        % 如果重要，就输入‘1’和正负标志到 Sn
                        Sn=[Sn,1];
                        if Sign
                            Sn=[Sn,1];
                        else
                            Sn=[Sn,0];
                        end
                        % 将该孩子添加到重要系数列表 LSP
                        LSP=[LSP;LIP(r,:)];
                    else
                        % 如果不重要，就输入‘0’到 Sn
                        Sn=[Sn,0];
                        % 将该孩子添加到不重要系数列表 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’型子孙树不重要，就输入‘0’到 Sn
                Sn=[Sn,0];
                % 将指针指向下一个 LIS 表项
                ls=ls+1;
            end
            % 更新当前 LIS 的表长，转入下一表项的扫描
            rlis=size(LIS,1);
        case 'L'
            % 对‘L’类表项，不需判断孩子的重要性
            % 读入该表项的坐标值
            rP=LIS(ls,1);
            cP=LIS(ls,2);
            % 生成‘L’型子孙树
            chL=coef_DOL(rP,cP,'L');
            % 判断该子孙树是否重要
            isImt=SnOut(chL,N);
            if isImt
                % 如果该子孙树重要，就输入‘1’到 Sn
                Sn=[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
                % 如果该子孙树是不重要的，就输入‘0’到 Sn
                Sn=[Sn,0];
                % 将指针指向下一个 LIS 表项
                ls=ls+1;
            end
            % 更新当前 LIS 的表长，转入下一表项的扫描
            rlis=size(LIS,1);
	end
end
% 对 LIS 的扫描结束，将本级阈值的指数减1，准备进入下一级编码
N=N-1;