palarspara.m

matlab写的lzw和huffman的编码

function px=palarspara(data,n,l)
%把图像按像素的平面分解,返回元胞类型的的矩阵

palar_num=linspace(0,255,n);

palar_num=round(palar_num);

tdata=data;

nn=size(data);

if nargin==2 
for k=1:n-1
 
if k==n
    td=palar_num(k+1);
else
    td=palar_num(k+1)-1;
end
d1=data>=palar_num(k) &data<=td;
px(k)={d1};
end

else
    
[a,b]=size(data);
c=zeros(a,b,n-1);
for i1=1:a
    for i2=1:b
        y=0;
        s=data(i1,i2);
        i=1;    
      while  s~=1 & s~=0
        y=mod(s,2);
        c(i1,i2,i)=y;
        s=fix(s/2);
        i=i+1;
      end
      c(i1,i2,i)=s;
  end  
end

for kk=1:n-1
px(kk)={c(:,:,kk)};
end

end