求助帮我修改下程序，这是个利用huffman编码来实现图像压缩的程序，请大虾们帮忙改改
%%%%%%%%数字图像Huffman压缩编码实现%%%%%%%%%%%%%%%%%%%clc; clear; J=imread('cameraman.tif');
imshow(J); title('原始图象'); I=double(J);
% load I;
=size(I);
%%%%%%将矩阵的不同数统计在数组c的第一列中%%%%%%%%%%%%
p1=1;s=m*n;
for k=1:m, for L=1:n
f=0;
for b=1:p1-1, if(c(b,1) == I(k,L)) f=1;break;end; end
if(f==0) c(p1,1)=I(k,L);p1=p1+1;end
end; end
%%%%%%%%%%%将相同的数占站整个数组总数的比例统计在数组p中%%%%%%%
for g=1:p1-1,p(g)=0; c(g,2)=0;
for k=1:m,for L=1:n
if(c(g,1) == I(k,L)) p(g)=p(g)+1;end
end; end
p(g)=p(g)/s;
end
p11=p;
%%%%%%%%找到最小的概率,相加直到等于1,把最小概率的序号存在tree第一列中,次小放在第二列,和放在p像素比例之后%%%%%%%
pn=0;po=1;
while(1)
if(pn >= 1.0)break;
else
=min(p(1:p1-1));p(p2)=1.1;
=min(p(1:p1-1));p(p3)=1.1;
pn=pm+pm2;p(p1)=pn;
tree(po,1)=p2;tree(po,2)=p3;
po=po+1;p1=p1+1;
end
end
%%%%%%%%%C数组第一维表示值,第二维表示代码数值大小,第三维表示代码的位数po=8 g=8
for k=1:po-1
tt=k;m1=1;
if(or(tree(k,1) <= g,tree(k,2) <= g))
if(tree(k,1) <= g)
c(tree(k,1),2)=c(tree(k,1),2)+m1;
m2=1;
while(tt < po-1)
m1=m1*2;
for L=tt:po-1
if(tree(L,1) == tt+g)
c(tree(k,1),2)=c(tree(k,1),2)+m1;
m2=m2+1;tt=L;break;
elseif(tree(L,2) == tt+g)
m2=m2+1;tt=L;break;
end
end
end
c(tree(k,1),3)=m2;
end
tt=k;m1=1;
if(tree(k,2) < g)
m2=1;
while(tt < po-1)
m1=m1*2;
for L=tt:po-1
if(tree(L,1) == tt+g)
c(tree(k,2),2)=c(tree(k,2),2)+m1;
m2=m2+1;tt=L;break;
elseif(tree(L,2) == tt+g)
m2=m2+1;tt=L;break;
end
end
end
c(tree(k,2),3)=m2;
end
end
end
%%%%%%%%%把概率小的值为1标识，概率大的值为0标识%%%%%%
=size(c);
disp('编码')
A1=dec2bin(c(1,2),c(1,3)) %%%说明：这里可以把编码存在高维数组或构架数组、元胞数组同时显示。
A2=dec2bin(c(2,2),c(2,3))
A3=dec2bin(c(3,2),c(3,3)), A4=dec2bin(c(4,2),c(4,3)), A5=dec2bin(c(5,2),c(5,3))
A6=dec2bin(c(6,2),c(6,3)), A7=dec2bin(c(7,2),c(7,3)), A8=dec2bin(c(8,2),c(8,3))
for m=1:M
if (p11(m)~=0) H(m)=-p11(m)*log2(p11(m)); end
end
disp('信源的熵'), H1=sum(H) %信源的熵
NN=0;
for i=1:M
NN=NN+p11(1,i)*c(i,3); %平均码长
end
disp('平均码长'), NN
disp('编码效率'). yita=H1/(NN*log2(2)) %效率
disp('冗余度'), Rd=1-yita %冗余度
[ 本帖最后由 ChaChing 于 2009-9-5 15:05 编辑 ] 请阅读 置顶贴：聚宝盆，然后重新整理你的问题，特别给出出错信息
[ 本帖最后由 ChaChing 于 2009-9-5 14:52 编辑 ]
页:
[1]