수다닷컴

질문 제목 : dct 코드 관련dct 코드를 분석중입니다. dct가 어떤 역할을 하는 모듈인지는 대략적으로 알고있는데,
분석 도중 이해가 되지 않는 부분이 있어서요 ㅜㅜ 막히네요.. 아침부터.. 이 함수만 붙잡고 있었어요 ㅜㅜ...
동영상 압축 프로세실 과정 중 dct 함수에 대하여 잘 아시는분 계시면 답변 부탁드려요.,질문 내용 :
int dct(short* inblock, short* outblock)
{
const short f0=32767;
const short f1=32138;
const short f2=30273;
const short f3=27245;
const short f4=23170;
const short f5=18204;
const short f6=12539;
const short f7=6392;
int j1, i, j;
short b[8];
short b1[8];
short d[8][8];
short *inblock_addr = inblock;
short *outblock_addr = outblock;
short *d_addr;

d_addr = &(d[0][0]);
for (i = 0; i 8; i++, d_addr+=8) {
for (j = 0; j 8; j++) {
b[j] = (*inblock_addr++);
}

b1[0] = b[0] + b[7];
b1[7] = b[0] - b[7];
b1[1] = b[1] + b[6];
b1[6] = b[1] - b[6];
b1[2] = b[2] + b[5];
b1[5] = b[2] - b[5];
b1[3] = b[3] + b[4];
b1[4] = b[3] - b[4];

b[0] = b1[0] + b1[3];
b[1] = b1[1] + b1[2];
b[2] = b1[1] - b1[2];
b[3] = b1[0] - b1[3];
b[4] = b1[4];

b[5] = ((b1[6] - b1[5]) * f0)16;
b[6] = ((b1[6] + b1[5]) * f0)16;
b[7] = b1[7];

*(d_addr) = ((b[0] + b[1]) * f4)18;
*(d_addr+4) = ((b[0] - b[1]) * f4)18;
*(d_addr+2) = (b[2] * f6 + b[3] * f2)18;
*(d_addr+6) = (b[3] * f6 - b[2] * f2)18;

b1[4] = b[4] + b[5];
b1[7] = b[7] + b[6];
b1[5] = b[4] - b[5];
b1[6] = b[7] - b[6];

*(d_addr+1) = (b1[4] * f7 + b1[7] * f1)18;
*(d_addr+5) = (b1[5] * f3 + b1[6] * f5)18;
*(d_addr+7) = (b1[7] * f7 - b1[4] * f1)18;
*(d_addr+3) = (b1[6] * f3 - b1[5] * f5)18;
}

진한 글자 부분의 코드가 왜 저렇게 작성되었는지 알려주시면 감사하겠습니다.
f0~f7이 들어오는 inblock을 계수행렬로 바꾸기 위한 변환계수의 역할을 하는것인지,
그리고 b[0]에서 b[7]까지 위에 작업 처럼 inblock 행렬이 좌우대칭이라 행렬을 또 줄이는 작업을 하는것인지...
shift 16과 shift 18은 왜 하는것이고, 최종적으로 배열b마다 f상수를 다 곱하는데, 정말 계수의 역할이라 저런 작업을 하는것인지..
dct 역할과 알고리즘 흐름정도만 알았지 ㅠㅠ 코드로 보니 왜 저렇게 짠건지 모르겠어요 ㅠㅜ

신청하기