treemat.m

来自「一种SPIHT的改进算法!!!mallat算法实现,还有待提高!多联系lhz_5」· M 代码 · 共 23 行

function [treePoint,chPoint]=treeMat(r,c)	
global rMat cMat
if (r==1)&&(c==1)
    tp1=[1,2;2,1;2,2];
    treePoint=[tp1;treeMat(1,2);treeMat(2,1);treeMat(2,2)];
elseif (r>rMat/2)||(c>cMat/2)
    treePoint=[];
else
    chPoint=[2*r-1,2*c-1;2*r-1,2*c;2*r,2*c-1;2*r,2*c];
    tm1=[];tm2=[];tm3=[];tm4=[];
    if (chPoint(4,1)<=rMat/2)&&(chPoint(4,2)<=cMat/2)
        [t1,c1]=treeMat(chPoint(1,1),chPoint(1,2));
        tm1=[tm1;t1];
        [t2,c2]=treeMat(chPoint(2,1),chPoint(2,2));
        tm2=[tm2;t2];
        [t3,c3]=treeMat(chPoint(3,1),chPoint(3,2));
        tm3=[tm3;t3];
        [t4,c4]=treeMat(chPoint(4,1),chPoint(4,2));
        tm4=[tm4;t4];
    end
    treePoint=[chPoint;tm1;tm2;tm3;tm4];
end