palard.m

matlab写的lzw和huffman的编码

function imgdata=palard(img,nn)
%行平面的割，img 为图像的路径输入，nn为要进行的分割层数
%返回imgdata的矩阵为3位数组，第三个下标为不同层的图像，前两个下标为不同层的图像数据，
tu=imread(img);
tu=rgb2gray(tu);
tu=double(tu);
[a,b]=size(tu);
c=zeros(a,b,nn);
for i1=1:a
    for i2=1:b
        y=0;
        s=tu(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:nn
%figure(kk),imshow(c(:,:,kk)) ,title(num2str(i-1)+'层')
%end

imgdata=c;