圖像處理實驗報告二

多媒體實驗報告二班級：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);//讀入數(shù)據(jù) 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文件實驗總結：本次實