標(biāo)題: 基于MATLAB，用哈弗曼編碼對(duì)圖像進(jìn)行編碼 [打印本頁(yè)]
作者: 51黑fan    時(shí)間: 2016-1-30 20:25
標(biāo)題: 基于MATLAB，用哈弗曼編碼對(duì)圖像進(jìn)行編碼


clear;
f0=imread('lena.jpg');
subplot(121)
imshow(uint8(f0));
xlabel('\fontsize{16}原始圖像');

f=abs(f0/4)-10;
[M,N]=size(f);

p=zeros(1,61);

for t=1:61
    count=0;
    for i=1:M
        for j=1:N
            if f(i,j)==t-1
                count=count+1;
            end
        end
    end
    p(t)=count;p0=p;
end

core=cell(61,1);
sign=zeros(61);

for hh=1:60
    re=M*N;
    for t=1:61
        if (p(t)<re)&(p(t)>0)
            re=p(t);
        end
    end
    t=1;
    while (p(t)~=re)&(t<61)
        t=t+1;
    end
    if sign(t,1)==0
        core{t}='0';
    else
        core{t}=['0',core{t}];
        i=1;
        while (sign(t,i)~=0)&(i<61)
            core{sign(t,i)}=['0',core{sign(t,i)}];
            i=i+1;
        end
    end
    p(t)=0;
    cou=t;

    re1=M*N;
    for t=1:61
        if (p(t)<re1)&(p(t)>0)
            re1=p(t);
        end
    end
    t=1;
    while (p(t)~=re1)&(t<61)
        t=t+1;
    end
    if sign(t,1)==0
        core{t}='1';
    else
        core{t}=['1',core{t}];
        i=1;
        while (sign(t,i)~=0)&(i<61)
            core{sign(t,i)}=['1',core{sign(t,i)}];
            i=i+1;
        end
    end
    p(t)=p(t)+re;
    cou1=t;
    i=1;
    while (sign(t,i)~=0)&(i<61)
        i=i+1;
    end
    sign(t,i)=cou;
    i=i+1;
    j=1;
    while (sign(cou,j)~=0)&(j<61)
        sign(t,i)=sign(cou,j);
        i=i+1;
        j=j+1;
    end
end  %產(chǎn)生huffman碼

fc=cell(M,N);
for i=1:M
    for j=1:N
        if f(i,j)<61
            fc{i,j}=core{f(i,j)+1};
        else
            fc{i,j}='0';
        end
    end
end %fc
imcore=char();
for i=1:M
    for j=1:N
        imcore=[imcore,fc{i,j}];
    end
end

save picture imcore core; %保存圖片碼流和編碼對(duì)應(yīng)表

  1. clear;
  2. load  picture.mat %載入圖片碼流和編碼對(duì)應(yīng)表
  3. Nc=size(core);
  4. Nic=size(imcore);
  5. flag=0;
  6. i=1;
  7. j=1;
  8. n=1;
  9. cz=char();
  10. f=zeros(128);
  11. for n=1:400930
  12.     if flag==0
  13.         cz=[cz,imcore(n)];
  14.     else
  15.         cz=imcore(n);
  16.         flag=0;
  17.     end
  18.     for t=1:61
  19.         if strcmp(cz,core{t})
  20.             f(j,i)=t;
  21.             flag=1;
  22.             if i>127;
  23.                 i=1;
  24.                 j=j+1;
  25.             else
  26.                 i=i+1;
  27.             end
  28.             break;
  29.         end
  30.     end
  31. end
  32. f=uint8(f*4+35);
  33. subplot(122)
  34. imshow(f);xlabel('\fontsize{16}解碼后的圖像');
復(fù)制代碼




tuxiang.rar (33.29 KB, 下載次數(shù): 7)







歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1