无损预测编码

%无损预测编码

f=imread('Aligned Matches.tif');
e=mat2lpc(f)%预测编码f
entropy(e)%计算e的一阶估计
c=mat2huff(e);%编码矩阵e
cr=imratio(f,c)%计算f和c的比率
[h,x]=hist(e(:)*512,512);%显示直方图
figure;bar(x,h,'k');
g=lpc2mat(huff2mat(c));%预测解码
compare(f,g)%比较误差