mat2huff.m

包含完整的程序代码

function y=mat2huff(x)
if ndims(x)~=2 | ~isreal(x) | (~isnumeric(x) & -islogical(x))
error('x必须是一个2D的实矩阵或是逻辑矩阵');
end
% 存储x的大小
y.size=uint32(size(x));
%找出x值的范围
x=round(double(x));
xmin=min(x(:));
xmax=max(x(:));
pmin=double(int16(xmin));
pmin=uint16(pmin+32768);
y.min=pmin;
%在xmin和xmax中计算输入的直方图
x=x(:)';
h=histc(x,xmin:xmax);
if max(h)>65535
h=65535*h/max(h);
end
h=uint16(h);
y.hist=h;
% 计算输入矩阵的编码并存储

  %生成huffman编码图
map=huffman(double(h));
hx=map(x(:)-xmin+1);
hx=char(hx);
hx=hx(:)';
hx(hx==' ')=[];
%计算编码大小
ysize=ceil(length(hx)/16);
hx16=repmat('0',1,ysize*16);
hx16(1:length(hx))=hx;
hx16=reshape(hx16,16,ysize);
%把二进制转化成十进制
hx16=hx16'-'0';
twos=pow2(15:-1:0);
y.code=uint16(sum(hx16.*twos(ones(ysize,1),:),2))';