dct.cpp

ARM汇编

#include "math.h"
#include <stdio.h>


#define PI (3.1415926536)
double Cm[64][64];
double doubleData[64][64];
unsigned char Data[64][64];

void FDCT(double p[][64])
{
	int i,j,k;
	double s[64][64],temp;
	for(i=0;i<64;i++)
		for(j=0;j<64;j++)
			{
				temp=0.0;
				for(k=0;k<64;k++)
					temp = temp + Cm[i][k]*p[k][j];
				s[i][j]=temp;
			}
	for(i=0;i<64;i++)
		for(j=0;j<64;j++)
			{
				temp=0.0;
				for(k=0;k<64;k++)
					temp+=s[i][k]*Cm[j][k];
				doubleData[i][j]=temp;
			}
}	

void Init_Matrix()
{
	int i,j;
	for(i=0;i<64;i++)
		Cm[0][i]=0.125;
	for(i=1;i<64;i++)
		for(j=0;j<64;j++)
			Cm[i][j]=cos(i*(2*j+1)*PI/128.0)/(sqrt(32.0));
}
int main(int argc, char* argv[])
{
	int i,j;
	Init_Matrix();//系数矩阵
	FILE *fp;
	if( (fp = fopen("test.raw" , "rb")) == NULL)
    {
        printf("Error Opening Input File\n");
        return 1;
    }
    else{
        fseek( fp, 0, SEEK_SET );
        fread(Data, 4096, 1, fp);  }
    fclose(fp);
    
    
    for(i = 0; i < 64; i++)
    {
        for(j = 0; j < 64; j++)
        {
            doubleData[i][j] = Data[i][j]; //doubleData[][]中存放原始图片数据
            printf("%f  ",doubleData[i][j]);
			 printf("\n");
        }
       
    }
    FDCT(doubleData);
    
    
    for(i = 0; i < 64; i++)
    {
        for(j = 0; j < 64; j++)
        {
            if(doubleData[i][j] >255){
                Data[i][j] = 255;}
            else if(doubleData[i][j] < 0){
                Data[i][j] = 0;}
            else{
                Data[i][j] = (unsigned char)doubleData[i][j];}
        }
    }
    if( (fp = fopen("new.raw" , "wb")) == NULL)
    {
        printf("Error Opening Output File\n");
        return 1;
    }
    else{
        fseek( fp, 0, SEEK_SET );
        fwrite(Data, 4096, 1, fp);}
    fclose(fp);
	return 0;
}