图像处理实验报告二

多媒体实验报告二班级：10011107学号：2011302542姓名：孔莹莹实验目的：使用2DDCT变换将testimg.txt转换到频域，并输出结果。要求以8x8为处理单元。实验要求：读入testimg.Txt;把testimg.Txt进行8x8的分块处理；对每一个分块进行2DDCT变换；把变换结果写入testdct.txt文件中；对testidct.txt进行分块处理对每一个分块进行IDCT变换；把变换结果写入testidct.txt文件中；实验流程图：开始开始读入读入testimg.txt结束是否读入False结束是否读入True把把testimg.Txt进行8x8的分块处理对每一个分块进行对每一个分块进行DCT变换把结果写入把结果写入testdct.txt中对对testdct.txt进行分块处理对每一个分块进行对每一个分块进行IDCT变换把变换结果写入把变换结果写入testidct.txt文件中结束结束实验代码：//DCT变换intDCT(intx,inty){ inti,j,v,u; doubles; doublec[8][8]; for(i=0;i<8;i++) { s=(i==0)?sqrt(0.125):0.5; for(j=0;j<8;j++) { c[i][j]=s*cos((PI/8)*i*(j+0.5)); } } doubletemp[8][8]; doubler[8][8]; for(i=0;i<8;i++)//一维DCT变换 { for(v=0;v<8;v++) { s=0.0; for(j=0;j<8;j++) s+=c[i][j]*code[j][v]; temp[i][v]=s; } } for(v=0;v<8;v++)//二维DCT变换 { for(u=0;u<8;u++) { s=0.0; for(i=0;i<8;i++) s+=temp[v][i]*c[u][i]; r[v][u]=s; } } s=r[x][y]; returns;}//IDCT变换intIDCT(intx,inty){ intp=0; inti=0,j=0; doublel=0,m,n,t; t=1/sqrt(2); for(i=0;i<8;i++) { for(j=0;j<8;j++) { m=cos(((2*x+1)*i*3.14159)/(16)); n=cos(((2*y+1)*j*3.14159)/(16)); if(i==0&&j==0) { l=l+m*n*0.5*code[i][j]; } elseif(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); returnp;}//对testimg进行分块并做dct变换处理 for(i=0;i<64;i++)//分成8*512的块 { for(j=0;j<8;j++) { for(k=0;k<512;k++) { fscanf(fp1,"%d",&node);//读入数据 code1[j][k]=node; } } for(l=0;l<64;l++) { for(m=0;m<8;m++) { for(n=0;n<8;n++) { code[m][n]=code1[m][l*8+n]; } } for(m=0;m<8;m++) { for(n=0;n<8;n++) { node=DCT(m,n);//对每一个点进行DCT变换 fprintf(fp2,"%d\t",node); } fprintf(fp2,"\n"); } } }//对testdct.txt进行分块处理并进行IDCT变换 for(i=0;i<64;i++) { for(h=0;h<64;h++) { for(j=0;j<8;j++) { for(k=0;k<8;k++) { fscanf(fp2,"%d",&node); code[j][k]=node; } } for(j=0;j<8;j++) { for(k=0;k<8;k++) { p=IDCT(j,k); code1[j][8*h+k]=p; } } } for(j=0;j<8;j++) { for(k=0;k<512;k++) { fprintf(fp3,"%d\t",code1[j][k]); } fprintf(fp3,"\n"); } }实验结果：未进行处理的testimg.txt文件进行了DCT变换的dct变换的testdct文件对testdct进行了IDCT变换的testidct文件实验总结：本次实