jdcenc.m

灰度图像的jpeg压缩处理

function b=jdcenc(x)
% Usage: b=jdcenc(x)
% Hoffman encoding of DC coefficients in JPEG
% copyright(c) 1997 by Yu Hen Hu
%
% category c = floor(log2(abs(x)))+1
% append code is the binary representation of abs(x) if x>0
% of the 1's complement of bianry rep of abs(x) if x < 0
% created: 11/25/97

% first figure out category number
if x ==0,
   %b=[0 1 0];
   b=[0 0];
   return % done
else
   c = floor(log2(abs(x)))+1; %c在这里相当于AC编码的组号，就是表示直流系数差分值所要用到的位数，而tab的每一行的第一个数，表示的是表示组号要用到的
                              %码长，最终编码的码长是C＋第一列中的选中的数
end


 
tbl=[2 0 0 0 0 0 0 0 0 0   
     3 0 1 0 0 0 0 0 0 0
     3 0 1 1 0 0 0 0 0 0
     3 1 0 0 0 0 0 0 0 0
     3 1 0 1 0 0 0 0 0 0
     3 1 1 0 0 0 0 0 0 0
     4 1 1 1 0 0 0 0 0 0
     5 1 1 1 1 0 0 0 0 0
     6 1 1 1 1 1 0 0 0 0
     7 1 1 1 1 1 1 0 0 0
     8 1 1 1 1 1 1 1 0 0
     9 1 1 1 1 1 1 1 1 0];

 % This is changed by Ravi Lakkundi

b=tbl(c+1,2:tbl(c+1,1)+1);    %这里的组号加1＝c+1，得到的是tab的行坐标。 tbl(c+1,1)是每一行的第一个数据，表示的是码字的长度，这里面得到的b就是查表得到
                              %得到的码字

tmp=int2bin(x,c);           % 这个函数在AC编码中已经详细的讲过了
% tmp is 1 by c+1 vector containing sign-mag
% representation of x, first bit is sign bit.

if tmp(1)==0, % if x > 0
   b=[b tmp(2:c+1)];
elseif tmp(1)==1, % if x < 0
   b=[b ones(1,c)-tmp(2:c+1)];
end