基于MATLAB，用哈弗曼编码对图像进行编码

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  %产生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; %保存图片码流和编码对应表

1. clear;
   
2. load  picture.mat %载入图片码流和编码对应表
   
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}解码后的图像');

复制代码

tuxiang.rar (33.29 KB, 下载次数: 7)

2016-1-30 20:24 上传

点击文件名下载附件
下载积分: 黑币 -5