수다닷컴

질문 제목 : dct계수 추출8x8 행렬로 dct계수만 추출하고 싶어요!질문 내용 : 아까 질문 올렸는데 dct를 계수만 추출하고 싶어요
8x8행렬로요. fopen사용해서 txt로 output까지 했어요.
이 소스코드로 dct idct 그림으로 출력은 모두 잘 되구요.
이제 계수만 정수값으로 얻어내면 되는데
선언할 때 unsigned char형인데 결과값 내야하는건 정수형이라 그런지
메모장에 255만 가득 떠요 ㅠ.ㅠ
이거 어떻게 해결할까요?ㅠㅠ
꼭 8x8로 만들어야해요..ㅠㅠ
아님 그냥 fopen사용하지 않고, printf 사용해서 하는 방법도 괜찮아요 ㅠㅠ소스코드 올려보아요
#include stdio.h
#include math.h
#define pi 3.141592void read_raw(file *);
void write_raw(file *);
void write_dc(file *);
void dct(int, int);
void idct(int, int);
void dct_diff(void);
unsigned char input[256][256];
unsigned char output[256][256];
double temp[256][256];
double output2[256][256];void main(void)
{
int i, j;
file *in, *out1, *out2, *out3;
file *fw; in=fopen(lena1.raw,rb);

if (in==null)
{
printf(file is not found!\n);
} out1=fopen(woman_dct.raw,wb);
out2=fopen(woman_idct.raw,wb);
out3=fopen(woman_diff.raw,wb);
fw=fopen(woman_dc.txt,w);

read_raw(in); for(i=0; i256; i+=8)
for(j=0; j256; j+=8)
dct(i,j);

write_raw(out1);
write_dc(fw);

for(i=0; i256; i+=8)
for(j=0; j256; j+=8)
idct(i,j);

write_raw(out2); dct_diff();

write_raw(out3); fclose(in);
fclose(out1);
fclose(out2);
fclose(out3);
fclose(fw);
}void read_raw(file *in)
{
int i, j; for(i=0;i256;i++)
for(j=0;j256;j++)
input[i][j]=fgetc(in);
}void write_raw(file *out)
{
int i, j; for(i=0;i256;i++)
for(j=0;j256;j++)
fputc(output[i][j],out);
}void write_dc(file *out)
{
int i, j; for(i=0;i256;i+=8)
for(j=0;j256;j+=8)
fprintf(out, %d\t, output[i][j]);
fprintf(out, \n);
}void dct(int i, int j)
{
int a, b, c, d;
double au, av;
double dct_output[8][8]={0,};
for(a=0; a8; a++)
for(b=0; b8; b++)
{
if(a==0) au=1/sqrt(2.0);
else au=1;
if(b==0) av=1/sqrt(2.0);
else av=1;

for(c=0; c8; c++)
for(d=0; d8; d++)
dct_output[a][b]+=input[i+c][j+d]*cos(((2*c+1)*a*pi)/16)*cos(((2*d+1)*b*pi)/16);
dct_output[a][b] *= (au*av)/4;
temp[i+a][j+b] = dct_output[a][b];

if(temp[i+a][j+b] 0)
output[i+a][j+b] = 0;
else if(temp[i+a][j+b] 255)
output[i+a][j+b] = 255;
else
output[i+a][j+b] = (int)dct_output[a][b];
}
}void idct(int i, int j)
{
int a, b, c, d;
double au, av;
double idct_output[8][8]={0,};
for(a=0; a8; a++)
for(b=0; b8; b++)
{
for(c=0; c8; c++)
for(d=0; d8; d++)
{
if(c==0) au=1/sqrt(2.0);
else au=1;
if(d==0) av=1/sqrt(2.0);
else av=1;
idct_output[a][b]+=au*av*temp[i+c][j+d]*cos(((2*a+1)*c*pi)/16)*cos(((2*b+1)*d*pi)/16);
}
idct_output[a][b] /= 4;
output2[output2[i+a][j+b] = idct_output[a][b];

if(output2[i+a][j+b] 0)
output[i+a][j+b] = 0;
else if(output2[i+a][j+b] 255)
output[i+a][j+b] = 255;
else
output[i+a][j+b] = (int)idct_output[a][b];
}
}void dct_diff(void)
{
int i, j;
double temp2;
for(i=0; i256; i++)
for(j=0; j256; j++)
{
temp2=(output2[i][j]-input[i][j])*100000;
if(temp2 0)
output[i][j] = 0;
else if(temp2 255)
output[i][j] = 255;
else
output[i][j]=(int)temp2;
}
}

신청하기