图像处理实验报告二.docxVIP

多媒体实验报告二班级 学号：2011302542姓名：孔莹莹实验目的： 使用2D DCT变换将testimg.txt转换到频域，并输出结果。要求以8x8为处理单元。实验要求：读入testimg.Txt;把testimg.Txt进行8x8的分块处理；对每一个分块进行2D DCT变换；把变换结果写入testdct.txt文件中；对testidct.txt进行分块处理对每一个分块进行IDCT变换；把变换结果写入testidct.txt文件中；实验流程图： 是否读入 False True 实验代码：//DCT变换int DCT(int x,int y){int i,j,v,u;double s;double c[8][8];for(i=0;i8;i++){s=(i==0)?sqrt(0.125):0.5;for(j=0;j8;j++){c[i][j]=s*cos((PI/8)*i*(j+0.5));}}double temp[8][8];double r[8][8];for(i=0;i8;i++)//一维DCT变换{for(v=0;v8;v++){s=0.0;for(j=0;j8;j++)s+=c[i][j]*code[j][v];temp[i][v]=s;}}for(v=0;v8;v++)//二维DCT变换{for(u=0;u8;u++){s=0.0;for(i=0;i8;i++)s+=temp[v][i]*c[u][i];r[v][u]=s;}}s=r[x][y];return s;}//IDCT变换int IDCT(int x,int y){int p=0;int i=0,j=0;double l=0,m,n,t;t=1/sqrt(2);for (i=0;i8;i++){for (j=0;j8;j++){m=cos(((2*x+1)*i*3.14159)/(16));n=cos(((2*y+1)*j*3.14159)/(16));if(i==0j==0){l=l+m*n*0.5*code[i][j];}else if(i==0||j==0){l=l+m*n*t*code[i][j];}else{l=l+m*n*code[i][j];}}}p=int(l*0.25);return p;}//对testimg进行分块并做dct变换处理for(i=0;i64;i++)//分成8*512的块{for(j=0;j8;j++){for(k=0;k512;k++){fscanf(fp1,%d,node);//读入数据code1[j][k]=node;}}for(l=0;l64;l++){for(m=0;m8;m++){for(n=0;n8;n++){code[m][n]=code1[m][l*8+n];}}for(m=0;m8;m++){for(n=0;n8;n++){node=DCT(m,n);//对每一个点进行DCT变换fprintf(fp2,%d\t,node);}fprintf(fp2,\n);}}}//对testdct.txt进行分块处理并进行IDCT变换for(i=0;i64;i++){for(h=0;h64;h++){for(j=0;j8;j++){for(k=0;k8;k++){fscanf(fp2,%d,node);code[j][k]=node;}}for(j=0;j8;j++){for(k=0;k8;k++){p=IDCT(j,k);code1[j][8*h+k]=p;}}}for(j=0;j8;j++){for(k=0;k512;k++){fprintf(fp3,%d\t,code1[j][k]);}fprintf(fp3,\n);}}实验结果：未进行处理的testimg.txt文件进行了DCT变换的dct变换的testdct文件对testdct进行了IDCT变换的testidct文件实验总结：本次实验中遇到的难点是如何对文件进行正确的分块处理，并把dct处理结果正确的输出，实验中我采用的方法是：对512x512的原文件先分成8x512的块，然后再把8x512的块分成8x8的块，此时对每一个块进行DCT变换，并把处理结果输出在testdct文件中。