实验三图像编码.docVIP

实验三 图像编码 一 实验目的 1．从编码原理的角度出发，了解图像压缩编码的分类方法。 2．掌握算术编码、预测编码和行程编码的图像压缩编码原理，编程实现各个算法，并比较原始图像和解码图像。 3．掌握变换编码的基本原理。 二 实验原理 三 实验内容 （编写MATLAB文件,要求调试通过，并给出程序清单及代码各个部分的说明，源代码中主要的变量、函数、循环结构等要有注释。） 1. 算术编码 clear all format long e; symbol =[abcd]; ps=[0.4 0.2 0.1 0.3]; inseq=(dacab); codeword=suanshubianma(symbol,ps,inseq) outseq=suanshujiema(symbol,ps,codeword,length(inseq)) function acode=suanshubianma(symbol,ps,inseq) high_range=[]; for k=1:length(ps) high_rabge=[high_range sum(ps(1:k))]; end low_range=[0 high_range(1:length(ps-1))]; sbidx=zeros(size(inseq)); for i=1:length(inseq) sbidx(i)=find(symbol==inseq(i)); end low=0; high=1; for i=1:length(inseq) range=high-low; high=low+range*high_range(sbidx(i)); low=low+range*low_range(sbidx(i)); end acode=low; function symbos=suanshujiema(symbol,ps,codeword,symlen) format long e high_range=[]; for k=1:length(ps) high_range=[high_range sum(ps(1:k))]; end low_range=[0 high_range(1:length(ps)-1)]; psmin=min(ps); symbos=[]; for i=1:symlen idx=max(find(low_range=codeword)); codeword=codeword-low_range(idx); if abs(codeword-ps(idx))0.01*psmin idx=idx+1; codeword=0; end symbos=[symbos symbol(idx)]; codeword=codeword/ps(idx); if abs(codeword)0.01*psmin i=symlen+1; end end 运行结果： 2. 预测编码 X=imread(peppers.bmp,bmp); imshow(X) title(原始图象) X=double(X) Y=Yucejiema(X) XX=Yucejiema(Y) fiqure,imshow(mat2gray(Y)); title(预测误差图象) e=double(X)-double(XX); [m,n]=size(e); erms=sqrt(sum(e(:).^2)/(m*n)); [h,x]=hist(X(:)); figure;bar(x,h,k); title(原图直方图) [h,x]=hist(Y(:)); figure:bar(x,h,k); title(预测误差直方图) XX=uint8(XX); figure,imshow(XX); title(解码图象) function Y=Yucebianma(x,f) error(nargchk(1,2,nargin)) if nargin2 f=1; end x=double(x); [m,n]=size(x); p=zeros(m,n); xs=x; zc=zeros(m,1); for j=1:length(f) xs=[zc xs(:,1:end-1)]; p=p+f(j)*xs; end Y=x-round(p); function x=Yucejiema(Y,f) error(nargchk(1,2,nargin)); if nargin2 f=1; end f=f(end:-1:1); [m,n]=siz