14020181051 杨宁 数据压缩作业3.docx

数据压缩作业3班级： 1402018 姓名：杨宁学号： 14020181051 一．摘要本次作业要求使用一种有损图像压缩方法，与上次无损压缩相比，压缩效果明显提升，我使用的是离散余弦变化(DCT)对图像进行压缩和解压缩处理。关键词：有损图像压缩，DCT，压缩，解压二．实验目的(1)编写DCT的压缩程序对图像(raw)进行压缩处理；(2)编写DCT的解压缩程序对压缩图像(raw)进行解压缩处理。三．实验内容(1)利用DCT变换压缩方法对256*256的图像(lena256x256.raw)进行压缩，压缩结果储存在lena256x256_code.raw文件中(2)编写相应的DCT解码程序对lena256x256_code.raw进行解压缩，恢复结果储存在lena256x256_decode.raw文件中。四．算法流程(1)读取256*256的灰度图像；(2)初始化8*8变换矩阵C；(3)将图像分成8*8的小像素块U；(4)进行DCT正变换，V=CU(5)将V中左上角10个dct系数保留；(6)将压缩编码依次输入至lena256x256_code.raw五．程序压缩性能评价源文件：压缩文件：压缩比例：压缩性能为：[(64-20)/64]·100％=68.75％对比无损压缩7%的压缩性能，有损压缩的压缩性能有了很大的提高，但其缺点是由压缩文件恢复的源文件有所损失。六．程序源代码(1)压缩程序DCT_Code.c#includestdio.h#includemath.hshort q[15],dct[8][8];unsigned char image[256][256];double temp,b[8],c[8][8],a[8][8];int i,j,u,v,t,l;FILE *fp;int r[8][8]={1,1,1,1,0,0,0,0, 1,1,1,0,0,0,0,0, 1,1,0,0,0,0,0,0, 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};int main(){fp=fopen(Lena.raw,rb);fread(image,sizeof(unsigned char),(256*256),fp);fclose(fp);fp=fopen(compress.raw,wb);for(i=0;i8;i++) for(j=0;j8;j++) c[i][j]=cos((2*i+1)*j*3.14159/16); b[0]=0 for(i=1;i8;i++) b[i]=0.5; for(t=0;t32;t++){for(l=0;l32;l++){for(i=0;i8;i++)for(j=0;j8;j++)a[i][j]=image[t*8+i][l*8+j];for(u=0;u8;u++) for(v=0;v8;v++){ temp=0.0; for(i=0;i8;i++) for(j=0;j8;j++) temp+=c[i][u]*c[j][v]*(a[i][j]-128); temp*=b[u]*b[v]; dct[u][v]=(short)temp;}u=0;for (i=0;i8;i++)for (j=0;j8;j++){if (r[i][j]==1){q[u]=dct[i][j];u++;}}fwrite(q,sizeof(short),10,fp);}}fclose(fp);return 0;}(2)解压程序DCT_Decode.c#includestdio.h#includemath.hint i,j,k,l,u,v;double c[8][8],b[8],temp;unsigned char image[256][256],a[8][8];int r[8][8]={1,1,1,1,0,0,0,0, 1,1,1,0,0,0,0,0, 1,1,0,0,0,0,0,0, 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};short q[15],dct[8][8]={0};FILE *fp1,*fp2;int main(){fp1=fopen(compress.raw,rb);fp2=fopen(restore.raw,wb);for(j=0;j8;j++) for(i=0;i8;i++) c[i][j]=cos((2*j+1)*i