实验四图像压缩编码程序研发设计.docVIP

实验名称：图像压缩编码程序设计压缩编码对原始图像数据和解压后的图像数据进行失真图像压缩对一幅数字图像进行压缩编码，计算出压缩率，并显示编码前图像2）对编码后的数据进行解压，并显示解压后图像3）改变算法参数对原始图像数据和解压后的图像数据进行失真clear I=imread(1.jpg); %读入原图像； R=I(:,:,1); R1=im2double(R); %将原图像转为双精度数据类型； 2、进行DCT矩阵变换（重复三次），并显示原图像大小、压缩图像大小、压缩比： T=dctmtx(8); %产生二维DCT变换矩阵 L=blkproc(R1,[8 8],P1*x*P2,T,T); %计算二维DCT，矩阵T及其转置T’是DCT函数P1*x*P2的参数 Mask=[ 1 0 1 0 0 0 0 0 1 1 0 0 1 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]; %二值掩膜，用来压缩DCT系数，只留下DCT系数中左上角的10个 L2=blkproc(L,[8 8],P1.*x,Mask); %只保留DCT变换的10个系数 [m n]=size(L2); J=[m n]; for i=1:m value=L2(i,1); num=1; for j=2:n if L2(i,j)==value num=num+1; else J=[J num value]; num=1; value=L2(i,j); end end J=[J num value ]; end disp(原图像大小) whos(L2); disp(压缩图像大小：) whos(J); disp(图像的压缩比：) disp(m*n/length(J))% 解压缩 t1=J(1); t2=J(2); K(1:t1,1:t2)=0; i1=1; j1=1; for i=3:2:length(J) c1=J(i); c2=J(i+1); for j=1:c1 K(i1,j1)=c2; j1=j1+1; if j1t2 i1=i1+1; j1=1; end end end R2= blkproc(K,[8,8],P1*x*P2,T,T); %逆DCT，重构图像 G=I(:,:,2); G1=im2double(G); %将原图像转为双精度数据类型； T=dctmtx(8); %产生二维DCT变换矩阵 L=blkproc(G1,[8 8],P1*x*P2,T,T); %计算二维DCT，矩阵T及其转置T’是DCT函数P1*x*P2的参数 Mask=[ 1 0 1 0 0 0 0 0 1 1 0 0 1 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]; %二值掩膜，用来压缩DCT系数，只留下DCT系数中左上角的10个 L2=blkproc(L,[8 8],P1.*x,Mask); %只保留DCT变换的10个系数 [m n]=size(L2); J=[m n]; for i=1:m value=L2(i,1); num=1; for j=2:n if L2(i,j)==value num=num+1; else J=[J num value]; num=1; value=L2(i,j); end end J=